Structural Space of the Duffy Antigen/Receptor for Chemokines' Intrinsically Disordered Ectodomain 1 Explored by Temperature Replica-Exchange Molecular Dynamics Simulations.

Plasmodium vivax malaria affects 14 million people each year. Its invasion requires interactions between the parasitic Duffy-binding protein (PvDBP) and the N-terminal extracellular domain (ECD1) of the host's Duffy antigen/receptor for chemokines (DARC). ECD1 is highly flexible and intrinsically disordered, therefore it can adopt different conformations. We computationally modeled the challenging ECD1 local structure. With T-REMD simulations, we sampled its dynamic behavior and collected its most representative conformations. Our results suggest that most of the DARC ECD1 domain remains in a disordered state during the simulated time. Globular local conformations are found in the analyzed local free-energy minima. These globular conformations share an α-helix spanning residues Ser18 to Ser29 and in many cases they comprise an antiparallel ß-sheet, whose ß-strands are formed around residues Leu10 and Ala49. The formation of a parallel ß-sheet is almost negligible. So far, progress in understanding the mechanisms forming the basis of the P. vivax malaria infection of reticulocytes has been hampered by experimental difficulties, along with a lack of DARC structural information. Our collection of the most probable ECD1 structural conformations will help to advance modeling of the DARC structure and to explore DARC-ECD1 interactions with a range of physiological and pathological ligands.

Metadynamics simulations of ligands binding to protein surfaces: a novel tool for rational drug design.

Structure-based drug design protocols may encounter difficulties to investigate poses when the biomolecular targets do not exhibit typical binding pockets. In this study, by providing two concrete examples from our labs, we suggest that the combination of metadynamics free energy methods (validated against affinity measurements), along with experimental structural information (by X-ray crystallography and NMR), can help to identify the poses of ligands on protein surfaces. The simulation workflow proposed here was implemented in a widely used code, namely GROMACS, and it could straightforwardly be applied to various drug-design campaigns targeting ligands' binding to protein surfaces.

A comparative study of receptor interactions between SARS-CoV and SARS-CoV-2 from molecular modeling.

The pandemic of COVID-19 severe acute respiratory syndrome, which was fatal for millions of people worldwide, triggered the race to understand in detail the molecular mechanisms of this disease. In this work, the differences of interactions between the SARS-CoV/SARS-CoV-2 Receptor binding domain (RBD) and the human Angiotensin Converting Enzyme 2 (ACE2) receptor were studied using in silico tools. Our results show that SARS-CoV-2 RBD is more stable and forms more interactions with ACE2 than SARS-CoV. At its interface, three stable binding patterns are observed and named red-K31, green-K353 and blue-M82 according to the central ACE2 binding residue. In SARS-CoV instead, only the first two binding patches are persistently formed during the MD simulation. Our MM/GBSA calculations indicate the binding free energy difference of about 2.5 kcal/mol between SARS-CoV-2 and SARS-CoV which is compatible with the experiments. The binding free energy decomposition points out that SARS-CoV-2 RBD-ACE2 interactions of the red-K31 ([Formula: see text]) and blue-M82 ([Formula: see text]) patterns contribute more to the binding affinity than in SARS-CoV ([Formula: see text] for red-K31), while the contribution of the green-K353 pattern is very similar in the two strains ([Formula: see text] and [Formula: see text] for SARS-CoV-2 and SARS-CoV, respectively). Five groups of mutations draw our attention at the RBD-ACE2 binding interface, among them, the mutation -PPA469-471/GVEG482-485 has the most important and favorable impact on SARS-CoV-2 binding to the ACE2 receptor. These results, highlighting the molecular differences in the binding between the two viruses, contribute to the common knowledge about the new corona virus and to the development of appropriate antiviral treatments, addressing the necessity of ongoing pandemics.

The Interplay of Cholesterol and Ligand Binding in hTSPO from Classical Molecular Dynamics Simulations.

