Copper(II) Complexes of l-Arginine as Netropsin Mimics Showing DNA Cleavage Activity in Red Light

Copper(II) complexes [Cu(L-arg)(2)](NO3)(2) (1) and [Cu(L-arg)(B)Cl]Cl (2-5), where B is a heterocyclic base, namely, 2,2'-bipyridine (bpy, 2), 1,10-phenanthroline (phen, 3), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq, 4), and dipyrido[3,2-a:2',3'-c)phenazine (dppz, 5), are prepared and their DNA binding and photoinduced DNA cleavage activity studied. Ternary complex 3, structurally characterized using X-ray crystallography, shows a square-pyramidal (4 + 1) coordination geometry in which the N,O-donor L-arginine and N,N-donor 1,10-phenanthroline form the basal plane with one chloride at the elongated axial site. The complex has a pendant cationic guanidinium moiety. The one-electron paramagnetic complexes display a metal-centered d-d band in the range of 590-690 nm in aqueous DMF They show quasireversible cyclic voltammetric response due to the Cu(II)/Cu(I) couple in the range of -0.1 to -0.3 V versus a saturated calomel electrode in a DMF-Tris HCl buffer (pH 7.2). The DNA binding propensity of the complexes is studied using various techniques. Copper(II) bis-arginate 1 mimics the minor groove binder netropsin by showing preferential binding to the AT-rich sequence of double-strand (ds) DNA. DNA binding study using calf thymus DNA gives an order: 5 (L-arg-dppz) >= 1 (biS-L-arg) > 4 (L-arg-dpq) > 3 (L-arg-phen) >> 2 (L-arg-bpy). Molecular docking calculations reveal that the complexes bind through extensive hydrogen bonding and electrostatic interactions with ds-DNA. The complexes cleave supercoiled pUC19 DNA in the presence of 3-mercaptopropionic acid as a reducing agent forming hydroxyl ((OH)-O-center dot) radicals. The complexes show oxidative photoinduced DNA cleavage activity in UV-A light of 365 nm and red light of 647.1 nm (Ar-Kr mixed-gas-ion laser) in a metal-assisted photoexcitation process forming singlet oxygen (O-1(2)) species in a type-II pathway. All of the complexes, barring complex 2, show efficient DNA photocleavage activity. Complexes 4 and 5 exhibit significant double-strand breaks of DNA in red light of 647.1 nm due to the presence of two photosensitizers, namely, L-arginine and dpq or dppz in the molecules.

Photo-induced double-strand DNA and site-specific protein cleavage activity of L-histidine (mu-oxo)diiron(III) complexes of heterocyclic bases

Three oxo-bridged diiron(III) complexes of L-histidine and heterocyclic bases [Fe-2(mu-O)(L-his)(2)(B)(2)](ClO4)(2) (1-3), where B is 2,2'-bipyridine (bpy),1,10-phenanthroline (phen), dipyrido[3,2-d:2',3'-f]quinoxaline (dpq), were prepared and characterized. The bpy complex 1 was structurally characterized by X-ray crystallography. The molecular structure showed a {Fe-2(mu-O)} core in which iron(III) in a FeN4O2 coordination is bound to tridentate monoanionic L-histidine and bidentate bpy ligands. The Fe center dot center dot center dot Fe distance is similar to 3.5 angstrom. The Fe-O-Fe unit is essentially linear, giving a bond angle of similar to 172 degrees. The complexes showed irreversible cyclic voltammetric cathodic response near -0.1 V vs. SCE in H2O-0.1 M KCl. The binuclear units displayed antiferromagnetic interaction between two high-spin (S = 5/2) iron(III) centers giving a -J value of -110 cm(-1). The complexes showed good DNA binding propensity giving a binding constant value of similar to 10(5) M-1. Isothermal titration calorimetric data indicated single binding mode to the DNA. The binding was found to be driven by negative free energy change and enthalpy. The dpq complex 3 showed oxidative double-strand DNA cleavage on exposure to UV-A and visible light. The phen complex 2 displayed single-strand photocleavage of DNA. The DNA double-strand breaks were rationalized from theoretical molecular docking calculations. Mechanistic investigations showed formation of hydroxyl radicals as the reactive species through photodecarboxylation of the L-histidine ligand. The complexes exhibited good binding propensity to bovine serum albumin (BSA) protein in Tris-HCl/NaCl buffer medium. The dpq complex 3 showed UV-A light-induced site-specific oxidative BSA cleavage forming fragments of similar to 45 kDa and similar to 20 kDa molecular weights via SOH pathway.

