Use of the FACTS solvation model for protein-ligand docking calculations. Application to EADock.

Protein-ligand docking has made important progress during the last decade and has become a powerful tool for drug development, opening the way to virtual high throughput screening and in silico structure-based ligand design. Despite the flattering picture that has been drawn, recent publications have shown that the docking problem is far from being solved, and that more developments are still needed to achieve high successful prediction rates and accuracy. Introducing an accurate description of the solvation effect upon binding is thought to be essential to achieve this goal. In particular, EADock uses the Generalized Born Molecular Volume 2 (GBMV2) solvent model, which has been shown to reproduce accurately the desolvation energies calculated by solving the Poisson equation. Here, the implementation of the Fast Analytical Continuum Treatment of Solvation (FACTS) as an implicit solvation model in small molecules docking calculations has been assessed using the EADock docking program. Our results strongly support the use of FACTS for docking. The success rates of EADock/FACTS and EADock/GBMV2 are similar, i.e. around 75% for local docking and 65% for blind docking. However, these results come at a much lower computational cost: FACTS is 10 times faster than GBMV2 in calculating the total electrostatic energy, and allows a speed up of EADock by a factor of 4. This study also supports the EADock development strategy relying on the CHARMM package for energy calculations, which enables straightforward implementation and testing of the latest developments in the field of Molecular Modeling.

Characterizing emission and breathing-zone concentrations following exposure cases to fluororesin-based waterproofing spray mists.

Measurements and simulations were performed to assess workers' exposure to solvent vapors and aerosols during the waterproofing of a tiled surface. This investigation followed two recent incidents in the same company where workers experienced acute respiratory illness after spraying a stain-repellent resin containing fluorinated polymers on stone-tiled walls and floors. Because the waterproofing activity had been done for years at the tile company without encountering any exposure problems prior to these cases, it was strongly suspected that the incidents were linked to a recent change in the composition of the coating mixture. Experimental measurements and simulations indicated that the emission rate of particles smaller than 10 microm may be estimated at 0.66 mg/sec (SD 0.10) for the old resin and at 0.37 mg/sec (SD 0.04) for the new one. The measurement of the solvent emission rate from surfaces coated with the two resins indicated that shortly after spraying, the emission was in the range of 18 to 20 mg/sec x m2 and was similar for both products. Solvent and overspray emission rates were introduced in a two-zone compartment model. The results obtained in the near-field indicate significant exposure to overspray mist (7 and 34 mg/m3 for new resin) and solvent vapors (80 to 350 ppm for the new resin). It was also shown that the introduction of the new resin tended to significantly decrease the levels of solvents and particulates in the workers' breathing zone. These results strongly suggest that cases of acute respiratory illness are related to the specific toxicity of the fluorinated polymer itself. The fact that the same polymer is used in various commercial products raises concern regarding other possible occupational and domestic exposures.

Fatal tolperisone poisoning: Autopsy and toxicology findings in three suicide cases.

Tolperisone (Mydocalm(®)) is a centrally acting muscle relaxant with few sedative side effects that is used for the treatment of chronic pain conditions. We describe three cases of suicidal tolperisone poisoning in three healthy young subjects in the years 2006, 2008 and 2009. In all cases, macroscopic and microscopic autopsy findings did not reveal the cause of death. Systematic toxicological analysis (STA) including immunological tests, screening for volatile substances and blood, urine and gastric content screening by GC-MS and HPLC-DAD demonstrated the presence of tolperisone in all cases. In addition to tolperisone, only the analgesics paracetamol (acetaminophen), ibuprofen and naproxen could be detected. The blood ethanol concentrations were all lower than 0.10g/kg. Tolperisone was extracted by liquid-liquid extraction using n-chlorobutane as the extraction solvent. The quantification was performed by GC-NPD analysis of blood, urine and gastric content. Tolperisone concentrations of 7.0mg/l, 14mg/l and 19mg/l were found in the blood of the deceased. In the absence of other autopsy findings, the deaths in these three cases were finally explained as a result of lethal tolperisone ingestion. To the best of our knowledge, these three cases are the first reported cases of suicidal tolperisone poisonings.

Neutron structure of type-III antifreeze protein allows the reconstruction of AFP-ice interface.

