Oxygen reactions with bacterial oxidases and globins: binding,reduction and regulation

Oxygen is favoured as terminal electron acceptor in aerobic and facultative microorganisms because of its appropriate physical state, satisfactory solubility and its desirable combinations of kinetic and thermodynamic properties. Oxygen is generally reduced by four electrons to yield oxygen, but there are important biological consequences of, and roles for, the partial reduction to superoxide and peroxide. Complex and multiple regulatory networks ensure (i) the utilization of oxygen in preference to other oxidants, (ii) the synthesis of oxygen-consuming enzymes with appropriate properties (particularly affinity for the ligand), and (iii) appropriate cellular protection in the event of oxidative stress. This contribution reviews the terminal respiratory oxidases of selected Gram-negative bacteria and microbial haemoglobin-like proteins.Recent studies of the cytochromebd-type oxidases ofEscherichia coli andAzotobacter vinelandii suggest that, despite probable similarity at the amino acid level, the reactivities of these oxidases with oxygen are strikingly different. The respiratory protection afforded to nitrogenase in the obligately aerobic diazotrophA. vinelandii by the cytochromebd complex appears to be accompanied by, and may be the result of, a low affinity for oxygen and a high V_max. The poorly characterized cytochromeo-containing oxidase in this bacterium is not required for respiratory protection. InE. coli, the cytochromebd-type oxidase has a remarkably high affinity for oxygen, consistent with the view that this is an oxygen-scavenging oxidase utilized under microaerobic conditions. The demonstration of substrate (i.e. oxygen) inhibition in this complex suggests a mechanism whereby wasteful electron flux through a non-proton-pumping oxidase is avoided at higher dissolved oxygen tensions. The demonstration of two ligandbinding sites (haemsd andb ₅₉₅) in oxidases of this type suggests plausible mechanisms for this phenomenon. InE. coli, assembly of the cytochromebd-type oxidase (and of periplasmic cytochromesb andc) requires the presence of an ABC transporter, which may serve to export haem or some assembly factor' to the periplasm.There is at least one additional oxygen-consuming protein inE. coli — the flavohaemoglobin encoded by thehmp gene. Globin-like proteins are also widely distributed in other bacteria, fungi and protozoa, but most have unknown functions. The function of HMP and the related chimaeric flavohaemoglobins in other bacteria and yeast is unknown; one of several possibilities for HMP is that its relatively low affinity for oxygen during turnover with NADH as substrate could enable it to function as a sensor of falling (or rising) cytoplasmic oxygen concentrations.(until October 1994: Section of Microbiology, Wing Hall, Cornell University, Ithaca, NY 14853-8101, USA)     

How different amino acid sequences determine similar protein structures: The structure and evolutionary dynamics of the globins

To determine how different amino acid sequences form similar protein structures, and how proteins adapt to mutations that change the volume of residues buried in their close-packed interiors, we have analysed and compared the atomic structures of nine different globins. The homology of the sequences in the two most distantly related molecules is only 16%.The principal determinants of three-dimensional structure of these proteins are the approximately 59 residues involved in helix to helix and helix to haem packings. Half of these residues are buried within the molecules. The observed variations in the sequence keep the side-chains of buried residues non-polar, but do not maintain their size: the mean variation of the volume among homologous amino acids is 56 ų.Changes in the volumes of buried residues are accompanied by changes in the geometry of the helix packings. The relative positions and orientations of homologous pairs of helices in the globins differ by rigid body shifts of up to 7 Å and 30 °. In order to retain functional activity these shifts are coupled so that the geometry of the residues forming the haem pocket is very similar in all the globins.We discuss the implications of these results for the mechanism of protein evolution.     

Polyclonal antibodies raised to phycocyanins contain components specific for the red-absorbing form of phytochrome

Polyclonal antibodies raised in rabbits to a mixture of sodium-dodecyl-sulphate-denatured C- and allo-phycocyanin, isolated from Anabaena cylindrica, cross-react with 124-kilodalton (kDa) phytochrome from etiolated oats, in enzyme-linked immunosorbent assays and on Western blots. The component(s) of the anti-phycocyanin serum that cross-reacts with phytochrome appears to be specific for the red-absorbing form of phytochrome (Pr). These antibodies can be detached from Pr by irradiation with red light, and thus show photoreversible binding. This property has been used to immunopurify the anti-phytochrome component from the antiserum using red light as the eluting agent. Competition assays and epitope-mapping studies indicate that the anti-phytochrome component may bind to a site located between 6 and 10 kDa from the amino-terminus of etiolated oat phytochrome.Abbreviations ELISA enzyme-linked immunosorbent assay - kDa kilodaton - FR far-red light - Pfr far-red-light-absorbing form of phytochrome - Pr red-light-absorbing form of phytochrome - R red light - SDS sodium dodecyl sulphate     

