Conditional mutual inclusive information enables accurate quantification of associations in gene regulatory networks

International audience

Automated morphometry provides accurate and reproducible virtual staging of liver fibrosis in chronic hepatitis C

International audience

A New Multi-Objective Approach for Molecular Docking Based on RMSD and Binding Energy

Ligand-protein docking is an optimization problem based on predicting the position of a ligand with the lowest binding energy in the active site of the receptor. Molecular docking problems are traditionally tackled with single-objective, as well as with multi-objective approaches, to minimize the binding energy. In this paper, we propose a novel multi-objective formulation that considers: the Root Mean Square Deviation (RMSD) difference in the coordinates of ligands and the binding (intermolecular) energy, as two objectives to evaluate the quality of the ligand-protein interactions. To determine the kind of Pareto front approximations that can be obtained, we have selected a set of representative multi-objective algorithms such as NSGA-II, SMPSO, GDE3, and MOEA/D. Their performances have been assessed by applying two main quality indicators intended to measure convergence and diversity of the fronts. In addition, a comparison with LGA, a reference single-objective evolutionary algorithm for molecular docking (AutoDock) is carried out. In general, SMPSO shows the best overall results in terms of energy and RMSD (value lower than 2A for successful docking results). This new multi-objective approach shows an improvement over the ligand-protein docking predictions that could be promising in in silico docking studies to select new anticancer compounds for therapeutic targets that are multidrug resistant.

Unraveling the performance of dispersion-corrected functionals for the accurate description of weakly bound natural polyphenols

Long-range non-covalent interactions play a key role in the chemistry of natural polyphenols. We have previously proposed a description of supramolecular polyphenol complexes by the B3P86 density functional coupled with some corrections for dispersion. We couple here the B3P86 functional with the D3 correction for dispersion, assessing systematically the accuracy of the new B3P86-D3 model using for that the well-known S66, HB23, NCCE31, and S12L datasets for non-covalent interactions. Furthermore, the association energies of these complexes were carefully compared to those obtained by other dispersion-corrected functionals, such as B(3)LYP-D3, BP86-D3 or B3P86-NL. Finally, this set of models were also applied to a database composed of seven non-covalent polyphenol complexes of the most interest.

Accurate Treatment of Large Supramolecular Complexes by Double-Hybrid Density Functionals Coupled with Nonlocal van der Waals Corrections

In this work, we present a thorough assessment of the performance of some representative double-hybrid density functionals (revPBE0-DH-NL and B2PLYP-NL) as well as their parent hybrid and GGA counterparts, in combination with the most modern version of the nonlocal (NL) van der Waals correction to describe very large weakly interacting molecular systems dominated by noncovalent interactions. Prior to the assessment, an accurate and homogeneous set of reference interaction energies was computed for the supramolecular complexes constituting the L7 and S12L data sets by using the novel, precise, and efficient DLPNO-CCSD(T) method at the complete basis set limit (CBS). The correction of the basis set superposition error and the inclusion of the deformation energies (for the S12L set) have been crucial for obtaining precise DLPNO-CCSD(T)/CBS interaction energies. Among the density functionals evaluated, the double-hybrid revPBE0-DH-NL and B2PLYP-NL with the three-body dispersion correction provide remarkably accurate association energies very close to the chemical accuracy. Overall, the NL van der Waals approach combined with proper density functionals can be seen as an accurate and affordable computational tool for the modeling of large weakly bonded supramolecular systems.

New cholinesterase inhibitors for Alzheimer's disease: Structure Activity Studies (SARs) and molecular docking of isoquinolone and azepanone derivatives

A library of isoquinolinone and azepanone derivatives were screened for both acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE) activity. The strategy adopted included (a) in vitro biological assays, against eel AChE (EeAChE) and equine serum BuChE (EqBuChE) in order to determine the compounds IC50 and their dose-response activity, consolidated by (b) molecular docking studies to evaluate the docking poses and interatomic interactions in the case of the hit compounds, validated by STD-NMR studies. Compound (1f) was identified as one of these hits with an IC50 of 89.5 mu M for EeAChE and 153.8 mu M for EqBuChE, (2a) was identified as a second hit with an IC50 of 108.4 mu M (EeAChE) and 277.8 mu M (EqBuChE). In order to gain insights into the binding mode and principle active site interactions of these molecules, (R)-(1f) along with 3 other analogues (also as the R-enantiomer) were docked into both RhAChE and hBuChE models. Galantamine was used as the benchmark. The docking study was validated by performing an STD-NMR study of (1f) with EeAChE using galantamine as the benchmark.

