Prediction of toluene/water partition coefficients in the SAMPL9 blind challenge: assessment of machine learning and IEF-PCM/MST continuum solvation models.

In recent years the use of partition systems other than the widely used biphasic n-octanol/water has received increased attention to gain insight into the molecular features that dictate the lipophilicity of compounds. Thus, the difference between n-octanol/water and toluene/water partition coefficients has proven to be a valuable descriptor to study the propensity of molecules to form intramolecular hydrogen bonds and exhibit chameleon-like properties that modulate solubility and permeability. In this context, this study reports the experimental toluene/water partition coefficients (log Ptol/w) for a series of 16 drugs that were selected as an external test set in the framework of the Statistical Assessment of the Modeling of Proteins and Ligands (SAMPL) blind challenge. This external set has been used by the computational community to calibrate their methods in the current edition (SAMPL9) of this contest. Furthermore, the study also investigates the performance of two computational strategies for the prediction of log Ptol/w. The first relies on the development of two machine learning (ML) models, which are built up by combining the selection of 11 molecular descriptors in conjunction with either the multiple linear regression (MLR) or the random forest regression (RFR) model to target a dataset of 252 experimental log Ptol/w values. The second consists of the parametrization of the IEF-PCM/MST continuum solvation model from B3LYP/6-31G(d) calculations to predict the solvation free energies of 163 compounds in toluene and benzene. The performance of the ML and IEF-PCM/MST models has been calibrated against external test sets, including the compounds that define the SAMPL9 log Ptol/w challenge. The results are used to discuss the merits and weaknesses of the two computational approaches.

Predicting the Relative Binding Affinity for Reversible Covalent Inhibitors by Free Energy Perturbation Calculations.

Covalent inhibitors are assuming central importance in drug discovery projects, especially in this pandemic scenario. Many research groups have focused their attention on inhibiting viral proteases or human proteases such as cathepsin L (hCatL). The inhibition of these critical enzymes may impair viral replication. However, molecular modeling of covalent ligands is challenging since covalent and noncovalent ligand-bound states must be considered in the binding process. In this work, we evaluated the suitability of free energy perturbation (FEP) calculations as a tool for predicting the binding affinity of reversible covalent inhibitors of hCatL. Our strategy relies on the relative free energy calculated for both covalent and noncovalent complexes and the free energy changes have been compared with experimental data for eight nitrile-based inhibitors, including three new inhibitors of hCatL. Our results demonstrate that the covalent complex can be employed to properly rank the inhibitors. Nevertheless, a comparison of the free energy changes in both noncovalent and covalent states is valuable to interpret the effect triggered by the formation of the covalent bond on the interactions played by functional groups distant from the warhead. Overall, FEP can be employed as a powerful predictor tool in developing and understanding the activity of reversible covalent inhibitors.

Dioxygen Binding and Sensing Proteins.

Oxygen binding proteins (O2BIP) have been actively investigated for the past five decades due to their rich redox chemistry and function as O2 carriers in blood cells, as well as their function as gasotransmitters and sensors that modulate cellular signaling. A series of meetings on the periodic advances in the knowledge gained in the field of globin structure and function are conducted typically on a biannual basis. In the fall of 2018, the XXth International Conference was conducted, and very important articles with breakthrough discoveries were presented and very enthusiastically discussed. This was yet another highly successful meeting in the series. Select articles from this meeting were recently reviewed, updated, and published over several issues of Antioxidants and Redox Signaling, as Forum articles communicating the latest advances in this important area of redox biology. This Forum editorial introduces these articles and highlights their scientific significance in advancing the field. Each of these articles grew out of lectures presented in the meeting, and appears either as an original contribution or a comprehensive review in the journal. Overall, the articles published in the Forum provide in-depth details on the recent developments in the field as well as point the way to future directions. These Forum articles thus serve as an important summary of progress and the ongoing direction of this field, and serve to highlight recent advances in our understanding of O2BIP.

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations.

Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O2 and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.