Effector protein structures: a tale of evolutionary relationship

Effector proteins are highly diverse, often lacking similarity in their protein sequences, making it challenging to determine their biological function. Using AlphaFold2 (AF2), Seong and Krasileva recently found that effector structures, but not sequences, share commonality. This helps further understanding of effector evolution across fungal species and reveals unique sequence-unrelated, structurally similar, effector families.     

Defined order of evolutionary adaptations: experimental evidence

Organisms often adapt to new conditions by means of beneficial mutations that become fixed in the population. Often, full adaptation requires several different mutations in the same cell, each of which may affect a different aspect of the behavior. Can one predict order in which these mutations become fixed? To address this, we experimentally studied evolution of Escherichia coli in a growth medium in which the effects of different adaptations can be easily classified as affecting growth rate or the lag‐phase duration. We find that adaptations are fixed in a defined and reproducible order: first reduction of lag phase, and then an increase of the exponential growth rate. A population genetics theory explains this order, and suggests growth conditions in which the order of adaptations is reversed. We experimentally find this order reversal under the predicted conditions. This study supports a view in which the evolutionary path to adaptation in a new environment can be captured by theory and experiment.     

Predicting protein interaction interfaces from protein sequences: case studies of subtilisin and phycocyanin

Identification of protein interaction interfaces is very important for understanding the molecular mechanisms underlying biological phenomena. Here, we present a novel method for predicting protein interaction interfaces from sequences by using PAM matrix (PIFPAM). Sequence alignments for interacting proteins were constructed and parsed into segments using sliding windows. By calculating distance matrix for each segment, the correlation coefficients between segments were estimated. The interaction interfaces were predicted by extracting highly correlated segment pairs from the correlation map. The predictions achieved an accuracy 0.41-0.71 for eight intraprotein interaction examples, and 0.07-0.60 for four interprotein interaction examples. Compared with three previously published methods, PIFPAM predicted more contacting site pairs for 11 out of the 12 example proteins, and predicted at least 34% more contacting site pairs for eight proteins of them. The factors affecting the predictions were also analyzed. Since PIFPAM uses only the alignments of the two interacting proteins as input, it is especially useful when no three-dimensional protein structure data are available.     

A method for simultaneous alignment of multiple protein structures

Here, we present MultiProt, a fully automated highly efficient technique to detect multiple structural alignments of protein structures. MultiProt finds the common geometrical cores between input molecules. To date, most methods for multiple alignment start from the pairwise alignment solutions. This may lead to a small overall alignment. In contrast, our method derives multiple alignments from simultaneous superpositions of input molecules. Further, our method does not require that all input molecules participate in the alignment. Actually, it efficiently detects high scoring partial multiple alignments for all possible number of molecules in the input. To demonstrate the power of MultiProt, we provide a number of case studies. First, we demonstrate known multiple alignments of protein structures to illustrate the performance of MultiProt. Next, we present various biological applications. These include: (1) a partial alignment of hinge-bent domains; (2) identification of functional groups of G-proteins; (3) analysis of binding sites; and (4) protein-protein interface alignment. Some applications preserve the sequence order of the residues in the alignment, whereas others are order-independent. It is their residue sequence order-independence that allows application of MultiProt to derive multiple alignments of binding sites and of protein-protein interfaces, making MultiProt an extremely useful structural tool.     

Strain distribution and linkage relationship of a mouse embryonic hemoglobin variant