Insights into (S)-rivastigmine inhibition of butyrylcholinesterase (BuChE): Molecular docking and saturation transfer difference NMR (STD-NMR)

Rivastigmine is a very important drug prescribed for the treatment of Alzheimer's disease (AD) symptoms. It is a dual inhibitor, in that it inhibits both acetylcholinesterase (AChE) and butyrylcholinesterase (BuChE). For our screening program on the discovery of new rivastigmine analogue hits for human butyrylcholinesterase (hBuChE) inhibition, we investigated the interaction of this inhibitor with BuChE using the complimentary approach of the biophysical method, saturation transfer difference (STD)-NMR and molecular docking. This allowed us to obtain essential information on the key binding interactions between the inhibitor and the enzyme to be used for screening of hit compounds. The main conclusions obtained from this integrated study was that the most dominant interactions were (a) H-bonding between the carbamate carbonyl of the inhibitor and the NH group of the imidazole unit of H434, (b) stacking of the aromatic unit of the inhibitor and the W82 aromatic unit in the choline binding pocket via pi-pi interactions and (c) possible CH/pi interactions between the benzylic methyl group and the N-methyl groups of the inhibitor and W82 of the enzyme.

Simulation of photoinduced processes in organic chromophores: mapping photochemically relevant states for accurate dynamics

The simulation of ultrafast photoinduced processes is a fundamental step towards the understanding of the underlying molecular mechanism and interpretation/prediction of experimental data. Performing a computer simulation of a complex photoinduced process is only possible introducing some approximations but, in order to obtain reliable results, the need to reduce the complexity must balance with the accuracy of the model, which should include all the relevant degrees of freedom and a quantitatively correct description of the electronic states involved in the process. This work presents new computational protocols and strategies for the parameterisation of accurate models for photochemical/photophysical processes based on state-of-the-art multiconfigurational wavefunction-based methods. The required ingredients for a dynamics simulation include potential energy surfaces (PESs) as well as electronic state couplings, which must be mapped across the wide range of geometries visited during the wavepacket/trajectory propagation. The developed procedures allow to obtain solid and extended databases reducing as much as possible the computational cost, thanks to, e.g., specific tuning of the level of theory for different PES regions and/or direct calculation of only the needed components of vectorial quantities (like gradients or nonadiabatic couplings). The presented approaches were applied to three case studies (azobenzene, pyrene, visual rhodopsin), all requiring an accurate parameterisation but for different reasons. The resulting models and simulations allowed to elucidate the mechanism and time scale of the internal conversion, reproducing or even predicting new transient experiments. The general applicability of the developed protocols to systems with different peculiarities and the possibility to parameterise different types of dynamics on an equal footing (classical vs purely quantum) prove that the developed procedures are flexible enough to be tailored for each specific system, and pave the way for exact quantum dynamics with multiple degrees of freedom.

Advancing Bag-of-visual-words Representations For Lesion Classification In Retinal Images.

Diabetic Retinopathy (DR) is a complication of diabetes that can lead to blindness if not readily discovered. Automated screening algorithms have the potential to improve identification of patients who need further medical attention. However, the identification of lesions must be accurate to be useful for clinical application. The bag-of-visual-words (BoVW) algorithm employs a maximum-margin classifier in a flexible framework that is able to detect the most common DR-related lesions such as microaneurysms, cotton-wool spots and hard exudates. BoVW allows to bypass the need for pre- and post-processing of the retinographic images, as well as the need of specific ad hoc techniques for identification of each type of lesion. An extensive evaluation of the BoVW model, using three large retinograph datasets (DR1, DR2 and Messidor) with different resolution and collected by different healthcare personnel, was performed. The results demonstrate that the BoVW classification approach can identify different lesions within an image without having to utilize different algorithms for each lesion reducing processing time and providing a more flexible diagnostic system. Our BoVW scheme is based on sparse low-level feature detection with a Speeded-Up Robust Features (SURF) local descriptor, and mid-level features based on semi-soft coding with max pooling. The best BoVW representation for retinal image classification was an area under the receiver operating characteristic curve (AUC-ROC) of 97.8% (exudates) and 93.5% (red lesions), applying a cross-dataset validation protocol. To assess the accuracy for detecting cases that require referral within one year, the sparse extraction technique associated with semi-soft coding and max pooling obtained an AUC of 94.2 ± 2.0%, outperforming current methods. Those results indicate that, for retinal image classification tasks in clinical practice, BoVW is equal and, in some instances, surpasses results obtained using dense detection (widely believed to be the best choice in many vision problems) for the low-level descriptors.

