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.

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 of distribution and evolution of mechanical dyssynchrony in a porcine model of myocardial infarction by cardiovascular magnetic resonance.

BACKGROUND: We sought to investigate the relationship between infarct and dyssynchrony post- myocardial infarct (MI), in a porcine model. Mechanical dyssynchrony post-MI is associated with left ventricular (LV) remodeling and increased mortality. METHODS: Cine, gadolinium-contrast, and tagged cardiovascular magnetic resonance (CMR) were performed pre-MI, 9 ± 2 days (early post-MI), and 33 ± 10 days (late post-MI) post-MI in 6 pigs to characterize cardiac morphology, location and extent of MI, and regional mechanics. LV mechanics were assessed by circumferential strain (eC). Electro-anatomic mapping (EAM) was performed within 24 hrs of CMR and prior to sacrifice. RESULTS: Mean infarct size was 21 ± 4% of LV volume with evidence of post-MI remodeling. Global eC significantly decreased post MI (-27 ± 1.6% vs. -18 ± 2.5% (early) and -17 ± 2.7% (late), p < 0.0001) with no significant change in peri-MI and MI segments between early and late time-points. Time to peak strain (TTP) was significantly longer in MI, compared to normal and peri-MI segments, both early (440 ± 40 ms vs. 329 ± 40 ms and 332 ± 36 ms, respectively; p = 0.0002) and late post-MI (442 ± 63 ms vs. 321 ± 40 ms and 355 ± 61 ms, respectively; p = 0.012). The standard deviation of TTP in 16 segments (SD16) significantly increased post-MI: 28 ± 7 ms to 50 ± 10 ms (early, p = 0.012) to 54 ± 19 ms (late, p = 0.004), with no change between early and late post-MI time-points (p = 0.56). TTP was not related to reduction of segmental contractility. EAM revealed late electrical activation and greatly diminished conduction velocity in the infarct (5.7 ± 2.4 cm/s), when compared to peri-infarct (18.7 ± 10.3 cm/s) and remote myocardium (39 ± 20.5 cm/s). CONCLUSIONS: Mechanical dyssynchrony occurs early after MI and is the result of delayed electrical and mechanical activation in the infarct.

Rufinamide Attenuates Mechanical Allodynia in a Model of Neuropathic Pain in the Mouse and Stabilizes Voltage-gated Sodium Channel Inactivated State.

BACKGROUND:: Voltage-gated sodium channels dysregulation is important for hyperexcitability leading to pain persistence. Sodium channel blockers currently used to treat neuropathic pain are poorly tolerated. Getting new molecules to clinical use is laborious. We here propose a drug already marketed as anticonvulsant, rufinamide. METHODS:: We compared the behavioral effect of rufinamide to amitriptyline using the Spared Nerve Injury neuropathic pain model in mice. We compared the effect of rufinamide on sodium currents using in vitro patch clamp in cells expressing the voltage-gated sodium channel Nav1.7 isoform and on dissociated dorsal root ganglion neurons to amitriptyline and mexiletine. RESULTS:: In naive mice, amitriptyline (20 mg/kg) increased withdrawal threshold to mechanical stimulation from 1.3 (0.6-1.9) (median [95% CI]) to 2.3 g (2.2-2.5) and latency of withdrawal to heat stimulation from 13.1 (10.4-15.5) to 30.0 s (21.8-31.9), whereas rufinamide had no effect. Rufinamide and amitriptyline alleviated injury-induced mechanical allodynia for 4 h (maximal effect: 0.10 ± 0.03 g (mean ± SD) to 1.99 ± 0.26 g for rufinamide and 0.25 ± 0.22 g to 1.92 ± 0.85 g for amitriptyline). All drugs reduced peak current and stabilized the inactivated state of voltage-gated sodium channel Nav1.7, with similar effects in dorsal root ganglion neurons. CONCLUSIONS:: At doses alleviating neuropathic pain, amitriptyline showed alteration of behavioral response possibly related to either alteration of basal pain sensitivity or sedative effect or both. Side-effects and drug tolerance/compliance are major problems with drugs such as amitriptyline. Rufinamide seems to have a better tolerability profile and could be a new alternative to explore for the treatment of neuropathic pain.

A comprehensive model for quantum noise characterization in digital mammography.