Antifreeze proteins (AFPs) inhibit ice growth at sub-zero temperatures. The prototypical type-III AFPs have been extensively studied, notably by X-ray crystallography, solid-state and solution NMR, and mutagenesis, leading to the identification of a compound ice-binding surface (IBS) composed of two adjacent ice-binding sections, each which binds to particular lattice planes of ice crystals, poisoning their growth. This surface, including many hydrophobic and some hydrophilic residues, has been extensively used to model the interaction of AFP with ice. Experimentally observed water molecules facing the IBS have been used in an attempt to validate these models. However, these trials have been hindered by the limited capability of X-ray crystallography to reliably identify all water molecules of the hydration layer. Due to the strong diffraction signal from both the oxygen and deuterium atoms, neutron diffraction provides a more effective way to determine the water molecule positions (as D(2) O). Here we report the successful structure determination at 293 K of fully perdeuterated type-III AFP by joint X-ray and neutron diffraction providing a very detailed description of the protein and its solvent structure. X-ray data were collected to a resolution of 1.05 Å, and neutron Laue data to a resolution of 1.85 Å with a "radically small" crystal volume of 0.13 mm(3). The identification of a tetrahedral water cluster in nuclear scattering density maps has allowed the reconstruction of the IBS-bound ice crystal primary prismatic face. Analysis of the interactions between the IBS and the bound ice crystal primary prismatic face indicates the role of the hydrophobic residues, which are found to bind inside the holes of the ice surface, thus explaining the specificity of AFPs for ice versus water.

Partition coefficient and molecular flexibility: the concept of lipophilicity space.

The rationale of this study was to investigate molecular flexibility and its influence on physicochemical properties with a view to uncovering additional information on the fuzzy concept of dynamic molecular structure. Indeed, it is now known that computed molecular interaction fields (MIFs) such as molecular electrostatic potentials (MEPs) and lipophilicity potentials (MLPs) are conformation-dependent, as are dipole moments. A database of 125 compounds was used whose conformational space was explored, while conformation-dependent parameters were computed for each non-redundant conformer found in the conformational space of the compounds. These parameters were the virtual log P (log P(MLP), calculated by a MLP approach), the apolar surface area (ASA), polar surface area (PSA), and solvent-accessible surface (SAS). For each compound, the range taken by each parameter (its property space) was divided by the number of rotors taken as an index of flexibility, yielding a parameter termed 'molecular sensitivity'. This parameter was poorly correlated with others (i.e., it contains novel information) and showed the compounds to fall into two broad classes. 'Sensitive' molecules are those whose computed property ranges are markedly sensitive to conformational effects, whereas 'insensitive' (in fact, less sensitive) molecules have property ranges which are comparatively less affected by conformational fluctuations. A pharmacokinetic application is presented.

Sex differences in urinary levels of several biological indicators of exposure: a human volunteer study

The aim of the study was to quantify the variability on biological indicators of exposure between men and women for three well known solvents: methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1-trichloroethane. Another purpose was to explore the effect of selected CYP2E1 polymorphisms on the toxicokinetic profile. Controlled human exposures were carried out in a 12m(3) exposure chamber for each solvent separately, during 6h and at half of the threshold limit value. The human volunteers groups were composed of ten young men and fifteen young women, including ten women using hormonal contraceptive. An analysis of variance mainly showed an effect on the urinary levels of several biomarkers of exposure among women due to the use of hormonal contraceptive, with an increase of more than 50% in metabolites concentrations and a decrease of up to 50% in unchanged substances concentrations, suggesting an increase in their metabolism rate. The results also showed a difference due to the genotype CYP2E1*6, when exposed to methyl ethyl ketone, with a tendency to increase CYP2E1 activity when volunteers were carriers of the mutant allele. Our study suggests that not only physiological differences between men and women but also differences due to sex hormones levels can have an impact on urinary concentrations of several biomarkers of exposure. The observed variability due to sex among biological exposure indices can lead to misinterpretation of biomonitoring results. This aspect should have its place in the approaches for setting limits of occupational exposure. [Authors]

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.

Chemodiversity and molecular plasticity: recognition processes as explored by property spaces.

