Fast computational mutation-response scanning of proteins

Studying the effect of perturbations on protein structure is a basic approach in protein research. Important problems, such as predicting pathological mutations and understanding patterns of structural evolution, have been addressed by computational simulations that model mutations using forces and predict the resulting deformations. In single mutation-response scanning simulations, a sensitivity matrix is obtained by averaging deformations over point mutations. In double mutation-response scanning simulations, a compensation matrix is obtained by minimizing deformations over pairs of mutations. These very useful simulation-based methods may be too slow to deal with large proteins, protein complexes, or large protein databases. To address this issue, I derived analytical closed formulas to calculate the sensitivity and compensation matrices directly, without simulations. Here, I present these derivations and show that the resulting analytical methods are much faster than their simulation counterparts.Fil: Echave, Julián. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Ciencias Físicas. - Universidad Nacional de San Martín. Instituto de Ciencias Físicas; Argentin

Relationship between protein thermodynamic constraints and variation of evolutionary rates among sites

Evolutionary-rate variation among sites within proteins depends on functional and biophysical properties that constrain protein evolution. It is generally accepted that proteins must be able to fold stably in order to function. However, the relationship between stability constraints and among-sites rate variation is not well understood. Here, we present a biophysical model that links the thermodynamic stability changes due to mutations at sites in proteins (ΔΔG) to the rate at which mutations accumulate at those sites over evolutionary time. We find that such a 'stability model' generally performs well, displaying correlations between predicted and empirically observed rates of up to 0.75 for some proteins. We further find that our model has comparable predictive power as does an alternative, recently proposed 'stress model' that explains evolutionary-rate variation among sites in terms of the excess energy needed for mutants to adopt the correct active structure (ΔΔG∗). The two models make distinct predictions, though, and for some proteins the stability model outperforms the stress model and vice versa. We conclude that both stability and stress constrain site-specific sequence evolution in proteins.Fil: Echave, Julián. Universidad Nacional de San Martín. Escuela de Ciencia y Tecnología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Jackson, Eleisha L.. University of Texas at Austin; Estados UnidosFil: Wilke, Claus O.. University of Texas at Austin; Estados Unido

The structurally constrained protein evolution model accounts for sequence patterns of the LβH superfamily

BACKGROUND: Structure conservation constrains evolutionary sequence divergence, resulting in observable sequence patterns. Most current models of protein evolution do not take structure into account explicitly, being unsuitable for investigating the effects of structure conservation on sequence divergence. To this end, we recently developed the Structurally Constrained Protein Evolution (SCPE) model. The model starts with the coding sequence of a protein with known three-dimensional structure. At each evolutionary time-step of an SCPE simulation, a trial sequence is generated by introducing a random point mutation in the current coding DNA sequence. Then, a "score" for the trial sequence is calculated and the mutation is accepted only if its score is under a given cutoff, λ. The SCPE score measures the distance between the trial sequence and a given reference sequence, given the structure. In our first brief report we used a "global score", in which the same reference sequence, the ancestral one, was used at each evolutionary step. Here, we introduce a new scoring function, the "local score", in which the sequence accepted at the previous evolutionary time-step is used as the reference. We assess the model on the UDP-N-acetylglucosamine acyltransferase (LPXA) family, as in our previous report, and we extend this study to all other members of the left-handed parallel beta helix fold (LβH) superfamily whose structure has been determined. RESULTS: We studied site-dependent entropies, amino acid probability distributions, and substitution matrices predicted by SCPE and compared with experimental data for several members of the LβH superfamily. We also evaluated structure conservation during simulations. Overall, SCPE outperforms JTT in the description of sequence patterns observed in structurally constrained sites. Maximum Likelihood calculations show that the local-score and global-score SCPE substitution matrices obtained for LPXA outperform the JTT model for the LPXA family and for the structurally constrained sites of class i of other members within the LβH superfamily. CONCLUSION: We extended the SCPE model by introducing a new scoring function, the local score. We performed a thorough assessment of the SCPE model on the LPXA family and extended it to all other members of known structure of the LβH superfamily

Functional Sites Induce Long-Range Evolutionary Constraints in Enzymes

Functional residues in proteins tend to be highly conserved over evolutionary time. However, to what extent functional sites impose evolutionary constraints on nearby or even more distant residues is not known. Here, we report pervasive conservation gradients toward catalytic residues in a dataset of 524 distinct enzymes: evolutionary conservation decreases approximately linearly with increasing distance to the nearest catalytic residue in the protein structure. This trend encompasses, on average, 80% of the residues in any enzyme, and it is independent of known structural constraints on protein evolution such as residue packing or solvent accessibility. Further, the trend exists in both monomeric and multimeric enzymes and irrespective of enzyme size and/or location of the active site in the enzyme structure. By contrast, sites in protein–protein interfaces, unlike catalytic residues, are only weakly conserved and induce only minor rate gradients. In aggregate, these observations show that functional sites, and in particular catalytic residues, induce long-range evolutionary constraints in enzymes.Fil: Jack, Benjamin R.. University of Texas at Austin; Estados UnidosFil: Meyer, Austin G.. University of Texas at Austin; Estados UnidosFil: Echave, Julián. Universidad Nacional de San Martín. Escuela de Ciencia y Tecnología; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Wilke, Claus O.. University of Texas at Austin; Estados Unido

