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.

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.

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.

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.

Immunotoxicity of nanoparticles: a computational study suggests that CNTs and C60 fullerenes might be recognized as pathogens by Toll-like receptors.

Over the last decade, a great deal of attention has been devoted to study the inflammatory response upon exposure to multi/single-walled carbon nanotubes (CNTs) and different fullerene derivatives. In particular, carbon nanoparticles are reported to provoke substantial inflammation in alveolar and bronchial epithelial cells, epidermal keratinocytes, cultured monocyte-macrophage cells, etc. We suggest a hypothetical model providing the potential mechanistic explanation for immune and inflammatory responses observed upon exposure to carbon nanoparticles. Specifically, we performed a theoretical study to analyze CNT and C60 fullerene interactions with the available X-ray structures of Toll-like receptors (TLRs) homo- and hetero-dimer extracellular domains. This assumption was based on the fact that similar to the known TLR ligands both CNTs and fullerenes induce, in cells, the secretion of certain inflammatory protein mediators, such as interleukins and chemokines. These proteins are observed within inflammation downstream processes resulted from the ligand molecule dependent inhibition or activation of TLR-induced signal transduction. Our computational studies have shown that the internal hydrophobic pockets of some TLRs might be capable of binding small-sized carbon nanostructures (5,5 armchair SWCNTs containing 11 carbon atom layers and C60 fullerene). High binding scores and minor structural alterations induced in TLR ectodomains upon binding C60 and CNTs further supported our hypothesis. Additionally, the proposed hypothesis is strengthened by the indirect experimental findings indicating that CNTs and fullerenes induce an excessive expression of specific cytokines and chemokines (i.e. IL-8 and MCP1).

Evaluation criteria for the quality of published experimental data on nanomaterials and their usefulness for QSAR modelling.

Nowadays nanotechnology is one of the most promising areas of science. The number and quantity of synthesized nanomaterials increase exponentially, therefore it is reasonable to expect that comprehensive risk assessment based only on empirical testing of all novel engineered nanoparticles (NPs) will very soon become impossible. Hence, the development of computational methods complementary to experimentation is very important. Quantitative structure-property relationship (QSPR) and quantitative structure-activity relationship (QSAR) models widely used in pharmaceutical chemistry and environmental science can also be modified and adopted for nanotechnology to predict physico-chemical properties and toxicity of empirically untested nanomaterials. All QSPR/QSAR modelling activities are based on experimentally derived data. It is important that, within a given data set, all values should be consistent, of high quality and measured according to a standardized protocol. Unfortunately, the amount of such data available for engineered nanoparticles in various data sources (i.e. databases and the literature) is very limited and seldom measured with a standardized protocol. Therefore, we have proposed a framework for collecting and evaluating the existing data, with the focus on possible applications for computational evaluation of properties and biological activities of nanomaterials.

CORAL: Monte Carlo Method as a Tool for the Prediction of the Bioconcentration Factor of Industrial Pollutants.

