A new protein-protein docking scoring function based on interface residue properties.

MOTIVATION: Protein-protein complexes are known to play key roles in many cellular processes. However, they are often not accessible to experimental study because of their low stability and difficulty to produce the proteins and assemble them in native conformation. Thus, docking algorithms have been developed to provide an in silico approach of the problem. A protein-protein docking procedure traditionally consists of two successive tasks: a search algorithm generates a large number of candidate solutions, and then a scoring function is used to rank them. RESULTS: To address the second step, we developed a scoring function based on a Voronoï tessellation of the protein three-dimensional structure. We showed that the Voronoï representation may be used to describe in a simplified but useful manner, the geometric and physico-chemical complementarities of two molecular surfaces. We measured a set of parameters on native protein-protein complexes and on decoys, and used them as attributes in several statistical learning procedures: a logistic function, Support Vector Machines (SVM), and a genetic algorithm. For the later, we used ROGER, a genetic algorithm designed to optimize the area under the receiver operating characteristics curve. To further test the scores derived with ROGER, we ranked models generated by two different docking algorithms on targets of a blind prediction experiment, improving in almost all cases the rank of native-like solutions. AVAILABILITY: http://genomics.eu.org/spip/-Bioinformatics-tools-

Encapsulation of mono- and oligo-nucleotides into aqueous-core nanocapsules in presence of various water-soluble polymers.

In a previous study, we have shown that cidofovir (CDV) and azidothymidine-triphosphate (AZT-TP) were poorly encapsulated in poly(iso-butylcyanoacrylate) (PIBCA) aqueous-core nanocapsules. This was attributed to the rapid leakage of these small and hydrophilic molecules through the thin polymer wall of the nanocapsules. In the present study, we have selected various water-soluble polymers as increasing Mw adjuvants and investigated their influence on the entrapment of mononucleotides (CDV, AZT-TP) as well as of oligonucleotides (ODN) into these PIBCA aqueous-core nanocapsules. We show here that the presence of cationic polymers (i.e. poly(ethyleneimine) (PEI) or chitosan) in the nanocapsule aqueous compartment allowed successful encapsulation of AZT-TP and ODN.

Encapsulation of antiviral nucleotide analogues azidothymidine-triphosphate and cidofovir in poly(iso-butylcyanoacrylate) nanocapsules.

Nucleoside analogues are widely used in the treatment of various viral infections. However, the poor in vivo conversion of the nucleoside analogues like azidothymidine (AZT) into their active triphosphate nucleotide counterpart limits their pharmacological efficacy. This could be overcome by the direct administration of azidothymidine triphosphate (AZT-TP), but it requires an appropriate drug delivery approach. Besides nucleoside analogues, nucleotide analogues like cidofovir (CDV) are also used in the treatment of viral infections. CDV has raised recent interest because of its promising activity against smallpox, but its use is limited by its poor bioavailability and nephrotoxicity. Here again, a proper drug delivery system should address these issues. In this study, we investigated the encapsulation of the nucleotide analogues AZT-TP and CDV into poly(iso-butylcyanoacrylate) aqueous core nanocapsules, known to efficiently entrap oligonucleotides. We show here that the encapsulation of these mono-nucleotides is less efficient than with oligonucleotides and that a rapid release of AZT-TP from the nanocapsules occurred in vitro. This highlights the importance of the molecular weight of the entrapped molecules which, if they are too small, are diffusing through the thin polymer membrane of the nanocapsules. On the other hand, a good protection of the encapsulated AZT-TP was observed.

Crystal structure of the S. cerevisiae D-ribose-5-phosphate isomerase: comparison with the archaeal and bacterial enzymes.

