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.

Endangerment and Likeability of Wildlife Species: How Important are they for Payments Proposed for Conservation?

This paper examines empirically the relative influence of the degree of endangerment of wildlife species and their stated likeability on individuals' allocation of funds for their conservation. To do this, it utilises data obtained from the IUCN Red List, and likeability and fund allocation data obtained from two serial surveys of a sample of the Australian public who were requested to assess 24 Australian wildlife species from three animal classes: mammals, birds and reptiles. Between the first and second survey, respondents were provided with extra information about the focal species. This information resulted in the dominance of endangerment as the major influence on the allocation of funding of respondents for the conservation of the focal wildlife species. Our results throw doubts on the proposition in the literature that the likeability of species is the dominant influence on willingness to pay for conservation of wildlife species. Furthermore, because the public's allocation of fund for conserving wildlife species seems to be more sensitive to information about the conservation status of species than to factors influencing their likeability, greater attention to providing accurate information about the former than the latter seems justified. Keywords: Conservation of wildlife species; Contingent valuation; Endangerment of species; Likeability of species; Willingness to pay

Effective Cell-Centred Time-Domain Maxwell's Equations Numerical Solvers

This research work analyses techniques for implementing a cell-centred finite-volume time-domain (ccFV-TD) computational methodology for the purpose of studying microwave heating. Various state-of-the-art spatial and temporal discretisation methods employed to solve Maxwell's equations on multidimensional structured grid networks are investigated, and the dispersive and dissipative errors inherent in those techniques examined. Both staggered and unstaggered grid approaches are considered. Upwind schemes using a Riemann solver and intensity vector splitting are studied and evaluated. Staggered and unstaggered Leapfrog and Runge-Kutta time integration methods are analysed in terms of phase and amplitude error to identify which method is the most accurate and efficient for simulating microwave heating processes. The implementation and migration of typical electromagnetic boundary conditions. from staggered in space to cell-centred approaches also is deliberated. In particular, an existing perfectly matched layer absorbing boundary methodology is adapted to formulate a new cell-centred boundary implementation for the ccFV-TD solvers. Finally for microwave heating purposes, a comparison of analytical and numerical results for standard case studies in rectangular waveguides allows the accuracy of the developed methods to be assessed.

A second order control-volume finite-element least-squares strategy for simulating diffusion in strongly anisotropic media

An unstructured mesh �nite volume discretisation method for simulating di�usion in anisotropic media in two-dimensional space is discussed. This technique is considered as an extension of the fully implicit hybrid control-volume �nite-element method and it retains the local continuity of the ux at the control volume faces. A least squares function recon- struction technique together with a new ux decomposition strategy is used to obtain an accurate ux approximation at the control volume face, ensuring that the overall accuracy of the spatial discretisation maintains second order. This paper highlights that the new technique coincides with the traditional shape function technique when the correction term is neglected and that it signi�cantly increases the accuracy of the previous linear scheme on coarse meshes when applied to media that exhibit very strong to extreme anisotropy ratios. It is concluded that the method can be used on both regular and irregular meshes, and appears independent of the mesh quality.

A Petrov-Galerkin method for a singularly perturbed ordinary differential equation with non-smooth data

In this paper, a singularly perturbed ordinary differential equation with non-smooth data is considered. The numerical method is generated by means of a Petrov-Galerkin finite element method with the piecewise-exponential test function and the piecewise-linear trial function. At the discontinuous point of the coefficient, a special technique is used. The method is shown to be first-order accurate and singular perturbation parameter uniform convergence. Finally, numerical results are presented, which are in agreement with theoretical results.

Numerical approximation of a fractional-in-space diffusion equation (II) - with nonhomogeneous boundary conditions

In this paper, a space fractional di®usion equation (SFDE) with non- homogeneous boundary conditions on a bounded domain is considered. A new matrix transfer technique (MTT) for solving the SFDE is proposed. The method is based on a matrix representation of the fractional-in-space operator and the novelty of this approach is that a standard discretisation of the operator leads to a system of linear ODEs with the matrix raised to the same fractional power. Analytic solutions of the SFDE are derived. Finally, some numerical results are given to demonstrate that the MTT is a computationally e±cient and accurate method for solving SFDE.

Using coloured petri nets to simulate DoS-resistant protocols

In this work, we examine unbalanced computation between an initiator and a responder that leads to resource exhaustion attacks in key exchange protocols. We construct models for two cryp-tographic protocols; one is the well-known Internet protocol named Secure Socket Layer (SSL) protocol, and the other one is the Host Identity Protocol (HIP) which has built-in DoS-resistant mechanisms. To examine such protocols, we develop a formal framework based on Timed Coloured Petri Nets (Timed CPNs) and use a simulation approach provided in CPN Tools to achieve a formal analysis. By adopting the key idea of Meadows' cost-based framework and re¯ning the de¯nition of operational costs during the protocol execution, our simulation provides an accurate cost estimate of protocol execution compar- ing among principals, as well as the percentage of successful connections from legitimate users, under four di®erent strategies of DoS attack.

Tracking people in 3D using position, size and shape

This paper presents a prototype tracking system for tracking people in enclosed indoor environments where there is a high rate of occlusions. The system uses a stereo camera for acquisition, and is capable of disambiguating occlusions using a combination of depth map analysis, a two step ellipse fitting people detection process, the use of motion models and Kalman filters and a novel fit metric, based on computationally simple object statistics. Testing shows that our fit metric outperforms commonly used position based metrics and histogram based metrics, resulting in more accurate tracking of people.