Conservation and diversity of ancient hemoglobins in Bacteria

A group of single-domain proteins in Bacteria similar to thermoglobin, an oxygen-avid hemoglobin representative of the ancestral form, reveals the primordial structure, function, and evolvability of the family. Conserved residues at specific positions function to bind ligand or participate in hydrophobic packing of the protein core during protein folding. A potential hydrogen bond network consisting of a tyrosine and glutamine residue in the distal ligand-binding site of most hemoglobins suggests that the ancestral protein bound oxygen avidly. Two divergent hemoglobins with mutations at generally conserved positions contain non-canonical ligand-binding sites, illustrating plasticity of the fold. One binds heme in a manner similar to cytochromes and may represent an evolutionary link to the precursor of the hemoglobin fold. Conservation suggests specific biochemical properties of the ancestral protein; diversity suggests an evolvability of this group of hemoglobins tolerant of mutations that perturb conserved biochemical properties for adaptation to novel functions.     

Evolutionary coupling of structural and functional sequence information in the intracellular lipid-binding protein family

We have mined the evolutionary record for the large family of intracellular lipid-binding proteins (iLBPs) by calculating the statistical coupling of residue variations in a multiple sequence alignment using methods developed by Ranganathan and coworkers (Lockless and Ranganathan, Science 1999:286;295-299). The 213 sequences analyzed have a wide range of ligand-binding functions as well as highly divergent phylogenetic origins, assuring broad sampling of sequence space. Emerging from this analysis were two major clusters of coupled residues, which when mapped onto the structure of a representative iLBP under study in our laboratory, cellular retinoic-acid binding protein I, are largely contiguous and provide useful points of comparison to available data for the folding of this protein. One cluster comprises a predominantly hydrophobic core away from the ligand-binding site and likely represents key structural information for the iLBP fold. The other cluster includes the portal region where ligand enters its binding site, regions of the ligand-binding cavity, and the region where the 10-stranded beta-barrel characteristic of this family closes (between strands 1' and 10). Linkages between these two clusters suggest that evolutionary pressures on this family constrain structural and functional sequence information in an interdependent fashion. The necessity of the structure to wrap around a hydrophobic ligand confounds the typical sequestration of hydrophobic side chains. Additionally, ligand entry and exit require these structures to have a capacity for specific conformational change during binding and release. We conclude that an essential and structurally apparent separation of local and global sequence information is conserved throughout the iLBP family.     

Signal Transduction at the Domain Interface of Prokaryotic Pentameric Ligand-Gated Ion Channels

Pentameric ligand-gated ion channels are activated by the binding of agonists to a site distant from the ion conduction path. These membrane proteins consist of distinct ligand-binding and pore domains that interact via an extended interface. Here, we have investigated the role of residues at this interface for channel activation to define critical interactions that couple conformational changes between the two structural units. By characterizing point mutants of the prokaryotic channels ELIC and GLIC by electrophysiology, X-ray crystallography and isothermal titration calorimetry, we have identified conserved residues that, upon mutation, apparently prevent activation but not ligand binding. The positions of nonactivating mutants cluster at a loop within the extracellular domain connecting β-strands 6 and 7 and at a loop joining the pore-forming helix M2 with M3 where they contribute to a densely packed core of the protein. An ionic interaction in the extracellular domain between the turn connecting β-strands 1 and 2 and a residue at the end of β-strand 10 stabilizes a state of the receptor with high affinity for agonists, whereas contacts of this turn to a conserved proline residue in the M2-M3 loop appear to be less important than previously anticipated. When mapping residues with strong functional phenotype on different channel structures, mutual distances are closer in conducting than in nonconducting conformations, consistent with a potential role of contacts in the stabilization of the open state. Our study has revealed a pattern of interactions that are crucial for the relay of conformational changes from the extracellular domain to the pore region of prokaryotic pentameric ligand-gated ion channels. Due to the strong conservation of the interface, these results are relevant for the entire family.     

Characterization of the ligand-binding site of the serotonin 5-HT3 receptor: the role of glutamate residues 97, 224, AND 235

Ligand-gated ion channels of the Cys loop family are receptors for small amine-containing neurotransmitters. Charged amino acids are strongly conserved in the ligand-binding domain of these receptor proteins. To investigate the role of particular residues in ligand binding of the serotonin 5-HT3AS receptor (5-HT3R), glutamate amino acid residues at three different positions, Glu97, Glu224, and Glu235, in the extracellular N-terminal domain were substituted with aspartate and glutamine using site-directed mutagenesis. Wild type and mutant receptor proteins were expressed in HEK293 cells and analyzed by electrophysiology, radioligand binding, fluorescence measurements, and immunochemistry. A structural model of the ligand-binding domain of the 5-HT3R based on the acetylcholine binding protein revealed the position of the mutated amino acids. Our results demonstrate that mutations of Glu97, distant from the ligand-binding site, had little effect on the receptor, whereas mutations Glu224 and Glu235, close to the predicted binding site, are indeed important for ligand binding. Mutations E224Q, E224D, and E235Q decreased EC50 and Kd values 5-20-fold, whereas E235D was functionally expressed at a low level and had a more than 100-fold increased EC50 value. Comparison of the fluorescence properties of a fluorescein-labeled antagonist upon binding to wild type 5-HT3R and E235Q, allowed us to localize Glu235 within a distance of 1 nm around the ligand-binding site, as proposed by our model.     

Evolutionary trace of G protein-coupled receptors reveals clusters of residues that determine global and class-specific functions

G protein-coupled receptor (GPCR) activation mediated by ligand-induced structural reorganization of its helices is poorly understood. To determine the universal elements of this conformational switch, we used evolutionary tracing (ET) to identify residue positions commonly important in diverse GPCRs. When mapped onto the rhodopsin structure, these trace residues cluster into a network of contacts from the retinal binding site to the G protein-coupling loops. Their roles in a generic transduction mechanism were verified by 211 of 239 published mutations that caused functional defects. When grouped according to the nature of the defects, these residues sub-divided into three striking sub-clusters: a trigger region, where mutations mostly affect ligand binding, a coupling region near the cytoplasmic interface to the G protein, where mutations affect G protein activation, and a linking core in between where mutations cause constitutive activity and other defects. Differential ET analysis of the opsin family revealed an additional set of opsin-specific residues, several of which form part of the retinal binding pocket, and are known to cause functional defects upon mutation. To test the predictive power of ET, we introduced novel mutations in bovine rhodopsin at a globally important position, Leu-79, and at an opsin-specific position, Trp-175. Both were functionally critical, causing constitutive G protein activation of the mutants and rapid loss of regeneration after photobleaching. These results define in GPCRs a canonical signal transduction mechanism where ligand binding induces conformational changes propagated through adjacent trigger, linking core, and coupling regions.     

Evaluation of site-directed spin labeling for characterizing protein-ligand complexes using simulated restraints

Simulation studies have been performed to evaluate the utility of site-directed spin labeling for determining the structures of protein-ligand complexes, given a known protein structure. Two protein-ligand complexes were used as model systems for these studies: a 1.9-A-resolution x-ray structure of a dihydrofolate reductase mutant complexed with methotrexate, and a 1.5-A-resolution x-ray structure of the V-Src tyrosine kinase SH2 domain complexed with a five-residue phosphopeptide. Nitroxide spin labels were modeled at five dihydrofolate reductase residue positions and at four SH2 domain residue positions. For both systems, after energy minimization, conformational ensembles of the spin-labeled residues were generated by simulated annealing while holding the remainder of the protein-ligand complex fixed. Effective distances, simulating those that could be obtained from (1)H-NMR relaxation measurements, were calculated between ligand protons and the spin labels. These were converted to restraints with several different levels of precision. Restrained simulated annealing calculations were then performed with the aim of reproducing target ligand-binding modes. The effects of incorporating a few supplementary short-range (< or =5.0 A) distance restraints were also examined. For the dihydrofolate reductase-methotrexate complex, the ligand-binding mode was reproduced reasonably well using relatively tight spin-label restraints, but methotrexate was poorly localized using loose spin-label restraints. Short-range and spin-label restraints proved to be complementary. For the SH2 domain-phosphopeptide complex without the short-range restraints, the peptide did not localize to the correct depth in the binding groove; nevertheless, the orientation and internal conformation of the peptide was reproduced moderately well. Use of the spin-label restraints in conjunction with the short-range restraints resulted in relatively well defined structural ensembles. These results indicate that restraints derived from site-directed spin labeling can contribute significantly to defining the orientations and conformations of bound ligands. Accurate ligand localization appears to require either a few supplementary short-range distance restraints, or relatively tight spin-label restraints, with at least one spin label positioned so that some of the restraints draw the ligand into the binding pocket in the latter case.     

Pre-existing soft modes of motion uniquely defined by native contact topology facilitate ligand binding to proteins

Modeling protein flexibility constitutes a major challenge in accurate prediction of protein-ligand and protein-protein interactions in docking simulations. The lack of a reliable method for predicting the conformational changes relevant to substrate binding prevents the productive application of computational docking to proteins that undergo large structural rearrangements. Here, we examine how coarse-grained normal mode analysis has been advantageously applied to modeling protein flexibility associated with ligand binding. First, we highlight recent studies that have shown that there is a close agreement between the large-scale collective motions of proteins predicted by elastic network models and the structural changes experimentally observed upon ligand binding. Then, we discuss studies that have exploited the predicted soft modes in docking simulations. Two general strategies are noted: pregeneration of conformational ensembles that are then utilized as input for standard fixed-backbone docking and protein structure deformation along normal modes concurrent to docking. These studies show that the structural changes apparently "induced" upon ligand binding occur selectively along the soft modes accessible to the protein prior to ligand binding. They further suggest that proteins offer suitable means of accommodating/facilitating the recognition and binding of their ligand, presumably acquired by evolutionary selection of the suitable three-dimensional structure.     

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.     

Solution structure and backbone dynamics of human Raf-1 kinase inhibitor protein

Human Raf-1 kinase inhibitor protein (hRKIP) is a small multi-functional protein of 187 residues. It contains a conserved pocket, which binds a wide range of ligands from various small molecules to distinct proteins. To provide a structural basis for the ligand diversity of RKIP, we herein determined the solution structure of hRKIP, and analyzed its structural dynamics. In solution, hRKIP mainly comprises two antiparallel β sheets, two α helices and two 3₁₀ helices. NMR dynamic analysis reveals that the overall structure of hRKIP is rigid, but its C-terminal helix which is close to the ligand-binding site is mobile. In addition, residues around the ligand-binding pocket exhibit significant conformational exchange on the μs–ms timescale. Conformational flexibility may allow the ligand-binding pocket and the C-terminal helix to adopt various conformations to interact with different substrates. This work may shed light on the underlying molecular mechanisms of how hRKIP recognizes and binds diverse substrate ligands.     

Ligand-Induced Protein Responses and Mechanical Signal Propagation Described by Linear Response Theories

In this study, a general linear response theory (LRT) is formulated to describe time-dependent and -independent protein conformational changes upon CO binding with myoglobin. Using the theory, we are able to monitor protein relaxation in two stages. The slower relaxation is found to occur from 4.4 to 81.2 picoseconds and the time constants characterized for a couple of aromatic residues agree with those observed by UV Resonance Raman (UVRR) spectrometry and time resolved x-ray crystallography. The faster “early responses”, triggered as early as 400 femtoseconds, can be best described by the theory when impulse forces are used. The newly formulated theory describes the mechanical propagation following ligand-binding as a function of time, space and types of the perturbation forces. The “disseminators”, defined as the residues that propagate signals throughout the molecule the fastest among all the residues in protein when perturbed, are found evolutionarily conserved and the mutations of which have been shown to largely change the CO rebinding kinetics in myoglobin.     

Activation of the luteinizing hormone receptor following substitution of Ser-277 with selective hydrophobic residues in the ectodomain hinge region

Glycoprotein hormone receptors are G protein-coupled receptors with ligand-binding ectodomains consisting of leucine-rich repeats. The ectodomain is connected by a conserved cysteine-rich hinge region to the seven transmembrane (TM) region. Gain-of-function mutants of luteinizing hormone (LH) and thyroid-stimulating hormone receptors found in patients allowed identification of residues important for receptor activation. Based on constitutively active mutations at Ser-281 in the hinge region of the thyroid-stimulating hormone receptor, we mutated the conserved serine in the LH (S277I) and follicle-stimulating hormone receptors (S273I) and observed increased basal cAMP production and ligand affinity by mutant receptors. For the LH receptor, conversion of Ser-277 to all natural amino acids led to varying degrees of receptor activation. Hydropathy index analysis indicated that substitution of neutral serine with selective nonpolar hydrophobic residues (Leu>Val>Met>Ile) confers constitutive receptor activation whereas serine deletion or substitution with charged Arg, Lys, or Asp led to defective receptor expression. Furthermore, mutation of the angular proline near Ser-273 to flexible Gly also led to receptor activation. The findings suggest the ectodomain of glycoprotein hormone receptors constrain the TM region. Point mutations in the hinge region of these proteins, or ligand binding to these receptors, could cause conformational changes in the TM region that result in G(s) activation.     

Mutational analysis of the human KDEL receptor: distinct structural requirements for Golgi retention, ligand binding and retrograde transport.

The KDEL receptor is a seven-transmembrane-domain protein that is responsible for the retrieval of endoplasmic reticulum (ER) proteins from the Golgi complex. It is a temporary resident of the Golgi apparatus: upon binding a KDEL-containing ligand, it moves to the ER, where the ligand is released. We have expressed mutant forms of the human receptor in COS cells and examined their intracellular locations and ligand-binding capacities. We show that ligand binding is dependent on charged residues within the transmembrane domains. Surprisingly, retrograde transport of occupied receptor is unaffected by most mutations in the cytoplasmic loops, but is critically dependent upon an aspartic acid residue in the seventh transmembrane domain. Retention in the Golgi apparatus requires neither ligand binding nor this aspartate residue, and thus is independent of receptor recycling. We suggest that movement of the receptor is controlled by conformational changes and intermolecular interactions within the membrane bilayer.     

Binding modes of Bcl-2 homology 3 (BH3) peptides with anti-apoptotic protein A1 and redesign of peptide inhibitors: a computational study

The interaction between protein and peptide ligand is a challenging problem in molecular biology and drug design. The binding of the Bcl-2 homology 3 (BH3) peptide to the anti-apoptotic protein A1 was revealed as a critical step in the regulation of apoptosis. These BH3 peptides hold high structural similarity, but are diverse in their regulation abilities. Based on molecular simulations and MM-P(G)BSA methods, this work presented a detailed analysis on binding mechanism of the BH3 peptides derived from PUMA and BMF. Residue-level energy decomposition showed that the core regions of BH3 peptides maintain in stable helical conformations and the four conserved hydrophobic residues together with an invariant aspartic acid contribute the major driving force for binding, whereas their two terminal segments exhibit obvious flexibility and various binding modes. Such kind of behavior was suggested as the reason for binding diversity and selectivity of BH3 peptides. As a further step, several BH3-mimetic peptides have been redesigned by computational mutation. Those new peptides showed not only stronger affinities when binding to protein A1, as well transferable binding patterns at some specific positions. A long-range coupling effect was disclosed for BH3 peptides, side-chain orientation and binding contribution of terminal residues were even affected by mutations at large sequence interval. Overall, this work reports that the binding modes of BH3 peptides are primarily dependent on its two terminal segments. The computational methods applied herein are also demonstrated to be of great assistance in the rational design of peptide inhibitors.     

Structural insights into the cellular retinaldehyde-binding protein (CRALBP)

Cellular retinaldehyde-binding protein (CRALBP) is an essential protein in the human visual cycle without a known three-dimensional structure. Previous studies associate retinal pathologies to specific mutations in the CRALBP protein. Here we use homology modeling and molecular dynamics methods to investigate the structural mechanisms by which CRALBP functions in the visual cycle. We have constructed two conformations of CRALBP representing two states in the process of ligand association and dissociation. Notably, our homology models map the pathology-associated mutations either directly in or adjacent to the putative ligand-binding cavity. Furthermore, six novel residues have been identified to be crucial for the hinge movement of the lipid-exchange loop in CRALBP. We conclude that the binding and release of retinoid involve large conformational changes in the lipid-exchange loop at the entrance of the ligand-binding cavity.     

Correlating efficacy and desensitization with GluK2 ligand-binding domain movements

Gating of AMPA- and kainate-selective ionotropic glutamate receptors can be defined in terms of ligand affinity, efficacy and the rate and extent of desensitization. Crucial insights into all three elements have come from structural studies of the ligand-binding domain (LBD). In particular, binding-cleft closure is associated with efficacy, whereas dissociation of the dimer formed by neighbouring LBDs is linked with desensitization. We have explored these relationships in the kainate-selective subunit GluK2 by studying the effects of mutating two residues (K531 and R775) that form key contacts within the LBD dimer interface, but whose truncation unexpectedly attenuates desensitization. One mutation (K531A) also switches the relative efficacies of glutamate and kainate. LBD crystal structures incorporating these mutations revealed several conformational changes that together explain their phenotypes. K531 truncation results in new dimer contacts, consistent with slower desensitization and sideways movement in the ligand-binding cleft correlating with efficacy. The tested mutants also disrupted anion binding; no chloride was detected in the dimer-interface site, including in R775A where absence of chloride was the only structural change evident. From this, we propose that the charge balance in the GluK2 LBD dimer interface maintains a degree of instability, necessary for rapid and complete desensitization.     

The Role of Protein-Ligand Contacts in Allosteric Regulation of the Escherichia coli Catabolite Activator Protein

Allostery is a fundamental process by which ligand binding to a protein alters its activity at a distant site. Both experimental and theoretical evidence demonstrate that allostery can be communicated through altered slow relaxation protein dynamics without conformational change. The catabolite activator protein (CAP) of Escherichia coli is an exemplar for the analysis of such entropically driven allostery. Negative allostery in CAP occurs between identical cAMP binding sites. Changes to the cAMP-binding pocket can therefore impact the allosteric properties of CAP. Here we demonstrate, through a combination of coarse-grained modeling, isothermal calorimetry, and structural analysis, that decreasing the affinity of CAP for cAMP enhances negative cooperativity through an entropic penalty for ligand binding. The use of variant cAMP ligands indicates the data are not explained by structural heterogeneity between protein mutants. We observe computationally that altered interaction strength between CAP and cAMP variously modifies the change in allosteric cooperativity due to second site CAP mutations. As the degree of correlated motion between the cAMP-contacting site and a second site on CAP increases, there is a tendency for computed double mutations at these sites to drive CAP toward noncooperativity. Naturally occurring pairs of covarying residues in CAP do not display this tendency, suggesting a selection pressure to fine tune allostery on changes to the CAP ligand-binding pocket without a drive to a noncooperative state. In general, we hypothesize an evolutionary selection pressure to retain slow relaxation dynamics-induced allostery in proteins in which evolution of the ligand-binding site is occurring.     

Extreme sequence divergence but conserved ligand-binding specificity in Streptococcus pyogenes M protein

Many pathogenic microorganisms evade host immunity through extensive sequence variability in a protein region targeted by protective antibodies. In spite of the sequence variability, a variable region commonly retains an important ligand-binding function, reflected in the presence of a highly conserved sequence motif. Here, we analyze the limits of sequence divergence in a ligand-binding region by characterizing the hypervariable region (HVR) of Streptococcus pyogenes M protein. Our studies were focused on HVRs that bind the human complement regulator C4b-binding protein (C4BP), a ligand that confers phagocytosis resistance. A previous comparison of C4BP-binding HVRs identified residue identities that could be part of a binding motif, but the extended analysis reported here shows that no residue identities remain when additional C4BP-binding HVRs are included. Characterization of the HVR in the M22 protein indicated that two relatively conserved Leu residues are essential for C4BP binding, but these residues are probably core residues in a coiled-coil, implying that they do not directly contribute to binding. In contrast, substitution of either of two relatively conserved Glu residues, predicted to be solvent-exposed, had no effect on C4BP binding, although each of these changes had a major effect on the antigenic properties of the HVR. Together, these findings show that HVRs of M proteins have an extraordinary capacity for sequence divergence and antigenic variability while retaining a specific ligand-binding function.