Ribose-5-phosphate isomerase A has an important role in sugar metabolism by interconverting ribose-5-phosphate and ribulose-5-phosphate. This enzyme is ubiquitous and highly conserved among the three kingdoms of life. We have solved the 2.1 A resolution crystal structure of the Saccharomyces cerevisiae enzyme by molecular replacement. This protein adopts the same fold as its archaeal and bacterial orthologs with two alpha/beta domains tightly packed together. Mapping of conserved residues at the surface of the protein reveals strong invariability of the active site pocket, suggesting a common ligand binding mode and a similar catalytic mechanism. The yeast enzyme associates as a homotetramer similarly to the archaeal protein. The effect of an inactivating mutation (Arg189 to Lys) is discussed in view of the information brought by this structure.

Adenosine phosphonoacetic acid is slowly metabolized by NDP kinase.

NDP kinase catalyzes the last step in the phosphorylation of nucleotides. It is also involved in the activation by cellular kinases of nucleoside analogs used in antiviral therapies. Adenosine phosphonoacetic acid, a close analog of ADP already proposed as an inhibitor of ribonucleotide reductase, was found to be a poor substrate for human NDP kinase, as well as a weak inhibitor with an equilibrium dissociation constant of 0.6 mM to be compared to 0.025 mM for ADP. The X-ray structure of a complex of adenosine phosphonoacetic acid and the NDP kinase from Dictyostelium was determined to 2.0 A resolution showing that the analog adopts a binding mode similar to ADP, but that no magnesium ion is present at the active site. As ACP may also interfere with other cellular kinases, its potential as a drug targeting NDP kinase or ribonucleotide reductase is likely to be limited due to strong side effects. The design of new molecules with a narrower specificity and a stronger affinity will benefit from the detailed knowledge of the complex ACP-NDP kinase.

Activation of anti-reverse transcriptase nucleotide analogs by nucleoside diphosphate kinase: improvement by alpha-boranophosphate substitution.

Nucleoside activation by nucleoside diphosphate kinase and inhibition of HIV-1 reverse transcriptase were studied comparatively for a new class of nucleoside analogs with a borano (BH3-) or a thio (SH) group on the alpha-phosphate. Both the alpha-Rp-borano derivatives of AZT and d4T improved phosphorylation by NDP kinase, inhibition of reverse transcription as well as stability of alpha-borano nonophosphate derivatives in terminated viral DNA chain.

High resolution crystal structures of T4 phage beta-glucosyltransferase: induced fit and effect of substrate and metal binding.

beta-Glucosyltransferase (BGT) is a DNA-modifying enzyme encoded by bacteriophage T4 that transfers glucose from uridine diphosphoglucose to 5-hydroxymethyl cytosine bases of phage T4 DNA. We report six X-ray structures of the substrate-free and the UDP-bound enzyme. Four also contain metal ions which activate the enzyme, including Mg(2+) in forms 1 and 2 and Mn(2+) or Ca(2+). The substrate-free BGT structure differs by a domain movement from one previously determined in another space group. Further domain movements are seen in the complex with UDP and the four UDP-metal complexes. Mg(2+), Mn(2+) and Ca(2+) bind near the beta-phosphate of the nucleotide, but they occupy slightly different positions and have different ligands depending on the metal and the crystal form. Whilst the metal site observed in these complexes with the product UDP is not compatible with a role in activating glucose transfer, it approximates the position of the positive charge in the oxocarbonium ion thought to form on the glucose moiety of the substrate during catalysis.

X-ray structure of HPr kinase: a bacterial protein kinase with a P-loop nucleotide-binding domain.

HPr kinase/phosphatase (HprK/P) is a key regulatory enzyme controlling carbon metabolism in Gram- positive bacteria. It catalyses the ATP-dependent phosphorylation of Ser46 in HPr, a protein of the phosphotransferase system, and also its dephosphorylation. HprK/P is unrelated to eukaryotic protein kinases, but contains the Walker motif A characteristic of nucleotide-binding proteins. We report here the X-ray structure of an active fragment of Lactobacillus casei HprK/P at 2.8 A resolution, solved by the multiwavelength anomalous dispersion method on a seleniated protein (PDB code 1jb1). The protein is a hexamer, with each subunit containing an ATP-binding domain similar to nucleoside/nucleotide kinases, and a putative HPr-binding domain unrelated to the substrate-binding domains of other kinases. The Walker motif A forms a typical P-loop which binds inorganic phosphate in the crystal. We modelled ATP binding by comparison with adenylate kinase, and designed a tentative model of the complex with HPr based on a docking simulation. The results confirm that HprK/P represents a new family of protein kinases, first identified in bacteria, but which may also have members in eukaryotes.

