Cytoglobin (Cygb) was recently discovered in the human genome and localized in different tissues. It was suggested to play tissue-specific protective roles, spanning from scavenging of reactive oxygen species in neurons to supplying oxygen to enzymes in fibroblasts. To shed light on the functioning of such versatile machinery, we have studied the processes supporting transport of gaseous heme ligands in Cygb. Carbon monoxide rebinding shows a complex kinetic pattern with several distinct reaction intermediates, reflecting rebinding from temporary docking sites, second order recombination, and formation (and dissociation) of a bis-histidyl heme hexacoordinated reaction intermediate. Ligand exit to the solvent occurs through distinct pathways, some of which exploit temporary docking sites. The remarkable change in energetic barriers, linked to heme bis-histidyl hexacoordination by HisE7, may be responsible for active regulation of the flux of reactants and products to and from the reaction site on the distal side of the heme. A substantial change in both protein dynamics and inner cavities is observed upon transition from the CO-liganded to the pentacoordinated and bis-histidyl hexacoordinated species, which could be exploited as a signalling state. These findings are consistent with the expected versatility of the molecular activity of this protein.

Background

Conjugation of human and animal hemoglobins with polyethylene glycol has been widely explored as a means to develop blood substitutes, a novel pharmaceutical class to be used in surgery or emergency medicine. However, PEGylation of human hemoglobin led to products with significantly different oxygen binding properties with respect to the unmodified tetramer and high NO dioxygenase reactivity, known causes of toxicity. These recent findings call for the biotechnological development of stable, low-affinity PEGylated hemoglobins with low NO dioxygenase reactivity.

Results

To investigate the effects of PEGylation on protein structure and function, we compared the PEGylation products of human hemoglobin and Trematomus bernacchii hemoglobin, a natural variant endowed with a remarkably low oxygen affinity and high tetramer stability. We show that extension arm facilitated PEGylation chemistry based on the reaction of T. bernacchii hemoglobin with 2-iminothiolane and maleimido-functionalyzed polyethylene glycol (MW 5000 Da) leads to a tetraPEGylated product, more homogeneous than the corresponding derivative of human hemoglobin. PEGylated T. bernacchii hemoglobin largely retains the low affinity of the unmodified tetramer, with a p50 50 times higher than PEGylated human hemoglobin. Moreover, it is still sensitive to protons and the allosteric effector ATP, indicating the retention of allosteric regulation. It is also 10-fold less reactive towards nitrogen monoxide than PEGylated human hemoglobin.

Conclusions

These results indicate that PEGylated hemoglobins, provided that a suitable starting hemoglobin variant is chosen, can cover a wide range of oxygen-binding properties, potentially meeting the functional requirements of blood substitutes in terms of oxygen affinity, tetramer stability and NO dioxygenase reactivity.

Background

There is a great interest in understanding and exploiting protein-protein associations as new routes for treating human disease. However, these associations are difficult to structurally characterize or model although the number of X-ray structures for protein-protein complexes is expanding. One feature of these complexes that has received little attention is the role of water molecules in the interfacial region.

Methodology

A data set of 4741 water molecules abstracted from 179 high-resolution (≤ 2.30 Å) X-ray crystal structures of protein-protein complexes was analyzed with a suite of modeling tools based on the HINT forcefield and hydrogen-bonding geometry. A metric termed Relevance was used to classify the general roles of the water molecules.

Results

The water molecules were found to be involved in: a) (bridging) interactions with both proteins (21%), b) favorable interactions with only one protein (53%), and c) no interactions with either protein (26%). This trend is shown to be independent of the crystallographic resolution. Interactions with residue backbones are consistent for all classes and account for 21.5% of all interactions. Interactions with polar residues are significantly more common for the first group and interactions with non-polar residues dominate the last group. Waters interacting with both proteins stabilize on average the proteins' interaction (−0.46 kcal mol−1), but the overall average contribution of a single water to the protein-protein interaction energy is unfavorable (+0.03 kcal mol−1). Analysis of the waters without favorable interactions with either protein suggests that this is a conserved phenomenon: 42% of these waters have SASA ≤ 10 Å2 and are thus largely buried, and 69% of these are within predominantly hydrophobic environments or “hydrophobic bubbles”. Such water molecules may have an important biological purpose in mediating protein-protein interactions.

