The spatial distribution of physical, chemical, energetic and conformational properties of amino acid residues in globular proteins

Natural Selection in Virulence Genes of Francisella tularensis

A fundamental tenet of evolution is that alleles that are under negative selection are often deleterious and confer no evolutionary advantage. Negatively selected alleles are removed from the gene pool and are eventually extinguished from the population. Conversely, alleles under positive selection do confer an evolutionary advantage and lead to an increase in the overall fitness of the organism. These alleles increase in frequency until they eventually become fixed in the population. Francisella tularensis is a zoonotic pathogen and a potential biothreat agent. The most virulent type of F. tularensis, Type A, is distributed across North America with Type A.I occurring mainly in the east and Type A.II appearing mainly in the west. F. tularensis is thought to be a genome in decay (losing genes) because of the relatively large number of pseudogenes present in its genome. We hypothesized that the observed frequency of gene loss/pseudogenes may be an artifact of evolution in response to a changing environment, and that genes involved in virulence should be under strong positive selection. To test this hypothesis, we sequenced and compared whole genomes of Type A.I and A.II isolates. We analyzed a subset of virulence and housekeeping genes from several F. tularensis subspecies genomes to ascertain the presence and extent of positive selection. Eleven previously identified virulence genes were screened for positive selection along with 10 housekeeping genes. Analyses of selection yielded one housekeeping gene and 7 virulence genes which showed significant evidence of positive selection at loci implicated in cell surface structures and membrane proteins, metabolism and biosynthesis, transcription, translation and cell separation, and substrate binding and transport. Our results suggest that while the loss of functional genes through disuse could be accelerated by negative selection, the genome decay in Francisella could also be the byproduct of adaptive evolution driven by complex interactions between host, pathogen, and thier environment, as evidenced by several of its virulence genes which are undergoing strong, positive selection.

Evidence of adaptive evolution of alpine pheasants to high-altitude environment from mitogenomic perspective

