Parallel‐ProBiS: Fast parallel algorithm for local structural comparison of protein structures and binding sites

The ProBiS algorithm performs a local structural comparison of the query protein surface against the nonredundant database of protein structures. It finds proteins that have binding sites in common with the query protein. Here, we present a new parallelized algorithm, Parallel‐ProBiS, for detecting similar binding sites on clusters of computers. The obtained speedups of the parallel ProBiS scale almost ideally with the number of computing cores up to about 64 computing cores. Scaling is better for larger than for smaller query proteins. For a protein with almost 600 amino acids, the maximum speedup of 180 was achieved on two interconnected clusters with 248 computing cores. Source code of Parallel‐ProBiS is available for download free for academic users at http://probis.cmm.ki.si/download . © 2012 Wiley Periodicals, Inc.     

Conformational analysis and the binding sites of nitrilotriacetamide: A computational study

The conformational analysis of nitrilotriacetamide has been carried out computationally, at both the semi‐empirical AM1 and density functional theory (DFT) (B3LYP/6‐31+G*) levels of theory. The lowest‐energy conformation predicted with the Monte Carlo search method, using the AM1 model, has two amide functionalities aligned on the same side; however, the DFT calculations at B3LYP/6‐31+G* predicted the global minimum with all three acetamide functionalities on the same side in the gas phase. In the aqueous phase, the DFT results predicted the orientations of amides similar to that of the reported crystal structure. The rotation barriers to transition to different low‐energy conformers of nitrilotriacetamide are lower in energy (5.0 kcal/mol) in water. The molecular electrostatic isopotentials (MESP) generated for the selected conformers at DFT level show that the nitrilotriacetamide could interact more effectively with the sodium chloride surface than that of its monomeric unit nitrilomonoacetamide. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

TargetATPsite: A template‐free method for ATP‐binding sites prediction with residue evolution image sparse representation and classifier ensemble

Understanding the interactions between proteins and ligands is critical for protein function annotations and drug discovery. We report a new sequence‐based template‐free predictor (TargetATPsite) to identify the Adenosine‐5′‐triphosphate (ATP) binding sites with machine‐learning approaches. Two steps are implemented in TargetATPsite: binding residues and pockets predictions, respectively. To predict the binding residues, a novel image sparse representation technique is proposed to encode residue evolution information treated as the input features. An ensemble classifier constructed based on support vector machines (SVM) from multiple random under‐samplings is used as the prediction model, which is effective for dealing with imbalance phenomenon between the positive and negative training samples. Compared with the existing ATP‐specific sequence‐based predictors, TargetATPsite is featured by the second step of possessing the capability of further identifying the binding pockets from the predicted binding residues through a spatial clustering algorithm. Experimental results on three benchmark datasets demonstrate the efficacy of TargetATPsite. © 2013 Wiley Periodicals, Inc.     

Insights into the N-terminal Cu(II) and Cu(I) binding sites of the human copper transporter CTR1

Abstract

Copper transporter 1 (CTR1) is the main copper transporter in the eukaryotic system. CTR1 has several important roles: It binds Cu(II) ions that are present in the blood; it reduces those Cu(II) ions to Cu(I); and it subsequently transfers Cu(I) to the cytoplasmic domain, where the ion is delivered to various cellular pathways. Here, we seek to identify CTR1 binding sites for Cu(II) and Cu(I) and to shed light on the Cu(II)-to-Cu(I) reduction process. We focus on the first 14 amino acids of CTR1. This N-terminal segment is rich with histidine and methionine residues, which are known to bind Cu(II) and Cu(I), respectively; thus, this region has been suggested to have an important function in recruiting Cu(II) and reducing it to Cu(I). We utilize electron paramagnetic resonance (EPR) spectroscopy together with nuclear magnetic resonance (NMR) and UV-VIS spectroscopy and alanine substitution to reveal Cu(II) and Cu(I) binding sites in the focal 14-amino-acid segment. We show that H5 and H6 directly coordinate to Cu(II), whereas M7, M9, and M12 are involved in Cu(I) binding. This research is another step on the way to a complete understanding of the cellular copper regulation mechanism in humans.     

Interaction of Cu(+) and Cu(2+) ions with alpha-alanine. A density functional study

