A coevolutionary residue network at the site of a functionally important conformational change in a phosphohexomutase enzyme family

Coevolution analyses identify residues that co-vary with each other during evolution, revealing sequence relationships unobservable from traditional multiple sequence alignments. Here we describe a coevolutionary analysis of phosphomannomutase/phosphoglucomutase (PMM/PGM), a widespread and diverse enzyme family involved in carbohydrate biosynthesis. Mutual information and graph theory were utilized to identify a network of highly connected residues with high significance. An examination of the most tightly connected regions of the coevolutionary network reveals that most of the involved residues are localized near an interdomain interface of this enzyme, known to be the site of a functionally important conformational change. The roles of four interface residues found in this network were examined via site-directed mutagenesis and kinetic characterization. For three of these residues, mutation to alanine reduces enzyme specificity to ～10% or less of wild-type, while the other has ～45% activity of wild-type enzyme. An additional mutant of an interface residue that is not densely connected in the coevolutionary network was also characterized, and shows no change in activity relative to wild-type enzyme. The results of these studies are interpreted in the context of structural and functional data on PMM/PGM. Together, they demonstrate that a network of coevolving residues links the highly conserved active site with the interdomain conformational change necessary for the multi-step catalytic reaction. This work adds to our understanding of the functional roles of coevolving residue networks, and has implications for the definition of catalytically important residues.     

Molecular Determinants Underlying Binding Specificities of the ABL Kinase Inhibitors: Combining Alanine Scanning of Binding Hot Spots with Network Analysis of Residue Interactions and Coevolution

Quantifying binding specificity and drug resistance of protein kinase inhibitors is of fundamental importance and remains highly challenging due to complex interplay of structural and thermodynamic factors. In this work, molecular simulations and computational alanine scanning are combined with the network-based approaches to characterize molecular determinants underlying binding specificities of the ABL kinase inhibitors. The proposed theoretical framework unveiled a relationship between ligand binding and inhibitor-mediated changes in the residue interaction networks. By using topological parameters, we have described the organization of the residue interaction networks and networks of coevolving residues in the ABL kinase structures. This analysis has shown that functionally critical regulatory residues can simultaneously embody strong coevolutionary signal and high network centrality with a propensity to be energetic hot spots for drug binding. We have found that selective (Nilotinib) and promiscuous (Bosutinib, Dasatinib) kinase inhibitors can use their energetic hot spots to differentially modulate stability of the residue interaction networks, thus inhibiting or promoting conformational equilibrium between inactive and active states. According to our results, Nilotinib binding may induce a significant network-bridging effect and enhance centrality of the hot spot residues that stabilize structural environment favored by the specific kinase form. In contrast, Bosutinib and Dasatinib can incur modest changes in the residue interaction network in which ligand binding is primarily coupled only with the identity of the gate-keeper residue. These factors may promote structural adaptability of the active kinase states in binding with these promiscuous inhibitors. Our results have related ligand-induced changes in the residue interaction networks with drug resistance effects, showing that network robustness may be compromised by targeted mutations of key mediating residues. This study has outlined mechanisms by which inhibitor binding could modulate resilience and efficiency of allosteric interactions in the kinase structures, while preserving structural topology required for catalytic activity and regulation.     

Integrated analysis of residue coevolution and protein structure in ABC transporters

Intraprotein side chain contacts can couple the evolutionary process of amino acid substitution at one position to that at another. This coupling, known as residue coevolution, may vary in strength. Conserved contacts thus not only define 3-dimensional protein structure, but also indicate which residue-residue interactions are crucial to a protein's function. Therefore, prediction of strongly coevolving residue-pairs helps clarify molecular mechanisms underlying function. Previously, various coevolution detectors have been employed separately to predict these pairs purely from multiple sequence alignments, while disregarding available structural information. This study introduces an integrative framework that improves the accuracy of such predictions, relative to previous approaches, by combining multiple coevolution detectors and incorporating structural contact information. This framework is applied to the ABC-B and ABC-C transporter families, which include the drug exporter P-glycoprotein involved in multidrug resistance of cancer cells, as well as the CFTR chloride channel linked to cystic fibrosis disease. The predicted coevolving pairs are further analyzed based on conformational changes inferred from outward- and inward-facing transporter structures. The analysis suggests that some pairs coevolved to directly regulate conformational changes of the alternating-access transport mechanism, while others to stabilize rigid-body-like components of the protein structure. Moreover, some identified pairs correspond to residues previously implicated in cystic fibrosis.     

Exploring a phylogenetic approach for the detection of correlated substitutions in proteins