Search for structural variants of three globin chains (x, y, z), synthesized only during mouse embryonic hematopoiesis, was carried out by electrophoretic analysis of blood from 12-day embryos, all with C57BL/6 mothers, and fathers from 115 inbred stocks selected for their diverse genetic origins. Structure of the -chains of adult hemoglobins differed among the tested strains, with 57 carrying the Hbb ^sallele, 56 the Hbb ^dallele, and two the Hbb ^pallele. The search revealed no x- or z-chain variants but confirmed and extended knowledge of a previously described y-chain variant. Blood of all embryos sired by males from the 57 Hbb ^sstrains contained only y¹-chains, while blood of all embryos sired by Hbb ^dor Hbb ^pmales contained y2-chains as well as the y¹-chains inherited from their C57 BL/6 mother. The locus controlling structure of the y-chain of mouse embryonic hemoglobins is thus extremely closely linked to the locus controlling structure of adult hemoglobin -chain, with maximum possible recombination frequency less than 0.019.This work was supported in part by Grants CA-01074 from the National Cancer Institute, USPHS, and GM 18684 from the National Institute of General Medical Sciences, in part by Grant ACS-VC58 from The American Cancer Society, in part by grants to the Jackson Laboratory from the Bushrod H. Campbell and Adah F. Hall Charity Fund and the Robert Sterling Clark Foundation, and in part by the Jackson Laboratory Endowment Fund. The Jackson Laboratory is fully accredited by the American Association for Accreditation of Laboratory Animal Care.     

A "Long Indel" model for evolutionary sequence alignment

We present a new probabilistic model of sequence evolution, allowing indels of arbitrary length, and give sequence alignment algorithms for our model. Previously implemented evolutionary models have allowed (at most) single-residue indels or have introduced artifacts such as the existence of indivisible "fragments." We compare our algorithm to these previous methods by applying it to the structural homology dataset HOMSTRAD, evaluating the accuracy of (1) alignments and (2) evolutionary time estimates. With our method, it is possible (for the first time) to integrate probabilistic sequence alignment, with reliability indicators and arbitrary gap penalties, in the same framework as phylogenetic reconstruction. Our alignment algorithm requires that we evaluate the likelihood of any specific path of mutation events in a continuous-time Markov model, with the event times integrated out. To this effect, we introduce a "trajectory likelihood" algorithm (Appendix A). We anticipate that this algorithm will be useful in more general contexts, such as Markov Chain Monte Carlo simulations.     

On the relationship between evolutionary and psychological definitions of altruism and selfishness

I examine the relationship between evolutionary definitions of altruism that are based on fitness effects and psychological definitions that are based on the motives of the actor. I show that evolutionary altruism can be motivated by proximate mechanisms that are psychologically either altruistic or selfish. I also show that evolutionary definitions do rely upon motives as a metaphor in which the outcome of natural selection is compared to the decisions of a psychologically selfish (or altruistic) individual. Ignoring the precise nature of both psychological and evolutionary definitions has obscured many important issues, including the biological roots of psychological altruism.     

Experimental evidence for friction-enhancing integumentary modifications of chameleons and associated functional and evolutionary implications

The striking morphological convergence of hair-like integumentary derivatives of lizards and arthropods (spiders and insects) demonstrates the importance of such features for enhancing purchase on the locomotor substrate. These pilose structures are responsible for the unique tractive abilities of these groups of animals, enabling them to move with seeming ease on overhanging and inverted surfaces, and to traverse inclined smooth substrates. Three groups of lizards are well known for bearing adhesion-promoting setae on their digits: geckos, anoles and skinks. Similar features are also found on the ventral subdigital and distal caudal skin of chameleons. These have only recently been described in any detail, and structurally and functionally are much less well understood than are the setae of geckos and anoles. The seta-like structures of chameleons are not branched (a characteristic of many geckos), nor do they terminate in spatulate tips (which is characteristic of geckos, anoles and skinks). They are densely packed and have attenuated blunt, globose tips or broad, blade-like shafts that are flattened for much of their length. Using a force transducer, we tested the hypothesis that these structures enhance friction and demonstrate that the pilose skin has a greater frictional coefficient than does the smooth skin of these animals. Our results are consistent with friction being generated as a result of side contact of the integumentary filaments. We discuss the evolutionary and functional implications of these seta-like structures in comparison with those typical of other lizard groups and with the properties of seta-mimicking synthetic structures.     

Comparison of the efficiency of evolutionary change-based and side chain orientation-based fold recognition potentials

The present article describes residue level knowledge based potential SORDIS. SORDIS incorporates the information on side-chain orientation in relation to hydrophobic core centres, distance of residue from the globule centre and secondary structure. SORDIS has been tested and compared with widespread evolutionary change-based substitution matrices (BLOSUM, PAM, GONNET, Johnson-Overington, BLAJ, HSDM, and STROMA) in fold recognition experiments within the zone of weak sequence similarity (<16%). The obtained results show that the lower is the amino acid similarity between homologous pairs the higher is the performance of SORDIS in comparison with the potentials, based on the information about the evolutionary changes. Therefore, we propose that the employment of SORDIS in fold recognition can be useful.     