Standard atomic volumes in double-stranded DNA and packing in protein--DNA interfaces.

Standard volumes for atoms in double-stranded B-DNA are derived using high resolution crystal structures from the Nucleic Acid Database (NDB) and compared with corresponding values derived from crystal structures of small organic compounds in the Cambridge Structural Database (CSD). Two different methods are used to compute these volumes: the classical Voronoi method, which does not depend on the size of atoms, and the related Radical Planes method which does. Results show that atomic groups buried in the interior of double-stranded DNA are, on average, more tightly packed than in related small molecules in the CSD. The packing efficiency of DNA atoms at the interfaces of 25 high resolution protein-DNA complexes is determined by computing the ratios between the volumes of interfacial DNA atoms and the corresponding standard volumes. These ratios are found to be close to unity, indicating that the DNA atoms at protein-DNA interfaces are as closely packed as in crystals of B-DNA. Analogous volume ratios, computed for buried protein atoms, are also near unity, confirming our earlier conclusions that the packing efficiency of these atoms is similar to that in the protein interior. In addition, we examine the number, volume and solvent occupation of cavities located at the protein-DNA interfaces and compared them with those in the protein interior. Cavities are found to be ubiquitous in the interfaces as well as inside the protein moieties. The frequency of solvent occupation of cavities is however higher in the interfaces, indicating that those are more hydrated than protein interiors. Lastly, we compare our results with those obtained using two different measures of shape complementarity of the analysed interfaces, and find that the correlation between our volume ratios and these measures, as well as between the measures themselves, is weak. Our results indicate that a tightly packed environment made up of DNA, protein and solvent atoms plays a significant role in protein-DNA recognition.

Silk fibroin: structural implications of a remarkable amino acid sequence.

The amino acid sequence of the heavy chain of Bombyx mori silk fibroin was derived from the gene sequence. The 5,263-residue (391-kDa) polypeptide chain comprises 12 low-complexity "crystalline" domains made up of Gly-X repeats and covering 94% of the sequence; X is Ala in 65%, Ser in 23%, and Tyr in 9% of the repeats. The remainder includes a nonrepetitive 151-residue header sequence, 11 nearly identical copies of a 43-residue spacer sequence, and a 58-residue C-terminal sequence. The header sequence is homologous to the N-terminal sequence of other fibroins with a completely different crystalline region. In Bombyx mori, each crystalline domain is made up of subdomains of approximately 70 residues, which in most cases begin with repeats of the GAGAGS hexapeptide and terminate with the GAAS tetrapeptide. Within the subdomains, the Gly-X alternance is strict, which strongly supports the classic Pauling-Corey model, in which beta-sheets pack on each other in alternating layers of Gly/Gly and X/X contacts. When fitting the actual sequence to that model, we propose that each subdomain forms a beta-strand and each crystalline domain a two-layered beta-sandwich, and we suggest that the beta-sheets may be parallel, rather than antiparallel, as has been assumed up to now.

Structure of the globular region of the prion protein Ure2 from the yeast Saccharomyces cerevisiae.