In the last few years, a need to account for molecular flexibility in drug-design methodologies has emerged, even if the dynamic behavior of molecular properties is seldom made explicit. For a flexible molecule, it is indeed possible to compute different values for a given conformation-dependent property and the ensemble of such values defines a property space that can be used to describe its molecular variability; a most representative case is the lipophilicity space. In this review, a number of applications of lipophilicity space and other property spaces are presented, showing that this concept can be fruitfully exploited: to investigate the constraints exerted by media of different levels of structural organization, to examine processes of molecular recognition and binding at an atomic level, to derive informative descriptors to be included in quantitative structure--activity relationships and to analyze protein simulations extracting the relevant information. Much molecular information is neglected in the descriptors used by medicinal chemists, while the concept of property space can fill this gap by accounting for the often-disregarded dynamic behavior of both small ligands and biomacromolecules. Property space also introduces some innovative concepts such as molecular sensitivity and plasticity, which appear best suited to explore the ability of a molecule to adapt itself to the environment variously modulating its property and conformational profiles. Globally, such concepts can enhance our understanding of biological phenomena providing fruitful descriptors in drug-design and pharmaceutical sciences.

Motility tests for drugs preventing PolyQ aggregation in recombinant Caenorhabditis elegans.

The pathological formation of proteinaceous aggregates that accumulate into the brain cells of patients are hallmarks of neurodegenerative diseases such as Alzheimer's disease, amyotrophic lateral sclerosis and the heterogeneous group of polyglutamine (polyQ) diseases. In the polyQ diseases, the most upstream events of the pathogenic cascade are the misfolding and aggregation of proteins, such as huntingtin in Huntington's disease, that contain expanded stretch of glutamine residues above 35--‐40 repeats. This expanded polyQ stretch triggers the misfolding and aggregation of cytotoxic polyQ proteins in the neurons that cause cell death through different processes, like apoptosis, excessive inflammation, formation of free radicals, eventually leading to neuronal loss and neurodegeneration. This study focuses on the cellular network of chaperone proteins that can prevent protein aggregation by binding misfolding intermediates and may, as in the case of HSP70, actively unfold misfolded proteins into refoldable non--‐toxic ones (Hinault et al., 2010; Sharma et al., 2011). The chaperones can also collaborate with the proteasome to convert stable harmful proteins into harmless amino acids. Thus, the chaperone proteins that are the most important cellular factors of prevention and curing of protein misfolding, are negatively affected by aging (Morley et al., 2002) and fail to act properly in the neurons of aged persons, which eventually may lead to neurodegenerative pathologies. The general aim of this research was to identify least toxic drugs that can upregulate the expression of chaperone genes in cells suffering from polyQ--‐ mediated protein aggregation and degeneration. The specific aim of this study was to observe the effect of ten drugs on polyQ aggregation in a recombinant nematode Caenorhabditis elegans expressing a chimeric protein containing a sequence of 35 glutamines (Q35) fused to the green fluorescent protein in muscle cells, which causes an age--‐ and temperature--‐ dependent phenotype of accelerated paralysis. The drugs were selected after having proven their causing the overexpression of chaperone proteins in a previous wide screening of 2000 drugs on the moss plant Physcomitrella patens. The screening that we performed in this study was on these ten drugs. It suggested that piroxicam and anisindione were good reducers of polyglutamine disease mediated paralysis. A hypothesis can be made that they may act as good enhancers of the heat shock response, which causes the overexpression of many HSP chaperones and thus reduce motility impairment of polyQ disease expressing nematodes. Piroxicam was found to have the greatest effect on reducing polyQ35 proteins aggregates mediated paralysis in a dose--‐dependent manner but was also found to either have a toxic effect on wild type C.elegans, either to change its natural motility behavior, eventually reducing its motility in both cases. Chloroform should be preferred over DMSO as a drug solvent as it appears to be less toxic to C.elegans.

Lipophilicity and its relationship with passive drug permeation.

In this review, we first summarize the structure and properties of biological membranes and the routes of passive drug transfer through physiological barriers. Lipophilicity is then introduced in terms of the intermolecular interactions it encodes. Finally, lipophilicity indices from isotropic solvent systems and from anisotropic membrane-like systems are discussed for their capacity to predict passive drug permeation across biological membranes such as the intestinal epithelium, the blood-brain barrier (BBB) or the skin. The broad evidence presented here shows that beyond the predictive power of lipophilicity parameters, the various intermolecular forces they encode allow a mechanistic interpretation of passive drug permeation.

Hand-disinfectant alcoholic vapors in incubators