Female choice in Drosophila: evidence from Hawaii and implications for evolutionary biology

Details of female choice of mate in Drosophila silvestris of Hawaii strikingly parallels epigamic behavioral systems in many other animals and may be common in other species of Drosophilidae. Females respond selectively to male circling, wing displays, songs and tactile stimulation with foreleg cilia, a quantitative character that is highly variable in some populations. I hypothesize that the female can exert choice based on these cues from individual males that differ genetically by quantitative trait loci. Laboratory tests show that one third of courting males are repeatedly rejected in favor of a minority of alpha males. This result imposes non-random mating at the local population level. Past multiple-choice lab tests, widely used to measure isolation between pairs of populations or species of Drosophila may be flawed, since random mating has been assumed in the interpretation of results. Pre-mating sexual selection is clearly a powerful intrapopulation force in population biology. This view creates difficulties for discerning any proposed simultaneous interpopulation selective events in the presence of strong female choice. The long-held theory assuming that there is significant selection for pre-mating isolation between groups is questionable.     

The evolutionary relationship of avian and mammalian myosin heavy-chain genes

Summary Sequence comparisons of avian and mammalian skeletal and cardiac myosin heavy-chain isoforms are used to examine the evolutionary relationships of sarcomeric myosin multigene families. Mammalian fast-myosin heavy-chain isoforms forms from different species, with comparable developmental expression, are more similar to each other than they are to other fast isoforms within the same genome. In contrast, the developmentally regulated chicken fast isoforms are more similar to each other than they are to myosin heavy-chain isoforms in other species. Extensive regions of nucleotide identity among the chicken fast myosin heavy chains and in the mouse and rat α- and β-cardiac myosin heavy-chain sequences suggest that geneconversion-like mechanisms have played a major role in the concerted evolution of these gene families. We also conclude that the chicken fast myosin heavy-chain multigene family has undergone recent expansion subsequent to the divergence of birds and mammals and that both the developmental regulation and the specialization of myosin isoforms have likely developed independently in birds and mammals.     

Taxonomic affinities and evolutionary history of the Early Pleistocene hominids of Java: dentognathic evidence

Temporal changes, within-group variation, and phylogenetic positions of the Early Pleistocene Javanese hominids remain unclear. Recent debate focused on the age of the oldest Javanese hominids, but the argument so far includes little morphological basis for the fossils. To approach these questions, we analyzed a comprehensive dentognathic sample from Sangiran, which includes most of the existing hominid mandibles and teeth from the Early Pleistocene of Java. The sample was divided into chronologically younger and older groups. We examined morphological differences between these chronological groups, and investigated their affinities with other hominid groups from Africa and Eurasia. The results indicated that 1) there are remarkable morphological differences between the chronologically younger and older groups of Java, 2) the chronologically younger group is morphologically advanced, showing a similar degree of dentognathic reduction to that of Middle Pleistocene Chinese H. erectus, and 3) the chronologically older group exhibits some features that are equally primitive as or more primitive than early H. erectus of Africa. These findings suggest that the evolutionary history of early Javanese H. erectus was more dynamic than previously thought. Coupled with recent discoveries of the earliest form of H. erectus from Dmanisi, Georgia, the primitive aspects of the oldest Javanese hominid remains suggest that hominid groups prior to the grade of ca. 1.8-1.5 Ma African early H. erectus dispersed into eastern Eurasia during the earlier Early Pleistocene, although the age of the Javanese hominids themselves is yet to be resolved. Subsequent periods of the Early Pleistocene witnessed remarkable changes in the Javanese hominid record, which are ascribed either to significant in situ evolution or replacement of populations.     

Sequence alignments and pair hidden Markov models using evolutionary history

This work presents a novel pairwise statistical alignment method based on an explicit evolutionary model of insertions and deletions (indels). Indel events of any length are possible according to a geometric distribution. The geometric distribution parameter, the indel rate, and the evolutionary time are all maximum likelihood estimated from the sequences being aligned. Probability calculations are done using a pair hidden Markov model (HMM) with transition probabilities calculated from the indel parameters. Equations for the transition probabilities make the pair HMM closely approximate the specified indel model. The method provides an optimal alignment, its likelihood, the likelihood of all possible alignments, and the reliability of individual alignment regions. Human alpha and beta-hemoglobin sequences are aligned, as an illustration of the potential utility of this pair HMM approach.     