Thermal Stability of Globins: Implications of Flexibility and Heme Coordination Studied by Molecular Dynamics Simulations.

Proteins are sensitive to temperature, and abrupt changes in the normal temperature conditions can have a profound impact on both structure and function, leading to protein unfolding. However, the adaptation of certain organisms to extreme conditions raises questions about the structural features that permit the structure and function of proteins to be preserved under these adverse conditions. To gain insight into the molecular basis of protein thermostability in the globin family, we have examined three representative examples: human neuroglobin, horse heart myoglobin, and Drosophila hemoglobin, which differ in their melting temperatures and coordination states of the heme iron in the absence of external ligands. In order to elucidate the possible mechanisms that govern the thermostability of these proteins, microsecond-scale classical molecular dynamics simulations were performed at different temperatures. Structural fluctuations and essential dynamics were analyzed, indicating that the flexibility of the CD region, which includes the two short C and D helixes and the connecting CD loop, is directly related to the thermostability. We observed that a larger inherent flexibility of the protein produces higher thermostability, probably concentrating the thermal fluctuations observed at high temperature in flexible regions, preventing unfolding. Globally, the results of this work improve our understanding of thermostability in the globin family.

Insertion of Isocyanides into N-Si Bonds: Multicomponent Reactions with Azines Leading to Potent Antiparasitic Compounds.

Trimethylsilyl chloride is an efficient activating agent for azines in isocyanide-based reactions, which then proceed through a key insertion of the isocyanide into a N-Si bond. The reaction is initiated by N activation of the azine, followed by nucleophilic attack of an isocyanide in a Reissert-type process. Finally, a second equivalent of the same or a different isocyanide inserts into the N-Si bond leading to the final adduct. The use of distinct nucleophiles leads to a variety of α-substituted dihydroazines after a selective cascade process. Based on computational studies, a mechanistic hypothesis for the course of these reactions was proposed. The resulting products exhibit significant activity against Trypanosoma brucei and T. cruzi, featuring favorable drug-like properties and safety profiles.

Combined experimental and computational investigation of the absorption spectra of E- and Z-cinnamic acids in solution: The peculiarity of Z-cinnamics.

Cinnamic acids are present in all kinds of plant tissues and hence in herbs and derived medicines, cosmetics and foods. The interest in their role in plants and their therapeutic applications has grown exponentially. Because of their molecular structure they can exist in E- and Z-forms, which are both found in plants. However, since only the E-forms are commercially available, very few in vitro and in vivo studies of the Z-form have been reported. In this work the physico-chemical properties of Z-cinnamic acids in solution have been examined by means of UV-absorption spectroscopy and high-level quantum mechanical computations. For each isomer similar absorption spectra were obtained in methanol and acetonitrile. However, distinct trends were found for Z- and E forms of cinnamic acids in water, where a higher hypsochromic shift of the Z-isomer relative to the E-form was observed. In general the wavelength of maximal absorption of the Z-form is dramatically blue shifted (-30 to -40 nm) to λ<280 nm, while a slightly blue shift of the absorption maxima for the corresponding E-form (+3 to -4 nm) was observed. This difference is associated with the non-planar, largely distorted, Z-structure and to the almost complete flat structure of the E-form. The results provide a basis for the study of functional and biotechnological roles of cinnamic acids and for the analysis of samples containing mixture of both geometric isomers.

Engineered chimeras reveal the structural basis of hexacoordination in globins: a case study of neuroglobin and myoglobin.