Objective: To analyze the atmosphere inside incubators regarding alcoholic solvent such as isopropanol or ethanol which are commonly used in hand disinfecting solutions. Design: Observational. Setting: The third level neonatal unit of the Centre Hospitalier Universitaire Vaudois, Lausanne, Switzerland. Patients: Nine neonates with median (range) gestational age of 29 4/7 (25 5/7-39 0/7) weeks and birth weight of 960 (550-3050) grams. All neonates were inside incubators. Interventions: Alcoholic vapors inside incubators were directly and cumulatively measured by photoionisation and gas chromatography respectively after absorption on a charcoal sampling tube. Results: Eleven studies (mean study time: 230 ± 19 minutes) were performed. Highly variable isopropanol/ethanol concentrations profiles were found inside incubators. Peak value for isopropanol was 1982 part per million and for ethanol was 906 part per million. Conclusions: Incubators' inner atmosphere can be highly polluted by alcohol vapors. To reduce them staff should respect long evaporation time between hands disinfection and manipulations inside incubators. The use of an ethanol-based disinfecting solution, because of its short evaporation time, could be favored. As alcohol vapor toxicity for neonate remains largely unknown, further studies could be welcome.

Biological monitoring of chemical exposure : toxicokinetic differences due to age and sex

Human biomonitoring is a widely used method in the assessment of occupational exposure to chemical substances and recommended biological limits are published periodically for interpretation and decision-making. However, it is increasingly recognized that a large variability is associated with biological monitoring, making interpretation less efficient than assumed. In order to improve the applicability of biological monitoring, specific factors responsible for this variability should be identified and their contribution quantified. Among these factors, age and sex are easily identifiable, and present knowledge about pharmaceutical chemicals suggests that they play an important role on the toxicokinetics of occupational chemical agents, and therefore on the biological monitoring results.The aim of the present research project was to assess the influence of age and sex on biological indicators corresponding to organic solvents. This has been done experimentally and by toxicokinetic computer simulation. Another purpose was to explore the effect of selected CYP2E1 polymorphisms on the toxicokinetic profile.Age differences were identified by numerical simulations using a general toxicokinetic model from a previous study which was applied to 14 chemicals, representing 21 specific biological entities, with, among others, toluene, phenol, lead and mercury. These models were runn with the modified parameters, indicating in some cases important differences due to age. The expected changes are mostly of the order of 10-20 %, but differences up to 50 % were observed in some cases. These differences appear to depend on the chemical and on the biological entity considered.Sex differences were quantified by controlled human exposures, which were carried out in a 12 m3 exposure chamber for three organic solvents separately: methyl ethyl ketone, 1-methoxy-2-propanol and 1,1,1-trichloroethane. The human volunteer groups were composed 12 of ten young men and fifteen young women, the latter subdivided into those with and without hormonal contraceptive. They were exposed during six hours at rest and at half of the threshold limit value. The kinetics of the parent compounds (organic volatiles) and their metabolite(s) were followed in blood, urine and expired air over time. Analyses of the solvent and their metabolites were performed by using headspace gas chromatography, CYP2E1 genotypes by using PCR-based RFLP methods. Experimental data were used to calibrate the toxicokinetic models developed for the three solvents. The results obtained for the different biomarkers of exposure mainly showed an effect on the urinary levels of several biomarkers among women due to the use of hormonal contraceptive, with an increase of about 50 % in the metabolism rate. The results also showed a difference due to the genotype CYP2E1*6, when exposed to methyl ethyl ketone, with a tendency to increase CYP2E1 activity when volunteers were carriers of the mutant allele. Simulations showed that it is possible to use simple toxicokinetic tools in order to predict internal exposure when exposed to organic solvents. Our study suggests that not only physiological differences but also exogenous sex hormones could influence CYP2E1 enzyme activity. The variability among the urinary biological indicators levels gives evidence of an interindividual susceptibility, an aspect that should have its place in the approaches for setting limits of occupational exposure.

Sublimation of new matrix candidates for high spatial resolution imaging mass spectrometry of lipids: enhanced information in both positive and negative polarities after 1,5-diaminonapthalene deposition.