NlmCategory="UNASSIGNED">A version of cascaded systems analysis was developed specifically with the aim of studying quantum noise propagation in x-ray detectors. Signal and quantum noise propagation was then modelled in four types of x-ray detectors used for digital mammography: four flat panel systems, one computed radiography and one slot-scan silicon wafer based photon counting device. As required inputs to the model, the two dimensional (2D) modulation transfer function (MTF), noise power spectra (NPS) and detective quantum efficiency (DQE) were measured for six mammography systems that utilized these different detectors. A new method to reconstruct anisotropic 2D presampling MTF matrices from 1D radial MTFs measured along different angular directions across the detector is described; an image of a sharp, circular disc was used for this purpose. The effective pixel fill factor for the FP systems was determined from the axial 1D presampling MTFs measured with a square sharp edge along the two orthogonal directions of the pixel lattice. Expectation MTFs were then calculated by averaging the radial MTFs over all possible phases and the 2D EMTF formed with the same reconstruction technique used for the 2D presampling MTF. The quantum NPS was then established by noise decomposition from homogenous images acquired as a function of detector air kerma. This was further decomposed into the correlated and uncorrelated quantum components by fitting the radially averaged quantum NPS with the radially averaged EMTF(2). This whole procedure allowed a detailed analysis of the influence of aliasing, signal and noise decorrelation, x-ray capture efficiency and global secondary gain on NPS and detector DQE. The influence of noise statistics, pixel fill factor and additional electronic and fixed pattern noises on the DQE was also studied. The 2D cascaded model and decompositions performed on the acquired images also enlightened the observed quantum NPS and DQE anisotropy.

Predicting trabecular bone microdamage initiation and accumulation using a non-linear perfect damage model

Studies evaluating the mechanical behavior of the trabecular microstructure play an important role in our understanding of pathologies such as osteoporosis, and in increasing our understanding of bone fracture and bone adaptation. Understanding of such behavior in bone is important for predicting and providing early treatment of fractures. The objective of this study is to present a numerical model for studying the initiation and accumulation of trabecular bone microdamage in both the pre- and post-yield regions. A sub-region of human vertebral trabecular bone was analyzed using a uniformly loaded anatomically accurate microstructural three-dimensional finite element model. The evolution of trabecular bone microdamage was governed using a non-linear, modulus reduction, perfect damage approach derived from a generalized plasticity stress-strain law. The model introduced in this paper establishes a history of microdamage evolution in both the pre- and post-yield regions

Determination of displacement, stress- and strain-distribution in the human heart: a FE-model on the basis of MR imaging.

The determination of characteristic cardiac parameters, such as displacement, stress and strain distribution are essential for an understanding of the mechanics of the heart. The calculation of these parameters has been limited until recently by the use of idealised mathematical representations of biventricular geometries and by applying simple material laws. On the basis of 20 short axis heart slices and in consideration of linear and nonlinear material behaviour we have developed a FE model with about 100,000 degrees of freedom. Marching Cubes and Phong's incremental shading technique were used to visualise the three dimensional geometry. In a quasistatic FE analysis continuous distribution of regional stress and strain corresponding to the endsystolic state were calculated. Substantial regional variation of the Von Mises stress and the total strain energy were observed at all levels of the heart model. The results of both the linear elastic model and the model with a nonlinear material description (Mooney-Rivlin) were compared. While the stress distribution and peak stress values were found to be comparable, the displacement vectors obtained with the nonlinear model were generally higher in comparison with the linear elastic case indicating the need to include nonlinear effects.

A musculoskeletal shoulder model based on pseudo-inverse and null-space optimization.

The goal of the present work was assess the feasibility of using a pseudo-inverse and null-space optimization approach in the modeling of the shoulder biomechanics. The method was applied to a simplified musculoskeletal shoulder model. The mechanical system consisted in the arm, and the external forces were the arm weight, 6 scapulo-humeral muscles and the reaction at the glenohumeral joint, which was considered as a spherical joint. The muscle wrapping was considered around the humeral head assumed spherical. The dynamical equations were solved in a Lagrangian approach. The mathematical redundancy of the mechanical system was solved in two steps: a pseudo-inverse optimization to minimize the square of the muscle stress and a null-space optimization to restrict the muscle force to physiological limits. Several movements were simulated. The mathematical and numerical aspects of the constrained redundancy problem were efficiently solved by the proposed method. The prediction of muscle moment arms was consistent with cadaveric measurements and the joint reaction force was consistent with in vivo measurements. This preliminary work demonstrated that the developed algorithm has a great potential for more complex musculoskeletal modeling of the shoulder joint. In particular it could be further applied to a non-spherical joint model, allowing for the natural translation of the humeral head in the glenoid fossa.

