Toward on-the-fly quantum mechanical/molecular mechanical (QM/MM) docking: development and benchmark of a scoring function.

We address the challenges of treating polarization and covalent interactions in docking by developing a hybrid quantum mechanical/molecular mechanical (QM/MM) scoring function based on the semiempirical self-consistent charge density functional tight-binding (SCC-DFTB) method and the CHARMM force field. To benchmark this scoring function within the EADock DSS docking algorithm, we created a publicly available dataset of high-quality X-ray structures of zinc metalloproteins ( http://www.molecular-modelling.ch/resources.php ). For zinc-bound ligands (226 complexes), the QM/MM scoring yielded a substantially improved success rate compared to the classical scoring function (77.0% vs 61.5%), while, for allosteric ligands (55 complexes), the success rate remained constant (49.1%). The QM/MM scoring significantly improved the detection of correct zinc-binding geometries and improved the docking success rate by more than 20% for several important drug targets. The performance of both the classical and the QM/MM scoring functions compare favorably to the performance of AutoDock4, AutoDock4Zn, and AutoDock Vina.

On some numerical aspects of an active strain model in cardiac mechanics

We are interested in the development, implementation and testing of an orthotropic model for cardiac contraction based on an active strain decomposition. Our model addresses the coupling of a transversely isotropic mechanical description at the cell level, with an orthotropic constitutive law for incompressible tissue at the macroscopic level. The main differences with the active stress model are addressed in detail, and a finite element discretization using Taylor-Hood and MINI elements is proposed and illustrated with numerical examples.

CNDOL: A fast and reliable method for the calculation of electronic properties of very large systems. Applications to retinal binding pocket in rhodopsin and gas phase porphine.

Very large molecular systems can be calculated with the so called CNDOL approximate Hamiltonians that have been developed by avoiding oversimplifications and only using a priori parameters and formulas from the simpler NDO methods. A new diagonal monoelectronic term named CNDOL/21 shows great consistency and easier SCF convergence when used together with an appropriate function for charge repulsion energies that is derived from traditional formulas. It is possible to obtain a priori molecular orbitals and electron excitation properties after the configuration interaction of single excited determinants with reliability, maintaining interpretative possibilities even being a simplified Hamiltonian. Tests with some unequivocal gas phase maxima of simple molecules (benzene, furfural, acetaldehyde, hexyl alcohol, methyl amine, 2,5 dimethyl 2,4 hexadiene, and ethyl sulfide) ratify the general quality of this approach in comparison with other methods. The calculation of large systems as porphine in gas phase and a model of the complete retinal binding pocket in rhodopsin with 622 basis functions on 280 atoms at the quantum mechanical level show reliability leading to a resulting first allowed transition in 483 nm, very similar to the known experimental value of 500 nm of "dark state." In this very important case, our model gives a central role in this excitation to a charge transfer from the neighboring Glu(-) counterion to the retinaldehyde polyene chain. Tests with gas phase maxima of some important molecules corroborate the reliability of CNDOL/2 Hamiltonians.

Assessment and analysis of mechanical allodynia-like behavior induced by spared nerve injury (SNI) in the mouse

Rrésumé: La première description dans une publication médicale des douleurs neuropathiques remonte à 1872, le Dr S.W. Mitchell les résumant ainsi [...]" la causalgie est la plus terrible des tortures qu'une lésion nerveuse puisse entraîner "[...]. Par définition, la douleur neuropathique est une douleur chronique faisant suite à une lésion ou dysfonction du système nerveux. Malgré les progrès faits dans la compréhension de ce syndrome, le détail des mécanismes impliqués nous échappe encore et son traitement reste insuffisant car moins de 50% des patients sont soulagés par les thérapies actuelles. Différents modèles expérimentaux ont été élaborés chez l'animal de laboratoire, en particulier des modèles de lésion de nerfs périphériques chez le rat, permettant des investigations tant moléculaires que fonctionnelles des mécanismes impliqués dans le développement de ces douleurs. En revanche, peu de modèles existent chez la souris, alors que cet animal, grâce à la transgénèse, est très fréquemment utilisé pour l'approche fonctionnelle ciblée sur un gène. Dans l'étude présentée ici, nous avons évalué chez la souris C57BL/6 l'adaptation d'un modèle neuropathique, proposé une nouvelle modalité de mesure de la sensibilité douloureuse adaptée à la souris et défini une méthode d'analyse performante des résultats. Ce modèle, dit de lésion avec épargne nerveuse (spared Werve injury, SNI), consiste en la lésion de deux des trois branches du nerf sciatique, soit les nerfs peronier commun et tibial. La troisième branche, le nerf sural est laissé intact et c'est dans le territoire cutané de ce dernier que la sensibilité douloureuse à des stimulations mécaniques est enregistrée. Des filaments calibrés de force croissante sont appliqués sur la surface de la patte impliquée et la fréquence relative de retrait de la patte a été modélisée mathématiquement et analysée par un modèle statistique intégrant tous les paramètres de l'expérience (mixed-effects model). Des variantes chirurgicales lésant séquentiellement les trois branches du nerf sciatique ainsi que la réponse en fonction du sexe de l'animal ont également été évaluées. La lésion SNI entraîne une hypersensibilité mécanique marquée comparativement aux souris avec chirurgie contrôle; cet effet est constant entre les animaux et persiste durant les quatre semaines de l'étude. De subtiles différences entre les variables, y compris une divergence de sensibilité mécanique entre les sexes, ont été démontrées. La nécessité de léser le nerf tibial pour le développement des symptômes a également été documentée par notre méthode d'évaluation et d'analyse. En conclusion, nous avons validé le modèle SNI chez la souris par l'apparition d'un symptôme reproductible et apparenté à l'allodynie mécanique décrite par les patients souffrant de douleurs neuropathiques. Nous avons développé des méthodes d'enregistrement et d'analyse de la sensibilité douloureuse sensibles qui permettent la mise en évidence de facteurs intrinsèques et extrinsèques de variation de la réponse. Le modèle SNI utilisé chez des souris génétiquement modifiées, de par sa précision et reproductibilité, pourra permettre la discrimination de facteurs génétiques et épigénétiques contribuant au développement et à la persistance de douleurs neuropathiques.