A rational approach to elucidate human monoamine oxidase molecular selectivity

Designing highly selective human monoamine oxidase (hMAO) inhibitors is a challenging goal on the road to a more effective treatment of depression and anxiety (inhibition of hMAO-A isoform) as well as neurodegenerative diseases (inhibition of hMAO-B isoform). To uncover the molecular rationale of hMAOs selectivity, two recently prepared 2H-chromene-2-ones, namely compounds 1 and 2, were herein chosen as molecular probes being highly selective toward hMAO-A and hMAO-B, respectively. We performed molecular dynamics (MD) studies on four different complexes, cross-simulating one at a time the two hMAO-isoforms (dimer embedded in a lipid bilayer) with the two considered probes. Our comparative analysis on the obtained 100 ns trajectories discloses a stable H-bond interaction between 1 and Gln215 as crucial for ligand selectivity toward hMAO-A whereas a water-mediated interaction might explain the observed hMAO-B selectivity of compound 2. Such hypotheses are further supported by binding free energy calculations carried out applying the molecular mechanics generalized Born surface area (MM-GBSA) method and allowing us to evaluate the contribution of each residue to the observed isoform selectivity. Taken as whole, this study represents the first attempt to explain at molecular level hMAO isoform selectivity and a valuable yardstick for better addressing the design of new and highly selective MAO inhibitors

Organic bioelectronics probing conformational changes in surface confined proteins

The study of proteins confined on a surface has attracted a great deal of attention due to its relevance in the development of bio-systems for laboratory and clinical settings. In this respect, organic bio-electronic platforms can be used as tools to achieve a deeper understanding of the processes involving protein interfaces. In this work, biotin-binding proteins have been integrated in two different organic thin-film transistor (TFT) configurations to separately address the changes occurring in the protein-ligand complex morphology and dipole moment. This has been achieved by decoupling the output current change upon binding, taken as the transducing signal, into its component figures of merit. In particular, the threshold voltage is related to the protein dipole moment, while the field-effect mobility is associated with conformational changes occurring in the proteins of the layer when ligand binding occurs. Molecular Dynamics simulations on the whole avidin tetramer in presence and absence of ligands were carried out, to evaluate how the tight interactions with the ligand affect the protein dipole moment and the conformation of the loops surrounding the binding pocket. These simulations allow assembling a rather complete picture of the studied interaction processes and support the interpretation of the experimental results

CATMoS: Collaborative Acute Toxicity Modeling Suite.