Primary Relapse Site Pattern In Women With Triple-negative Breast Cancer.

Despite the remarkable improvements in breast cancer (BC) characterization, accurate prediction of BC clinical behavior is often still difficult to achieve. Some studies have investigated the association between the molecular subtype, namely the basal-like BC and the pattern of relapse, however only few investigated the association between relapse pattern and immunohistochemical defined triple-negative breast cancers (TNBCs). The aim of this study was to evaluate the pattern of relapse in patients with TNBC, namely the primary distant relapse site. One-hundred twenty nine (129) invasive breast carcinomas with follow-up information were classified according to the molecular subtype using immunohistochemistry for ER, PgR and Her2. The association between TNBC and distant relapse primary site was analyzed by logistic regression. Using multivariate logistic regression analysis patients with TNBC displayed only 0.09 (95% CI: 0.00-0.74; p=0.02) the odds of the non-TNBC patients of developing bone primary relapse. Regarding visceral and lymph-node relapse, no differences between in this cohort were found. Though classically regarded as aggressive tumors, TNBCs rarely development primary relapse in bone when compared to non-TNBC, a clinical relevant fact when investigating a metastasis of an occult or non-sampled primary BC.

Synthesis And Antiproliferative Activity Of 8-hydroxyquinoline Derivatives Containing A 1,2,3-triazole Moiety.

Twelve novel 8-hydroxyquinoline derivatives were synthesized with good yields by performing copper-catalyzed Huisgen 1,3-dipolar cycloaddition (click reaction) between an 8-O-alkylated-quinoline containing a terminal alkyne and various aromatic or protected sugar azides. These compounds were evaluated in vitro for their antiproliferative activity on various cancer cell types. Protected sugar derivative 16 was the most active compound in the series, exhibiting potent antiproliferative activity and high selectivity toward ovarian cancer cells (OVCAR-03, GI50 < 0.25 μg mL(-1)); this derivative was more active than the reference drug doxorubicin (OVCAR-03, GI50 = 0.43 μg mL(-1)). In structure-activity relationship (SAR) studies, the physico-chemical parameters of the compounds were evaluated and docking calculations were performed for the α-glucosidase active site to predict the possible mechanism of action of this series of compounds.

Developmental Rate Of Immatures Of Two Fly Species Of Forensic Importance: Sarcophaga (liopygia) Ruficornis And Microcerella Halli (diptera: Sarcophagidae).

Since insect species are poikilothermic organisms, they generally exhibit different growth patterns depending on the temperature at which they develop. This factor is important in forensic entomology, especially for estimating postmortem interval (PMI) when it is based on the developmental time of the insects reared in decomposing bodies. This study aimed to estimate the rates of development, viability, and survival of immatures of Sarcophaga (Liopygia) ruficornis (Fabricius 1794) and Microcerella halli (Engel 1931) (Diptera: Sarcophagidae) reared in different temperatures: 10, 15, 20, 25, 30, and 35 ± 1 °C. Bovine raw ground meat was offered as food for all experimental groups, each consisting of four replicates, in the proportion of 2 g/larva. To measure the evolution of growth, ten specimens of each group were randomly chosen and weighed every 12 h, from initial feeding larva to pupae, and then discarded. Considering the records of weight gain, survival rates, and stability of growth rates, the range of optimum temperature for the development of S. (L.) ruficornis is between 20 and 35 °C, and that of M. halli is between 20 and 25 °C. For both species, the longest times of development were in the lowest temperatures. The survival rate at extreme temperatures (10 and 35 °C) was lower in both species. Biological data such as the ones obtained in this study are of great importance to achieve a more accurate estimate of the PMI.

Propylthiouracil Quantification In Human Plasma By High-performance Liquid Chromatography Coupled With Electrospray Tandem Mass Spectrometry: Application In A Bioequivalence Study.