The translocator protein (TSPO) is a 18kDa transmembrane protein, ubiquitously present in human mitochondria. It is overexpressed in tumor cells and at the sites of neuroinflammation, thus representing an important biomarker, as well as a promising drug target. In mammalian TSPO, there are cholesterol-binding motifs, as well as a binding cavity able to accommodate different chemical compounds. Given the lack of structural information for the human protein, we built a model of human (h) TSPO in the apo state and in complex with PK11195, a molecule routinely used in positron emission tomography (PET) for imaging of neuroinflammatory sites. To better understand the interactions of PK11195 and cholesterol with this pharmacologically relevant protein, we ran molecular dynamics simulations of the apo and holo proteins embedded in a model membrane. We found that: (i) PK11195 stabilizes hTSPO structural fold; (ii) PK11195 might enter in the binding site through transmembrane helices I and II of hTSPO; (iii) PK11195 reduces the frequency of cholesterol binding to the lower, N-terminal part of hTSPO in the inner membrane leaflet, while this impact is less pronounced for the upper, C-terminal part in the outer membrane leaflet, where the ligand binding site is located; (iv) very interestingly, cholesterol most frequently binds simultaneously to the so-called CRAC and CARC regions in TM V in the free form (residues L150-X-Y152-X(3)-R156 and R135-X(2)-Y138-X(2)-L141, respectively). However, when the protein is in complex with PK11195, cholesterol binds equally frequently to the CRAC-resembling motif that we observed in TM I (residues L17-X(2)-F20-X(3)-R24) and to CRAC in TM V. We expect that the CRAC-like motif in TM I will be of interest in future experimental investigations. Thus, our MD simulations provide insight into the structural features of hTSPO and the previously unknown interplay between PK11195 and cholesterol interactions with this pharmacologically relevant protein.

Data set of intrinsically disordered proteins analysed at a local protein conformation level.

Intrinsic Disorder Proteins (IDPs) have become a hot topic since their characterisation in the 90s. The data presented in this article are related to our research entitled "A structural entropy index to analyse local conformations in Intrinsically Disordered Proteins" published in Journal of Structural Biology [1]. In this study, we quantified, for the first time, continuum from rigidity to flexibility and finally disorder. Non-disordered regions were also highlighted in the ensemble of disordered proteins. This work was done using the Protein Ensemble Database (PED), which is a useful database collecting series of protein structures considered as IDPs. The data set consists of a collection of cleaned protein files in classical pdb format that can be readily used as an input with most automatic analysis software. The accompanying data include the coding of all structural information in terms of a structural alphabet, namely Protein Blocks (PBs). An entropy index derived from PBs that allows apprehending the continuum between protein rigidity to flexibility to disorder is included, with information from secondary structure assignment, protein accessibility and prediction of disorder from the sequences. The data may be used for further structural bioinformatics studies of IDPs. It can also be used as a benchmark for evaluating disorder prediction methods.

Identification of Anti-Mycobacterium and Anti-Legionella Compounds With Potential Distinctive Structural Scaffolds From an HD-PBL Using Phenotypic Screens in Amoebae Host Models.

Tubercular Mycobacteria and Legionella pneumophila are the causative agents of potentially fatal respiratory diseases due to their intrinsic pathogenesis but also due to the emergence of antibiotic resistance that limits treatment options. The aim of our study was to explore the antimicrobial activity of a small ligand-based chemical library of 1255 structurally diverse compounds. These compounds were screened in a combination of three assays, two monitoring the intracellular growth of the pathogenic bacteria, Mycobacterium marinum and L. pneumophila, and one assessing virulence of M. marinum. We set up these assays using two amoeba strains, the genetically tractable social amoeba Dictyostelium discoideum and the free-living amoeba Acanthamoeba castellanii. In summary, 64 (5.1%) compounds showed anti-infective/anti-virulence activity in at least one of the three assays. The intracellular assays hit rate varied between 1.7% (n = 22) for M. marinum and 2.8% (n = 35) for L. pneumophila with seven compounds in common for both pathogens. In parallel, 1.2% (n = 15) of the tested compounds were able to restore D. discoideum growth in the presence of M. marinum spiked in a lawn of food bacteria. We also validated the generality of the hits identified in the A. castellanii-M. marinum anti-infective screen using the D. discoideum-M. marinum host-pathogen model. The characterization of anti-infective and antibacterial hits in the latter infection model revealed compounds able to reduce intracellular growth more than 50% at 30 µM. Moreover, the chemical space and physico-chemical properties of the anti-M. marinum hits were compared to standard and candidate Mycobacterium tuberculosis (Mtb) drugs using ChemGPS-NP. A principle component analysis identified separate clusters for anti-M. marinum and anti-L. pneumophila hits unveiling the potentially new physico-chemical properties of these hits compared to standard and candidate M. tuberculosis drugs. Our studies underscore the relevance of using a combination of low-cost and low-complexity assays with full 3R compliance in concert with a rationalized focused library of compounds to identify new chemical scaffolds and to dissect some of their properties prior to taking further steps toward compound development.

