Computational Analysis of CRTh2 receptor antagonist: A Ligand-based CoMFA and CoMSIA approach

CRTh2 receptor is an important mediator of inflammatory effects and has attracted much attention as a therapeutic target for the treatment of conditions such as asthma, COPD, allergic rhinitis and atopic dermatitis. In pursuit of better CRTh2 receptor antagonist agents, 3D-QSAR studies were performed on a series of 2-(2-(benzylthio)-1H-benzo[d]imidazol-1-yl) acetic acids. There is no crystal structure information available on this protein; hence in this work, ligand-based comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were performed by atom by atom matching alignment using systematic search and simulated annealing methods. The 3D-QSAR models were generated with 10 different combinations of test and training set molecules, since the robustness and predictive ability of the model is very important. We have generated 20 models for CoMFA and 100 models for CoMSIA based on two different alignments. Each model was validated with statistical cut off values such as q² > 0.4, r² > 0.5 and r²_pred > 0.5. Based on better q² and r²_pred values, the best predictions were obtained for the CoMFA (model 5 q² = 0.488, r²_pred = 0.732), and CoMSIA (model 45 q² = 0.525, r²_pred = 0.883) from systematic search conformation alignment. The high correlation between the cross-validated/predicted and experimental activities of a test set revealed that the CoMFA and CoMSIA models were robust. Statistical parameters from the generated QSAR models indicated the data is well fitted and have high predictive ability. The generated models suggest that steric, electrostatic, hydrophobic, H-bond donor and acceptor parameters are important for activity. Our study serves as a guide for further experimental investigations on the synthesis of new CRTh2 antagonist.     

Ligand-based virtual screening and in silico design of new antimalarial compounds using nonstochastic and stochastic total and atom-type quadratic maps

Malaria has been one of the most significant public health problems for centuries. It affects many tropical and subtropical regions of the world. The increasing resistance of Plasmodium spp. to existing therapies has heightened alarms about malaria in the international health community. Nowadays, there is a pressing need for identifying and developing new drug-based antimalarial therapies. In an effort to overcome this problem, the main purpose of this study is to develop simple linear discriminant-based quantitative structure-activity relationship (QSAR) models for the classification and prediction of antimalarial activity using some of the TOMOCOMD-CARDD (TOpological MOlecular COMputer Design-Computer Aided "Rational" Drug Design) fingerprints, so as to enable computational screening from virtual combinatorial datasets. In this sense, a database of 1562 organic chemicals having great structural variability, 597 of them antimalarial agents and 965 compounds having other clinical uses, was analyzed and presented as a helpful tool, not only for theoretical chemists but also for other researchers in this area. This series of compounds was processed by a k-means cluster analysis in order to design training and predicting sets. Afterward, two linear classification functions were derived in order to discriminate between antimalarial and nonantimalarial compounds. The models (including nonstochastic and stochastic indices) correctly classify more than 93% of the compound set, in both training and external prediction datasets. They showed high Matthews' correlation coefficients, 0.889 and 0.866 for the training set and 0.855 and 0.857 for the test one. The models' predictivity was also assessed and validated by the random removal of 10% of the compounds to form a new test set, for which predictions were made using the models. The overall means of the correct classification for this process (leave group 10% full-out cross validation) using the equations with nonstochastic and stochastic atom-based quadratic fingerprints were 93.93% and 92.77%, respectively. The quadratic maps-based TOMOCOMD-CARDD approach implemented in this work was successfully compared with four of the most useful models for antimalarials selection reported to date. The developed models were then used in a simulation of a virtual search for Ras FTase (FTase = farnesyltransferase) inhibitors with antimalarial activity; 70% and 100% of the 10 inhibitors used in this virtual search were correctly classified, showing the ability of the models to identify new lead antimalarials. Finally, these two QSAR models were used in the identification of previously unknown antimalarials. In this sense, three synthetic intermediaries of quinolinic compounds were evaluated as active/inactive ones using the developed models. The synthesis and biological evaluation of these chemicals against two malaria strains, using chloroquine as a reference, was performed. An accuracy of 100% with the theoretical predictions was observed. Compound 3 showed antimalarial activity, being the first report of an arylaminomethylenemalonate having such behavior. This result opens a door to a virtual study considering a higher variability of the structural core already evaluated, as well as of other chemicals not included in this study. We conclude that the approach described here seems to be a promising QSAR tool for the molecular discovery of novel classes of antimalarial drugs, which may meet the dual challenges posed by drug-resistant parasites and the rapid progression of malaria illnesses.     

