Distribution of amino acids in functional sites of proteins with high melting temperature

The stability of proteins in its native state has an important implication on its function and evolution. The functional site analysis may lead to better understanding of how these amino acid distributions influence the melting temperature of proteins. It has been reported that increasing the fraction of hydrophobic contacts in a protein tends to raise melting temperature; increasing the fraction of repulsive charge contacts decrease the melting temperature and consistent with a destabilizing effect. The role of amino acid distribution as hydrophobic, charged and polar residues in proteins and mainly in its functional sites has been studied. Due to limited data availability, redundancy check and controlled environment parameters, the study was carried out with ten single chain-wild proteins having melting temperature above 80°C at pH 7. The analysis depicts that, the entire protein, hydrophobic residues are more frequent in single chain proteins and charged residues are more frequent in multi-chains proteins. In functional sites of these proteins, hydrophobic and charged residues are equally frequent in single chain proteins and charged residues are very high in multi-chains proteins. But, the polar residue distribution remains same for both single chain and multi-chain proteins and its functional sites.     

Spatial proximity statistics suggest a regulatory role of protein phosphorylation on compound binding

Phosphorylation is an important post‐translational modification that regulates protein function by the attachment of negatively charged phosphate groups to phosphorylatable amino acid residues. As a mode of action, an influence of phosphorylation on the binding of compounds to proteins has been discussed and described for a number of proteins in the literature. However, a systematic statistical survey probing for enriched phosphorylation sites close to compound binding sites in support of this notion and with properly chosen random reference distributions has not been presented yet. Using high‐resolution protein structures from the Protein Data Bank including their co‐crystallized non‐covalently bound compounds and experimentally determined phosphorylation sites, we analyzed the pairwise distance distributions of phosphorylation and compound binding sites on protein surfaces. We found that phosphorylation sites are indeed located at significantly closer distances to compounds than expected by chance holding true specifically also for the subset of compound binding sites serving as catalytic sites of metabolic reactions. This tendency was particularly evident when treating phosphorylation sites as collective sets supporting the relevance of phosphorylation hotspots. Interestingly, phosphorylation sites were found to be closer to negatively charged than to positively charged compounds suggesting a stronger modulation of the binding of negatively charged compounds in dependence on phosphorylation status than on positively charged compounds. The enrichment of phosphorylation sites near compound binding sites confirms a regulatory role of phosphorylation in compound binding and provides a solid statistical basis for the literature‐reported selected events. Proteins 2016; 84:565–579. © 2016 Wiley Periodicals, Inc.     

Effective activation energy of enzymatic and nonenzymatic reactions. Evolution-imposed requirements to enzyme structure

The difference of the activation energies in a protein globule and water has been treated in terms of the theory of an elementary act of charge transfer reaction with regards to the energy spent on the transfer of charged reactants from water into the protein. The protein was treated as a structureless dielectric with a given optical and static dielectric constants surrounded by the aqueous phase. Reactions of different types (charge exchange between reactants, charge separation, neutralization, etc.) have been analyzed both under prevalence of purely electrostatic effects and under considerable nonelectrostatic contributions to the activation energies. It is shown that for all one-electron and most multi-electron reactions involving two reaction centres the energy spent for charged reactant transfer from water into protein is greater than the concomitant activation energy gain. The same effect takes place in a number of cases for multi-centre processes as well. To overcome the entropy hindrances, the reactants and catalysts must combine into multiparticle complexes, i.e. form microscopic regions of low dielectric constant. This results in increased effective activation energy as compared to reactions in water. It has been hypothesized that in order to make up for this loss the evolution has selected the proteins which are characterized by considerable intraglobular permanent electric fields. The presence in proteins of high concentrations of strongly polar peptide groups renders them advantageous in this respect over other polymers that are less polar.     

Spatial sign-alternating charge clusters in globular proteins