The CORAL software (http://www.insilico.eu/coral/) has been evaluated for application in QSAR modeling of the bioconcentration factor in fish (logBCF). The data used include 237 organic substances (industrial pollutants). Six random splits of the data into sub-training (30-50 %), calibration (20-30 %), test (13-30 %), and validation sets (7-25 %) have been carried out. The following numbers display the average statistical characteristics of the models for the external validation set: correlation coefficient r(2) =0.880±0.017 and standard error of estimation s=0.559±0.131. The best models were obtained with a combined representation of the molecular structure by SMILES together with hydrogen suppressed graph.

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.

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.

CORAL: classification model for predictions of anti-sarcoma activity.

A modified version of the CORAL software (http://www.insilico.eu/coral) allows building up the classification model for the case of the Yes/No data on the anti-sarcoma activity of organic compounds. Three random splits into the sub-training, calibration, and test sets of the data for 3017 compounds were examined. The performance of the proposed approach is satisfactory. The average values of the statistical characteristics for external test set on three random splits are as follows: n=1173-1234, sensitivity = 0.8903±0.0390, specificity = 0.9869±0.0013, and accuracy = 0.9759±0.0043. Mechanistic interpretation of the suggested model is discussed.

QSAR analysis of the toxicity of nitroaromatics in Tetrahymena pyriformis: structural factors and possible modes of action.

The Hierarchical Technology for Quantitative Structure-Activity Relationships (HiT QSAR) was applied to 95 diverse nitroaromatic compounds (including some widely known explosives) tested for their toxicity (50% inhibition growth concentration, IGC50) against the ciliate Tetrahymena pyriformis. The dataset was divided into subsets according to putative mechanisms of toxicity. The Classification and Regression Trees (CART) approach implemented within HiT QSAR has been used for prediction of mechanism of toxicity for new compounds. The resulting models were shown to have ~80% accuracy for external datasets indicating that the mechanistic dataset division was sensible. The Partial Least Squares (PLS) statistical approach was then used to develop 2D QSAR models. Validated PLS models were explored to: (1) elucidate the effects of different substituents in nitroaromatic compounds on toxicity; (2) differentiate compounds by probable mechanisms of toxicity based on their structural descriptors; and (3) analyse the role of various physical-chemical factors responsible for compounds' toxicity. Models were interpreted in terms of molecular fragments promoting or interfering with toxicity. It was also shown that mutual influence of substituents in benzene ring plays the determining role in toxicity variation. Although chemical mechanism based models were statistically significant and externally predictive (r²(ext) = 0.64 for the external set of 63 nitroaromatics identified after all calculations have been completed), they were also shown to have limited coverage (57% for modelling and 76% for external set).

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.

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.

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.

Aconitum and Delphinium alkaloids of curare-like activity. QSAR analysis and molecular docking of alkaloids into AChBP.

Early studies have shown that some of diterpenoid alkaloids, found in highly toxic plants of the genera Aconitum and Delphinium, act at neuronal nicotinic acetylcholine receptors (nAChRs) and exhibit potent N-cholinolytic activity. In the current study, GA-MLRA and GA-PLS approaches have been used to build QSAR models to predict N-cholinolytic activity measured in vivo (blockade of neuromuscular conductivity, BNMC and third eyelid relaxing activity, TYRA) and in vitro (suppression of frog's abdominal straight muscles on acetylcholine, SAM) for a series of diterpenoid alkaloids. Random splitting of a data set (five trials in total) produced QSAR models of a good level of correlation between experimental in vitro/in vivo and calculated N-cholinolytic activity expressed as log(1/ED(50)) with following average statistical parameters: log BNMC (r(2) = 0.87, s = 0.14, q(2) = 0.82), log TYRA (r(2) = 0.80, s = 0.29, q(2) = 0.67), log SAM (r(2) = 0.84, s = 29, q(2) = 0.64). QSAR results suggest descriptors accounting for H-bond capability of molecules influence all three type of N-cholinolytic activity with additional contribution of steric and reactivity features as identified for TYRA and SAM data, respectively. The alkaloid-receptor complexes were further analyzed by means of AutoDock Vina docking program using the binding site of MLA complexed with AChBP (homolog of the ligand binding domain of nAChRs) as template. All compounds were shown to be well fitted in the binding pocket of native MLA with good correlation exhibited between their ED(50) and AutoDock Vina binding free energy. An analysis of the possible factors significant for the ligand recognition has been enhanced by comparative docking studies performed for structurally related lycoctonine-type alkaloids (lappaconitine and aconitine) that are known to bind to voltage-gated Na(+) channel, but not to nAChRs.

Structural variability and the nature of intermolecular interactions in Watson-Crick B-DNA base pairs.

A set of nearly 100 crystallographic structures was analyzed using ab initio methods in order to verify the effect of the conformational variability of Watson-Crick guanine-cytosine and adenine-thymine base pairs on the intermolecular interaction energy and its components. Furthermore, for the representative structures, a potential energy scan of the structural parameters describing mutual orientation of the base pairs was carried out. The results were obtained using the hybrid variational-perturbational interaction energy decomposition scheme. The electron correlation effects were estimated by means of the second-order Møller-Plesset perturbation theory and coupled clusters with singles and doubles method adopting AUG-cc-pVDZ basis set. Moreover, the characteristics of hydrogen bonds in complexes, mimicking those appearing in B-DNA, were evaluated using topological analysis of the electron density. Although the first-order electrostatic energy is usually the largest stabilizing component, it is canceled out by the associated exchange repulsion in majority of the studied crystallographic structures. Therefore, the analyzed complexes of the nucleic acid bases appeared to be stabilized mainly by the delocalization component of the intermolecular interaction energy which, in terms of symmetry adapted perturbation theory, encompasses the second- and higher-order induction and exchange-induction terms. Furthermore, it was found that the dispersion contribution, albeit much smaller in terms of magnitude, is also a vital stabilizing factor. It was also revealed that the intermolecular interaction energy and its components are strongly influenced by four (out of six) structural parameters describing mutual orientation of bases in Watson-Crick pairs, namely shear, stagger, stretch, and opening. Finally, as a part of a model study, much of the effort was devoted to an extensive testing of the UBDB databank. It was shown that the databank quite successfully reproduces the electrostatic energy determined with the aid of ab initio methods.

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.

Structure and reactivity of TNT and related species: application of spectroscopic approaches and quantum-chemical approximations toward understanding transformation mechanisms.

This paper presents our latest findings regarding the structure and reactivity of the nitroaromatics, TNT and selected derivatives, within their environmental context. We also demonstrate the useful and proactive role of combined computational chemistry and spectroscopy tools in studying competing transformation mechanisms, particularly those with toxic potential. TNT and selected derivatives were reacted via alkaline hydrolysis as well as via free radical initiators through monochromatic irradiation and through Fenton reactions in complex competing transformation mechanisms. Only alkaline hydrolysis produced consistent and effective transformation intermediate and final products in this research. However, irradiation of the product generated by alkaline hydrolysis at 450 nm (wavelength of maximum absorption) caused complete disappearance of the spectra.