The Cu(+) and Cu(2+) preferred binding sites on alpha-alanine and their affinity values for this amino acid were determined at the density functional level using three different hybrid exchange correlation potentials and the 6-311++G** basis set. The results demonstrated that the two ions both give stable complexes with alpha-alanine but the stability order of the metalated species and the coordination sites are different depending on the nature of the cation. In particular, the Cu(+)-alpha-alanine ground-state structure is characterized by an eta(2)-N,O coordination with the nitrogen and oxygen atoms belonging to the amino and carbonyl groups, respectively. In contrast, the most stable complex of the Cu(2+)-alpha-alanine system has an eta(2)-O,O coordination with the cation bonded to the -CO(2) (-) moiety of the zwitterionic form of the amino acid. Comparison with the Cu(+) and Cu(2+) affinity values for glycine, computed at the same levels of theory, demonstrated that the relative values do not change significantly as different hybrid functionals are used, although the absolute affinities are strongly influenced by the choice of the hybrid potential.     

N-glycosylation of High Mobility Group Box 1 protein (HMGB1) modulates the interaction with glycyrrhizin: A molecular modeling study

High Mobility Group Box 1 protein (HMGB1) is an abundant protein with multiple functions in cells, acting as a DNA chaperone and damage-associated molecular pattern molecule. It represents an attractive target for the treatment of inflammatory diseases and cancers. The plant natural product glycyrrhizin (GLR) is a well-characterized ligand of HMGB1 and a drug used to treat diverse liver and skin diseases. The drug is known to bind to each of the two adjacent HMG boxes of the non-glycosylated protein. In cells, HMGB1 is N-glycosylated at three asparagine residues located in boxes A and B, and these N-glycans are essential for the nucleocytoplasmic transport of the protein. But the impact of the N-glycans on drug binding is unknown. Here we have investigated the effect of the N-glycosylation of HMGB1 on its interaction with GLR using molecular modelling, after incorporation of three N-glycans on a Human HMGB1 structure (PDB code 2YRQ). Sialylated bi-antennary N-glycans were introduced on the protein and exposed in a folded or an extended conformation for the drug binding study. The docking of the drug was performed using both 18α- and 18β-epimers of GLR and the conformations and potential energy of interaction (ΔE) of the different drug-protein complexes were compared. The N-glycans do not shield the drug binding sites on boxes A and B but can modulate the drug-protein interaction, via both direct and indirect effects. The calculations indicate that binding of 18α/β-GLR to the HMG box is generally reduced when the protein is N-glycosylated vs. the non-glycosylated protein. In particular, the N-glycans in an extended configuration significantly weaken the binding of GLR to box-B. The effects of the N-glycans are mostly indirect, but in one case a direct contact with the drug, via a carbohydrate-carbohydrate interaction, was observed with 18β-GLR bound to Box-B of glycosylated HMGB1. For the first time, it is shown (at least in silico) that N-glycosylation, one of the many post-translational modifications of HMGB1, can affect drug binding.     

A computational model of the 5-HT3 receptor extracellular domain: search for ligand binding sites

A three-dimensional model of the 5-HT₃ receptor extarcellular domain has been derived on the basis of the nicotinic acetylcholine receptor model recently published by Tsigelny et al. Maximum complementarity between the position and characteristics of mutated residues putatively involved in ligand interaction and the pharmacophoric elements derived by the indirect approach applied on several series of 5-HT₃ ligands have been exploited to gain insights into the ligand binding modalities and to speculate on the mechanistic role of the structural components. The analysis of the three-dimensional model allows one to distinguish among amino acids that exert key roles in ligand interactions, subunit architecture, receptor assembly and receptor dynamics. For some of these, alternative roles with respect to the ones hypothesized by experimentalists are assigned. Different binding modalities for agonists and antagonists are highlighted, and residues which probably play a role in the transduction of binding into a change in conformational state of the receptor are suggested. Received: 27 July 2000 / Accepted: 15 September 2000 / Published online: 21 December 2000     

N-glycosylation and ubiquitinylation of PD-L1 do not restrict interaction with BMS-202: A molecular modeling study

