Structural and dynamic roles of permanent water molecules in ligand molecular recognition by chicken liver bile acid binding protein

An Efficient Implementation of the Nwat-MMGBSA Method to Rescore Docking Results in Medium-Throughput Virtual Screenings

Nwat-MMGBSA is a variant of MM-PB/GBSA based on the inclusion of a number of explicit water molecules that are the closest to the ligand in each frame of a molecular dynamics trajectory. This method demonstrated improved correlations between calculated and experimental binding energies in both protein-protein interactions and ligand-receptor complexes, in comparison to the standard MM-GBSA. A protocol optimization, aimed to maximize efficacy and efficiency, is discussed here considering penicillopepsin, HIV1-protease, and BCL-XL as test cases. Calculations were performed in triplicates on both classic HPC environments and on standard workstations equipped by a GPU card, evidencing no statistical differences in the results. No relevant differences in correlation to experiments were also observed when performing Nwat-MMGBSA calculations on 4 or 1 ns long trajectories. A fully automatic workflow for structure-based virtual screening, performing from library set-up to docking and Nwat-MMGBSA rescoring, has then been developed. The protocol has been tested against no rescoring or standard MM-GBSA rescoring within a retrospective virtual screening of inhibitors of AmpC β-lactamase and of the Rac1-Tiam1 protein-protein interaction. In both cases, Nwat-MMGBSA rescoring provided a statistically significant increase in the ROC AUCs of between 20 and 30%, compared to docking scoring or to standard MM-GBSA rescoring.

All-Purpose Containers? Lipid-Binding Protein - Drug Interactions

The combined use of in vitro (¹⁹F-NMR) and in silico (molecular docking) procedures demonstrates the affinity of a number of human calycins (lipid-binding proteins from ileum, liver, heart, adipose tissue and epidermis, and retinol-binding protein from intestine) for different drugs (mainly steroids and vastatins). Comparative evaluations on the complexes outline some of the features relevant for interaction (non-polar character of the drugs; amino acids and water molecules in the protein calyx most often involved in binding). Dissociation constants (K_i) for drugs typically lie in the same range as K_i for natural ligands; in most instances (different proteins and docking conditions), vastatins are the strongest interactors, with atorvastatin ranking top in half of the cases. The affinity of some calycins for some of the vastatins is in the order of magnitude of the drug C_max after systemic administration in humans. The possible biological implications of this feature are discussed in connection with drug delivery parameters (route of administration, binding to carrier proteins, distribution to, and accumulation in, human tissues).

A conserved loop in polynucleotide phosphorylase (PNPase) essential for both RNA and ADP/phosphate binding

Polynucleotide phosphorylase (PNPase) reversibly catalyzes RNA phosphorolysis and polymerization of nucleoside diphosphates. Its homotrimeric structure forms a central channel where RNA is accommodated. Each protomer core is formed by two paralogous RNase PH domains: PNPase1, whose function is largely unknown, hosts a conserved FFRR loop interacting with RNA, whereas PNPase2 bears the putative catalytic site, ∼20 Å away from the FFRR loop. To date, little is known regarding PNPase catalytic mechanism. We analyzed the kinetic properties of two Escherichia coli PNPase mutants in the FFRR loop (R79A and R80A), which exhibited a dramatic increase in Km for ADP/Pi binding, but not for poly(A), suggesting that the two residues may be essential for binding ADP and Pi. However, both mutants were severely impaired in shifting RNA electrophoretic mobility, implying that the two arginines contribute also to RNA binding. Additional interactions between RNA and other PNPase domains (such as KH and S1) may preserve the enzymatic activity in R79A and R80A mutants. Inspection of enzyme structure showed that PNPase has evolved a long-range acting hydrogen bonding network that connects the FFRR loop with the catalytic site via the F380 residue. This hypothesis was supported by mutation analysis. Phylogenetic analysis of PNPase domains and RNase PH suggests that such network is a unique feature of PNPase1 domain, which coevolved with the paralogous PNPase2 domain.

The simulation approach to lipid-protein interactions

The interactions between lipids and proteins are crucial for a range of biological processes, from the folding and stability of membrane proteins to signaling and metabolism facilitated by lipid-binding proteins. However, high-resolution structural details concerning functional lipid/protein interactions are scarce due to barriers in both experimental isolation of native lipid-bound complexes and subsequent biophysical characterization. The molecular dynamics (MD) simulation approach provides a means to complement available structural data, yielding dynamic, structural, and thermodynamic data for a protein embedded within a physiologically realistic, modelled lipid environment. In this chapter, we provide a guide to current methods for setting up and running simulations of membrane proteins and soluble, lipid-binding proteins, using standard atomistically detailed representations, as well as simplified, coarse-grained models. In addition, we outline recent studies that illustrate the power of the simulation approach in the context of biologically relevant lipid/protein interactions.

In silico identification of new ligands for GPR17: a promising therapeutic target for neurodegenerative diseases

GPR17, a previously orphan receptor responding to both uracil nucleotides and cysteinyl-leukotrienes, has been proposed as a novel promising target for human neurodegenerative diseases. Here, in order to specifically identify novel potent ligands of GPR17, we first modeled in silico the receptor by using a multiple template approach, in which extracellular loops of the receptor, quite complex to treat, were modeled making reference to the most similar parts of all the class-A GPCRs crystallized so far. A high-throughput virtual screening exploration of GPR17 binding site with more than 130,000 lead-like compounds was then applied, followed by the wet functional and pharmacological validation of the top-scoring chemical structures. This approach revealed successful for the proposed aim, and allowed us to identify five agonists or partial agonists with very diverse chemical structure. None of these compounds could have been expected 'a priori' to act on a GPCR, and all of them behaved as much more potent ligands than GPR17 endogenous activators.

Insights from modeling the 3D structure of New Delhi metallo-β-lactamse and its binding interactions with antibiotic drugs

New Delhi metallo-beta-lactamase (NDM-1) is an enzyme that makes bacteria resistant to a broad range of beta-lactam antibiotic drugs. This is because it can inactivate most beta-lactam antibiotic drugs by hydrolyzing them. For in-depth understanding of the hydrolysis mechanism, the three-dimensional structure of NDM-1 was developed. With such a structural frame, two enzyme-ligand complexes were derived by respectively docking Imipenem and Meropenem (two typical beta-lactam antibiotic drugs) to the NDM-1 receptor. It was revealed from the NDM-1/Imipenem complex that the antibiotic drug was hydrolyzed while sitting in a binding pocket of NDM-1 formed by nine residues. And for the case of NDM-1/Meropenem complex, the antibiotic drug was hydrolyzed in a binding pocket formed by twelve residues. All these constituent residues of the two binding pockets were explicitly defined and graphically labeled. It is anticipated that the findings reported here may provide useful insights for developing new antibiotic drugs to overcome the resistance problem.

Insights from investigating the interactions of adamantane-based drugs with the M2 proton channel from the H1N1 swine virus

The M2 proton channel is one of indispensable components for the influenza A virus that plays a vital role in its life cycle and hence is an important target for drug design against the virus. In view of this, the three-dimensional structure of the H1N1-M2 channel was developed based on the primary sequence taken from a patient recently infected by the H1N1 (swine flu) virus. With an explicit water-membrane environment, molecular docking studies were performed for amantadine and rimantadine, the two commercial drugs generally used to treat influenza A infection. It was found that their binding affinity to the H1N1-M2 channel is significantly lower than that to the H5N1-M2 channel, fully consistent with the recent report that the H1N1 swine virus was resistant to the two drugs. The findings and the relevant analysis reported here might provide useful structural insights for developing effective drugs against the new swine flu virus.