Matrix sublimation has demonstrated to be a powerful approach for high-resolution matrix-assisted laser desorption ionization (MALDI) imaging of lipids, providing very homogeneous solvent-free deposition. This work presents a comprehensive study aiming to evaluate current and novel matrix candidates for high spatial resolution MALDI imaging mass spectrometry of lipids from tissue section after deposition by sublimation. For this purpose, 12 matrices including 2,5-dihydroxybenzoic acid (DHB), sinapinic acid (SA), α-cyano-4-hydroxycinnamic acid (CHCA), 2,6-dihydroxyacetphenone (DHA), 2',4',6'-trihydroxyacetophenone (THAP), 3-hydroxypicolinic acid (3-HPA), 1,8-bis(dimethylamino)naphthalene (DMAN), 1,8,9-anthracentriol (DIT), 1,5-diaminonapthalene (DAN), p-nitroaniline (NIT), 9-aminoacridine (9-AA), and 2-mercaptobenzothiazole (MBT) were investigated for lipid detection efficiency in both positive and negative ionization modes, matrix interferences, and stability under vacuum. For the most relevant matrices, ion maps of the different lipid species were obtained from tissue sections at high spatial resolution and the detected peaks were characterized by matrix-assisted laser desorption ionization time-of-flight/time-of-flight (MALDI-TOF/TOF) mass spectrometry. First proposed for imaging mass spectrometry (IMS) after sublimation, DAN has demonstrated to be of high efficiency providing rich lipid signatures in both positive and negative polarities with high vacuum stability and sub-20 μm resolution capacity. Ion images from adult mouse brain were generated with a 10 μm scanning resolution. Furthermore, ion images from adult mouse brain and whole-body fish tissue sections were also acquired in both polarity modes from the same tissue section at 100 μm spatial resolution. Sublimation of DAN represents an interesting approach to improve information with respect to currently employed matrices providing a deeper analysis of the lipidome by IMS.

Intravitreous injection of PLGA microspheres encapsulating GDNF promotes the survival of photoreceptors in the rd1/rd1 mouse.

PURPOSE: To evaluate the potential delay of the retinal degeneration in rd1/rd1 mice using recombinant human glial cell line-derived neurotrophic factor (rhGDNF) encapsulated in poly(D,L-lactide-co-glycolide) (PLGA) microspheres. METHODS: rhGDNF-loaded PLGA microspheres were prepared using a water in oil in water (w/o/w) emulsion solvent extraction-evaporation process. In vitro, the rhGDNF release profile was assessed using radiolabeled factor. In vivo, rhGDNF microspheres, blank microspheres, or microspheres loaded with inactivated rhGDNF were injected into the vitreous of rd1/rd1 mice at postnatal day 11 (PN11). The extent of retinal degeneration was examined at PN28 using rhodopsin immunohistochemistry on whole flat-mount retinas, outer nuclear layer (ONL) cell counting on histology sections, and electroretinogram tracings. Immunohistochemical reactions for glial fibrillary acidic protein (GFAP), F4/80, and rhodopsin were performed on cryosections. RESULTS: Significant delay of rod photoreceptors degeneration was observed in mice receiving the rhGDNF-loaded microspheres compared to either untreated mice or to mice receiving blank or inactivated rhGDNF microspheres. The degeneration delay in the eyes receiving the rhGDNF microspheres was illustrated by the increased rhodopsin positive signals, the preservation of significantly higher number of cell nuclei within the ONL, and significant b-wave increase. A reduction of the subretinal glial proliferation was also observed in these treated eyes. No significant intraocular inflammatory reaction was observed after the intravitreous injection of the various microspheres. CONCLUSIONS: A single intravitreous injection of rhGDNF-loaded microspheres slows the retinal degeneration processes in rd1/rd1 mice. The use of injectable, biodegradable polymeric systems in the vitreous enables the efficient delivery of therapeutic proteins for the treatment of retinal diseases.

Contribution of electrostatic interactions, compactness and quaternary structure to protein thermostability: lessons from structural genomics of Thermotoga maritima.

Studies of the structural basis of protein thermostability have produced a confusing picture. Small sets of proteins have been analyzed from a variety of thermophilic species, suggesting different structural features as responsible for protein thermostability. Taking advantage of the recent advances in structural genomics, we have compiled a relatively large protein structure dataset, which was constructed very carefully and selectively; that is, the dataset contains only experimentally determined structures of proteins from one specific organism, the hyperthermophilic bacterium Thermotoga maritima, and those of close homologs from mesophilic bacteria. In contrast to the conclusions of previous studies, our analyses show that oligomerization order, hydrogen bonds, and secondary structure play minor roles in adaptation to hyperthermophily in bacteria. On the other hand, the data exhibit very significant increases in the density of salt-bridges and in compactness for proteins from T.maritima. The latter effect can be measured by contact order or solvent accessibility, and network analysis shows a specific increase in highly connected residues in this thermophile. These features account for changes in 96% of the protein pairs studied. Our results provide a clear picture of protein thermostability in one species, and a framework for future studies of thermal adaptation.