A structural entropy index to analyse local conformations in intrinsically disordered proteins.

Sequence - structure - function paradigm has been revolutionized by the discovery of disordered regions and disordered proteins more than two decades ago. While the definition of rigidity is simple with X-ray structures, the notion of flexibility is linked to high experimental B-factors. The definition of disordered regions is more complex as in these same X-ray structures; it is associated to the position of missing residues. Thus a continuum so seems to exist between rigidity, flexibility and disorder. However, it had not been precisely described. In this study, we used an ensemble of disordered proteins (or regions) and, we applied a structural alphabet to analyse their local conformation. This structural alphabet, namely Protein Blocks, had been efficiently used to highlight rigid local domains within flexible regions and so discriminates deformability and mobility concepts. Using an entropy index derived from this structural alphabet, we underlined its interest to measure these local dynamics, and to quantify, for the first time, continuum states from rigidity to flexibility and finally disorder. We also highlight non-disordered regions in the ensemble of disordered proteins in our study.

In Silico Driven Design and Synthesis of Rhodanine Derivatives as Novel Antibacterials Targeting the Enoyl Reductase InhA.

Here, we report on the design, synthesis, and biological evaluation of 4-thiazolidinone (rhodanine) derivatives targeting Mycobacterial tuberculosis (Mtb) trans-2-enoyl-acyl carrier protein reductase (InhA). Compounds having bulky aromatic substituents at position 5 and a tryptophan residue at position N-3 of the rhodanine ring were the most active against InhA, with IC50 values ranging from 2.7 to 30 µM. The experimental data showed consistent correlations with computational studies. Their antimicrobial activity was assessed against Mycobacterium marinum (Mm) (a model for Mtb), Pseudomonas aeruginosa (Pa), Legionella pneumophila (Lp), and Enterococcus faecalis (Ef) by using anti-infective, antivirulence, and antibiotic assays. Nineteen out of 34 compounds reduced Mm virulence at 10 µM. 33 exhibited promising antibiotic activity against Mm with a MIC of 0.21 µM and showed up to 89% reduction of Lp growth in an anti-infective assay at 30 µM. 32 showed high antibiotic activity against Ef, with a MIC of 0.57 µM.

HIV-1 Protease Dimerization Dynamics Reveals a Transient Druggable Binding Pocket at the Interface.

The binding mechanism of HIV-1 protease monomers leading to the catalytically competent dimeric enzyme has been investigated by means of state-of-the-art atomistic simulations. The emerging picture allows a deeper understanding of experimental observations and reveals that water molecules trapped at the interface have an important role in slowing down the kinetics of the association process. Unexpectedly, a cryptic binding pocket is identified at the interface of the complex, corresponding to a partially bound dimer that lacks enzymatic function. The pocket has a transient nature with a lifetime longer than 1 µs, and it displays very favorable druggability features. Docking as well as MM-GBSA free-energy calculations further support the possibility to target the new binding site by means of inhibitors able to prevent the complete dimerization by capturing the inactive conformation. This discovery could open the way to the rational design of a new class of anti-HIV drugs.

Adrenergic antagonists restrict replication of Legionella.

Legionella pneumophila is a facultative intracellular bacterium, which upon inhalation can cause a potentially fatal pneumonia termed Legionnaires' disease. The opportunistic pathogen grows in environmental amoebae and mammalian macrophages within a unique membrane-bound compartment, the 'Legionella-containing vacuole'. Bacteria are exposed to many environmental cues including small signalling molecules from eukaryotic cells. A number of pathogenic bacteria sense and respond to catecholamine hormones, such as adrenalin and noradrenalin, a process mediated via the QseBC two-component system in some bacteria. In this study, we examined the effect of adrenergic compounds on L. pneumophila, and discovered that the adrenergic receptor antagonists benoxathian, naftopidil, propranolol and labetalol, as well as the QseC sensor kinase inhibitor LED209, reduced the growth of L. pneumophila in broth or amoebae, while replication in macrophages was enhanced. Growth restriction was common to members of the genus Legionella and Mycobacterium, and was observed for L. pneumophila in the replicative but not stationary phase of the biphasic life cycle. Deletion of the L. pneumophila qseBC genes indicated that growth inhibition by adrenergics or LED209 is mediated only to a minor extent by this two-component system, implying the presence of other adrenergic sensing systems. This study identifies adrenergic molecules as novel inhibitors of extra- and intracellular growth of Legionella and reveals LED209 as a potential lead compound to combat infections with Legionella or Mycobacterium spp.