Structural distance and evolutionary relationship of networks

Exploring common features and universal qualities shared by a particular class of networks in biological and other domains is one of the important aspects of evolutionary study. In an evolving system, evolutionary mechanism can cause functional changes that forces the system to adapt to new configurations of interaction pattern between the components of that system (e.g. gene duplication and mutation play a vital role for changing the connectivity structure in many biological networks. The evolutionary relation between two systems can be retraced by their structural differences). The eigenvalues of the normalized graph Laplacian not only capture the global properties of a network, but also local structures that are produced by graph evolutions (like motif duplication or joining). The spectrum of this operator carries many qualitative aspects of a graph. Given two networks of different sizes, we propose a method to quantify the topological distance between them based on the contrasting spectrum of normalized graph Laplacian. We find that network architectures are more similar within the same class compared to between classes. We also show that the evolutionary relationships can be retraced by the structural differences using our method. We analyze 43 metabolic networks from different species and mark the prominent separation of three groups: Bacteria, Archaea and Eukarya. This phenomenon is well captured in our findings that support the other cladistic results based on gene content and ribosomal RNA sequences. Our measure to quantify the structural distance between two networks is useful to elucidate evolutionary relationships.     

Homology versus analogy: possible evolutionary relationship of immunoglobulins, cupredoxins, and Cu,Zn-superoxide dismutase

The 'immunoglobulin-like' fold is one of most common structural motifs observed in proteins. This topology is found in more than 80 superfamilies of proteins, including Cu,Zn-superoxide dismutase (SOD) and cupredoxin. Evolutionary relationships have not been identified, but may exist. The challenge remains, therefore, of resolving the issue of whether the diverse distribution of the fold is accounted for by divergent evolution of function or convergent evolution of structure following multiple independent origins of function. Since the early studies that revealed conformational similarity of immunoglobulins and other proteins, the number of primary structures available for comparison has dramatically increased and new computational approaches for analysis of sequences have been developed. It now appears that a hypothesis of a common evolutionary origin for cupredoxins, Cu,Zn-SOD, and immunoglobulins may be credible. The distinction between protein homology and protein analogy is fundamental. The immunoglobulin-like fold may represent a robust system within which to examine again the issue of protein homology versus analogy.     

An evolutionary analysis of the relationship between spite and altruism

We investigate the selective pressures on a social trait when evolution occurs in a population of constant size. We show that any social trait that is spiteful simultaneously qualifies as altruistic. In other words, any trait that reduces the fitness of less related individuals necessarily increases that of related ones. Our analysis demonstrates that the distinction between "Hamiltonian spite" and "Wilsonian spite" is not justified on the basis of fitness effects. We illustrate this general result with an explicit model for the evolution of a social act that reduces the recipient's survival ("harming trait"). This model shows that the evolution of harming is favoured if local demes are of small size and migration is low (philopatry). Further, deme size and migration rate determine whether harming evolves as a selfish strategy by increasing the fitness of the actor, or as a spiteful/altruistic strategy through its positive effect on the fitness of close kin.     

Hematopoietic cytokines: similarities and differences in the structures, with implications for receptor binding.

Crystal and NMR structures of helical cytokines--interleukin-4 (IL-4), granulocyte-macrophage colony-stimulating factor (GM-CSF), and interleukin-2 (IL-2)--have been compared. Root mean square deviations in the C alpha coordinates for the conserved regions of the helices were 1-2 A between different cytokines, about twice the differences observed for independently determined crystal and solution structures of IL-4. Considerable similarity in amino acid sequence in the areas expected to interact with the receptors was detected, and the available mutagenesis data for these cytokines were correlated with structure conservation. Models of cytokine-receptor interactions were postulated for IL-4 based on its structure as well as on the published structure of human growth hormone interacting with its receptors (de Vos, A.M., Ultsch, M., & Kossiakoff, A.A., 1992, Science 255, 306-312). Patches of positively charged residues on the surfaces of helices C and D of IL-4 may be responsible for the interactions with the negatively charged residues found in the complementary parts of the IL-4 receptors.