BACKGROUND: Myoglobin (Mb) and neuroglobin (Ngb) are representative members of pentacoordinated and bis-histidyl, hexacoordinated globins. In spite of their low sequence identity, they show surprisingly similar three-dimensional folds. The ability of Ngb to form a hexacoordinated bis-histidyl complex with the distal HisE7 has a strong impact on ligand affinity. The factors governing such different behaviors have not been completely understood yet, even though they are extremely relevant to establish structure-function relationships within the globin superfamily. METHODS: In this work we generated chimeric proteins by swapping a previously identified regulatory segment - the CD region - and evaluated comparatively the structural and functional properties of the resulting proteins by molecular dynamics simulations, and spectroscopic and kinetic investigations. RESULTS: Our results show that chimeric proteins display heme coordination properties displaced towards those expected for the corresponding CD region. In particular, in the absence of exogenous ligands, chimeric Mb is found as a partially hexacoordinated bis-histidyl species, whereas chimeric Ngb shows a lower equilibrium constant for forming a hexacoordinated bis-histidyl species. CONCLUSIONS: While these results confirm the regulatory role of the CD region for bis-histidyl hexacoordination, they also suggest that additional sources contribute to fine tune the equilibrium. General significance Globins constitute a ubiquitous group of heme proteins widely found in all kingdoms of life. These findings raise challenging questions regarding the structure-function relationships in these proteins, as bis-histidyl hexacoordination emerges as a novel regulatory mechanism of the physiological function of globins.

Role of the distal hydrogen-bonding network in regulating oxygen affinity in the truncated hemoglobin III from Campylobacter jejuni.

Oxygen affinity in heme-containing proteins is determined by a number of factors, such as the nature and conformation of the distal residues that stabilize the heme bound-oxygen via hydrogen-bonding interactions. The truncated hemoglobin III from Campylobacter jejuni (Ctb) contains three potential hydrogen-bond donors in the distal site: TyrB10, TrpG8, and HisE7. Previous studies suggested that Ctb exhibits an extremely slow oxygen dissociation rate due to an interlaced hydrogen-bonding network involving the three distal residues. Here we have studied the structural and kinetic properties of the G8(WF) mutant of Ctb and employed state-of-the-art computer simulation methods to investigate the properties of the O(2) adduct of the G8(WF) mutant, with respect to those of the wild-type protein and the previously studied E7(HL) and/or B10(YF) mutants. Our data indicate that the unique oxygen binding properties of Ctb are determined by the interplay of hydrogen-bonding interactions between the heme-bound ligand and the surrounding TyrB10, TrpG8, and HisE7 residues.

Molecular modeling of class I and II alleles of the major histocompatibility complex in Salmo salar.

Knowledge of the 3D structure of the binding groove of major histocompatibility (MHC) molecules, which play a central role in the immune response, is crucial to shed light into the details of peptide recognition and polymorphism. This work reports molecular modeling studies aimed at providing 3D models for two class I and two class II MHC alleles from Salmo salar (Sasa), as the lack of experimental structures of fish MHC molecules represents a serious limitation to understand the specific preferences for peptide binding. The reliability of the structural models built up using bioinformatic tools was explored by means of molecular dynamics simulations of their complexes with representative peptides, and the energetics of the MHC-peptide interaction was determined by combining molecular mechanics interaction energies and implicit continuum solvation calculations. The structural models revealed the occurrence of notable differences in the nature of residues at specific positions in the binding groove not only between human and Sasa MHC proteins, but also between different Sasa alleles. Those differences lead to distinct trends in the structural features that mediate the binding of peptides to both class I and II MHC molecules, which are qualitatively reflected in the relative binding affinities. Overall, the structural models presented here are a valuable starting point to explore the interactions between MHC receptors and pathogen-specific interactions and to design vaccines against viral pathogens.

Role of heme distortion on oxygen affinity in heme proteins: the protoglobin case.

The chemical properties of heme proteins largely reflect the electronic properties of their heme group. Often, the porphyrin ring of the heme exhibits significant distortions from its isolated structure, but the impact of these distortions on the chemical properties of the heme is yet uncertain. A systematic study focused on the effects of the distortion of the macrocycle on the binding affinity for oxygen is presented. The results show that out-of-plane distortions decrease the binding affinity, while in-plane distortions can increase or decrease it. Among in-plane distortions, only the breathing mode, which involves the symmetric compression-expansion of the porphyrin ring, strongly modulates the binding affinity. These findings shed light into the peculiar binding affinity of Methanosarcina acetivorans protoglobin, a protein that contains a highly distorted heme. Overall, the results highlight that in-plane distortions might be exploited by certain classes of heme proteins to modulate the ligand affinity.