The remarkable conservation of protein structure, compared with that of sequences, suggests that in the course of evolution, residue substitutions which tend to destabilize a particular structure must be compensated by other substitutions that confer greater stability on that structure. Several approaches have been designed to detect correlated changes in a set of homologous sequences. However, most of them do not take into account the phylogeny of the sequences, and it has been shown that their detection power is weak. It remains unclear whether coevolution could be a general process at the level of amino acids of proteins. In the present study, we analyze the phylogenetic reconstruction of 15 sets of homologous proteins to assess, under different conditions, whether a significant amount of coevolving sites can be detected. Two criteria are used to detect significantly cosubstituting sites. One criterion corresponds to that of Shindyalov, Kolchanov, and Sander. The second one is based on intensive simulations of evolution of protein sequences along a phylogeny to estimate the significance of the number of observed cosubstitutions for pairs of sites. Our results show an important sensitivity of the detection of cosubstituting sites to the model used for the phylogenetic reconstruction. Not considering the uncertainty associated with the reconstructed data might lead to detecting numerous false-positive pairs of sites. Finally, significant amounts of coevolving pairs could be found only when substitutions affecting the physicochemical properties of the amino acids were considered. Such results suggest evidence of a cosubstitution mechanism in protein evolution. However, the identification of nonambiguous coevolving sites is still unresolved.     

Large-Scale Conformational Transitions and Dimerization Are Encoded in the Amino-Acid Sequences of Hsp70 Chaperones

Hsp70s are a class of ubiquitous and highly conserved molecular chaperones playing a central role in the regulation of proteostasis in the cell. Hsp70s assist a myriad of cellular processes by binding unfolded or misfolded substrates during a complex biochemical cycle involving large-scale structural rearrangements. Here we show that an analysis of coevolution at the residue level fully captures the characteristic large-scale conformational transitions of this protein family, and predicts an evolutionary conserved–and thus functional–homo-dimeric arrangement. Furthermore, we highlight that the features encoding the Hsp70 dimer are more conserved in bacterial than in eukaryotic sequences, suggesting that the known Hsp70/Hsp110 hetero-dimer is a eukaryotic specialization built on a pre-existing template.     

Structure of the prion Ure2p in protein fibrils assembled in vitro

The Ure2 protein from the yeast Saccharomyces cerevisiae has prion properties. In vitro and at neutral pH, soluble Ure2p spontaneously forms long, straight, insoluble protein fibrils. Two models have been proposed to account for the assembly of Ure2p into protein fibrils. The "amyloid backbone" model postulates that a segment ranging from 40 to 70 amino acids in the flexible N-terminal domain from different Ure2p molecules forms a parallel superpleated beta-structure running along the fibrils. The second model hypothesizes that assembly of full-length Ure2p is driven by limited conformational rearrangements and non-native inter- and/or intramolecular interactions between Ure2p monomers. Here, we performed a cysteine scan on residues located in the N- and C-terminal parts of Ure2p to determine whether these domains interact. Amino acid sequences centered around residue 6 in the N-terminal domain of Ure2p and residue 137 in the C-terminal moiety interacted at least transiently via intramolecular interactions. We documented the assembly properties of a Ure2p variant in which a disulfide bond was established between the N- and C-terminal domains and showed that it possesses assembly properties indistinguishable from those of wild-type Ure2p. We probed the structure of Ure2pC6C137 within the fibrils and demonstrate that the polypeptide is in a conformation similar to that of its soluble assembly-competent state. Our results constitute the first structural characterization of the N-terminal domain of Ure2p in both its soluble assembly-competent and fibrillar forms. Our data indicate that the flexibility of the N-terminal domain and conformational changes within this domain are essential for fibril formation and provide new insight into the conformational rearrangements that lead to the assembly of Ure2p into fibrils and the propagation of the [URE3] phenotype in yeast.     

Immobilized prion protein undergoes spontaneous rearrangement to a conformation having features in common with the infectious form

It is hypothesized that infectious prions are generated as the cellular form of the prion protein (PrP(C)) undergoes pronounced conformational change under the direction of an infectious PrP(Sc) template. Conversion to the infectious conformer is particularly associated with major structural rearrangement in the central portion of the protein (residues 90-120), which has an extended flexible structure in the PrP(C) isoform. Using a panel of recombinant antibodies reactive with different parts of PrP, we show that equivalent major structural rearrangements occur spontaneously in this region of PrP immobilized on a surface. In contrast, regions more towards the termini of the protein remain relatively unaltered. The rearrangements occur even under conditions where individual PrP molecules should not contact one another. The propensity of specific unstructured regions of PrP to spontaneously undergo large and potentially deleterious conformational changes may have important implications for prion biology.     

