For Guldberg and Waage, with love and cratic entropy.
Explore the source record for details and available documents.
Entropy↗
Biomedical subjects
J Janin
Publications and source records attributed to J Janin.
Molecular systematics in 3D: the classification of protein structures. Paris, April 9-11, 1996.
Explore the source record for details and available documents.
Databases, Factual↗
Quantifying biological specificity: the statistical mechanics of molecular recognition.
The Random Energy Model of statistical physics is applied to the problem of the specificity of recognition between two biological (macro)molecules forming a non-covalent complex. In this model, the native mode of association is separated by an energy gap from a large body of non-native modes. Whereas the native mode is unique, the non-native modes form an energy spectrum which is approximated by a gaussian distribution. Specificity can then be estimated by writing the partition function and calculating the ratio r of non-native to native modes at thermodynamic equilibrium. We examine three situations: (i) recognition in the absence of a competitor; (ii) recognition in the presence of a competing ligand; (iii) recognition in a heterogeneous mixture. We derive the dependence of the ratio r on temperature and on the concentration of competing ligands, and we estimate the effect of a local perturbation such as can result from a point mutation. Cases (i) and (iii) are modeled by docking experiments in the computer. In case (iii), which is representative of a wide variety of biological situations, we show that increasing the heterogeneity of a mixture affects the specificity of recognition, even when the concentration of competing species is kept constant.
Binding Sites↗
Molecular docking programs successfully predict the binding of a beta-lactamase inhibitory protein to TEM-1 beta-lactamase.
Crystallization of the 1:1 molecular complex between the beta-lactamase TEM-1 and the beta-lactamase inhibitory protein BLIP has provided an opportunity to put a stringent test on current protein-docking algorithms. Prior to the successful determination of the structure of the complex, nine laboratory groups were given the refined atomic coordinates of each of the native molecules. Other than the fact that BLIP is an effective inhibitor of a number of beta-lactamase enzymes (KI for TEM-1 approximately 100 pM) no other biochemical or structural data were available to assist the practitioners in their molecular docking. In addition, it was not known whether the molecules underwent conformational changes upon association or whether the inhibition was competitive or non-competitive. All six of the groups that accepted the challenge correctly predicted the general mode of association of BLIP and TEM-1.
Amino Acid Sequence↗
X-ray structure of human nucleoside diphosphate kinase B complexed with GDP at 2 A resolution.
BACKGROUND: Nucleoside diphosphate (NDP) kinases provide precursors for DNA and RNA synthesis. In mammals, these enzymes are also involved in cell regulations. Human NDP kinase B, product of the human nm23-H2 gene, is both an enzyme and a transcription factor. It activates transcription of the c-myc oncogene independently of its catalytic function, by binding to its promoter DNA. How do the two functions coexist? RESULTS: Recombinant human NDP kinase B was co-crystallized with GDP. The X-ray structure was solved at 2.0 A resolution by molecular replacement from the homologous Drosophila Awd protein. Both enzymes are homo-hexamers with a characteristic beta alpha beta beta alpha beta fold. GDP binds near the active site His118. The guanine base is in a surface cleft and interacts with the C terminus of another subunit. CONCLUSIONS: The beta alpha beta beta alpha beta fold, also present in the 'palm' domain of Escherichia coli DNA polymerase I and HIV reverse transcriptase, is both a mononucleotide- and a polynucleotide-binding fold. If NDP kinase B binds DNA in the same way as the polymerases, the enzyme must undergo a conformation change in order to carry out gene activation.
Amino Acid Sequence↗
X-ray structure and catalytic mechanism of lobster enolase.
Enolase prepared from lobster tail muscle yielded trigonal crystals with one 47 kDa subunit per asymmetric unit. X-ray data were collected on the apoenzyme at 2.4 A resolution and on a complex with Mn2+ and the inhibitor phosphoglycolate at 2.2 A resolution. The corresponding cDNA was amplified from a library of lobster muscle cDNA, and a sequence corresponding to residues 27-398 was determined. It is highly homologous to other enolases, including yeast enolase for which an X-ray structure is available. Yeast enolase was used as a starting point for crystallographic refinement, which led to models of lobster enolase having R-factors below 22% and good stereochemistry. These models are very similar to yeast enolase; they have the same fold with a beta 3 alpha 4 N-terminal domain followed by an atypical alpha/beta barrel. Lobster apoenolase and the ternary complex differ only in the position of three mobile loops. In the complex, a single Mn2+ ion is seen ligated to three carboxylates and three water molecules. Phosphoglycolate binds near, but not directly to, the metal. His 157, which belongs to one of the mobile loops, is in contact with the C2 atom of the ligand. A water molecule hydrogen-bonds to the carboxylate of the ligand and to those of Glu 166 and Glu 209. We suggest that His 157 is the base that abstracts the C2H proton, whereas the water molecule is part of a proton relay system keeping the substrate in the carboxylic acid form where the pKa of the C2H group is low enough for proton transfer to His 157. The resulting catalytic mechanism is different from those proposed on the basis of the yeast enzyme X-ray structures, but it fits with earlier biochemical and spectroscopic data.
Amino Acid Sequence↗
Mechanism of phosphate transfer by nucleoside diphosphate kinase: X-ray structures of the phosphohistidine intermediate of the enzymes from Drosophila and Dictyostelium.
Nucleoside diphosphate kinase (NDP kinase) has a ping-pong mechanism with a phosphohistidine intermediate. Crystals of the enzymes from Dictyostelium discoideum and from Drosophila melanogaster were treated with phosphoramidate, and their X-ray structures were determined at 2.1 and 2.2 A resolution, respectively. The atomic models, refined to R factors below 20%, show no conformation change relative to the free proteins. In both enzymes, the active site histidine was phosphorylated on N delta, and it was the only site of phosphorylation. The phosphate group interacts with the hydroxyl group of Tyr56 and with protein-bound water molecules. Its environment is compared with that of phosphohistidines in succinyl-CoA synthetase and in phosphocarrier proteins. The X-ray structures of phosphorylated NDP kinase and of previously determined complexes with nucleoside diphosphates provide a basis for modeling the Michaelis complex with a nucleoside triphosphate, that of the phosphorylated protein with a nucleoside diphosphate, and the transition state of the phosphate transfer reaction in which the gamma-phosphate is pentacoordinated.
Animals↗
Phasing with mercury at 1 A wavelength.
Synchrotron sources provide a continuously tunable X-ray beam which makes it possible to optimize the anomalous contribution to phase determination using heavy-atom replacement. This method was used to solve two protein structures, those of Dictyostelium discoideum nucleoside diphosphate kinase and of lobster enolase. The first had 17 kDa of protein in the asymmetric unit, the second, 47 kDa. In both cases, a single mercury derivative yielded single isomorphous replacement with anomalous-scattering phases from which an interpretable electron-density map was derived by solvent flattening. The efficient solution of the X-ray structure was largely due to the large anomalous scattering of mercury at a wavelength shorter than the L(III) absorption edge.
Journal Article↗
Thermodynamics of the temperature-induced unfolding of globular proteins.
The heat capacity, enthalpy, entropy, and Gibbs energy changes for the temperature-induced unfolding of 11 globular proteins of known three-dimensional structure have been obtained by microcalorimetric measurements. Their experimental values are compared to those we calculate from the change in solvent-accessible surface area between the native proteins and the extended polypeptide chain. We use proportionality coefficients for the transfer (hydration) of aliphatic, aromatic, and polar groups from gas phase to aqueous solution, we estimate vibrational effects, and we discuss the temperature dependence of each constituent of the thermodynamic functions. At 25 degrees C, stabilization of the native state of a globular protein is largely due to two favorable terms: the entropy of non-polar group hydration and the enthalpy of interactions within the protein. They compensate the unfavorable entropy change associated with these interactions (conformational entropy) and with vibrational effects. Due to the large heat capacity of nonpolar group hydration, its stabilizing contribution decreases quickly at higher temperatures, and the two unfavorable entropy terms take over, leading to temperature-induced unfolding.
Calorimetry↗
Elusive affinities.
The affinity of two molecules for each other and its temperature dependence are determined by the change in enthalpy, free enthalpy, entropy, and heat capacity upon dissociation. As we know the forces that stabilize protein-protein or protein-DNA association and the three-dimensional structures of the complex, we can in principle derive values for each one of these parameters. The calculation is done first in gas phase by molecular mechanics, then in solution with the help of hydration parameters calibrated on small molecules. However, estimates of enthalpy and entropy changes in gas phase have excessively large error bars even under the approximation that the components of the complex associate as rigid bodies. No reliable result can be expected at the end. The fit to experimental values derived from binding and calorimetric measurements is poor, except for the dissociation heat capacity. This parameter can be attributed mostly to the hydration step and it correlates with the size of the interface. Many protein-protein complexes have interface areas in the range 1200-2000 A2 and only small conformation changes, so the rigid body approximation applies. It is less generally valid in protein-DNA complexes, which have interfaces covering 2200-3100 A2, large dissociation heat capacities, and affect both the conformation and the dynamics of their components.
DNA↗
Protein-protein interaction at crystal contacts.
Packing contacts are crystal artifacts, yet they make use of the same forces that govern specific recognition in protein-protein complexes and oligomeric proteins. They provide examples of a nonspecific protein-protein interaction which can be compared to biologically relevant ones. We evaluate the number and size of pairwise interfaces in 152 crystal forms where the asymmetric unit contains a monomeric protein. In those crystal forms that have no element of 2-fold symmetry, we find that molecules form 8 to 10 pairwise interfaces. The total area of the surface buried on each molecule is large, up to 4400 A2. Pairwise interfaces bury 200-1200 A2, like interfaces generated at random in a computer simulation, and less than interfaces in protease-inhibitor or antigen-antibody complexes, which bury 1500 A2 or more. Thus, specific contacts occurring in such complexes extend over a larger surface than nonspecific ones. In crystal forms with 2-fold symmetry, pairwise interfaces are fewer and larger on average than in the absence of 2-fold symmetry. Some bury 1500-2500 A2, like interfaces in oligomeric proteins, and create "crystal oligomers" which may have formed in the solution before crystallizing.
Adenylate Kinase↗
Principles of protein-protein recognition from structure to thermodynamics.
Specific recognition is illustrated by X-ray structures of protease-inhibitor, antigen-antibody and other high affinity complexes including five electron transfer complexes. We attempt to give a physical definition to affinity and specificity on the basis of these data. In a protein-protein complex, specific recognition results from the assembly of complementary surfaces into well-packed interfaces that cover about 1500 A2 and contain about ten hydrogen bonds. These interfaces are larger than between molecules in protein crystals, and smaller than between subunits in oligomeric proteins. We relate the size and chemical nature of interfaces in complexes to the thermodynamical parameters that characterize affinity: the heat capacity and free enthalpy (Gibbs energy) of dissociation at equilibrium, the activation free enthalpy for the dissociation reaction. The same structural and thermodynamical parameters are inadequate for representing the specificity of recognition. We propose instead to describe specificity with the help of statistical physics, and we illustrate the application of the random energy model to antigen-antibody recognition by analyzing results of computer simulations by docking.
Chemical Phenomena↗
Protein-protein recognition.
Explore the source record for details and available documents.
Bacterial Proteins↗
Refined X-ray structure of Dictyostelium discoideum nucleoside diphosphate kinase at 1.8 A resolution.
The X-ray structure of the nucleoside diphosphate kinase (NDP kinase) from Dictyostelium discoideum has been refined at 1.8 A resolution from a hexagonal crystal form with a 17 kDa monomer in its asymmetric unit. The atomic model was derived from the previously determined structure of a point mutant of the protein. It contains 150 amino acid residues out of 155, and 95 solvent molecules. The R-factor is 0.196 and the estimated accuracy of the average atomic position, 0.25 A. The Dictyostelium structure is described in detail and compared to those of Drosophila and Myxococcus xanthus NDP kinases. The protein is a hexamer with D3 symmetry. Residues 8 to 138 of each subunit form a globular alpha/beta domain. The four-stranded beta-sheet is antiparallel; its topology is different from other phosphate transfer enzymes, and also from the HPr protein which, like NDP kinase, carries a phosphorylated histidine. The same topology is nevertheless found in several other proteins that bind mononucleotides, RNA or DNA. Strand connections in NDP kinase involve alpha-helices and a 20-residue segment called the Kpn loop. The beta-sheet is regular except for a beta-bulge in edge strand beta 2 and a gamma-turn at residue Ile120 just preceding strand beta 4. The latter may induce strain in the main chain near the active site His122. The alpha 1 beta 2 motif participates in forming dimers within the hexamer, helices alpha 1 and alpha 3, the Kpn loop and C terminus, in forming trimers. The subunit fold and dimer interactions found in Dictyostelium are conserved in other NDP kinases. Trimer interactions probably occur in all eukaryotic enzymes. They are absent in the bacterial Myxococcus xanthus enzyme which is a tetramer, even though the subunit structure is very similar. In Dictyostelium, contacts between Kpn loops near the 3-fold axis block access to a central cavity lined with polar residues and filled with well-defined solvent molecules. Biochemical data on point mutants highlight the contribution of the Kpn loop to protein stability. In Myxococcus, the Kpn loops are on the tetramer surface and their sequence is poorly conserved. Yet, their conformation is maintained and they make a similar contribution to the substrate binding site.
Amino Acid Sequence↗
X-ray structure of nucleoside diphosphate kinase complexed with thymidine diphosphate and Mg2+ at 2-A resolution.
We report the crystal structure of nucleoside diphosphate kinase (NDP kinase) from Dictyostelium discoideum with thymidine diphosphate (dTDP) and Mg2+ bound at the active site. The structure has been refined to an R-factor of 18.3% at 2-A resolution. The base stacks on the aromatic ring of Phe 64 near the protein surface and is wedged between the side chains of Phe 64 and Val 116. The sugar and the pyrophosphate are deeper inside the protein and make numerous H-bonds with protein side chains. There is no backbone interaction with the nucleotide. A Mg2+ ion bridges the alpha- and beta-phosphates and interacts with the protein via water molecules. NDP kinase shows little specificity toward ribonucleotides and deoxyribonucleotides. This property, required by the enzyme biological function, can now be analyzed by comparing the crystal structures of free, ADP-ligated, and dTDP-ligated enzymes. The most significant differences are located in residues 60-64, which adapt their conformation to allow Phe 64 to stack on both types of bases. Nonspecific binding is achieved by the absence of polar interaction between the base and protein atoms. The ribose of ADP and the deoxyribose of dTDP occupy similar positions, their hydroxyl groups interacting with Lys 16 and Asn 119. The H-bond between Lys 16 and the O2' hydroxyl of ADP is replaced by a similar interaction with a water molecule in the dTDP complex. The beta-phosphate position is the same for ADP and dTDP, suggesting that the mechanism of phosphate transfer is the same for all substrates ofNDP kinase.(ABSTRACT TRUNCATED AT 250 WORDS)
Amino Acid Sequence↗
Proteins with a ring.
Proteins come in all sizes and shapes. Those which fold into a ring with a large hole in the middle may act as a clamp on DNA, a polysaccharide or another protein.
Bacterial Proteins↗
Adenosine 5'-diphosphate binding and the active site of nucleoside diphosphate kinase.
The X-ray structure of nucleoside diphosphate kinase (NDP kinase) from the slime mold Dictyostelium discoideum has been determined to 2.2-A resolution and refined to an R-factor of 0.19 with and without bound ADP-Mg2+. The nucleotide binds near His 122, a residue which becomes phosphorylated during the catalytic cycle. The mode of binding is different from that observed in other phosphokinases, and it involves no glycine-rich sequence. The adenine base makes only nonpolar contacts with the protein. It points outside, explaining the lack of specificity of NDP kinase toward the base. The ribose 2'- and 3'-hydroxyls and the pyrophosphate moiety are H-bonded to polar side chains. A Mg2+ ion bridges the alpha- to the beta-phosphate which approaches the imidazole group of His 122 from the N delta side. The geometry at the active site in the ADP-Mg2+ complex suggests a mechanism for catalysis whereby the gamma-phosphate of a nucleoside triphosphate can be transferred onto His 122 with a minimum of atomic motion.
Adenosine Diphosphate↗
Rigid-body docking with mutant constraints of influenza hemagglutinin with antibody HC19.
An automatic docking algorithm has been applied to the modeling of the complex between hemagglutinin from influenza virus and the Fab fragment of a monoclonal antibody raised against this antigen. We have introduced here the use of biochemical information provided by mutants of hemagglutinin. The docking procedure finds a small number of candidate solutions where three sites of escape mutations are buried and form hydrogen bonds in the interface. The localization of the epitope is improved by additional biochemical data about mutants that do not affect antibody binding. Five candidate solutions with low energy, reasonably well-packed interfaces, and six to ten hydrogen bonds are compatible with mutant information. One of the five stands out as generally better than the others from these points of views.
Algorithms↗