Unraveling the molecular basis for ligand binding in truncated hemoglobins: the trHbO Bacillus subtilis case.

Truncated hemoglobins (trHbs) are heme proteins present in bacteria, unicellular eukaryotes, and higher plants. Their tertiary structure consists in a 2-over-2 helical sandwich, which display typically an inner tunnel/cavity system for ligand migration and/or storage. The microorganism Bacillus subtilis contains a peculiar trHb, which does not show an evident tunnel/cavity system connecting the protein active site with the solvent, and exhibits anyway a very high oxygen association rate. Moreover, resonant Raman results of CO bound protein, showed that a complex hydrogen bond network exists in the distal cavity, making it difficult to assign unambiguously the residues involved in the stabilization of the bound ligand. To understand these experimental results with atomistic detail, we performed classical molecular dynamics simulations of the oxy, carboxy, and deoxy proteins. The free energy profiles for ligand migration suggest that there is a key residue, GlnE11, that presents an alternate conformation, in which a wide ligand migration tunnel is formed, consistently with the kinetic data. This tunnel is topologically related to the one found in group I trHbs. On the other hand, the results for the CO and O(2) bound protein show that GlnE11 is directly involved in the stabilization of the cordinated ligand, playing a similar role as TyrB10 and TrpG8 in other trHbs. Our results not only reconcile the structural data with the kinetic information, but also provide additional insight into the general behaviour of trHbs. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

A microscopic study of the deoxyhemoglobin-catalyzed generation of nitric oxide from nitrite anion.

There is recent evidence suggesting that nitrite anion (NO 2 (-)) represents the major intravascular NO storage molecule whose transduction to NO is facilitated by a reduction mechanism catalyzed by deoxygenated hemoglobin (deoxy-Hb). In this work, we provide a detailed microscopic study of deoxy-Hb nitrite reductase (NIR) activity by combining classical molecular dynamics and hybrid quantum mechanical-molecular mechanical simulations. Our results point out that two alternative mechanisms could be operative and suggest that the most energetic barriers should stem from either reprotonation of the distal histidine or NO dissociation from the ferric heme. In the first proposed mechanism, which is similar to that proposed for bacterial NIRs, nitrite anion or nitrous acid coordinates to the heme through the N atom. This pathway involves HisE7 in a one or two proton transfer process, depending on whether the active species is nitrite anion or nitrous acid, to yield an intermediate Fe(III)NO species which eventually dissociates leading to NO and methemoglobin. In the second mechanism, the nitrite anion coordinates to the heme through the O atom. This pathway requires only one proton transfer from HisE7 and leads directly to the formation of a hydroxo Fe(III) complex and NO.

Structural determinants of ligand migration in Mycobacterium tuberculosis truncated hemoglobin O.

Mycobacterium tuberculosis is the causative agent of human tuberculosis, one of the most prevalent infectious diseases in the world. Its genome hosts the glbN and glbO genes coding for two proteins, truncated hemoglobin N (trHbN) and truncated hemoglobin O (trHbO), that belong to different groups (I and II, respectively) of the recently discovered trHb family of hemeproteins. The different expression pattern and kinetics rates constants for ligand association and NO oxidation rate suggest different functions for these proteins. Previous experimental and theoretical studies showed that, in trHbs, ligand migration along the internal tunnel cavity system is a key issue in determining the ligand-binding characteristics. The X-ray structure of trHbO has been solved and shows several internal cavities and secondary-docking sites. In this work, we present an extensive investigation of the tunnel/cavity system ofM. tuberculosis trHbO by means of computer-simulation techniques. We have computed the free-energy profiles for ligand migration along three found tunnels in the oxy and deoxy w.t. and mutant trHbO proteins. Our results show that multiple-ligand migration paths are possible and that several conserved residues such as TrpG8 play a key role in the ligand-migration regulation.