Photo-CIDNP nuclear magnetic resonance as a probe for conformational changes in epidermal growth factor

Photo-chemically induced dynamic nuclear polarization (CIDNP)-NMR spectroscopy at 360 MHz has been used to investigate pH-induced conformational transitions in mouse epidermal growth factor. At about pH 9, all five tyrosine residues and both tryptophan residues are, to various extents, solvent-exposed, while the His-22 residue is buried in the protein matrix. Tyr-13 is the least exposed of the tyrosine residues and also the most immobilized. As the pH is decreased to 5.9, the tryptophan residues gradually become less exposed, while the Tyr-13 residue becomes internalized in the protein. These data suggest that the C-terminus and part of the N-terminal structural domain are affected by a conformational transition in mouse epidermal growth factor occurring between pH 6 and 8 via breakage of the His-22 inter-residue linkage. Above pH 9, a decreased photo-CIDNP effect is evident for both tryptophans and for Tyr-10 and Try-13; this information suggests that a second conformational change takes place at basic pH, which may simply be incipient denaturation.     

Side-chain flexibility in proteins upon ligand binding

Ligand binding may involve a wide range of structural changes in the receptor protein, from hinge movement of entire domains to small side-chain rearrangements in the binding pocket residues. The analysis of side chain flexibility gives insights valuable to improve docking algorithms and can provide an index of amino-acid side-chain flexibility potentially useful in molecular biology and protein engineering studies. In this study we analyzed side-chain rearrangements upon ligand binding. We constructed two non-redundant databases (980 and 353 entries) of "paired" protein structures in complexed (holo-protein) and uncomplexed (apo-protein) forms from the PDB macromolecular structural database. The number and identity of binding pocket residues that undergo side-chain conformational changes were determined. We show that, in general, only a small number of residues in the pocket undergo such changes (e.g., approximately 85% of cases show changes in three residues or less). The flexibility scale has the following order: Lys > Arg, Gln, Met > Glu, Ile, Leu > Asn, Thr, Val, Tyr, Ser, His, Asp > Cys, Trp, Phe; thus, Lys side chains in binding pockets flex 25 times more often then do the Phe side chains. Normalizing for the number of flexible dihedral bonds in each amino acid attenuates the scale somewhat, however, the clear trend of large, polar amino acids being more flexible in the pocket than aromatic ones remains. We found no correlation between backbone movement of a residue upon ligand binding and the flexibility of its side chain. These results are relevant to 1. Reduction of search space in docking algorithms by inclusion of side-chain flexibility for a limited number of binding pocket residues; and 2. Utilization of the amino acid flexibility scale in protein engineering studies to alter the flexibility of binding pockets.     

Integrated Analysis of Residue Coevolution and Protein Structures Capture Key Protein Sectors in HIV-1 Proteins

HIV type 1 (HIV-1) is characterized by its rapid genetic evolution, leading to challenges in anti-HIV therapy. However, the sequence variations in HIV-1 proteins are not randomly distributed due to a combination of functional constraints and genetic drift. In this study, we examined patterns of sequence variability for evidence of linked sequence changes (termed as coevolution or covariation) in 15 HIV-1 proteins. It shows that the percentage of charged residues in the coevolving residues is significantly higher than that in all the HIV-1 proteins. Most of the coevolving residues are spatially proximal in the protein structures and tend to form relatively compact and independent units in the tertiary structures, termed as “protein sectors”. These protein sectors are closely associated with anti-HIV drug resistance, T cell epitopes, and antibody binding sites. Finally, we explored candidate peptide inhibitors based on the protein sectors. Our results can establish an association between the coevolving residues and molecular functions of HIV-1 proteins, and then provide us with valuable knowledge of pathology of HIV-1 and therapeutics development.     

The Structural Determinants of Intra-Protein Compensatory Substitutions