BACKGROUND: The [URE3] non-Mendelian element of the yeast S. cerevisiae is due to the propagation of a transmissible form of the protein Ure2. The infectivity of Ure2p is thought to originate from a conformational change of the normal form of the prion protein. This conformational change generates a form of Ure2p that assembles into amyloid fibrils. Hence, knowledge of the three-dimensional structure of prion proteins such as Ure2p should help in understanding the mechanism of amyloid formation associated with a number of neurodegenerative diseases. RESULTS: Here we report the three-dimensional crystal structure of the globular region of Ure2p (residues 95--354), also called the functional region, solved at 2.5 A resolution by the MAD method. The structure of Ure2p 95--354 shows a two-domain protein forming a globular dimer. The N-terminal domain is composed of a central 4 strand beta sheet flanked by four alpha helices, two on each side. In contrast, the C-terminal domain is entirely alpha-helical. The fold of Ure2p 95--354 resembles that of the beta class glutathione S-transferases (GST), in line with a weak similarity in the amino acid sequence that exists between these proteins. Ure2p dimerizes as GST does and possesses a potential ligand binding site, although it lacks GST activity. CONCLUSIONS: The structure of the functional region of Ure2p is the first crystal structure of a prion protein. Structure comparisons between Ure2p 95--354 and GST identified a 32 amino acid residues cap region in Ure2p exposed to the solvent. The cap region is highly flexible and may interact with the N-terminal region of the partner subunit in the dimer. The implication of this interaction in the assembly of Ure2p into amyloid fibrils is discussed.

Binding of nucleotides to nucleoside diphosphate kinase: a calorimetric study.

The source of affinity for substrates of human nucleoside diphosphate (NDP) kinases is particularly important in that its knowledge could be used to design more effective antiviral nucleoside drugs (e.g., AZT). We carried out a microcalorimetric study of the binding of enzymes from two organisms to various nucleotides. Isothermal titration calorimetry has been used to characterize the binding in terms of Delta G degrees, Delta H degrees and Delta S degrees. Thermodynamic parameters of the interaction of ADP with the hexameric NDP kinase from Dictyostelium discoideum and with the tetrameric enzyme from Myxococcus xanthus, at 20 degrees C, were similar and, in both cases, binding was enthalpy-driven. The interactions of ADP, 2'deoxyADP, GDP, and IDP with the eukaryotic enzyme differed in enthalpic and entropic terms, whereas the Delta G degrees values obtained were similar due to enthalpy--entropy compensation. The binding of the enzyme to nonphysiological nucleotides, such as AMP--PNP, 3'deoxyADP, and 3'-deoxy-3'-amino-ADP, appears to differ in several respects. Crystallography of the protein bound to 3'-deoxy-3'-amino-ADP showed that the drug was in a distorted position, and was unable to interact correctly with active site side chains. The interaction of pyrimidine nucleoside diphosphates with the hexameric enzyme is characterized by a lower affinity than that with purine nucleotides. Titration showed the stoichiometry of the interaction to be abnormal, with 9--12 binding sites/hexamer. The presence of supplementary binding sites might have physiological implications.

Mechanism of phosphoryl transfer by nucleoside diphosphate kinase pH dependence and role of the active site Lys16 and Tyr56 residues.

Nucleoside diphosphate (NDP) kinase phosphorylates nucleoside diphosphates with little specificity for the base and the sugar. Although nucleotide analogues used in antiviral therapies are also metabolized to their triphosphate form by NDP kinase, their lack of the 3'-hydroxyl of the ribose, which allows them to be DNA chain terminators, severely impairs the catalytic efficiency of NDP kinase. We have analyzed the kinetics parameters of several mutant NDP kinases modified on residues (Lys16, Tyr56, Asn119) interacting with the gamma-phosphate and/or the 3'-OH of the Mg2+-ATP substrate. We compared the relative contributions of the active-site residues and the substrate 3'-OH for point mutations on Lys16, Tyr56 and Asn119. Analysis of additional data from pH profiles identify the ionization state of these residues in the enzyme active form. X-ray structure of K16A mutant NDP kinase shows no detectable rearrangement of the residues of the active site.

Chemical rescue of phosphoryl transfer in a cavity mutant: a cautionary tale for site-directed mutagenesis.