The inhibition of cysteine biosynthesis in prokaryotes and protozoa has been proposed to be relevant for the development of antibiotics. Haemophilus influenzae O-acetylserine sulfhydrylase (OASS), catalyzing L-cysteine formation, is inhibited by the insertion of the C-terminal pentapeptide (MNLNI) of serine acetyltransferase into the active site. 400 MNXXI pentapeptides were generated, docked into OASS active site using GOLD and scored with HINT. The terminal P5 Ile accounts for about 50% of the binding energy. Glu or Asp at position P4, and to a lesser extent, at position P3, also significantly contribute to the binding interaction. The predicted affinity of 14 selected pentapeptides correlated well with the experimentally determined dissociation constants. The X-ray structure of three high affinity pentapeptides-OASS complexes were compared with the docked poses. These results, combined with a GRID analysis of the active site, allowed us to define a pharmacophoric scaffold for the design of peptidomimetic inhibitors.

The oligomerization/co-localization of protein complexes and their cooperative regulation in protein function is a key feature in many biological systems. The synergistic regulation in different subunits often enhances the functional properties of the multi-enzyme complex. The present study used molecular dynamics and Brownian dynamics simulations to study the effects of allostery, oligomerization and intermediate channeling on enhancing the protein function of tryptophan synthase (TRPS). TRPS uses a set of α/β–dimeric units to catalyze the last two steps of L-tryptophan biosynthesis, and the rate is remarkably slower in the isolated monomers. Our work shows that without their binding partner, the isolated monomers are stable and more rigid. The substrates can form fairly stable interactions with the protein in both forms when the protein reaches the final ligand–bound conformations. Our simulations also revealed that the α/β–dimeric unit stabilizes the substrate–protein conformation in the ligand binding process, which lowers the conformation transition barrier and helps the protein conformations shift from an open/inactive form to a closed/active form. Brownian dynamics simulations with a coarse-grained model illustrate how protein conformations affect substrate channeling. The results highlight the complex roles of protein oligomerization and the fine balance between rigidity and dynamics in protein function.

Author Summary

Conformational changes of enzymes are often related to regulating and creating an optimal environment for efficient chemistry. An increasing number of evidences also indicate that oligomerization/co-localization of proteins contributes to the efficiency of metabolic pathways. Although static structures have been available for many multi-enzyme complexes, their efficiency is also governed by the synergistic regulation between the multi-units. Our study applies molecular dynamics and Brownian dynamics simulations to the model system, the tryptophan synthase complex. The multi-enzyme complex is a bienzyme nanomachine and its catalytic activity is intimately related to the allosteric signaling and the metabolite transfer between its α– and β–subunits connected by a 25-Å long channel. Our studies suggest that the binding partner is crucial for the ligand binding processes. Although the isolated monomers are stable in the ligand–free state and can form stable interaction if the substrate is in the final bound conformation, it has higher energy barrier when ligand binds to the active site. We also show that the channel does not always exist, but it may be blocked before the enzyme reaches its final bound conformation. The results highlight the importance of forming protein complexes and the cooperative changes during different states.

In the search for xenoestrogens within food additives, we have analyzed the Joint FAO-WHO expert committee database, containing 1500 compounds, using an integrated in silico and in vitro approach. This analysis identified 31 potential estrogen receptor α ligands that were reduced to 13 upon applying a stringent filter based on ligand volume and binding mode. Among the 13 potential xenoestrogens, four were already known to exhibit an estrogenic activity, and the other nine were assayed in vitro, determining the binding affinity to the receptor and biological effects. Propyl gallate was found to act as an antagonist, and 4-hexylresorcinol was found to act as a potent transactivator; both ligands were active at nanomolar concentrations, as predicted by the in silico analysis. Some caution should be issued for the use of propyl gallate and 4-hexylresorcinol as food additives.

Background

To understand the energetics of the interaction between protein and DNA we analyzed 39 crystallographically characterized complexes with the HINT (Hydropathic INTeractions) computational model. HINT is an empirical free energy force field based on solvent partitioning of small molecules between water and 1-octanol. Our previous studies on protein-ligand complexes demonstrated that free energy predictions were significantly improved by taking into account the energetic contribution of water molecules that form at least one hydrogen bond with each interacting species.

Results

An initial correlation between the calculated HINT scores and the experimentally determined binding free energies in the protein-DNA system exhibited a relatively poor r2 of 0.21 and standard error of ± 1.71 kcal mol-1. However, the inclusion of 261 waters that bridge protein and DNA improved the HINT score-free energy correlation to an r2 of 0.56 and standard error of ± 1.28 kcal mol-1. Analysis of the water role and energy contributions indicate that 46% of the bridging waters act as linkers between amino acids and nucleotide bases at the protein-DNA interface, while the remaining 54% are largely involved in screening unfavorable electrostatic contacts.

Conclusion

This study quantifies the key energetic role of bridging waters in protein-DNA associations. In addition, the relevant role of hydrophobic interactions and entropy in driving protein-DNA association is indicated by analyses of interaction character showing that, together, the favorable polar and unfavorable polar/hydrophobic-polar interactions (i.e., desolvation) mostly cancel.