An in silico approach to elucidate structure based functional evolution of oxacillinase

Bacterial Oxacillinases (OXAs), genetically being extremely diverse and highly versatile in hydrolyzing antibiotics of different classes, holds utmost significant clinical importance. Hence, to analyze functional evolution of this enzyme, plausible changes in drug profile, affinity and binding stability of different subclasses of OXA with their preferred drugs, viz. penicillin, ceftazidime, imipenem/meropenem were investigated. Maximum-Likelihood dendrogram was constructed and based on tree topology, the least and most divergent variants of each clade were selected. Pocket characterization, enzyme structural stability and mutational effect were analyzed in silico. Modes of interaction of selected OXA variants with respective antibiotics were analyzed by Autodock4.0 and LIGPLOT. Comparative mobility profiling and subsequent ΔG° and Km calculations of representative OXA variants revealed that after RSBL evolution, perhaps, two competitive strategies evolved among the OXA variants. Either loops flanking helix5 gets stabilized or it becomes more flexible. Therefore, while OXA variants (e.g. OXA-2, OXA-32, OXA-23, OXA-133, OXA-24, OXA-25, OXA-51 and OXA-75) with highly stabilized loops flanking helix5 exhibited improved binding stability and affinity towards carbapenems, especially meropenem, OXA variants (e.g. OXA-10, OXA-251, OXA-48 and OXA-247) possessing highly flexibile loops flanking helix5 revealed their catalytic proficiency towards ceftazidime. Moreover, LIGPLOT and PROMALS3D jointly identified ten consensuses/conserved residues, viz. P68, A69, F72, K73, W105, V120, W164, L169, K216 and G218 to be critical for drug hydrolysis. Hence, novel inhibitors could be designed to target these sites.     

A statistical approach to delta layer depth profiling

We present a simple statistical model describing the removal and relocation of material during a sputter depth profiling experiment. All input parameters are determined from low‐fluence molecular dynamics simulations, making the model de facto parameter free. The model can be used to extrapolate data from the molecular dynamics simulations to projectile fluences relevant to sputter depth profiling experiments. As a result, the erosion of the surface is calculated in terms of fluence‐dependent filling factors of different sample layers. Using input data determined for the 20‐keV C₆₀ cluster bombardment of silicon, it is found that a steady‐state erosion profile is reached after removal of approximately 20 monolayer equivalents of material. Plotting the contribution of particles from a specific layer to the instantaneous sputtered flux, one can directly determine the delta layer response function predicted from such a model. It is shown that this function can be parameterized by the semiempirical Dowsett response function, and the resulting fitting parameters are compared with published depth profile data. The model is then used to study the role of different processes influencing the observed depth resolution. We find that the statistical nature of the sputtering process suffices to explain many features of experimentally measured delta layer depth profiles. Copyright © 2012 John Wiley & Sons, Ltd.     

A global metabolite profiling approach to identify protein-metabolite interactions

Understanding the biochemical functions of proteins is an important factor in elucidating their cellular and physiological functions. Due to the predominance of biopolymer interactions in biology, many methods have been designed to interrogate and identify biologically relevant interactions that proteins make to DNA, RNA, and other proteins. Complementary approaches that can elucidate binding interactions between proteins and small molecule metabolites will impact the understanding of protein-metabolite interactions and fill a need that is outside the scope of current methods. Here, we demonstrate the ability to identify natural protein-metabolite interactions from complex metabolite mixtures by combining a protein-mediated small molecule enrichment step with a global metabolite profiling platform.     

A new approach to express ToF SIMS depth profiling