Adaptive evolutions to high-altitude adaptation have been intensively studied in mammals. However, considering the additional vertebrate groups, new perception regarding selection challenged by high-altitude stress on mitochondrial genome can be gained. To test this hypothesis, we compiled and analyzed the mitochondrial genomes of 5 alpine pheasants and 12 low-altitude species in Phasianidae. The results that evolutionary rates of ATP6 and ND6 showing significant fluctuation among branches when involved with five alpine pheasants revealed both genes might have implications with adapting to highland environment. The radical physico-chemical property changes identified by the modified MM01 model, including composition (C) and equilibrium constant (ionization of COOH) (Pk') in ATP6 and beta-structure tendencies (Pβ), Pk', and long-range non-bonded energy (El) in ND6, suggested that minor overall adjustments in size, protein conformation and relative orientation of reaction interfaces have been optimized to provide the ideal environments for electron transfer, proton translocation and generation of adenosine triphosphate (ATP). Additionally, three unique substitution sites were identified under selection in ND6, which could be potentially important adaptive changes contributing to cellular energy production. Our findings suggested that adaptive evolution may occur in alpine pheasants, which are an important complement to the knowledge of genetic mechanisms against the high-altitude environment in non-mammal animals.

Directional Darwinian Selection in proteins

Background

Molecular evolution is a very active field of research, with several complementary approaches, including dN/dS, HON90, MM01, and others. Each has documented strengths and weaknesses, and no one approach provides a clear picture of how natural selection works at the molecular level. The purpose of this work is to present a simple new method that uses quantitative amino acid properties to identify and characterize directional selection in proteins.

Methods

Inferred amino acid replacements are viewed through the prism of a single physicochemical property to determine the amount and direction of change caused by each replacement. This allows the calculation of the probability that the mean change in the single property associated with the amino acid replacements is equal to zero (H₀: μ = 0; i.e., no net change) using a simple two-tailed t-test.

Results

Example data from calanoid and cyclopoid copepod cytochrome oxidase subunit I sequence pairs are presented to demonstrate how directional selection may be linked to major shifts in adaptive zones, and that convergent evolution at the whole organism level may be the result of convergent protein adaptations.

Conclusions

Rather than replace previous methods, this new method further complements existing methods to provide a holistic glimpse of how natural selection shapes protein structure and function over evolutionary time.

Divergent evolution and molecular adaptation in the Drosophila odorant-binding protein family: inferences from sequence variation at the OS-E and OS-F genes

BACKGROUND

The Drosophila Odorant-Binding Protein (Obp) genes constitute a multigene family with moderate gene number variation across species. The OS-E and OS-F genes are the two phylogenetically closest members of this family in the D. melanogaster genome. In this species, these genes are arranged in the same genomic cluster and likely arose by tandem gene duplication, the major mechanism proposed for the origin of new members in this olfactory-system family.

RESULTS

We have analyzed the genomic cluster encompassing OS-E and OS-F genes (Obp83 genomic region) to determine the role of the functional divergence and molecular adaptation on the Obp family size evolution. We compared nucleotide and amino acid variation across 18 Drosophila and 4 mosquito species applying a phylogenetic-based maximum likelihood approach complemented with information of the OBP three-dimensional structure and function. We show that, in spite the OS-E and OS-F genes are currently subject to similar and strong selective constraints, they likely underwent divergent evolution. Positive selection was likely involved in the functional diversification of new copies in the early stages after the gene duplication event; moreover, it might have shaped nucleotide variation of the OS-E gene concomitantly with the loss of functionally related members. Besides, molecular adaptation likely affecting the functional OBP conformational changes was supported by the analysis of the evolution of physicochemical properties of the OS-E protein and the location of the putative positive selected amino acids on the OBP three-dimensional structure.

CONCLUSION

Our results support that positive selection was likely involved in the functional differentiation of new copies of the OBP multigene family in the early stages after their birth by gene duplication; likewise, it might shape variation of some members of the family concomitantly with the loss of functionally related genes. Thus, the stochastic gene gain/loss process coupled with the impact of natural selection would influence the observed OBP family size.

Physicochemical Evolution and Molecular Adaptation of the Cetacean and Artiodactyl Cytochrome b Proteins

Cetaceans have most likely experienced metabolic shifts since evolutionarily diverging from their terrestrial ancestors, shifts that may be reflected in the proteins such as cytochrome b that are responsible for metabolic efficiency. However, accepted statistical methods for detecting molecular adaptation are largely biased against even moderately conservative proteins because the primary criterion involves a comparison of nonsynonymous and synonymous substitution rates (dN/dS); they do not allow for the possibility that adaptation may come in the form of very few amino acid changes. We apply the MM01 model to the possible molecular adaptation of cytochrome b among cetaceans because it does not rely on a dN/dS ratio, instead evaluating positive selection in terms of the amino acid properties that comprise protein phenotypes that selection at the molecular level may act upon. We also apply the codon-degeneracy model (CDM), which focuses on evaluating overall patterns of nucleotide substitution in terms of base exchange, codon position, and synonymy to estimate the overall effect of selection. Using these relatively new models, we characterize the molecular adaptation that has occurred in the cetacean cytochrome b protein by comparing revealed amino acid replacement patterns to those found among artiodactyls, the modern terrestrial mammals found to be most closely related to cetaceans. Our findings suggest that several regions of the cetacean cytochrome b protein have experienced molecular adaptation. Also, these adaptations are spatially associated with domain structure, protein function, and the structure and function of the cytochrome bc(1) complex and its constituents. We also have found a general correlation between the results of the analytical software programs TreeSAAP (which implements the MM01 model) and CDM (which implements the codon-degeneracy model).

Factors influencing the stability of alpha-helices and beta-strands in thermophilic ribonuclease H

Understanding the influence of structural parameters is crucial to enhance the thermal stability of proteins. In this work, the stability (deltaG) of residues in different secondary structures of Ribonuclease H (RNase H) has been analyzed with 48 amino acid properties. The properties reflecting hydrophobicity show a good correlation with stability. Further, the linear distribution of surrounding hydrophobicity in alpha-helices, obtained from the three dimensional structure of thermophilic RNase H, agrees well with experimental deltaG values. Moreover, the stability parameters correlate better in alpha-helices than those did in beta-strand segments. Multiple regression analysis, incorporating combinations of three properties from among all possible combinations of the 48 properties, increased the correlation coefficient to 0.77.

Important amino acid properties for enhanced thermostability from mesophilic to thermophilic proteins

Understanding the role of various interactions in enhancing the thermostability of proteins is important not only for clarifying the mechanism of protein stability but also for designing stable proteins. In this work, we have analyzed the thermostability of 16 different families by comparing mesophilic and thermophilic proteins with 48 various physicochemical, energetic and conformational properties. We found that the increase in shape, s (location of branch point in side chain) increases the thermostability, whereas, an opposite trend is observed for Gibbs free energy change of hydration for native proteins, GhN, in 14 families. A good correlation is observed between these two properties and the simultaneous increases of -GhN and s is necessary to enhance the thermostability from mesophile to thermophile. The increase in shape, which tends to increase with increasing number of carbon atoms both for polar and non-polar residues, may generate more packing and compactness, and the position of beta and higher order branches may be important for better packing. On the other hand, the increase in -GhN in thermophilic proteins increases the solubility of the proteins. This tendency counterbalances the increases in insolubility and unfolding heat capacity change due to the increase in the number of carbon atoms. Thus, the present results suggest that the stability of thermophilic proteins may be achieved by a balance between better packing and solubility.

Role of structural and sequence information in the prediction of protein stability changes: comparison between buried and partially buried mutations

Predicting mutation-induced changes in protein stability is one of the greatest challenges in molecular biology. In this work, we analyzed the correlation between stability changes caused by buried and partially buried mutations and changes in 48 physicochemical, energetic and conformational properties. We found that properties reflecting hydrophobicity strongly correlated with stability of buried mutations, and there was a direct relation between the property values and the number of carbon atoms. Classification of mutations based on their location within helix, strand, turn or coil segments improved the correlation of mutations with stability. Buried mutations within beta-strand segments correlated better than did those in alpha-helical segments, suggesting stronger hydrophobicity of the beta-strands. The stability changes caused by partially buried mutations in ordered structures (helix, strand and turn) correlated most strongly and were mainly governed by hydrophobicity. Due to the disordered nature of coils, the mechanism underlying their stability differed from that of the other secondary structures: the stability changes due to mutations within the coil were mainly influenced by the effects of entropy. Further classification of mutations within coils, based on their hydrogen-bond forming capability, led to much stronger correlations. Hydrophobicity was the major factor in determining the stability of buried mutations, whereas hydrogen bonds, other polar interactions and hydrophobic interactions were all important determinants of the stability of partially buried mutations. Information about local sequence and structural effects were more important for the prediction of stability changes caused by partially buried mutations than for buried mutations; they strengthened correlations by an average of 27% among all data sets.

A new spatial representation of amino acid residues in globular proteins

For the globular proteins with known three-dimensional structures, an ellipsoid model of each protein was constructed with least volume and its dimensions were derived. The spatial arrangements were made for the C alpha and side chain atoms of that protein within that ellipsoid. This new spatial representation shows the residue position from the centroid, as well as the depth from the surface. The average spatial parameters were then calculated. The correlations between these new spatial parameters and the existing parameters of the amino acid residues were then derived.

Hydrophobic distribution and spatial arrangement of amino acid residues in membrane proteins

The analysis of known three-dimensional structures of membrane proteins provides an opportunity to understand their structure and stability. In this article we analyse the hydrophobic variation of amino acid residues at various ranges in membrane and aqueous parts of membrane proteins. The numerical indices for several properties of amino acid residues in membrane proteins, such as surrounding hydrophobicity, gain in surrounding hydrophobicity, hydrophobic gain ratio, accessible surface area, preference of amino acid residues in the interior and surface parts, solvent accessible reduction ratio and buriedness, were set up. The relative preference of amino acid residues at various positions of membrane proteins were obtained in a very realistic approach.

Somatic hypermutagenesis in immunoglobulin genes. I. Correlation between somatic mutations and repeats. Somatic mutation properties and clonal selection

Based on the analysis of some immunoglobulin V-gene sequences, somatic mutations are assumed to occur by correction of complementary violations in heteroduplexes formed by direct or inverted repeats. Correlation between somatic mutations and repeats is investigated by a statistical weights method in 323 somatic substitutions in 14 V-genes. Assuming absence of correlation, the probability of observing data in the sample would be very low (0.00004). This result supports the idea that somatic mutations may arise from heteroduplex repair. The high frequency of these mutations in complementarity-determining regions (CDRs) of V-genes may be due to a high concentration of repeats in these regions. Analysis of somatic substitutions has revealed that stabilizing selection seems to provide conservation of framework regions (FRs) (which leads to preservation of the protein's three-dimensional structure). Positive selection may be provided by B-lymphocyte proliferation with large changes in CDRs.

Variation of amino acid properties in protein secondary structures, alpha-helices and beta-strands

A study was made on the physical, chemical, energetic, conformational, geometric, and dynamic property potentials of amino acid residues in protein secondary structures: alpha-helix and beta-strand. Property patterns were obtained by computing the average property values for specified residue units partitioned longitudinally and transversely about the chain. It was found that in alpha-helices with not more than 15 residues, there exist longitudinally opposing portions, one characteristically higher in average property potentials than the other. The helical chain, in general, acquires either an increasing or decreasing average potential in the N-terminal to C-terminal direction. The sequence-wise and surface-wise variations of property potentials in the elements of beta-structure also revealed such general patterns. Possible wrong predictions in statistical methods of one secondary structural class over the other are pointed out.

Genetic code preferentially conserves long-range interactions among the amino acids

The physical properties of amino acids were investigated in order to evaluate their possible relationship to the assignment of codons for amino acids in the genetic code. A comparison of the interconversion probability between amino acids and the distances between the amino acids for individual physical properties revealed a striking hierarchy among the physical properties. Surprisingly, it is the long-range/solvent interactions and not the short-range/stereochemical properties which are preferentially conserved in the genetic code.

On the appearance and role of a spacer group in the protein amino acids

Of the 20 protein amino acids, 16 have a methylene group at the beta position, and a further three bear a methine group. No aromatic, carboxamido, carboxylic carbon, or hetero atoms are attached directly to the alpha carbon, but they are separated by this methylene or occasionally by a longer n-alkylene spacer group. Therefore, the structure of the protein amino acids should rather be formulated as H2N-CH((CH2)n-R')-COOH instead of the generally accepted H2N-CH(R)-COOH. The appearance of and the role played by the spacer group are discussed in an evolutionary context. It is suggested that the spacer group appeared as a result of prebiotic selection, based on the relative abundance, racemization rate, and suitability for thermal polymerization of the protein amino acids and their homologs with various spacer group lengths. At the biotic level of evolution the requirements for ribosomal polymerization, as well as the abilities of polypeptides to maintain a stable and flexible three-dimensional structure and to bind ligands are considered and are proposed to have been responsible for the possible exclusion of longer spacer groups. It is concluded that the general role of the spacer group is to ensure the uniformity of the constant regions H2N-CH(-)-COOH and the individuality of the R' contact groups by spatially separating them.

Analysis of the factors and implications of an empirical method for estimating the stability of mutant human haemoglobins

An empirical method for estimating the effects of single amino acid substitutions on structural stability of proteins with known spatial structure is developed. Twenty physical and chemical properties of amino acids and characteristics of protein tertiary structure were analysed to determine those most involved in producing instability. We employed data on 330 mutant variants of the alpha- and beta-subunits of human haemoglobin in choice of the parameters of the method developed which yielded a 81% of prediction accuracy of stability estimates for human mutant haemoglobins.

Differential equation model to study dynamic behaviour of globular proteins

By regarding the globular proteins as spheroidal shaped bodies of uniform density, a differential equation model to study their low amplitude fluctuations was developed. It was then applied to the crystal structures of pancreatic trypsin inhibitor and ferrocytochrome c. The results were tested by comparing them with those of the dynamic simulation and temperature factor studies on the same proteins.

Spatial constraints and group behaviour in globular proteins

A comparative analysis has been attempted on the spatial placement of amino acid residues derived from radial, ellipsoidal and exposure arrangements. The group behaviour of residues and their restraining influence in protein folding have been brought out. A study is also made on the geometry of proteins, the exposure arrangement of residues and the spatial distribution of the physical properties of the residues in globular proteins. It has been shown that the group constraints along with the information on the shape of the globular proteins would be highly useful in assigning the spatial and exposure arrangements of residues in globular proteins.

Conservation of physico-chemical amino acid properties during the evolution of proteins

Based on a similarity ring constructed from a substitution probability matrix, we have analyzed the conservation of some amino acid properties in the evolution of proteins. Refractive index and bulkiness are highly conserved, hydrophobicity and polarity are fairly well conserved while optical rotation appears to be a less relevant property. On the other hand, the analysis of the correspondence between phenotype and genotype shows that the most frequent amino acid substitutions in proteins do not always correspond to the most feasible codon changes. The apparent disagreement between amino acid substitutions in modern proteins and the primordial amino acid-codon assignment is discussed.

Code dependent conservation of the physico-chemical properties in amino acid substitutions

The frequency of amino acid replacements in families of typical proteins has been elegantly analyzed by Argyle (1980) showing that the most frequent replacements involve a conservation of the amino acid chemical properties. The cyclic arrangement of the twenty amino acids resulting from the most frequent replacements has been described as an amino acid chemical ring. In this work, a novel amino acid replacement frequency ring is proposed, for which a conservation of over 90% of the most general physico-chemical properties can be deduced. The amino acid chemical similarity ring is also analyzed in terms of the genetic code base probability changes, showing that the discrepancy that exists between the standard deviation value of the amino acid replacement frequency matrix and its respective ideal value is almost equal to that deduced from the corresponding base codon replacement probability matrices. These differences are finally evaluated and discussed in terms of the restrictions imposed by the structure of the genetic code and the physico-chemical dissimilarities between some codons of amino acids which are chemically similar.

Application of pattern recognition techniques to the analysis of protein crystal structure data. I. Characteristics of per-residue side chain contact frequency distributions

A statistical study of amino acid side chain contact interactions was carried out using a data set based on 36 protein structures. For each type of amino acid, a distribution of per-residue inter-side-chain contacts was obtained, over the observed span of zero to 11 contacts per residue. Significant observations included the following: 1) The mean number of inter-side-chain contacts is proportional to side chain surface area with the exception of Lys and Arg. 2) The mean number of contacts was greater for amino acids in beta-sheet relative to alpha-helical regions. 3) The more polar or surface-loving amino acids exhibited non-normal distributions, whereas distributions for the non-polar or interior-loving amino acids fell within accepted limits of normality.

A neutral theory of biogenesis

The selective Darwinian theory of chemical evolution is critically reviewed and the tentative conclusion is reached that neither the theoretical analyses nor the experiments with phages can really prove it. An alternative proposal is put forth which considers the possibility that the biogenetic process has been driven by stochastic forces, e.g. it took place in the absence of Darwinian selection which, in turn, started only when the first protocells came into existence. The dynamics of the early self-organization of living structures should be understood in terms of self-assembly. The complexification of living matter is thus not represented as a gradual phenomenon but as a series of abrupt and relatively fast transitions consisting in the aggregation of pre-systems which had evolved by their own. The shift towards new and variegated states proposed by the bifurcation theory are not considered particularly relevant for reasons reported in the test, nor is it believed that dissipation can entirely account for the order observed in living cells.

Hydrophobic packing and spatial arrangement of amino acid residues in globular proteins

Amino acid residues acquire characteristic hydrophobic environments in globular proteins. Using the crystal data on 21 proteins, a new scale of hydrophobic indices for the residues is set up. This scale provides valuable information with regard to hydrophobic domains, nucleation sites, surface domains, loop sites and the spatial positions of residues in protein molecules.