Large sign-alternating charge clusters formed by the charged side groups of amino acid residues and N- and C-terminal groups were found in the majority of considered globular proteins, namely 235 in a total of 274 protein structures, i.e. 85.8%. The clusters were determined by the criteria proposed earlier: charged groups were included in the cluster if their charged N and O atoms were located at distances between 2.4 and 7.0 A. The set of selected proteins consisted of known non-homologous protein structures from the Protein Data Bank with a resolution less than or equal to 2.5 A and pair sequence similarity less than 25%. Molecular masses of the proteins were from 5.5 to 91.5 kDa and protein chain length from 50 to 830 residues. The distribution of charged groups on the protein surface between isolated charged groups, small clusters with two and three groups, and large clusters with four or more groups were found to be approximately similar making 33, 35 and 32% of the total amount of protein charged groups, respectively. The large sign-alternating charge clusters with four or more charged groups were studied in greater detail. The amount of such clusters depends on the protein chain length. The small proteins contain 1-3 clusters while the large proteins display 4-6 or more clusters. On average, 1.5 clusters per each 100 residues were observed. In contrast with this, the size of a cluster, i.e. the number of charged groups inside a cluster, does not depend on the protein molecular mass, and large clusters are observed for proteins from a range of molecular masses. Clusters consisting of four to six charged groups occur most frequently, although extra large clusters are also often revealed. We can conclude that sign-alternating charge clusters are a common feature of the protein surface of globular protein. They are suggested to play a general functional role as a local polar factor of protein surface.     

Extreme differences in charge changes during protein evolution

Summary The maintenance of a proper distribution of charged amino acid residues might be expected to be an important factor in protein evolution. We therefore compared the inferred changes in charge during the evolution of 43 protein families with the changes expected on the basis of random base substitutions. It was found that certain proteins, like the eye lens crystallins and most histones, display an extreme avoidance of changes in charge. Other proteins, like phospholipase A2 and ferredoxin, apparently have sustained more charged replacements than expected, suggesting a positive selection for changes in charge. Depending on function and structure of a protein, charged residues apparently can be important targets for selective forces in protein evolution. It appears that actual biased codon usage tends to decrease the proportion of charged amino acid replacements. The influence of nonrandomness of mutations is more equivocal. Genes that use the mitochondrial instead of the universal code lower the probability that charge changes will occur in the encoded proteins.     

Intermediate filament structure: 1. Analysis of IF protein sequence data

Amino acid sequence data for intermediate filament proteins have been analysed with a view to identifying structurally invariant segments and determining their likely secondary structure. The sequences in these segments have also been analysed for periodic distributions of particular types of residue. The results support the classification of intermediate filament proteins into three main groups and also reinforce the concept of a molecular structure with a central domain of coiled-coil segments, together with essentially non-helical N-terminal and C-terminal domains of variable size and composition. Regions exhibiting the greatest homology between the three types of IF chain are identified and significant variation in charged residue disposition along the length of individual chains is noted. The conservation in all IF protein chains of specific sites of coiled-coil rope interruption are discussed in terms of the probable molecular structure. Stabilizing ionic interactions between coiled-coil chain segments have been investigated quantitatively as a function of the relative chain stagger. In all cases and calculations favour ropes in which the constituent chains are in-register and parallel rather than antiparallel.     

The structure of interfaces between subunits of dimeric and tetrameric proteins

The structures of the interfaces of nine dimeric and nine tetrameric proteins have been analyzed and have been seen to follow general principles. These interfaces are combinations of four structural motifs, which resemble features of monomeric proteins. These are: (i) extended beta sheet; (ii) helix-helix packing; (iii) sheet-sheet packing; and (iv) loop interactions. Other common structural features in the interfaces studied are two-fold symmetry, charged hydrogen bonds and channel formation (found only in tetramers). Monomer-monomer interfaces are intermediate in hydrophobicity and charge between the interfaces between secondary structures of monomeric proteins and the exteriors of monomeric proteins. A typical interface has one of the first three of the structural motifs at its centre and loop interactions around the outside, where most of the charge resides.     

Electrostatic Potential of B-DNA: Effect of Interionic Correlations