Cross Talk between Receptor Guanylyl Cyclase C and c-src Tyrosine Kinase Regulates Colon Cancer Cell Cytostasis

Increased activation of c-src seen in colorectal cancer is an indicator of a poor clinical prognosis, suggesting that identification of downstream effectors of c-src may lead to new avenues of therapy. Guanylyl cyclase C (GC-C) is a receptor for the gastrointestinal hormones guanylin and uroguanylin and the bacterial heat-stable enterotoxin. Though activation of GC-C by its ligands elevates intracellular cyclic GMP (cGMP) levels and inhibits cell proliferation, its persistent expression in colorectal carcinomas and occult metastases makes it a marker for malignancy. We show here that GC-C is a substrate for inhibitory phosphorylation by c-src, resulting in reduced ligand-mediated cGMP production. Consequently, active c-src in colonic cells can overcome GC-C-mediated control of the cell cycle. Furthermore, docking of the c-src SH2 domain to phosphorylated GC-C results in colocalization and further activation of c-src. We therefore propose a novel feed-forward mechanism of activation of c-src that is induced by cross talk between a receptor GC and a tyrosine kinase. Our findings have important implications in understanding the molecular mechanisms involved in the progression and treatment of colorectal cancer.

Toward the Discovery of Vaccine Adjuvants: Coupling In Silico Screening and In Vitro Analysis of Antagonist Binding to Human and Mouse CCR4 Receptors

Background: Adjuvants enhance or modify an immune response that is made to an antigen. An antagonist of the chemokine CCR4 receptor can display adjuvant-like properties by diminishing the ability of CD4+CD25+ regulatory T cells (Tregs) to down-regulate immune responses. Methodology: Here, we have used protein modelling to create a plausible chemokine receptor model with the aim of using virtual screening to identify potential small molecule chemokine antagonists. A combination of homology modelling and molecular docking was used to create a model of the CCR4 receptor in order to investigate potential lead compounds that display antagonistic properties. Three-dimensional structure-based virtual screening of the CCR4 receptor identified 116 small molecules that were calculated to have a high affinity for the receptor; these were tested experimentally for CCR4 antagonism. Fifteen of these small molecules were shown to inhibit specifically CCR4-mediated cellmigration, including that of CCR4(+) Tregs. Significance: Our CCR4 antagonists act as adjuvants augmenting human T cell proliferation in an in vitro immune response model and compound SP50 increases T cell and antibody responses in vivo when combined with vaccine antigens of Mycobacterium tuberculosis and Plasmodium yoelii in mice.

The Role of UPF0157 in the Folding of M-tuberculosis Dephosphocoenzyme A Kinase and the Regulation of the Latter by CTP

Background: Targeting the biosynthetic pathway of Coenzyme A (CoA) for drug development will compromise multiple cellular functions of the tubercular pathogen simultaneously. Structural divergence in the organization of the penultimate and final enzymes of CoA biosynthesis in the host and pathogen and the differences in their regulation mark out the final enzyme, dephosphocoenzyme A kinase (CoaE) as a potential drug target. Methodology/Principal Findings: We report here a complete biochemical and biophysical characterization of the M. tuberculosis CoaE, an enzyme essential for the pathogen's survival, elucidating for the first time the interactions of a dephosphocoenzyme A kinase with its substrates, dephosphocoenzyme A and ATP; its product, CoA and an intrinsic yet novel inhibitor, CTP, which helps modulate the enzyme's kinetic capabilities providing interesting insights into the regulation of CoaE activity. We show that the mycobacterial enzyme is almost 21 times more catalytically proficient than its counterparts in other prokaryotes. ITC measurements illustrate that the enzyme follows an ordered mechanism of substrate addition with DCoA as the leading substrate and ATP following in tow. Kinetic and ITC experiments demonstrate that though CTP binds strongly to the enzyme, it is unable to participate in DCoA phosphorylation. We report that CTP actually inhibits the enzyme by decreasing its Vmax. Not surprisingly, a structural homology search for the modeled mycobacterial CoaE picks up cytidylmonophosphate kinases, deoxycytidine kinases, and cytidylate kinases as close homologs. Docking of DCoA and CTP to CoaE shows that both ligands bind at the same site, their interactions being stabilized by 26 and 28 hydrogen bonds respectively. We have also assigned a role for the universal Unknown Protein Family 0157 (UPF0157) domain in the mycobacterial CoaE in the proper folding of the full length enzyme. Conclusions/Significance: In view of the evidence presented, it is imperative to assign a greater role to the last enzyme of Coenzyme A biosynthesis in metabolite flow regulation through this critical biosynthetic pathway.

Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis

Background: Fatty acids are indispensable constituents of mycolic acids that impart toughness & permeability barrier to the cell envelope of M. tuberculosis. Biotin is an essential co-factor for acetyl-CoA carboxylase (ACC) the enzyme involved in the synthesis of malonyl-CoA, a committed precursor, needed for fatty acid synthesis. Biotin carboxyl carrier protein (BCCP) provides the co-factor for catalytic activity of ACC. Methodology/Principal Findings: BPL/BirA (Biotin Protein Ligase), and its substrate, biotin carboxyl carrier protein (BCCP) of Mycobacterium tuberculosis (Mt) were cloned and expressed in E. coli BL21. In contrast to EcBirA and PhBPL, the similar to 29.5 kDa MtBPL exists as a monomer in native, biotin and bio-5'AMP liganded forms. This was confirmed by molecular weigt profiling by gel filtration on Superdex S-200 and Dynamic Light Scattering (DLS). Computational docking of biotin and bio-5'AMP to MtBPL show that adenylation alters the contact residues for biotin. MtBPL forms 11 H-bonds with biotin, relative to 35 with bio-5'AMP. Docking simulations also suggest that bio-5'AMP hydrogen bonds to the conserved `GRGRRG' sequence but not biotin. The enzyme catalyzed transfer of biotin to BCCP was confirmed by incorporation of radioactive biotin and by Avidin blot. The K-m for BCCP was similar to 5.2 mu M and similar to 420 nM for biotin. MtBPL has low affinity (K-b = 1.06 x 10(-6) M) for biotin relative to EcBirA but their K-m are almost comparable suggesting that while the major function of MtBPL is biotinylation of BCCP, tight binding of biotin/bio-5'AMP by EcBirA is channeled for its repressor activity. Conclusions/Significance: These studies thus open up avenues for understanding the unique features of MtBPL and the role it plays in biotin utilization in M. tuberculosis.

Design, Synthesis, and DNA Binding Properties of Photoisomerizable Azobenzene−Distamycin Conjugates: An Experimental and Computational Study

Here, we present the synthesis, photochemical, and DNA binding properties of three photoisomerizable azobenzene−distamycin conjugates in which two distamycin units were linked via electron-rich alkoxy or electron-withdrawing carboxamido moieties with the azobenzene core. Like parent distamycin A, these molecules also demonstrated AT-specific DNA binding. Duplex DNA binding abilities of these conjugates were found to depend upon the nature and length of the spacer, the location of protonatable residues, and the isomeric state of the conjugate. The changes in the duplex DNA binding efficiency of the individual conjugates in the dark and with their respective photoirradiated forms were examined by circular dichroism, thermal denaturation of DNA, and Hoechst displacement assay with poly[d(A-T).d(T-A)] DNA in 150 mM NaCl buffer. Computational structural analyses of the uncomplexed ligands using ab initio HF and MP2 theory and molecular docking studies involving the conjugates with duplex d[(GC(AT)10CG)]2 DNA were performed to rationalize the nature of binding of these conjugates.

DNA cleavage in red light promoted by copper(II) complexes of -amino acids and photoactive phenanthroline bases

Ternary copper(II) complexes [Cu(L-trp)(B)(H2O)](NO3) ( 1–3) and [Cu(L-phe)(B)(H2O)](NO3) ( 4–6) of L-tryptophan (L-trp) and L-phenylalanine (L-phe) having phenanthroline bases (B), viz. 1,10-phenanthroline (phen, 1 and 4), dipyrido[3,2-d:2,3-f]quinoxaline (dpq, 2 and 5) and dipyrido[3,2-a:2,3-c]phenazine (dppz, 3 and 6), were prepared and characterized by physico-chemical techniques. Complexes 3 and 6 were structurally characterized by X-ray crystallography and show the presence of a square pyramidal (4 + 1) CuN3O2 coordination geometry in which the N,O-donor amino acid (L-trp or L-phe) and N,N-donor phenanthroline base bind at the equatorial plane with an aqua ligand coordinated at the elongated axial site. Complex 3 shows significant distortion from the square pyramidal geometry and a strong intramolecular – stacking interaction between the pendant indole ring of L-trp and the planar dppz aromatic moiety. All the complexes display good binding propensity to the calf thymus DNA giving an order: 3, 6 (dppz) > 2, 5 (dpq) > 1, 4 (phen). The binding constant (Kb) values are in the range of 2.1 × 104–1.1 × 106 mol-1 with the binding site size (s) values of 0.17–0.63. The phen and dpq complexes are minor groove binders while the dppz analogues bind at the DNA major groove. Theoretical DNA docking studies on 2 and 3 show the close proximity of two photosensitizers, viz. the indole moiety of L-trp and the quinoxaline/phenazine of the dpq/dppz bases, to the complementary DNA strands. Complexes 2 and 3 show oxidative DNA double strand breaks (dsb) of supercoiled (SC) DNA forming a significant quantity of linear DNA along with the nicked circular (NC) form on photoexposure to UV-A light of 365 nm and red light of 647.1 nm (Ar–Kr laser). Complexes 1, 5 and 6 show only single strand breaks (ssb) forming NC DNA. The red light induced DNA cleavage involves metal-assisted photosensitization of L-trp and dpq/dppz base resulting in the formation of a reactive singlet oxygen (1O2) species.

Functional interaction of diphenols with polyphenol oxidase - Molecular determinants of substrate/inhibitor specificity SO FEBS JOURNAL

Polyphenol oxidase (PPO) catalyzes the oxidation of o-diphenols to their respective quinones. The quinones autopolymerize to form dark pigments, an undesired effect. PPO is therefore the target for the development of antibrowning and antimelanization agents. A series of phenolic compounds experimentally evaluated for their binding affinity and inhibition constants were computationally docked to the active site of catechol oxidase. Docking studies suggested two distinct modes of binding, dividing the docked ligands into two groups. Remarkably, the first group corresponds to ligands determined to be substrates and the second group corresponds to reversible inhibitors. Analyses of the complexes provide structural explanations for correlating subtle changes in the position and nature of the substitutions on o-diphenols to their functional properties as substrates and inhibitors. Higher reaction rates and binding are reckoned by additional interactions of the substrates with key residues that line the hydrophobic cavity. The docking results suggest that inhibition of oxidation stems from an interaction between the aromatic carboxylic acid group and the apical His 109 of the four coordinates of the trigonal pyramidal coordination polyhedron of CuA. The spatial orientation of the hydroxyl in relation to the carboxylic group either allows a perfect fit in the substrate cavity, leading to inhibition, or because of a steric clash flips the molecule vertically, facilitating oxidation. This is the first study to explain, at the molecular level, the determinants Of substrate and inhibitor specificity of a catechol oxidase, thereby providing a platform for the design of selective inhibitors useful to both the food and pharmaceutical industries.

Parallel implementation of AutoDock

Computational docking of ligands to protein structures is a key step in structure-based drug design. Currently, the time required for each docking run is high and thus limits the use of docking in a high-throughput manner, warranting parallelization of docking algorithms. AutoDock, a widely used tool, has been chosen for parallelization. Near-linear increases in speed were observed with 96 processors, reducing the time required for docking ligands to HIV-protease from 81 min, as an example, on a single IBM Power-5 processor ( 1.65 GHz), to about 1 min on an IBM cluster, with 96 such processors. This implementation would make it feasible to perform virtual ligand screening using AutoDock.

Protein sequence design on the basis of topology optimization techniques -Using continuous modeling of discrete amino acid types

The notion of optimization is inherent in protein design. A long linear chain of twenty types of amino acid residues are known to fold to a 3-D conformation that minimizes the combined inter-residue energy interactions. There are two distinct protein design problems, viz. predicting the folded structure from a given sequence of amino acid monomers (folding problem) and determining a sequence for a given folded structure (inverse folding problem). These two problems have much similarity to engineering structural analysis and structural optimization problems respectively. In the folding problem, a protein chain with a given sequence folds to a conformation, called a native state, which has a unique global minimum energy value when compared to all other unfolded conformations. This involves a search in the conformation space. This is somewhat akin to the principle of minimum potential energy that determines the deformed static equilibrium configuration of an elastic structure of given topology, shape, and size that is subjected to certain boundary conditions. In the inverse-folding problem, one has to design a sequence with some objectives (having a specific feature of the folded structure, docking with another protein, etc.) and constraints (sequence being fixed in some portion, a particular composition of amino acid types, etc.) while obtaining a sequence that would fold to the desired conformation satisfying the criteria of folding. This requires a search in the sequence space. This is similar to structural optimization in the design-variable space wherein a certain feature of structural response is optimized subject to some constraints while satisfying the governing static or dynamic equilibrium equations. Based on this similarity, in this work we apply the topology optimization methods to protein design, discuss modeling issues and present some initial results.

Studies on active site mutants of P-falciparum adenylosuccinate synthetase: Insights into enzyme catalysis and activation

Adenylosuccinate synthetase catalyzes a reversible reaction utilizing IMP, GTP and aspartate in the presence of Mg2+ to form adenylosuccinate, GDP and inorganic phosphate. Comparison of similarly liganded complexes of Plasmodium falciparum, mouse and Escherichia coil AdSS reveals H-bonding interactions involving nonconserved catalytic loop residues (Asn429, Lys62 and Thr307) that are unique to the parasite enzyme. Site-directed mutagenesis has been used to examine the role of these interactions in catalysis and structural organization of P. falciparum adenylosuccinate synthetase (PfAdSS). Mutation of Asn429 to Val, Lys62 to Leu and Thr307 to Val resulted in an increase in K-m values for IMP, GTP and aspartate, respectively along with a 5 fold drop in the k(cat) value for N429V mutant suggesting the role of these residues in ligand binding and/or catalysis. We have earlier shown that the glycolytic intermediate, fructose 1,6 bisphosphate, which is an inhibitor of mammalian AdSS is an activator of the parasite enzyme. Enzyme kinetics along with molecular docking suggests a mechanism for activation wherein F16BP seems to be binding to the Asp loop and inducing a conformation that facilitates aspartate binding to the enzyme active site. Like in other AdSS, a conserved arginine residue (Arg155) is involved in dimer crosstalk and interacts with IMP in the active site of the symmetry related subunit of PfAdSS. We also report on the iochemical characterization of the arginine mutants (R155L, R155K and R155A) which suggests that unlike in E. coil AdSS, Arg155 in PfAdSS influences both ligand binding and catalysis. (C) 2010 Elsevier B.V. All rights reserved.

Crystal structures, dynamics and functional implications of molybdenum-cofactor biosynthesis protein MogA from two thermophilic organisms

Molybdenum-cofactor (Moco) biosynthesis is an evolutionarily conserved pathway in almost all kingdoms of life, including humans. Two proteins, MogA and MoeA, catalyze the last step of this pathway in bacteria, whereas a single two-domain protein carries out catalysis in eukaryotes. Here, three crystal structures of the Moco-biosynthesis protein MogA from the two thermophilic organisms Thermus thermophilus (TtMogA; 1.64 angstrom resolution, space group P2(1)) and Aquifex aeolicus (AaMogA; 1.70 angstrom resolution, space group P2(1) and 1.90 angstrom resolution, space group P1) have been determined. The functional roles and the residues involved in oligomerization of the protein molecules have been identified based on a comparative analysis of these structures with those of homologous proteins. Furthermore, functional roles have been proposed for the N- and C-terminal residues. In addition, a possible protein-protein complex of MogA and MoeA has been proposed and the residues involved in protein-protein interactions are discussed. Several invariant water molecules and those present at the subunit interfaces have been identified and their possible structural and/or functional roles are described in brief. In addition, molecular-dynamics and docking studies with several small molecules (including the substrate and the product) have been carried out in order to estimate their binding affinities towards AaMogA and TtMogA. The results obtained are further compared with those obtained for homologous eukaryotic proteins.

Modeling of angiogenin-3 '-NMP complex

Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the li,ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.

DEPTH: a web server to compute depth and predict small-molecule binding cavities in proteins

Depth measures the extent of atom/residue burial within a protein. It correlates with properties such as protein stability, hydrogen exchange rate, protein-protein interaction hot spots, post-translational modification sites and sequence variability. Our server, DEPTH, accurately computes depth and solvent-accessible surface area (SASA) values. We show that depth can be used to predict small molecule ligand binding cavities in proteins. Often, some of the residues lining a ligand binding cavity are both deep and solvent exposed. Using the depth-SASA pair values for a residue, its likelihood to form part of a small molecule binding cavity is estimated. The parameters of the method were calibrated over a training set of 900 high-resolution X-ray crystal structures of single-domain proteins bound to small molecules (molecular weight < 1.5 KDa). The prediction accuracy of DEPTH is comparable to that of other geometry-based prediction methods including LIGSITE, SURFNET and Pocket-Finder (all with Matthew's correlation coefficient of similar to 0.4) over a testing set of 225 single and multi-chain protein structures. Users have the option of tuning several parameters to detect cavities of different sizes, for example, geometrically flat binding sites. The input to the server is a protein 3D structure in PDB format. The users have the option of tuning the values of four parameters associated with the computation of residue depth and the prediction of binding cavities. The computed depths, SASA and binding cavity predictions are displayed in 2D plots and mapped onto 3D representations of the protein structure using Jmol. Links are provided to download the outputs. Our server is useful for all structural analysis based on residue depth and SASA, such as guiding site-directed mutagenesis experiments and small molecule docking exercises, in the context of protein functional annotation and drug discovery.