Variations in the Golgi complex of mouse somatotrops at different times in a 24-hour period

El aparato de Golgi de las celulas somatotropas de la pars distalis de la hipofisis del raton, presenta importantes diferencias ultraestructurales en distintos momentos de un periodo de 24 horas. Es pequeno a mediodia y marcadamente hipertrofico a mediodia y 6 de la tarde.Facultad de Ciencias Médica

Sulfide-binding hemoglobins: Effects of mutations on active-site flexibility

The dynamics of Hemoglobin I (HbI) from the clam Lucina pectinata, from wild-type sperm whale (SW) myoglobin, and from the L29F/H64Q/V68F triple mutant of SW, both unligated and bound to hydrogen sulfide (H2S), have been studied in molecular dynamics simulations. Features that account for differences in H2S affinity among the three have been examined. Our results verify the existence of an unusual heme rocking motion in unligated HbI that can promote the entrance of large ligands such as H2S. The FQF-mutant partially reproduces the amplitude and relative orientation of the motion of HbI's heme group. Therefore, besides introducing favorable electrostatic interactions with H2S, the three mutations in the distal pocket change the dynamic properties of the heme group. The active-site residues Gln-64(E7), Phe-43(CD1), and His-93(F8) are also shown to be more flexible in unligated HbI than in FQF-mutant and SW. Further contributions to H2S affinity come from differences in hydrogen bonding between the heme propionate groups and nearby amino acid residues.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

Sulfide-binding hemoglobins: Effects of mutations on active-site flexibility

The dynamics of Hemoglobin I (HbI) from the clam Lucina pectinata, from wild-type sperm whale (SW) myoglobin, and from the L29F/H64Q/V68F triple mutant of SW, both unligated and bound to hydrogen sulfide (H2S), have been studied in molecular dynamics simulations. Features that account for differences in H2S affinity among the three have been examined. Our results verify the existence of an unusual heme rocking motion in unligated HbI that can promote the entrance of large ligands such as H2S. The FQF-mutant partially reproduces the amplitude and relative orientation of the motion of HbI's heme group. Therefore, besides introducing favorable electrostatic interactions with H2S, the three mutations in the distal pocket change the dynamic properties of the heme group. The active-site residues Gln-64(E7), Phe-43(CD1), and His-93(F8) are also shown to be more flexible in unligated HbI than in FQF-mutant and SW. Further contributions to H2S affinity come from differences in hydrogen bonding between the heme propionate groups and nearby amino acid residues.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada

The pars distalis of the pituitary as a source of liver growth stimulating factors during liver regeneration

Plasma and pituitary extracts from hepatectomized mice as well as pure growth hormone, stimulate DNA synthesis in hepatocyte population of adult intact mouse liver. This effect is not observed with plasma and pituitary extracts from intact or sham operated mice. The corresponding effects on the litoral cell population of the liver follow the same pattern, but the results are not so clearcut as in the hepatocyte population. The possibility that growth hormone, released by STH cells of the pituitary, is one of the factors responsible for the stimulating effect of plasma during liver regeneration, is discussed on the basis of the present results and of those already reported in the literature.Facultad de Ciencias Médica

Evolutionary Conservation of Protein Backbone Flexibility

Internal protein dynamics is essential for biological function. During evolution, protein divergence is functionally constrained: properties more relevant for function vary more slowly than less important properties. Thus, if protein dynamics is relevant for function, it should be evolutionary conserved. In contrast with the well-studied evolution of protein structure, the evolutionary divergence of protein dynamics has not been addressed systematically before, apart from a few case studies. X-Ray diffraction analysis gives information not only on protein structure but also on B-factors, which characterize the flexibility that results from protein dynamics. Here we study the evolutionary divergence of protein backbone dynamics by comparing the Cα flexibility (B-factor) profiles for a large dataset of homologous proteins classified into families and superfamilies. We show that Cα flexibility profiles diverge slowly, so that they are conserved at family and superfamily levels, even for pairs of proteins with nonsignificant sequence similarity. We also analyze and discuss the correlations among the divergences of flexibility, sequence, and structure.Instituto de Investigaciones Fisicoquímicas Teóricas y Aplicada