Compensatory substitutions happen when one mutation is advantageously selected because it restores the loss of fitness induced by a previous deleterious mutation. How frequent such mutations occur in evolution and what is the structural and functional context permitting their emergence remain open questions. We built an atlas of intra-protein compensatory substitutions using a phylogenetic approach and a dataset of 1,630 bacterial protein families for which high-quality sequence alignments and experimentally derived protein structures were available. We identified more than 51,000 positions coevolving by the mean of predicted compensatory mutations. Using the evolutionary and structural properties of the analyzed positions, we demonstrate that compensatory mutations are scarce (typically only a few in the protein history) but widespread (the majority of proteins experienced at least one). Typical coevolving residues are evolving slowly, are located in the protein core outside secondary structure motifs, and are more often in contact than expected by chance, even after accounting for their evolutionary rate and solvent exposure. An exception to this general scheme is residues coevolving for charge compensation, which are evolving faster than noncoevolving sites, in contradiction with predictions from simple coevolutionary models, but similar to stem pairs in RNA. While sites with a significant pattern of coevolution by compensatory mutations are rare, the comparative analysis of hundreds of structures ultimately permits a better understanding of the link between the three-dimensional structure of a protein and its fitness landscape.     

Characterization of the reversible conformational equilibrium of the cytoplasmic domain of erythrocyte membrane band 3

The cytoplasmic domain of the erythrocyte membrane protein, band 3, contains binding sites for hemoglobin, several glycolytic enzymes, and ankyrin, the linkage to the cytoskeleton. In an earlier study, we found evidence which suggested that band 3 might undergo a native conformational change. We demonstrate here that the cytoplasmic domain of band 3 does exist in a reversible, pH-dependent conformational equilibrium among 3 native states. At physiological salt concentrations this equilibrium is characterized by apparent pKa values of 7.2 and 9.2; however, these apparent pKa values change if the domain's sulfhydryl groups are modified. A major component of the structural change appears to involve the pivoting of two subdomains of the cytoplasmic domain at a central hinge, as evidenced by both hydrodynamic and fluorescence energy transfer measurements. The probable site of this hinge is between residues 176 and 191, a region highly accessible to proteases and also rich in proline. These structural rearrangements also apparently extend to the cluster of tryptophan residues near the N terminus, since the domain's intrinsic fluorescence more than doubles between pH 6.5 and 9.5. No measurable change in band 3 secondary or quaternary structure could be detected during the conformational transitions. A structural model of the cytoplasmic domain of band 3 is presented to show the possible spatial relationships between the regions of conformational change and the sites of peripheral protein binding.     

The contribution of coevolving residues to the stability of KDO8P synthase

Background

The evolutionary tree of 3-deoxy-D-manno-octulosonate 8-phosphate (KDO8P) synthase (KDO8PS), a bacterial enzyme that catalyzes a key step in the biosynthesis of bacterial endotoxin, is evenly divided between metal and non-metal forms, both having similar structures, but diverging in various degrees in amino acid sequence. Mutagenesis, crystallographic and computational studies have established that only a few residues determine whether or not KDO8PS requires a metal for function. The remaining divergence in the amino acid sequence of KDO8PSs is apparently unrelated to the underlying catalytic mechanism.

Methodology/Principal Findings

The multiple alignment of all known KDO8PS sequences reveals that several residue pairs coevolved, an indication of their possible linkage to a structural constraint. In this study we investigated by computational means the contribution of coevolving residues to the stability of KDO8PS. We found that about 1/4 of all strongly coevolving pairs probably originated from cycles of mutation (decreasing stability) and suppression (restoring it), while the remaining pairs are best explained by a succession of neutral or nearly neutral covarions.

Conclusions/Significance

Both sequence conservation and coevolution are involved in the preservation of the core structure of KDO8PS, but the contribution of coevolving residues is, in proportion, smaller. This is because small stability gains or losses associated with selection of certain residues in some regions of the stability landscape of KDO8PS are easily offset by a large number of possible changes in other regions. While this effect increases the tolerance of KDO8PS to deleterious mutations, it also decreases the probability that specific pairs of residues could have a strong contribution to the thermodynamic stability of the protein.     

pH-induced conformational changes in the soluble manganese-stabilizing protein of photosystem II

In this paper, we analyzed the pH-induced changes in the conformational states of the manganese-stabilizing protein (MSP) of photosystem II. Distinct conformational states of MSP were identified using fluorescence spectra, far-UV circular dichroism, and pressure-induced unfolding at varying suspension pH values, and four different conformational states of MSP were clearly distinguished using the center of fluorescence spectra mass when suspension pH was altered from 2 to 12. MSP was completely unfolded at a suspension pH above 11 and partly unfolded below a pH of 3. Analysis of the center of fluorescence spectral mass showed that the MSP structure appears stably folded around pH 6 and 4. The conformational state of MSP at pH 4 seems more stable than that at pH 6. Studies of peak positions of tryptophan fluorescence and MSP-bound 1-anilinonaphthalene-8-sulfonic acid fluorescence spectra supported this observation. A decrease in the suspension pH to 2 resulted in significant alterations in the MSP structure possibly because of protonation of unprotonated residues at lower pH, suggesting the existence of a large number of unprotonated amino acid residues at neutral pH possibly useful for proton transport in oxygen evolution. The acidic pH-induced conformational changes of MSP were reversible upon increase of pH to neutral pH; however, N-bromosuccinimide modification of tryptophan (Trp241) blocks the recovery of pH-induced conformational changes in MSP, implying that Trp241 is a key residue for the unfolded protein to form a functional structure. Thus, pH-induced structural changes of stable MSP (pH 6-4) may be utilized to analyze its functionality as a cofactor for oxygen evolution.     