Modified Poisson-Boltzmann (MPB) equations have been numerically solved to study ionic distributions and mean electrostatic potentials around a macromolecule of arbitrarily complex shape and charge distribution. Results for DNA are compared with those obtained by classical Poisson-Boltzmann (PB) calculations. The comparisons were made for 1:1 and 2:1 electrolytes at ionic strengths up to 1 M. It is found that ion-image charge interactions and interionic correlations, which are neglected by the PB equation, have relatively weak effects on the electrostatic potential at charged groups of the DNA. The PB equation predicts errors in the long-range electrostatic part of the free energy that are only ∼1.5 kJ/mol per nucleotide even in the case of an asymmetrical electrolyte. In contrast, the spatial correlations between ions drastically affect the electrostatic potential at significant separations from the macromolecule leading to a clearly predicted effect of charge overneutralization.     

FG nucleoporins feature unique patterns that distinguish them from other IDPs

FG nucleoporins (FG Nups) are intrinsically disordered proteins and are the putative regulators of nucleocytoplasmic transport. They allow fast, yet selective, transport of molecules through the nuclear pore complex, but the underlying mechanism of nucleocytoplasmic transport is not yet fully discovered. As a result, FG Nups have been the subject of extensive research in the past two decades. Although most studies have been focused on analyzing the conformation and function of FG Nups from a biophysical standpoint, some recent studies have investigated the sequence-function relationship of FG Nups, with a few investigating amino acid sequences of a large number of FG Nups to understand common characteristics that might enable their function. Previously, we identified an evolutionarily conserved feature in FG Nup sequences, which are extended subsequences with low charge density, containing only positive charges, and located toward the N-terminus of FG Nups. We named these patterns longest positive like charge regions (lpLCRs). These patterns are specific to positively charged residues, and negatively charged residues do not demonstrate such a pattern. In this study, we compare FG Nups with other disordered proteins obtained from the DisProt and UniProt database in terms of presence of lpLCRs. Our results show that the lpLCRs are virtually exclusive to FG Nups and are not observed in other disordered proteins. Also, lpLCRs are what differentiate FG Nups from DisProt proteins in terms of charge distribution, meaning that excluding lpLCRs from the sequences of FG Nups make them similar to DisProt proteins in terms of charge distribution. We also previously showed the biophysical effect of lpLCRs in conformation of FG Nups. The results of this study are in line with our previous findings and imply that lpLCRs are virtually exclusive and functionally significant characteristics of FG Nups and nucleocytoplasmic transport.     

Nature of plasmalemmal functional “hemichannels”

The molecular identity of the protein forming “hemichannels” at non-junctional membranes is disputed. The family of gap junction proteins, innexins, connexins, and pannexins share several common features, including permeability characteristics and sensitivity to blocking agents. Such overlap in properties renders the identification of which of these protein species actually establishes the non-junctional membrane conductance and permeability quite complicated, especially because in vertebrates pannexins and connexins have largely overlapping distributions in tissues. Recently, attempts to establish criteria to identify events that are “hemichannel” mediated and those to allow the distinction between connexin- from pannexin-mediated events have been proposed.Here, I present an update on that topic and discuss the most recent findings related to the nature of functional “hemichannels” focusing on connexin43 and pannexin1. This article is part of a Special Issue entitled: The Communicating junctions, composition, structure and characteristics.     

The origin and impact of bound water around intrinsically disordered proteins

Proteins and water couple dynamically over a wide range of time scales. Motivated by their central role in protein function, protein-water dynamics and thermodynamics have been extensively studied for structured proteins, where correspondence to structural features has been made. However, properties controlling intrinsically disordered protein (IDP)-water dynamics are not yet known. We report results of megahertz-to-terahertz dielectric spectroscopy and molecular dynamics simulations of a group of IDPs with varying charge content along with structured proteins of similar size. Hydration water around IDPs is found to exhibit more heterogeneous rotational and translational dynamics compared with water around structured proteins of similar size, yielding on average more restricted dynamics around individual residues of IDPs, charged or neutral, compared with structured proteins. The on-average slower water dynamics is found to arise from excess tightly bound water in the first hydration layer, which is related to greater exposure to charged groups. The more tightly bound water to IDPs correlates with the smaller hydration shell found experimentally, and affects entropy associated with protein-water interactions, the contribution of which we estimate based on the dielectric measurements and simulations. Water-IDP dynamic coupling at terahertz frequencies is characterized by the dielectric measurements and simulations.     