The Programmed cell Death protein-1/Ligand 1 (PD-1/L1) checkpoint is a major target in oncology. Monoclonal antibodies targeting PD-1 or PD-L1 are used to treat different types of solid tumors and lymphoma. PD-L1-binding small molecules are also actively searched. The lead compound is the biphenyl drug BMS-202 which stabilizes PD-L1 protein dimers and displays a potent antitumor activity in experimental models. Here we have investigated the effect of N-glycosylation (at N35, N192, N200 and N219) and mono-ubiquitination (at K178) of PD-L1 on the interaction with BMS-202 by molecular modeling. Two complementary tridimensional models of PD-L1, based on available crystallographic structures, were constructed with BMS-202 bound. The structures were glycosylated, with a fucosylated bi-antennary N-glycan and ubiquitinated. Model 1 refers to glycoPD-L1 bearing 16 N-glycans, with or without 4 ubiquitin residues. Model 2 presents 8 N-glycans and 2 ubiquitin residues. In both cases, BMS-202 was bound to the protein interface, stabilizing a PD-L1 dimer. The incorporation of the N-glycans or the ubiquitins did not significantly alter the drug-protein recognition. The interface of the drug-stabilized protein dimer is unaffected by the glycosylation or ubiquitination. Calculations of the binding energies indicated that the glycosylation slightly reduces the stability of the drug-protein complexes but does not prevent the drug binding process. Our modeling study suggests that the drug can target efficiently the different forms of PD-L1 in cells, glycosylated, ubiquitinated or not. These models of N-glycosylated and ubiquitinated PD-L1 will be useful to study other PD-L1 protein complexes.     

Preliminary study on application of impregnated synthetic peptide TLC stationary phases for the pre-screening of 5-HT1A ligands

Chromatographic parameters (deltaR(f)), defined as a difference in the migration of tested compound on the control and peptide impregnated silica gel TLC plates, were determined for 42 arylpiperazine derivatives. An amino acid sequence of the peptide used for impregnation was derived from the III transmembrane segment of the 5-HT(1A) receptor in the close vicinity of aspartic acid (Asp 166) residue. It was found that the deltaR(f) values obtained in a model employing tetrapeptide P4LA (ADVL), as well as the calculated logP correlate with 5-HT(1A) receptor affinity of the studied compounds.     

Functional assessment of “in vivo” and “in silico” mutations in the guanine binding site of RNase T1: A DFT study

Experimentally, the functional assessment of amino acid side chains in proteins is carried out by comparing parameters such as binding constants for the wild‐type protein and a mutant protein in which the considered side chain is deleted. In the present study, we apply a density functional theory (DFT) methodology to obtain changes in binding energy upon mutations in the enzyme ribonuclease T₁. Mutant structures were either taken directly from crystallographic data (“in vivo”) allowing for conformational changes upon mutation, or derived from the wild‐type (“in silico”). Excluding entropic contributions, the computed interaction energy changes upon mutation in vivo correlate qualitatively well with experimental binding free energy changes. In contrast, the in silico approach does not perform as well, especially for residues that contribute largely to binding. Subsequently, we assessed the applicability of the in vivo approach by analyzing the functional cooperativity between pairs of side chains. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Improved 3D-QSAR prediction by multiple-conformational alignment: A case study on PTP1B inhibitors

Three-dimension quantitative structure activity relationship (3D-QSAR) was one of the major statistical techniques to investigate the correlation of biological activity with structural properties of candidate molecules, and the accuracy of statistic greatly depended on molecular alignment methodology. Exhaustive conformational search and successful conformational superposition could extremely improve the predictive accuracy of QSAR modeling. In this work, we proposed a solution to optimize QSAR prediction by multiple-conformational alignment methods, with a set of 40 flexible PTP1B inhibitors as case study. Three different molecular alignment methods were used for the development of 3D-QSAR models listed as following: (1) docking-based alignment (DBA); (2) pharmacophore-based alignment (PBA) and (3) co-crystallized conformer-based alignment (CCBA). Among these three alignments, it was indicated that the CCBA was the best and the fastest strategy in 3D-QSAR development, with the square correlation coefficient (r²) and cross-validated squared correlation coefficient (q²) of comparative molecular field analysis (CoMFA) were 0.992 and 0.694; the r² and q² of comparative molecular similarity indices analysis (CoMSIA) were 0.972 and 0.603, respectively. The alignment methodologies used here not only generated a robust QSAR model with useful molecular field contour maps for designing novel PTP1B inhibitors, but also provided a solution for constructing accurate 3D-QSAR model for various disease targets. Undoubtedly, such attempt in QSAR analysis would greatly help us to understand essential structural features of inhibitors required by its target, and so as to discover more promising chemical derivatives.     

Exploring the molecular basis of selectivity in A1 adenosine receptors agonists: a case study

Adenosine is a naturally occurring purine nucleoside that has a wide variety of well-documented regulatory functions and physiological roles. Selective activation of the adenosine A₁ receptor has drawn attention in drug discovery for the therapeutic effects on neural and cardiovascular disorders. We have developed a model of the human A₁ adenosine receptor using bovine rhodopsin as a template. A flexible docking approach has been subsequently carried out for evaluating the molecular interactions of twenty-one selective A₁ agonists with the receptor model. The results of these studies are consistent with mutational and biochemical data. In particular, they highlight a wide hydrogen-bonding network between the nucleoside portion of the ligands and the A₁ receptor as well as key amino acids for hydrophobic interactions with the different N⁶-groups of the agonists. The models presented here provide a detailed molecular map for the selective stimulation of the adenosine A₁ receptor subtype and a steady basis for the rational design of new A₁ selective ligands.     