Nitric oxide reactivity with globins as investigated through computer simulation.

This chapter reviews the application of classical and quantum-mechanical atomistic simulation tools used in the investigation of several relevant issues in nitric oxide reactivity with globins and presents different simulation strategies based on classical force fields: standard molecular dynamics, essential dynamics, umbrella sampling, multiple steering molecular dynamics, and a novel technique for exploring the protein energy landscape. It also presents hybrid quantum-classical schemes as a tool to obtain relevant information regarding binding energies and chemical reactivity of globins. As illustrative examples, investigations of the structural flexibility, ligand migration profiles, oxygen affinity, and reactivity toward nitric oxide of truncated hemoglobin N of Mycobacterium tuberculosis are presented.

Mechanism of product release in NO detoxification from Mycobacterium tuberculosis truncated hemoglobin N.

The capability of Mycobacterium tuberculosis to rest in latency in the infected organism appears to be related to the disposal of detoxification mechanisms, which converts the nitric oxide (NO) produced by macrophages during the initial growth infection stage into a nitrate anion. Such a reaction appears to be associated with the truncated hemoglobin N (trHbN). Even though previous experimental and theoretical studies have examined the pathways used by NO and O2 to access the heme cavity, the eggression pathway of the nitrate anion is still a challenging question. In this work we present results obtained by means of classical and quantum chemistry simulations that show that trHbN is able to release rapidly the nitrate anion using an eggression pathway other than those used for the entry of both O2 and NO and that its release is promoted by hydration of the heme cavity. These results provide a detailed understanding of the molecular basis of the NO detoxification mechanism used by trHbN to guarantee an efficient NO detoxification and thus warrant survival of the microorganism under stress conditions.

Theoretical study of the truncated hemoglobin HbN: exploring the molecular basis of the NO detoxification mechanism.

Mycobacterium tuberculosis is the causative agent of human tuberculosis. The nitric oxide reaction with oxy-truncated hemoglobin N (trHbN) has been proposed to be responsible for the resistance mechanism by which this microorganism can evade the toxic effects of NO. In this work, we explore the molecular basis of the NO detoxification mechanism using a combination of classical and hybrid quantum-classical (QM-MM) simulation techniques. We have investigated the structural flexibility of the protein, the ligand affinity properties, and the nitric oxide reaction with coordinated O2. The analysis of the classical MD trajectory allowed us to identify Phe62 as the gate of the main channel for ligand diffusion to the active site. Moreover, the opening of the channel stems from the interplay between collective backbone motions and local rearrangements in the side chains of the residues that form the bottleneck of the tunnel. Even though the protein environment is not found to make a significant contribution to the heme moiety catalyzed reaction, the binding site influences the physiological function of the enzyme at three different levels. First, by isolating the intermediates formed in the reaction, it prevents nondesired reactions from proceeding. Second, it modulates the ligand (O2, NO) affinity of the protein, which can be ascribed to both distal and proximal effects. Finally, the stabilization of the Tyr33-Gln58 pair upon O2 binding might alter the essential dynamics of the protein, leading in turn to a mechanism for ligand-induced regulation.

Acute intestinal anisakiasis in Spain: a fourth-stage Anisakis simplex larva.

A case of acute intestinal anisakiasis has been reported; a nematode larva being found in the submucosa of the ileum of a woman in Jaén (Spain). The source of infection was the ingestion of raw Engraulis encrasicholus. On the basis of its morphology, the worm has been identified as a fourth-stage larva of Anisakis simplex. In Spain, this is the ninth report of human anisakiasis and also probably the first case of anisakiasis caused by a fourth-stage larva of A. simplex.