We have explored the ability of a nucleoside diphosphate kinase (NDPK) mutant in which the nucleophilic histidine has been replaced by glycine (H122G) to transfer phosphate from ATP to alcohols of varying pK(a), size, shape, and polarity. This cavity mutant does indeed act as a primitive alcohol kinase. The rate of its phosphoryl transfer to alcohols varies considerably, with values spanning a DeltaDeltaG(double dagger) range of 4 kcal/mol, whereas the alcohols have very similar intrinsic reactivities. Analysis of these results suggests that the ability to carry out phosphoryl transfer within the cavity is not a simple function of being small enough to enter the cavity, but rather is a complex function of steric, solvation, entropic, van der Waals packing, and electrostatic properties of the alcohol. In addition, large differences are observed between the reactivities of alcohols within the nucleophile cavity of H122G and the reactivities of the same alcohols within the nucleophile cavity of H122A, a mutant NDPK that differs from H122G by a single methyl group within the cavity. The crystal structures of the two cavity mutants are very similar to one another and to wild-type NDPK, providing no evidence for a structurally perturbed active site. The differences in reactivity between the two mutant proteins illustrate a fundamental limitation of energetic analysis from site-directed mutagenesis: although removal of a side chain is generally considered to be a conservative change, the energetic effects of any given mutation are inextricably linked to the molecular properties of the created cavity and the surrounding protein environment.

Structural basis for activation of alpha-boranophosphate nucleotide analogues targeting drug-resistant reverse transcriptase.

AIDS chemotherapy is limited by inadequate intracellular concentrations of the active triphosphate form of nucleoside analogues, leading to incomplete inhibition of viral replication and the appearance of drug-resistant virus. Drug activation by nucleoside diphosphate kinase and inhibition of HIV-1 reverse transcriptase were studied comparatively. We synthesized analogues with a borano (BH(3)(-)) group on the alpha-phosphate, and found that they are substrates for both enzymes. X-ray structures of complexes with nucleotide diphosphate kinase provided a structural basis for their activation. The complex with d4T triphosphate displayed an intramolecular CH.O bond contributing to catalysis, and the R(p) diastereoisomer of thymidine alpha-boranotriphosphate bound like a normal substrate. Using alpha-(R(p))-boranophosphate derivatives of the clinically relevant compounds AZT and d4T, the presence of the alpha-borano group improved both phosphorylation by nucleotide diphosphate kinase and inhibition of reverse transcription. Moreover, repair of blocked DNA chains by pyrophosphorolysis was reduced significantly in variant reverse transcriptases bearing substitutions found in drug-resistant viruses. Thus, the alpha-borano modification of analogues targeting reverse transcriptase may be of generic value in fighting viral drug resistance.

Fine organization of Bombyx mori fibroin heavy chain gene.

The complete sequence of the Bombyx mori fibroin gene has been determined by means of combining a shotgun sequencing strategy with physical map-based sequencing procedures. It consists of two exons (67 and 15 750 bp, respectively) and one intron (971 bp). The fibroin coding sequence presents a spectacular organization, with a highly repetitive and G-rich (approximately 45%) core flanked by non-repetitive 5' and 3' ends. This repetitive core is composed of alternate arrays of 12 repetitive and 11 amorphous domains. The sequences of the amorphous domains are evolutionarily conserved and the repetitive domains differ from each other in length by a variety of tandem repeats of subdomains of approximately 208 bp which are reminiscent of the repetitive nucleosome organization. A typical composition of a subdomain is a cluster of repetitive units, Ua, followed by a cluster of units, Ub, (with a Ua:Ub ratio of 2:1) flanked by conserved boundary elements at the 3' end. Moreover some repeats are also perfectly conserved at the peptide level indicating that the evolutionary pressure is not identical along the sequence. A tentative model for the constitution and evolution of this unusual gene is discussed.

The human nm23-H4 gene product is a mitochondrial nucleoside diphosphate kinase.