Ab initio quantum mechanical study of the binding energies of human estrogen receptor alpha with its ligands: an application of fragment molecular orbital method

We have theoretically examined the relative binding affinities (RBA) of typical ligands, 17beta-estradiol (EST), 17alpha-estradiol (ESTA), genistein (GEN), raloxifene (RAL), 4-hydroxytamoxifen (OHT), tamoxifen (TAM), clomifene (CLO), 4-hydroxyclomifene (OHC), diethylstilbestrol (DES), bisphenol A (BISA), and bisphenol F (BISF), to the alpha-subtype of the human estrogen receptor ligand-binding domain (hERalpha LBD), by calculating their binding energies. The ab initio fragment molecular orbital (FMO) method, which we have recently proposed for the calculations of macromolecules such as proteins, was applied at the HF/STO-3G level. The receptor protein was primarily modeled by 50 amino acid residues surrounding the ligand. The number of atoms in these model complexes is about 850, including hydrogen atoms. For the complexes with EST, RAL, OHT, and DES, the binding energies were calculated again with the entire ERalphaLBD consisting of 241 residues or about 4000 atoms. No significant difference was found in the calculated binding energies between the model and the real protein complexes. This indicates that the binding between the protein and its ligands is well characterized by the model protein with the 50 residues. The calculated binding energies relative to EST were very well correlated with the experimental RBA (the correlation coefficient r=0.837) for the ligands studied in this work. We also found that the charge transfer between ER and ligands is significant on ER-ligand binding. To our knowledge, this is the first achievement of ab initio quantum mechanical calculations of large molecules such as the entire ERalphaLBD protein.     

A 1H NMR study of 1,4-bis(N-hexadecyl-N, N-dimethylammonium)butane dibromide/sodium anthranilate system: spherical to rod-shaped transition

The effect of addition of sodium anthranilate to 5 mM micellar solutions of gemini surfactant 1,4-bis(N-hexadecyl-N,N-dimethylammonium)butane dibromide is investigated by ¹H NMR. The solubilization site of anthranilate anion near the micellar surface is inferred. In the micelles, the An⁻ ions intercalate among the surfactant headgroups producing morphological changes.     

1-(2′-Pyridylazo)-2-naphtholate complexes of ruthenium: Synthesis,characterization, and DNA binding properties

Reaction of 1-(2′-pyridylazo)-2-naphthol (Hpan) with [Ru(dmso)₄Cl₂] (dmso = dimethylsulfoxide), [Ru(trpy)Cl₃] (trpy = 2,2′,2″-terpyridine), [Ru(bpy)Cl₃] (bpy = 2,2′-bipyridine) and [Ru(PPh₃)₃Cl₂] in refluxing ethanol in the presence of a base (NEt₃) affords, respectively, the [Ru(pan)₂], [Ru(trpy)(pan)]⁺ (isolated as perchlorate salt), [Ru(bpy)(pan)Cl] and [Ru(PPh₃)₂(pan)Cl] complexes. Structures of these four complexes have been determined by X-ray crystallography. In each of these complexes, the pan ligand is coordinated to the metal center as a monoanionic tridentate N,N,O-donor. Reaction of the [Ru(bpy)(pan)Cl] complex with pyridine (py) and 4-picoline (pic) in the presence of silver ion has yielded the [Ru(bpy)(pan)(py)]⁺ and [Ru(bpy)(pan)(pic)]⁺ complexes (isolated as perchlorate salts), respectively. All the complexes are diamagnetic (low-spin d⁶, S = 0) and show characteristic ¹H NMR signals and intense MLCT transitions in the visible region. Cyclic voltammetry on all the complexes shows a Ru(II)–Ru(III) oxidation on the positive side of SCE. Except in the [Ru(pan)₂] complex, a second oxidative response has been observed in the other five complexes. Reductions of the coordinated ligands have also been observed on the negative side of SCE. The [Ru(trpy)(pan)]ClO₄, [Ru(bpy)(pan)(py)]ClO₄ and [Ru(bpy)(pan)(pic)]ClO₄ complexes have been observed to bind to DNA, but they have not been able to cleave super-coiled DNA on UV irradiation.