Amoebae-Based Screening Reveals a Novel Family of Compounds Restricting Intracellular Legionella pneumophila.

The causative agent of Legionnaires' disease, Legionella pneumophila, grows in environmental amoebae and mammalian macrophages within a distinct compartment, the 'Legionella-containing vacuole' (LCV). Intracellular bacteria are protected from many antibiotics, and thus are notoriously difficult to eradicate. To identify novel compounds that restrict intracellular bacterial replication, we previously developed an assay based on a coculture of amoebae and GFP-producing L. pneumophila. This assay was used to screen a pathway-based, highly diverse chemical library, referred to as the Sinergia library. In this work, we chose to focus on a group of 11 hit compounds, the majority of which originated from the query molecule CN585, a compound that targets the protein phosphatase calcineurin. Further studies on 78 related compound variants revealed crucial structural attributes, namely a triple-ring scaffold with a central triazine moiety, substituted in positions 3 and 5 by two piperidine or pyrrolidine rings, and in position 1 by an amine group bearing a single aliphatic chain moiety. The most effective compound, ZINC00615682, inhibited intracellular replication of L. pneumophila with an IC50 of approximately 20 nM in Acanthamoeba castellanii and slightly less efficiently in Dictyostelium discoideum or macrophages. Pharmacological and genetic attempts to implicate calcineurin in the intracellular replication of L. pneumophila failed. Taken together, these results show that the amoebae-based screen and structure-activity relationship analysis is suitable for the identification of novel inhibitors of the intracellular replication of L. pneumophila. The most potent compound identified in this study targets (an) as yet unidentified host factor(s).

Modulation of alpha-synuclein aggregation by dopamine analogs.

The action of dopamine on the aggregation of the unstructured alpha-synuclein (alpha-syn) protein may be linked to the pathogenesis of Parkinson's disease. Dopamine and its oxidation derivatives may inhibit alpha-syn aggregation by non-covalent binding. Exploiting this fact, we applied an integrated computational and experimental approach to find alternative ligands that might modulate the fibrillization of alpha-syn. Ligands structurally and electrostatically similar to dopamine were screened from an established library. Five analogs were selected for in vitro experimentation from the similarity ranked list of analogs. Molecular dynamics simulations showed they were, like dopamine, binding non-covalently to alpha-syn and, although much weaker than dopamine, they shared some of its binding properties. In vitro fibrillization assays were performed on these five dopamine analogs. Consistent with our predictions, analyses by atomic force and transmission electron microscopy revealed that all of the selected ligands affected the aggregation process, albeit to a varying and lesser extent than dopamine, used as the control ligand. The in silico/in vitro approach presented here emerges as a possible strategy for identifying ligands interfering with such a complex process as the fibrillization of an unstructured protein.

Regulation of bestrophins by Ca2+: a theoretical and experimental study.

Bestrophins are a recently discovered family of Cl(-) channels, for which no structural information is available. Some family members are activated by increased intracellular Ca2+ concentration. Bestrophins feature a well conserved Asp-rich tract in their COOH terminus (Asp-rich domain), which is homologous to Ca2+-binding motifs in human thrombospondins and in human big-conductance Ca2+- and voltage-gated K+ channels (BK(Ca)). Consequently, the Asp-rich domain is also a candidate for Ca2+ binding in bestrophins. Based on these considerations, we constructed homology models of human bestrophin-1 (Best1) Asp-rich domain using human thrombospondin-1 X-ray structure as a template. Molecular dynamics simulations were used to identify Asp and Glu residues binding Ca2+ and to predict the effects of their mutations to alanine. We then proceeded to test selected mutations in the Asp-rich domain of the highly homologous mouse bestrophin-2. The mutants expressed in HEK-293 cells were investigated by electrophysiological experiments using the whole-cell voltage-clamp technique. Based on our molecular modeling results, we predicted that Asp-rich domain has two defined binding sites and that D301A and D304A mutations may impact the binding of the metal ions. The experiments confirmed that these mutations do actually affect the function of the protein causing a large decrease in the Ca2+-activated Cl(-) current, fully consistent with our predictions. In addition, other studied mutations (E306A, D312A) did not decrease Ca2+-activated Cl(-) current in agreement with modeling results.