Real-Time Nanopore-Based Recognition of Protein Translocation Success

A growing number of new technologies are supported by a single- or multi-nanopore architecture for capture, sensing, and delivery of polymeric biomolecules. Nanopore-based single-molecule DNA sequencing is the premier example. This method relies on the uniform linear charge density of DNA, so that each DNA strand is overwhelmingly likely to pass through the nanopore and across the separating membrane. For disordered peptides, folded proteins, or block copolymers with heterogeneous charge densities, by contrast, translocation is not assured, and additional strategies to monitor the progress of the polymer molecule through a nanopore are required. Here, we demonstrate a single-molecule method for direct, model-free, real-time monitoring of the translocation of a disordered, heterogeneously charged polypeptide through a nanopore. The crucial elements are two “selectivity tags”—regions of different but uniform charge density—at the ends of the polypeptide. These affect the selectivity of the nanopore differently and enable discrimination between polypeptide translocation and retraction. Our results demonstrate exquisite sensitivity of polypeptide translocation to applied transmembrane potential and prove the principle that nanopore selectivity reports on biopolymer substructure. We anticipate that the selectivity tag technique will be broadly applicable to nanopore-based protein detection, analysis, and separation technologies, and to the elucidation of protein translocation processes in normal cellular function and in disease.     

Buried charged surface in proteins

BACKGROUND: The traditional picture of charged amino acids in globular proteins is that they are almost exclusively on the outside exposed to the solvent. Buried charges, when they do occur, are assumed to play an essential role in catalysis and ligand binding, or in stabilizing structure as, for instance, helix caps. RESULTS: By analyzing the amount and distribution of buried charged surface and charges in proteins over a broad range of protein sizes, we show that buried charge is much more common than is generally believed. We also show that the amount of buried charge rises with protein size in a manner which differs from other types of surfaces, especially aromatic and polar uncharged surfaces. In large proteins such as hemocyanin, 35% of all charges are greater than 75% buried. Furthermore, at all sizes few charged groups are fully exposed. As an experimental test, we show that replacement of the buried D178 of muconate lactonizing enzyme by N stabilizes the enzyme by 4.2 degrees C without any change in crystallographic structure. In addition, free energy calculations of stability support the experimental results. CONCLUSIONS: Nature may use charge burial to reduce protein stability; not all buried charges are fully stabilized by a prearranged protein environment. Consistent with this view, thermophilic proteins often have less buried charge. Modifying the amount of buried charge at carefully chosen sites may thus provide a general route for changing the thermophilicity or psychrophilicity of proteins.     

Charge-Rich Regions Modulate the Anti-Aggregation Activity of Hsp90

Protein aggregation can have dramatic effects on cellular function and plays a causative role in many human diseases. In all cells, molecular chaperones bind to aggregation-prone proteins and hinder aggregation. The ability of a protein to resist aggregation and remain soluble in aqueous solution is linked to the physical properties of the protein. Numerous physical studies demonstrate that charged atoms favor solubility. We note that many molecular chaperones possess a substantial negative charge that may allow them to impart solubility on aggregation-prone proteins. Hsp90 is one such negatively charged molecular chaperone. The charge on Hsp90 is largely concentrated in two highly acidic regions. To investigate the relationship between chaperone charge and protein solubility, we deleted these charge-rich regions and analyzed the resulting Hsp90 constructs for anti-aggregation activity. We found that deletion of both charge-rich regions dramatically impaired Hsp90 anti-aggregation activity. The anti-aggregation role of the deleted charge-rich regions could be due to net charge or sequence-specific features. To distinguish these possibilities, we attached an acid-rich region with a distinct amino acid sequence to our double-deleted Hsp90 construct. This charge rescue construct displayed effective anti-aggregation activity indicating that the net charge of Hsp90 contributes to its anti-aggregation activity.     