We propose a new approach to express SIMS depth profiling on a TOF.SIMS‐5 time‐of‐flight mass spectrometer. The approach is based on the instrument capability to independently perform raster scans of sputter and probe ion beams. The probed area can be much smaller than the diameter of a sputter ion beam, like in the AES depth profiling method. This circumstance alleviates limitations on the sputter beam–raster size relation, which are critical in other types of SIMS, and enables analysis on a curved‐bottomed sputter crater. By considerably reducing the raster size, it is possible to increase the depth profiling speed by an order of magnitude without radically degrading the depth resolution. A technique is proposed for successive improvement of depth resolution through profile recovery with account for the developing curvature of the sputtered crater bottom in the probed area. Experimental study of the crater bottom form resulted in implementing a method to include contribution of the instrumental artifacts in a nonstationary depth resolution function within the Hofmann's mixing–roughness–information depth model. The real‐structure experiment has shown that the analysis technique combining reduction of a raster size with a successive nonstationary recovery ensures high speed of profiling at ~100 µm/h while maintaining the depth resolution of about 30 nm at a 5 µm depth. Copyright © 2015 John Wiley & Sons, Ltd.     

An ensemble approach for in silico prediction of Ames mutagenicity

In this paper, we evaluate three learning algorithms based on supervised projections for molecular activity prediction. Using an approach based on supervised projections of the input space to construct ensembles of classifiers, three algorithms were tested. We constructed the projections by considering only instances that were misclassified by a previous classifier using the hidden layer of an Artificial Neural Network. We applied a supervised linear projection of the input space using a Nonparametric Discriminant Analysis method. Finally, we projected onto a subspace that minimizes the weighted error for each step. Using these three methods to construct ensembles of classifiers for the in silico prediction of Ames mutagenicity, we demonstrated the improved behavior of our proposal compared to classical methods.     

In silico and in vitro assessment of androgen receptor antagonists

There is a growing concern for male reproductive health as studies suggest that there is a sharp increase in prostate cancer and other fertility related problems. Apart from lifestyle, pollutants are also known to negatively affect the reproductive system. In addition to many other compounds that have been shown to alter androgen signaling, several environmental pollutants are known to disrupt androgen signaling via binding to androgen receptor (AR) or indirectly affecting the androgen synthesis. We analyzed here the molecular mechanism of the interaction between the human AR Ligand Binding Domain (hAR-LBD) and two environmental pollutants, linuron (a herbicide) and procymidone (a pesticide), and compared with the steroid agonist dihydrotestosterone (DHT) and well-known hAR antagonists bicalutamide and enzalutamide. Using molecular docking and dynamics simulations, we showed that the co-activator interaction site of the hAR-LBD is disrupted in different ways by different ligands. Binding free energies of the ligands were also ordered in increasing order as follows: linuron, procymidone, DHT, bicalutamide, and enzalutamide. These data were confirmed by in vitro assays. Reporter assay with MDA-kb2 cells showed that linuron, procymidone, bicalutamide and enzalutamide can inhibit androgen mediated activation of luciferase activity. Gene expression analysis further showed that these compounds can inhibit the expression of prostate specific antigen (PSA) and microseminoprotein beta (MSMB) in prostate cell line LNCaP. Comparative analysis showed that procymidone is more potent than linuron in inhibiting AR activity. Furthermore, procymidone at 10 μM dose showed equivalent and higher activity to AR inhibitor enzalutamide and bicalutamide respectively.     

The group theoretical approach to nuclear quadrupole resonance

The phenomenon of nuclear quadrupole resonance is presented from the grouptheoretical point of view. This approach allows one to discuss the splitting of the nuclear spin states without prior wave-mechanical calculation. Nuclear quadrupole moments are also compared with moecular quadrupole moments and the differences discussed. Finally the selection rules are deduced.     

An integrative in silico approach to the structure of Omp33-36 in Acinetobacter baumannii

Omp33-36 in A. baumannii, a bacterium causing serious nosocomial infections, is a virulence factor associated with the pathogen metabolic fitness as well as its adherence and invasion to human epithelial cells. This protein is also involved in interaction of the bacteria with host cells by binding to fibronectin. Moreover, Omp33-36 renders cytotoxicity to A. baumannii in addition to inducing apoptosis and modulation of autophagy. In the present study, an integrated strategy is launched to pierce into the 3D structure of Omp33-36 protein. The signal peptide within the sequence was determined, then, topology as well as secondary and tertiary structures of the protein were predicted. The mature protein assigned as a 14-stranded barrel in which residues 1–19 is removed as signal peptide. The obtained 3D models were evaluated in terms of quality; and then, served as queries to find similar protein structures. The hits were analyzed regarding topology among which 14-stranded were considered. The most qualified model was refined and then its sequence aligned to its counterpart similar structure protein (CymA from Klebsiella oxytoca). The determined structure of Omp33-36 could justify its porin function and carbapenem-resistance associated with the loss of this protein.     