Neuroscience, quantum indeterminism and the Cartesian soul.

Quantum indeterminism is frequently invoked as a solution to the problem of how a disembodied soul might interact with the brain (as Descartes proposed), and is sometimes invoked in theories of libertarian free will even when they do not involve dualistic assumptions. Taking as example the Eccles-Beck model of interaction between self (or soul) and brain at the level of synaptic exocytosis, I here evaluate the plausibility of these approaches. I conclude that Heisenbergian uncertainty is too small to affect synaptic function, and that amplification by chaos or by other means does not provide a solution to this problem. Furthermore, even if Heisenbergian effects did modify brain functioning, the changes would be swamped by those due to thermal noise. Cells and neural circuits have powerful noise-resistance mechanisms, that are adequate protection against thermal noise and must therefore be more than sufficient to buffer against Heisenbergian effects. Other forms of quantum indeterminism must be considered, because these can be much greater than Heisenbergian uncertainty, but these have not so far been shown to play a role in the brain.

First Report of a Newborn Rat Ventilation Model for Bronchopulmonary Dysplasia Permitting Evaluation of Long-Term Outcome

Background: Bronchopulmonary dysplasia (BPD) remains the leading cause of chronic pulmonary morbidity among preterm neonates. However, the exact pathophysiology is still unknown. Here we present the first results from a new model inteAbstracts, 25th International Workshop on Surfactant Replacement 400 Neonatology 2010;97:395-400 grating the most common risk factors for BPD (lung immaturity, inflammation, mechanical ventilation (MV), oxygen), which allows long-term outcome evaluation due to a non-traumatic intubation procedure. Objectives: To test the feasibility of a new rat model by investigating effects of MV, inflammation and oxygen applied to immature lungs after a ventilation-free interval. Methods: On day 4, 5, or 6 newborn rats were given an intraperitoneal injection of lipopolysaccharides to induce a systemic inflammation. 24 h later they were anesthetized, endotracheally intubated and ventilated for 8 h with 60% oxygen. After weaning of anesthesia and MV the newborn rats were extubated and returned to their mothers. Two days later they were killed and outcome measurements were performed (histology, quantitative RT-PCR) and compared to animals investigated directly after MV. Results: Directly after MV, histological signs of ventilator-induced lung injury were found. After 48 h, the first signs of early BPD were seen with delayed alveolar formation. Expression of inflammatory genes was only transiently increased. After 48 h genes involved in alveolarization, such as matrix metalloproteinase-9 and tropoelastin, showed a significant change of their expression. Conclusion: For the first time we can evaluate in a newborn rat model the effects of MV after a ventilation-free interval. This allows discrimination between immediate response genes and delayed changes of expression of more structural genes involved in alveolarization.

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.

Biometal muscle to restore atrial transport function in a permanent atrial fibrillation animal model: a potential tool in the treatment of end-stage heart failure.

BACKGROUND: Half of the patients with end-stage heart failure suffer from persistent atrial fibrillation (AF). Atrial kick (AK) accounts for 10-15% of the ejection fraction. A device restoring AK should significantly improve cardiac output (CO) and possibly delay ventricular assist device (VAD) implantation. This study has been designed to assess the mechanical effects of a motorless pump on the right chambers of the heart in an animal model. METHODS: Atripump is a dome-shaped biometal actuator electrically driven by a pacemaker-like control unit. In eight sheep, the device was sutured onto the right atrium (RA). AF was simulated with rapid atrial pacing. RA ejection fraction (EF) was assessed with intracardiac ultrasound (ICUS) in baseline, AF and assisted-AF status. In two animals, the pump was left in place for 4 weeks and then explanted. Histology examination was carried out. The mean values for single measurement per animal with +/-SD were analysed. RESULTS: The contraction rate of the device was 60 per min. RA EF was 41% in baseline, 7% in AF and 21% in assisted-AF conditions. CO was 7+/-0.5 l min(-1) in baseline, 6.2+/-0.5 l min(-1) in AF and 6.7+/-0.5 l min(-1) in assisted-AF status (p<0.01). Histology of the atrium in the chronic group showed chronic tissue inflammation and no sign of tissue necrosis. CONCLUSIONS: The artificial muscle restores the AK and improves CO. In patients with end-stage cardiac failure and permanent AF, if implanted on both sides, it would improve CO and possibly delay or even avoid complex surgical treatment such as VAD implantation.