A rapid, sensitive and specific method for quantifying propylthiouracil in human plasma using methylthiouracil as the internal standard (IS) is described. The analyte and the IS were extracted from plasma by liquid-liquid extraction using an organic solvent (ethyl acetate). The extracts were analyzed by high performance liquid chromatography coupled with electrospray tandem mass spectrometry (HPLC-MS/MS) in negative mode (ES-). Chromatography was performed using a Phenomenex Gemini C18 5μm analytical column (4.6mm×150mm i.d.) and a mobile phase consisting of methanol/water/acetonitrile (40/40/20, v/v/v)+0.1% of formic acid. For propylthiouracil and I.S., the optimized parameters of the declustering potential, collision energy and collision exit potential were -60 (V), -26 (eV) and -5 (V), respectively. The method had a chromatographic run time of 2.5min and a linear calibration curve over the range 20-5000ng/mL. The limit of quantification was 20ng/mL. The stability tests indicated no significant degradation. This HPLC-MS/MS procedure was used to assess the bioequivalence of two propylthiouracil 100mg tablet formulations in healthy volunteers of both sexes in fasted and fed state. The geometric mean and 90% confidence interval CI of Test/Reference percent ratios were, without and with food, respectively: 109.28% (103.63-115.25%) and 115.60% (109.03-122.58%) for Cmax, 103.31% (100.74-105.96%) and 103.40% (101.03-105.84) for AUClast. This method offers advantages over those previously reported, in terms of both a simple liquid-liquid extraction without clean-up procedures, as well as a faster run time (2.5min). The LOQ of 20ng/mL is well suited for pharmacokinetic studies. The assay performance results indicate that the method is precise and accurate enough for the routine determination of the propylthiouracil in human plasma. The test formulation with and without food was bioequivalent to reference formulation. Food administration increased the Tmax and decreased the bioavailability (Cmax and AUC).

Atomic Charge Transfer-counter Polarization Effects Determine Infrared Ch Intensities Of Hydrocarbons: A Quantum Theory Of Atoms In Molecules Model.

Atomic charge transfer-counter polarization effects determine most of the infrared fundamental CH intensities of simple hydrocarbons, methane, ethylene, ethane, propyne, cyclopropane and allene. The quantum theory of atoms in molecules/charge-charge flux-dipole flux model predicted the values of 30 CH intensities ranging from 0 to 123 km mol(-1) with a root mean square (rms) error of only 4.2 km mol(-1) without including a specific equilibrium atomic charge term. Sums of the contributions from terms involving charge flux and/or dipole flux averaged 20.3 km mol(-1), about ten times larger than the average charge contribution of 2.0 km mol(-1). The only notable exceptions are the CH stretching and bending intensities of acetylene and two of the propyne vibrations for hydrogens bound to sp hybridized carbon atoms. Calculations were carried out at four quantum levels, MP2/6-311++G(3d,3p), MP2/cc-pVTZ, QCISD/6-311++G(3d,3p) and QCISD/cc-pVTZ. The results calculated at the QCISD level are the most accurate among the four with root mean square errors of 4.7 and 5.0 km mol(-1) for the 6-311++G(3d,3p) and cc-pVTZ basis sets. These values are close to the estimated aggregate experimental error of the hydrocarbon intensities, 4.0 km mol(-1). The atomic charge transfer-counter polarization effect is much larger than the charge effect for the results of all four quantum levels. Charge transfer-counter polarization effects are expected to also be important in vibrations of more polar molecules for which equilibrium charge contributions can be large.

Selection And Validation Of Reference Genes For Functional Studies In The Calliphoridae Family.

The genera Cochliomyia and Chrysomya contain both obligate and saprophagous flies, which allows the comparison of different feeding habits between closely related species. Among the different strategies for comparing these habits is the use of qPCR to investigate the expression levels of candidate genes involved in feeding behavior. To ensure an accurate measure of the levels of gene expression, it is necessary to normalize the amount of the target gene with the amount of a reference gene having a stable expression across the compared species. Since there is no universal gene that can be used as a reference in functional studies, candidate genes for qPCR data normalization were selected and validated in three Calliphoridae (Diptera) species, Cochliomyia hominivorax Coquerel, Cochliomyia macellaria Fabricius, and Chrysomya albiceps Wiedemann . The expression stability of six genes ( Actin, Gapdh, Rp49, Rps17, α -tubulin, and GstD1) was evaluated among species within the same life stage and between life stages within each species. The expression levels of Actin, Gapdh, and Rp49 were the most stable among the selected genes. These genes can be used as reliable reference genes for functional studies in Calliphoridae using similar experimental settings.