Phenotypic Interactions among Mutations in a Thermus thermophilus 16S rRNA Gene Detected with Genetic Selections and Experimental Evolution

During protein synthesis, the ribosome undergoes conformational transitions between functional states, requiring communication between distant structural elements of the ribosome. Despite advances in ribosome structural biology, identifying the protein and rRNA residues governing these transitions remains a significant challenge. Such residues can potentially be identified genetically, given the predicted deleterious effects of mutations stabilizing the ribosome in discrete conformations and the expected ameliorating effects of second-site compensatory mutations. In this study, we employed genetic selections and experimental evolution to identify interacting mutations in the ribosome of the thermophilic bacterium Thermus thermophilus. By direct genetic selections, we identified mutations in 16S rRNA conferring a streptomycin dependence phenotype and from these derived second-site suppressor mutations relieving dependence. Using experimental evolution of streptomycin-independent pseudorevertants, we identified additional compensating mutations. Similar mutations could be evolved from slow-growing streptomycin-resistant mutants. While some mutations arose close to the site of the original mutation in the three-dimensional structure of the 30S ribosomal subunit and probably act directly by compensating for local structural distortions, the locations of others are consistent with long-range communication between specific structural elements within the ribosome.     

Chasing long-range evolutionary couplings in the AlphaFold era

Coevolution between protein residues is normally interpreted as direct contact. However, the evolutionary record of a protein sequence contains rich information that may include long-range functional couplings, couplings that report on homo-oligomeric states or even conformational changes. Due to the complexity of the sequence space and the lack of structural information on various members of a protein family, it has been difficult to effectively mine the additional information encoded in a multiple sequence alignment (MSA). Here, taking advantage of the recent release of the AlphaFold (AF) database we attempt to identify coevolutionary couplings that cannot be explained simply by spatial proximity. We propose a simple computational method that performs direct coupling analysis on a MSA and searches for couplings that are not satisfied in any of the AF models of members of the identified protein family. Application of this method on 2012 protein families suggests that ~12% of the total identified coevolving residue pairs are spatially distant and more likely to be disordered than their contacting counterparts. We expect that this analysis will help improve the quality of coevolutionary distance restraints used for structure determination and will be useful in identifying potentially functional/allosteric cross-talk between distant residues.     

The cytoplasmic isoenzyme of horse liver aldehyde dehydrogenase. Relationship to the corresponding human isoenzyme

The structural divergence between the cytoplasmic isoenzymes of aldehyde dehydrogenase from different species was investigated by analysis of peptides from the horse protein, and correlation of the results with the complete primary structure of the human isoenzyme. The amino acid sequences of these two proteins show a high degree of homology (91% of residues compared are identical). The differences observed are spread over the entire polypeptide chains, with only one cluster, which is close to a reactive cysteine residue and also adjacent to the most conserved region (covering 68 residues) in the primary structures of the whole enzymes. The secondary structure predicted for the human isoenzyme is mainly unaffected by the residue differences in the horse isoenzyme, although limited conformational changes might be compatible with an unexpected overrepresentation of differences involving isoleucine (12 of 43 exchanges represent a loss of Ile in the horse protein). Two cysteine residues that correlate with catalytic activity are identically positioned in the enzyme from the two species.     

Integrated prediction of one-dimensional structural features and their relationships with conformational flexibility in helical membrane proteins

Background  

Many structural properties such as solvent accessibility, dihedral angles and helix-helix contacts can be assigned to each residue in a membrane protein. Independent studies exist on the analysis and sequence-based prediction of some of these so-called one-dimensional features. However, there is little explanation of why certain residues are predicted in a wrong structural class or with large errors in the absolute values of these features. On the other hand, membrane proteins undergo conformational changes to allow transport as well as ligand binding. These conformational changes often occur via residues that are inherently flexible and hence, predicting fluctuations in residue positions is of great significance.