Aquaporin-4 (AQP4) is the predominant water channel in different organs and tissues. An alteration of its physiological functioning is responsible for several disorders of water regulation and, thus, is considered an attractive target with a promising therapeutic and diagnostic potential. Molecular dynamics (MD) simulations performed on the AQP4 tetramer embedded in a bilayer of lipid molecules allowed us to analyze the role of spontaneous fluctuations occurring inside the pore. Following the approach by Hashido et al. [Hashido M, Kidera A, Ikeguchi M (2007) Biophys J 93: 373-385], our analysis on 200ns trajectory discloses three domains inside the pore as key elements for water permeation. Herein, we describe the gating mechanism associated with the well-known selectivity filter on the extracellular side of the pore and the crucial regulation ensured by the NPA motifs (asparagine, proline, alanine). Notably, on the cytoplasmic side, we find a putative gate formed by two residues, namely, a cysteine belonging to the loop D (C178) and a histidine from loop B (H95). We observed that the spontaneous reorientation of the imidazole ring of H95 acts as a molecular switch enabling H-bond interaction with C178. The occurrence of such local interaction seems to be responsible for the narrowing of the pore and thus of a remarkable decrease in water flux rate. Our results are in agreement with recent experimental observations and may represent a promising starting point to pave the way for the discovery of chemical modulators of AQP4 water permeability.

We explore the relation between the morphological and the charge transport properties of poly(3-hexylthiophene) (P3HT) and poly(2,5-bis(3-alkylthiophen-2-yl)thieno[3,2-b]thiophene) (PBTTT) semiconductor polymers in both amorphous and crystalline phases. Using molecular dynamics to simulate bulk supercells and the Marcus theory to analyze the transport properties we found that amorphous systems display a reduced hole mobility due to the loss of nematic order and π-π stacking leading to a reduction in the electronic coupling between two chains. In the crystal phase, PBTTT displays a larger charge mobility than P3HT due to the interdigitation of the side chains enhancing the stability of the conjugated rings on the backbones. This more stable π-π stacking reduces the energetic disorder with respect to P3HT and increases the electronic coupling. In contrast, in the amorphous phase, PBTTT displays a reduced charge mobility with respect to P3HT due to the absence of side chains attached to the thienothiophenes, which increases their fluctuations and the energetic disorder. In addition, we show that it is possible to calculate the reorganization energy neglecting the side chains of the polymers and thus saving computational time. Within this approximation, we obtained mobility values matching the experimental measurements, thus confirming that the side chains are crucial to shape the morphology of the polymeric systems but are not involved in the charge transport process.

Organic thin film transistors (OTFT) are metal-insulator-semiconductor field-effect transistors in which the semiconductor is a conjugated organic material. They are the subject of intense industrial research because their fabrication process is less expensive when compared with inorganic TFTs. Among the others, the organic material mostly employed in their construction consists of two semiconductor polymers, namely poly(3-hexylthiophene) (P3HT) and poly(2,5-bis(3-alkylthiophen-2-yl)thieno[3,2-b]thiophene) (PBTTT). Despite the large amount of experimental efforts in the characterization of the electronic properties of these devices, several questions regarding their morphological arrangement in bulk and at interfaces remain wide open. Here, we report results obtained by classical molecular dynamics simulations of P3HT and PBTTT inspired by OTFT fabrication techniques. In particular, we investigate how the annealing fabrication process and the presence of residual solvent molecules left over after spin coating might modify the morphology and the dynamics of the amorphous phase of these two polymers. Simulations of both polymer deposits at 300 K after annealing show an increase in the number of interdigitation events between the alkyl chains of two polymeric macromolecules; moreover, we find that the increased stability of the pi-pi stacking is caused by an improved layering of the films, which may account for the better charge transport properties reported in experiments. Our results strongly suggest that thin semiconductor films are required to boost the performances of the devices and that a minimal presence of residual solvent does not alter dramatically the microscopic structure and stability of the polymeric films

We present MuSSeL, a multifingerprint similarity search algorithm, able to predict putative drug targets for a given query small molecule as well as to return a quantitative assessment of its bioactivity in terms of K i or IC 50 values. Predictions are automatically made exploiting a large collection of high quality experimental bioactivity data available from ChEMBL (version 22.1) combining, in a consensus-like approach, predictions resulting from a similarity search performed using 13 different fingerprint definitions. Importantly, the herein proposed algorithm is also effective in detecting and handling activity cliffs. A calibration set including small molecules present in the last updated version of ChEMBL (version 23) was employed to properly tune the algorithm parameters. Three randomly built external sets were instead challenged for model performances. The potential use of MuSSeL was also challenged by a prospective exercise for the prediction of five bioactive compounds taken from articles published in the Journal of Medicinal Chemistry just few months ago. The paper emphasizes the importance of implementing multifingerprint consensus strategies to increase the confidence in prediction of similarity search algorithms and provides a fast and easy-to-run tool for drug target and bioactivity prediction.

Key pointsr Loss-of-function mutations of the skeletal muscle ClC-1 channel cause myotonia congenita with variable phenotypes.r Using patch clamp we show that F484L, located in the conducting pore, probably induces mild dominant myotonia by right-shifting the slow gating of ClC-1 channel, without exerting a dominant-negative effect on the wild-type (WT) subunit.r Molecular dynamics simulations suggest that F484L affects the slow gate by increasing the frequency and the stability of H-bond formation between E232 in helix F and Y578 in helix R.r Three other myotonic ClC-1 mutations are shown to produce distinct effects on channel function: L198P shifts the slow gate to positive potentials, V640G reduces channel activity, while L628P displays a WT-like behaviour (electrophysiology data only).r Our results provide novel insight into the molecular mechanisms underlying normal and altered ClC-1 function.Abstract Myotonia congenita is an inherited disease caused by loss-of-function mutations of the skeletal muscle ClC-1 chloride channel, characterized by impaired muscle relaxation after contraction and stiffness. In the present study, we provided an in-depth characterization of F484L, a mutation previously identified in dominant myotonia, in order to define the genotype-phenotype correlation, and to elucidate the contribution of this pore residue to the mechanisms of ClC-1 gating. Patch-clamp recordings showed that F484L reduced chloride currents at every tested potential and dramatically right-shifted the voltage dependence of slow gating, thus contributing to the mild clinical phenotype of affected heterozygote carriers. Unlike dominant mutations located at the dimer interface, no dominant-negative effect was observed when F484L mutant subunits were co-expressed with wild type. Molecular dynamics simulations further revealed that F484L affected the slow gate by increasing the frequency and stability of the H-bond formation between the pore residue E232 and the R helix residue Y578. In addition, using patch-clamp electrophysiology, we characterized three other myotonic ClC-1 mutations. We proved that the dominant L198P mutation in the channel pore also right-shifted the voltage dependence of slow gating, recapitulating mild myotonia. The recessive V640G mutant drastically reduced channel function, which probably accounts for myotonia. In contrast, the recessive L628P mutant produced currents very similar to wild type, suggesting that the occurrence of the compound truncating mutation (Q812X) or other muscle-specific mechanisms accounted for the severe symptoms observed in this family. Our results provide novel insight into the molecular mechanisms underlying normal and altered ClC-1 function.