Automated carbohydrate profiling by capillary electrophoresis: a bioindustrial approach

Automated, high-resolution, quantitative, high-throughput analysis of mono- and oligosaccharides, produced by enzymatic digestion of cellohexaose (model substrate) and lignocellulosic biomass, is demonstrated using high-performance capillary electrophoresis in conjunction with a single-step fluorophore labeling strategy for sensitive laser-induced fluorescence detection. Unattended batch sample processing from 96-well plates enabled reliable industrial-scale carbohydrate analysis. Excellent resolution of mono- and oligosaccharides was achieved under suppressed electroosmotic flow conditions, using either covalently or dynamically coated fused-silica capillary columns. The proposed approach readily supports the demands of bioindustrial operation environments with respect to high-throughput carbohydrate profiling.     

Molecular and structural determinants of adamantyl susceptibility to HLA-DRs allelic variants: an in silico approach to understand the mechanism of MLEs

Class II major histocompatibility complex (MHC II) molecules as expressed by antigen-presenting cells are heterodimeric cell-surface glycoprotein receptors that are fundamental in initiating and propagating an immune response by presenting tumor-associated antigenic peptides to CD4⁺/T_H cells. The loading efficiency of such peptides can be improved by small organic compounds (MHC Loading Enhancers—MLEs), that convert the non-receptive peptide conformation of MHC II to a peptide-receptive conformation. In a reversible reaction, these compounds open up the binding site of MHC II molecules by specific interactions with a yet undefined pocket. Here, we performed molecular docking and molecular dynamics simulation studies of adamantyl compounds on the predicted cavity around the P1 pocket of 2 allelic variants of HLA-DRs. The purpose was to investigate the suitability of adamantyl compounds as MLEs at the dimorphic β86 position. Docking studies revealed that besides numerous molecular interactions formed by the adamantyl compounds, Asnβ82, Tyrβ83, and Thrβ90 are the crucial amino acid residues that are characterized as the “sensors” of peptide loading. Molecular dynamics simulation studies exposed the dynamical structural changes that HLA-DRs adopted as a response to binding of 3-(1-adamantyl)-5-hydrazidocarbonyl-1H-pyrazole (AdCaPy). The conformations of AdCaPy complexed with the Glyβ86 HLA-DR allelic variant are well correlated with the stabilized form of peptide-loaded HLA-DRs, further confirming the role of AdCaPy as a MLE. Hydrogen bonding interaction analysis clearly demonstrated that after making suitable contacts with AdCaPy, HLA-DR changes its local conformation. However, AdCaPy complexed with HLA-DR having Valβ86 at the dimorphic position did not accommodate AdCaPy as MLE due to steric hindrance caused by the valine.     

Ligand-based molecular modeling study on a chemically diverse series of cholecystokinin-B/gastrin receptor antagonists: generation of predictive model

Pharmacophore hypotheses were developed for six structurally diverse series of cholecystokinin-B/gastrin receptor (CCK-BR) antagonists. A training set consisting of 33 compounds was carefully selected. The activity spread of the training set molecules was from 0.1 to 2100 nM. The most predictive pharmacophore model (hypothesis 1), consisting of four features, namely, two hydrogen bond donors, one hydrophobic aliphatic, and one hydrophobic aromatic feature, had a correlation (r) of 0.884 and a root-mean-square deviation of 1.1526, and the cost difference between null cost and fixed cost was 81.5 bits. The model was validated on a test set consisting of six different series of 27 structurally diverse compounds and performed well in classifying active and inactive molecules correctly. This validation approach provides confidence in the utility of the predictive pharmacophore model developed in this work as a 3D query tool in the virtual screening of drug-like molecules to retrieve new chemical entities as potent CCK-BR antagonists. The model can also be used to predict the biological activities of compounds prior to their costly and time-consuming synthesis.     

Maltase-glucoamylase inhibition potency and cytotoxicity of pyrimidine-fused compounds: An in silico and in vitro approach