BACKGROUND: Humans are exposed to tens of thousands of chemical substances that need to be assessed for their potential toxicity. Acute systemic toxicity testing serves as the basis for regulatory hazard classification, labeling, and risk management. However, it is cost- and time-prohibitive to evaluate all new and existing chemicals using traditional rodent acute toxicity tests. In silico models built using existing data facilitate rapid acute toxicity predictions without using animals. OBJECTIVES: The U.S. Interagency Coordinating Committee on the Validation of Alternative Methods (ICCVAM) Acute Toxicity Workgroup organized an international collaboration to develop in silico models for predicting acute oral toxicity based on five different end points: Lethal Dose 50 (LD50 value, U.S. Environmental Protection Agency hazard (four) categories, Globally Harmonized System for Classification and Labeling hazard (five) categories, very toxic chemicals [LD50 (LD50≤50mg/kg)], and nontoxic chemicals (LD50>2,000mg/kg). METHODS: An acute oral toxicity data inventory for 11,992 chemicals was compiled, split into training and evaluation sets, and made available to 35 participating international research groups that submitted a total of 139 predictive models. Predictions that fell within the applicability domains of the submitted models were evaluated using external validation sets. These were then combined into consensus models to leverage strengths of individual approaches. RESULTS: The resulting consensus predictions, which leverage the collective strengths of each individual model, form the Collaborative Acute Toxicity Modeling Suite (CATMoS). CATMoS demonstrated high performance in terms of accuracy and robustness when compared with in vivo results. DISCUSSION: CATMoS is being evaluated by regulatory agencies for its utility and applicability as a potential replacement for in vivo rat acute oral toxicity studies. CATMoS predictions for more than 800,000 chemicals have been made available via the National Toxicology Program's Integrated Chemical Environment tools and data sets (ice.ntp.niehs.nih.gov). The models are also implemented in a free, standalone, open-source tool, OPERA, which allows predictions of new and untested chemicals to be made. https://doi.org/10.1289/EHP8495

P3HT and PBTTT: a comparative molecular dynamics study inspired by OTFT fabrication techniques

P3HT and PBTTT: a comparative molecular dynamics study inspired by OTFT fabrication techniques: 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. The most used organic material for their construction consists in 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 and the interface with biologic material (phospholipid bilayer) remain wide open. We have recently approached this field of research by applying Molecular Dynamics simulations to investigate bulk and interface properties of these two polymers. Our simulations were inspired by OTFT fabrication techniques: in particular, we have characterized the effects produced by the annealing process, widely used to enhance the charge transport properties of these devices, and by the presence of residual solvent molecules after the spin coating deposition procedure. We will present the results obtained in bulk and at the interfaces polymer/vacuum, polymer/SiO2 and polymer/phospholipidic membrane, discuss their relevance in the phenomenological context and highlight a roadmap for future plans of research

Permeability Coefficients of Lipophilic Compounds Estimated by Computer Simulations

The ability of a drug to cross the intestine-blood barrier is a key quantity for drug design and employment and is normally quantified by the permeability coefficient P, often evaluated in the so-called Caco-2 assay. This assay is based on measuring the initial growth rate of the concentration of the drug beyond the cellular barrier but not its steady-state flux through the membrane. This might lead to confusion since, in the case of lipophilic drugs, the initial slope is strongly affected by the retention of the drug in the membrane. This effect is well known but seldom considered in the assay. Here, we exploit all-atoms molecular dynamics and bias exchange metadynamics to calculate the concentration of two lipophilic drugs across a model membrane as a function of time. This allows estimating both the steady-state flux and the initial slope of the concentration growth and comparing Caco-2 and steady-state estimates of P. We show that our computational procedure is able to reproduce the experimental values, although these may differ from the permeability coefficients by orders of magnitude. Our findings are generalized by a simplified one-dimensional model of the permeation process that may act as a roadmap to assess which measure of membrane permeability would be more appropriate and, consequently, whether retention corrections should be included in estimates based on Caco-2 assays

Prediction of Acute Oral Systemic Toxicity Using a Multifingerprint Similarity Approach

The implementation of nonanimal approaches is of particular importance to regulatory agencies for the prediction of potential hazards associated with acute exposures to chemicals. This work was carried out in the framework of an international modeling initiative organized by the Acute Toxicity Workgroup (ATWG) of the Interagency Coordinating Committee on the Validation of Alternative Methods (ICCVAM) with the participation of 32 international groups across government, industry, and academia. Our contribution was to develop a multifingerprints similarity approach for predicting five relevant toxicology endpoints related to the acute oral systemic toxicity that are: The median lethal dose (LD 50) point prediction, the "nontoxic" (LD 50 > 2000 mg/kg) and "very toxic" (LD 50 <50 mg/kg) binary classification, and the multiclass categorization of chemicals based on the United States Environmental Protection Agency and Globally Harmonized System of Classification and Labeling of Chemicals schemes. Provided by the ICCVAM's ATWG, the training set used to develop the models consisted of 8944 chemicals having high-quality rat acute oral lethality data. The proposed approach integrates the results coming from a similarity search based on 19 different fingerprint definitions to return a consensus prediction value. Moreover, the herein described algorithm is tailored to properly tackling the so-called toxicity cliffs alerting that a large gap in LD 50 values exists despite a high structural similarity for a given molecular pair. An external validation set made available by ICCVAM and consisting in 2896 chemicals was employed to further evaluate the selected models. This work returned high-Accuracy predictions based on the evaluations conducted by ICCVAM's ATWG

Human Aquaporin-4 and Molecular Modeling: Historical Perspective and View to the Future

Among the different aquaporins (AQPs), human aquaporin-4 (hAQP4) has attracted the greatest interest in recent years as a new promising therapeutic target. Such a membrane protein is, in fact, involved in a multiple sclerosis-like immunopathology called Neuromyelitis Optica (NMO) and in several disorders resulting from imbalanced water homeostasis such as deafness and cerebral edema. The gap of knowledge in its functioning and dynamics at the atomistic level of detail has hindered the development of rational strategies for designing hAQP4 modulators. The application, lately, of molecular modeling has proved able to fill this gap providing a breeding ground to rationally address compounds targeting hAQP4. In this review, we give an overview of the important advances obtained in this field through the application of Molecular Dynamics (MD) and other complementary modeling techniques. The case studies presented herein are discussed with the aim of providing important clues for computational chemists and biophysicists interested in this field and looking for new challenges

Effects of Different Self-Assembled Monolayers on Thin-Film Morphology: A Combined DFT/MD Simulation Protocol

Organic thin film transistors (OTFTs) are multilayer field-effect transistors that employ an organic conjugated material as semiconductor. Several experimental groups have recently demonstrated that the insertion of an organic self-assembled monolayer (SAM) between the dielectric and the semiconductive layer is responsible for a sensible improvement of the OTFT performances in terms of an increased charge carrier mobility caused by a higher degree of order in the organic semiconductor layer. Here, we describe a combined periodic density functional theory (DFT) and classical molecular dynamics (MD) protocol applied to four different SAMs and a pentacene monolayer deposited onto their surfaces. In particular, we investigate the morphology and the surface of the four SAMs and the translational, orientational, and nematic order of the monolayer through the calculation of several distribution functions and order parameters pointing out the differences among the systems and relating them to known experimental results. Our calculations also suggest that small differences in the SAM molecular design will produce remarkable differences in the SAM surface and monolayer order. In particular, our simulations explain how a SAM with a bulky terminal group results in an irregular and rough surface that determines the deposition of a disordered semiconductive monolayer. On the contrary, SAMs with a small terminal group generate smooth surfaces with uninterrupted periodicity, thus favoring the formation of an ordered pentacene monolayer that increases the mobility of charge carriers and improves the overall performances of the OTFT devices. Our results clearly point out that the in silico procedure presented here might be of help in tuning the design of SAMs in order to improve the quality of OTFT devices