Development of antibody fragments for immunotherapy of prion diseases.

Prions are infectious proteins responsible for a group of fatal neurodegenerative diseases called TSEs (transmissible spongiform encephalopathies) or prion diseases. In mammals, prions reproduce themselves by recruiting the normal cellular protein PrP(C) and inducing its conversion into the disease-causing isoform denominated PrP(Sc). Recently, anti-prion antibodies have been shown to permanently cure prion-infected cells. However, the inability of full-length antibodies and proteins to cross the BBB (blood-brain barrier) hampers their use in the therapy of TSEs in vivo. Alternatively, brain delivery of prion-specific scFv (single-chain variable fragment) by AAV (adeno-associated virus) transfer delays the onset of the disease in infected mice, although protection is not complete. We investigated the anti-prion effects of a recombinant anti-PrP (D18) scFv by direct addition to scrapie-infected cell cultures or by infection with both lentivirus and AAV-transducing vectors. We show that recombinant anti-PrP scFv is able to reduce proteinase K-resistant PrP content in infected cells. In addition, we demonstrate that lentiviruses are more efficient than AAV in gene transfer of the anti-PrP scFv gene and in reducing PrP(Sc) content in infected neuronal cell lines. Finally, we have used a bioinformatic approach to construct a structural model of the D18scFv-PrP(C) complex. Interestingly, according to the docking results, Arg(PrP)(151) (Arg(151) from prion protein) is the key residue for the interactions with D18scFv, anchoring the PrP(C) to the cavity of the antibody. Taken together, these results indicate that combined passive and active immunotherapy targeting PrP might be promising strategies for therapeutic intervention in prion diseases.

Docking Ligands on Protein Surfaces: The Case Study of Prion Protein.

Molecular docking of ligands targeting proteins undergoing fibrillization in neurodegenerative diseases is difficult because of the lack of deep binding sites. Here we extend standard docking methods with free energy simulations in explicit solvent to address this issue in the context of the prion protein surface. We focus on a specific ligand (2-pyrrolidin-1-yl-N-[4-[4-(2-pyrrolidin-1-yl-acetylamino)-benzyl]-phenyl]-acetamide), which binds to the structured part of the protein as shown by NMR (Kuwata, K. et al. Proc Natl Acad Sci U.S.A. 2007, 104, 11921-11926). The calculated free energy of dissociation (7.8 ± 0.9 kcal/mol) is in good agreement with the value derived by the experimental dissociation constant (Kd = 3.9 µM, corresponding to ΔG(0) = -7.5 kcal/mol). Several binding poses are predicted, including the one reported previously. Our prediction is fully consistent with the presence of multiple binding sites, emerging from NMR measurements. Our molecular simulation-based approach emerges, therefore, as a useful tool to predict poses and affinities of ligand binding to protein surfaces.

Carcinogenicity of acrylamide: a computational study.

This paper reports a series of ab initio, density functional theory (DFT), and semiempirical molecular orbital (MO) calculations concerning the reaction between the ultimate carcinogen of acrylamide and guanine. Acrylamide--a product of the Maillard reaction--is present in a variety of fried and oven-cooked food. After intake, it is epoxidized by cytochrome P450 2E1 to yield the ultimate carcinogen--glycidamide. Effects of solvation were considered using the Langevin dipoles (LD) model of Florian and Warshel and the solvent reaction field (SCRF) model of Tomasi and co-workers. In silico activation free energies are in very good agreement with the experimental value of 22.8 kcal/mol. This agreement presents strong evidence in favor of the validity of the proposed S N2 reaction mechanism and points to the applicability of quantum chemical methods to studies of reactions associated with carcinogenesis. In addition, insignificant stereoselectivity of the studied reaction was predicted. Finally, the competing reaction of glycidamide with adenine was simulated, and the experimentally observed regioselectivity was successfully reproduced.

Guanine alkylation by ethylene oxide: calculation of chemical reactivity.

In this paper we report on calculations of the activation free energy for a chemical reaction between ethylene oxide and guanine. Ethylene oxide is biologically relevant per se and is also a model compound for numerous ultimate carcinogens. Calculations were performed on the medium-high ab initio, DFT, and semiempirical MO levels. Effects of solvation were considered using the Langevine dipole method and solvent reaction field method of Tomasi and co-workers. The calculated activation free energies are in reasonable agreement with the experimental value.