Charges on proteins and distances of electron transfer in metalloprotein redox reactions

A modified form of the Debye-Marcus equation relating electron transfer rate constants to charges on proteins and distances of electron transfer has been applied to the reaction of chemically modified cytochrome f, in which positively charged amino groups are replaced with negatively charged carboxyl groups. The rate of electron transfer from reduced cytochrome f to ferricyanide decreased with increasing ionic strength when the native and singly substituted cytochrome f were used, although a sharp decrease was observed in the former case. When doubly or more than triply substituted cytochrome f was used, the rate of electron transfer was almost constant or increased with increasing ionic strength, respectively. The kinetic-ionic strength effects on this reaction can be well explained by the Debye-Marcus equation in which the charge and radius of the protein are treated as variable parameters. The results show the importance of local positive charges of about 2.0 on native cytochrome f and effective radius of about 11 A of cytochrome f for the electron transfer to ferricyanide. Since the net charge on the native cytochrome f is negative and the calculated radius of the protein is 22.8 A, the above results indicate that positive charges on the electron transfer site control the electrostatic interactions in this reaction. Previously reported data which had been analyzed by using the total net charge and full radius of the protein, were also well explained by the local charge and effective radius of the protein.(ABSTRACT TRUNCATED AT 250 WORDS)     

Moments of the distribution of the amino acid residues in tertiary structures of globular proteins

Moments of the distributions of the C^α and “side-chain atoms” and associated properties were examined in 22 globular proteins, considered as statistical aggregates of atoms. Although the distributions are generally anisotropic, the densities of the evaluated distributions are highly uniform in the interior of a single protein, as well as among the proteins investigated. The tertiary structure of proteins is characterized by a compact and uniform distribution of amino acids, independent of their molecular weight and hydrophobic character, and by an isotropic distribution of the virtual bond directions in the polypeptide folding. While the general uniformity of the density of distributions in the bulk of proteins can be justified by the architectural requirements of high thermodynamic stability, significant differences in the distribution of the C^α with respect to the “side-chain atoms” suggest a plausible explanation of the general anisotropic morphology of the proteins. The invariance of the density of distributions allows easy recognition of proteinlike domains in more complex proteins and suggests a practical way to predict the following path in single proteins.     

Intrinsically disordered regions have specific functions in mitochondrial and nuclear proteins

Proteins in general consist not only of globular structural domains (SDs), but also of intrinsically disordered regions (IDRs), i.e. those that do not assume unique three-dimensional structures by themselves. Although IDRs are especially prevalent in eukaryotic proteins, the functions are mostly unknown. To elucidate the functions of IDRs, we first divided eukaryotic proteins into subcellular localizations, identified IDRs by the DICHOT system that accurately divides entire proteins into SDs and IDRs, and examined charge and hydropathy characteristics. On average, mitochondrial proteins have IDRs more positively charged than SDs. Comparison of mitochondrial proteins with orthologous prokaryotic proteins showed that mitochondrial proteins tend to have segments attached at both N and C termini, high fractions of which are IDRs. Segments added to the N-terminus of mitochondrial proteins contain not only signal sequences but also mature proteins and exhibit a positive charge gradient, with the magnitude increasing toward the N-terminus. This finding is consistent with the notion that positively charged residues are added to the N-terminus of proteobacterial proteins so that the extended proteins can be chromosomally encoded and efficiently transported to mitochondria after translation. By contrast, nuclear proteins generally have positively charged SDs and negatively charged IDRs. Among nuclear proteins, DNA-binding proteins have enhanced charge tendencies. We propose that SDs in nuclear proteins tend to be positively charged because of the need to bind to negatively charged nucleotides, while IDRs tend to be negatively charged to interact with other proteins or other regions of the same proteins to avoid premature proteasomal degradation.