We demonstrate here the catalytic activity and subcellular localization of the Nm23-H4 protein, product of nm23-H4, a new member of the human nm23/nucleoside diphosphate (NDP) kinase gene family (Milon, L., Rousseau-Merck, M., Munier, A., Erent, M., Lascu, I., Capeau, J., and Lacombe, M. L. (1997) Hum. Genet. 99, 550-557). Nm3-H4 was synthesized in escherichia coli as the full-length protein and as a truncated form missing the N-terminal extension characteristic of mitochondrial targeting. The truncated form possesses NDP kinase activity, whereas the full-length protein is inactive, suggesting that the extension prevents enzyme folding and/or activity. X-ray crystallographic analysis was performed on active truncated Nm23-H4. Like other eukaryotic NDP kinases, it is a hexamer. Nm23-H4 naturally possesses a serine residue at position 129, equivalent to the K-pn mutation of the Drosophila NDP kinase. The x-ray structure shows that the presence of Ser(129) has local structural effects that weaken subunit interactions. Site-directed mutagenesis shows that the serine is responsible for the lability of Nm23-H4 to heat and urea treatment, because the S129P mutant is greatly stabilized. Examination of human embryonic kidney 293 cells transfected with green fluorescent protein fusions by confocal microscopy shows a specific mitochondrial localization of Nm23-H4 that was also demonstrated by Western blot analysis of subcellular fractions of these cells. Import into mitochondria is accompanied by cleavage of the N-terminal extension that results in NDP kinase activity. Submitochondrial fractionation indicates that Nm23-H4 is associated with mitochondrial membranes, possibly to the contact sites between the outer and inner membranes.

Three-dimensional structure of nucleoside diphosphate kinase.

Three-dimensional structures are known from X-ray studies of the nucleoside diphosphate (NDP) kinase of many organisms from bacteria to human. All NDP kinases have subunits of about 150 residues with a very similar fold based on the alphabeta sandwich or ferredoxin fold. This fold is found in many nucleotide or polynucleotide-binding proteins with no sequence relationship to NDP kinase. This common fold is augmented here with specific features: a surface alpha-helix hairpin, the Kpn loop, and the C-terminal extension. The alpha-helix hairpin and Kpn loop make up the nucleotide binding site, which is unique to NDP kinase and different from that of other kinases or ATPases. The Kpn loop and the C-terminal extension are also involved in the quaternary structure. Whereas all known eukaryotic NDP kinases, including mitochondral enzymes, are hexamers, some bacterial enzymes are tetramers. However, hexameric and tetrameric NDP kinases are built from the same dimer. The structural environment of the active histidine is identical in all. The nucleotide binding site is also fully conserved, except for a feature implicating C-terminal residues in the hexamer, but not in the tetramer. Structural data on the native and phosphorylated enzyme, complexes with substrates, inhibitor, and a transition state analog, give a solid basis to a mechanism of phosphate transfer in which the largest contributors to catalysis are the 3'-OH of the sugar and the bound Mg2+ in the nucleotide substrate. In contrast, we still lack structural data relating to DNA binding and other functions of NDP kinases.

p55-hGRF, a short natural form of the Ras-GDP exchange factor high yield production and characterization.

p55-hGRF, a natural short form of the guanine-nucleotide-releasing factor for p21-Ras from human brain, was expressed at high level in Escherichia coli as well as an engineered truncated form, p39-hGRF. A T7 polymerase expression system was used, resulting in the formation of insoluble cytoplasmic protein aggregates. The recombinant products were resolubilized, renatured and purified to homogeneity. The exchange activity of the refolded hGRF samples on H-Ras was comparable with that published for the soluble catalytic domain of the mouse counterpart, CDC25 Mm. Both p55-hGRF and p39-hGRF form dimers. We established a procedure to prepare and purify the complex with Ras. The results of the characterization study are consistent with a stoichiometry of 1:1 and an equilibrium between dimeric and monomeric forms of the complex.