The prevalence of diabetes mellitus has been incremented in the current century and the need for novel therapeutic compounds to treat this disease has been significantly increased. One of the most promising approaches is to inhibit intestinal alpha glucosidases. Based on our previous studies, four pyrimidine-fused heterocycles (PFH) were selected as they revealed satisfactory inhibitory action against mammalian α-glucosidase. The interaction of these compounds with both active domains of human maltase-glucoamylase (MGAM) and their effect on human Caco-2 cell line were investigated. The docking assessments suggested that binding properties of these ligands were almost similar to that of acarbose by establishing hydrogen bonds especially with Tyr1251 and Arg526 in both C-terminal and N-terminal MGAM, respectively. Also, these compounds indicated a stronger affinity for C-terminal of MGAM. L2 and L4 made tightly complexes with both terminals of MGAM which in turn revealed the importance of introducing pyrimidine scaffold and its hinge compartment. The results of molecular dynamics simulation analyses confirmed the docking data and showed deep penetration of L2 and L4 into the active site of MGAM. Based on cell cytotoxicity assessments, no significant cell death induction was observed. Hence, these functional MGAM inhibitors might be considered as new potential therapeutic compounds in treatment of diabetes and its complications.     

Expanded metabolomics approach to profiling endogenous carbohydrates in the serum of ovarian cancer patients

We applied hydrophilic interaction liquid chromatography coupled with tandem mass spectrometry to the quantitative analysis of serum from 58 women, including ovarian cancer patients, ovarian benign tumor patients, and healthy controls. All of these ovarian cancer and ovarian benign tumor patients have elevated cancer antigen 125, which makes them clinically difficult to differentiate the malignant from the benign. All of the 16 endogenous carbohydrates were quantitatively detected in the human sera, of which, eight endogenous carbohydrates were significantly different (P‐value < 0.05) between the ovarian cancer and healthy control. According to the receiver operating characteristic curve analysis, arabitol was the most potentially specific biomarker for discriminating ovarian cancer from healthy control, having an area under the curve of 0.911. A panel of metabolite markers composed of maltose, maltotriose, raffinose, and mannitol was selected, which was able to discriminate the ovarian cancer from the benign ovarian tumor counterparts, with an area under concentration‐time curve value of 0.832. Endogenous carbohydrates in the expanded metabolomics approach after the global metabolic profiling are characterized and are potential biomarkers for the early diagnosis of ovarian cancer.     

Noise reduction method for molecular interaction energy: application to in silico drug screening and in silico target protein screening

We developed a new method to improve the accuracy of molecular interaction data using a molecular interaction matrix. This method was applied to enhance the database enrichment of in silico drug screening and in silico target protein screening using a protein-compound affinity matrix calculated by a protein-compound docking software. Our assumption was that the protein-compound binding free energy of a compound could be improved by a linear combination of its docking scores with many different proteins. We proposed two approaches to determine the coefficients of the linear combination. The first approach is based on similarity among the proteins, and the second is a machine-learning approach based on the known active compounds. These methods were applied to in silico screening of the active compounds of several target proteins and in silico target protein screening.     

Degradation kinetics of fluvoxamine in buffer solutions: In silico ADMET profiling and identification of degradation products by LC-MS/ESI

The kinetics of hydrolysis of fluvoxamine maleate (FLV) has been investigated over the pH range 1.0–12.0 at 40, 60 and 80 °C. FLV degradation follows pseudo-first-order kinetics which is consistent with the kinetics of drugs that are not readily dissolved in aqueous medium. The hydrolytic degradation rate constant (k_obs) range from 0.92 (pH 6.0) to 13.8 × 10⁻⁴ min⁻¹ (pH 1.0). The k_obs represents the sum of six different degradation rate constants; the k_H has been found to be higher than k_OH. The FLV exhibits a typical rate- pH profile with a flat bottom over the pH range 3.0–6.0 which indicates its maximum stability at pH 6.0. Ten FLV degradants have been predicted by Zeneth software and among them four degradation products (D1, D2, D3 and D4) have been identified in degraded samples. The in-silico pharmacokinetics and toxicity of degradation products have been determined using Swiss ADME and admetSAR software. The toxicity profile reveals that D2 is both AMES toxic and carcinogenic while the rest of the products are non-AMES toxic and non-carcinogenic. All of the degradation products are high in causing fish toxicity thus their presence in pharmaceutical waste is alarming for environmental safety.