Interplay of buried histidine protonation and protein stability in prion misfolding

A role for pH dynamics regulating transcription factor DNA binding selectivity

Intracellular pH (pHi) dynamics regulates diverse cell processes such as proliferation, dysplasia, and differentiation, often mediated by the protonation state of a functionally critical histidine residue in endogenous pH sensing proteins. How pHi dynamics can directly regulate gene expression and whether transcription factors can function as pH sensors has received limited attention. We tested the prediction that transcription factors with a histidine in their DNA binding domain (DBD) that forms hydrogen bonds with nucleotides can have pH-regulated activity, which is relevant to more than 85 transcription factors in distinct families, including FOX, KLF, SOX and MITF/Myc. Focusing on FOX family transcription factors, we used unbiased SELEX-seq to identify pH-dependent DNA binding motif preferences, then confirm pH-regulated binding affinities for FOXC2, FOXM1, and FOXN1 to a canonical FkhP DNA motif that are 2.5 to 7.5 greater at pH 7.0 compared with pH 7.5. For FOXC2, we also find greater activity for an FkhP motif at lower pHi in cells and that pH-regulated binding and activity are dependent on a conserved histidine (His122) in the DBD. RNA-seq with FOXC2 also reveals pH-dependent differences in enriched promoter motifs. Our findings identify pH-regulated transcription factor-DNA binding selectivity with relevance to how pHi dynamics can regulate gene expression for myriad cell behaviours.

E-Protein Protonation Titration-Induced Single-Particle Chemical Force Spectroscopy for Microscopic Understanding and pI Estimation of Infectious DENV

The ionization state of amino acids on the outer surface of a virus regulates its physicochemical properties toward the sorbent surface. Serologically different strains of the dengue virus (DENV) show different extents of infectivity depending upon their interactions with a receptor on the host cell. To understand the structural dependence of E-protein protonation over its sequence dependence, we have followed E-protein titration kinetics both experimentally and theoretically for two differentially infected dengue serotypes, namely, DENV-2 and DENV-4. We have performed E-protein protonation titration-induced single-particle chemical force spectroscopy using an atomic force microscope (AFM) to measure the surface chemistry of DENV in physiological aqueous solutions not only to understand the charge distribution dynamics on the virus surface but also to estimate the isoelectric point (pI) accurately for infectious dengue viruses. Cryo-EM structure-based theoretical pI calculations of the DENV-2 surface protein were shown to be consistent with the evaluated pI value from force spectroscopy measurements. We also highlighted here the role of the microenvironment around the titrable residues (in the 3D-folded structure of the protein) in altering the pKa. This is a comprehensive study to understand how the cumulative charge distribution on the outer surface of a specific serotype of DENV regulates a prominent role of infectivity over minute changes at the genetic level.

Met/Val129 polymorphism of the full-length human prion protein dictates distinct pathways of amyloid formation

Methionine/valine polymorphism at position 129 of the human prion protein, huPrP, is tightly associated with the pathogenic phenotype, disease progress, and age of onset of neurodegenerative diseases such as Creutzfeldt–Jakob disease or Fatal Familial Insomnia. This raises the question of whether and how the amino acid type at position 129 influences the structural properties of huPrP, affecting its folding, stability, and amyloid formation behavior. Here, our detailed biophysical characterization of the 129M and 129V variants of recombinant full-length huPrP(23–230) by amyloid formation kinetics, CD spectroscopy, molecular dynamics simulations, and sedimentation velocity analysis reveals differences in their aggregation propensity and oligomer content, leading to deviating pathways for the conversion into amyloid at acidic pH. We determined that the 129M variant exhibits less secondary structure content before amyloid formation and higher resistance to thermal denaturation compared to the 129V variant, whereas the amyloid conformation of both variants shows similar thermal stability. Additionally, our molecular dynamics simulations and rigidity analyses at the atomistic level identify intramolecular interactions responsible for the enhanced monomer stability of the 129M variant, involving more frequent minimum distances between E196 and R156, forming a salt bridge. Removal of the N-terminal half of the 129M full-length variant diminishes its differences compared to the 129V full-length variant and highlights the relevance of the flexible N terminus in huPrP. Taken together, our findings provide insight into structural properties of huPrP and the effects of the amino acid identity at position 129 on amyloid formation behavior.

The Physical Basis for pH Sensitivity in Biomolecular Structure and Function, With Application to the Spike Protein of SARS-CoV-2

Since pH sensitivity has a fundamental role in biology, much effort has been committed to establishing physical models to rationalize and predict pH dependence from molecular structures. Two of the key challenges are to accurately calculate ionizable group solvation and hydration and then to apply this modeling to all conformations relevant to the process in question. Explicit solvent methods coupled to molecular dynamics simulation are increasingly complementing lower resolution implicit solvent techniques, but equally, the scale of biological data acquisition leaves a role for high-throughput modeling. Additionally, determination of ranges of structures for a system allows sampling of key stages in solvation. In a review of the area, it is emphasized that pH sensors in biology beyond the most obvious candidate (histidine side chain, with an unshifted pK a near neutral pH) should be considered; that modeling can benefit from other concepts in bioinformatics, in particular modulation of interactions and function in families of homologs; and that it can also be beneficial to incorporate as many experimental structures as possible, to mitigate against small variations in conformation and to analyze larger, functional, conformational changes. These aspects are then demonstrated with new work on the spike protein of SARS-CoV-2, looking at the pH dependence of variants, including prediction of a change in the balance of locked, closed, and open forms at neutral pH for the Omicron variant spike protein.

Deuterium-Enhanced Raman Spectroscopy for Histidine pKa Determination in a pH-Responsive Hydrogel

We report here a method for the determination of the pK_a of histidine in complex or heterogeneous systems amenable to neither solid-state nor solution NMR spectroscopy. Careful synthesis of a fluorenylmethyloxycarbonyl- and trityl-protected, C2-deuterated histidine produces a vibrational-probe-equipped amino acid that can readily be incorporated into any peptide accessible by standard solid-phase methods. The frequency of the unique, Raman-active stretching vibration of this C2-D probe is a clear reporter of the protonation state of histidine. We investigate here a pH-sensitive peptide that self-assembles to form a hydrogel at neutral pH. The pK_a of the lone histidine residue in the peptide, which is likely responsible for this pH-dependent behavior, cannot be investigated by NMR spectroscopy because of the supramolecular, soft nature of the gel. However, after synthesizing a C2-deuterated-histidine-containing peptide, we were able to follow the protonation state of histidine throughout a pH titration using Raman difference spectroscopy, thereby precisely determining the pK_a of interest.

Examination of pH dependency and orientation differences of membrane spanning alpha helices carrying a single or pair of buried histidine residues

We have employed the peptide framework of GWALP23 (acetyl-GGALWLALALALALALALWLAGA-amide) to examine the orientation, dynamics and pH dependence of peptides having buried single or pairs of histidine residues. When residue L8 is substituted to yield GWALP23-H⁸, acetyl-GGALWLAH⁸ALALALALALWLAGA-amide, the deuterium NMR spectra of ²H-labeled core alanine residues reveal a helix that occupies a single transmembrane orientation in DLPC, or in DMPC at low pH, yet shows multiple states at higher pH or in bilayers of DOPC. Moreover, a single histidine at position 8 or 16 in the GWALP23 framework is sensitive to pH. Titration points are observed near pH 3.5 for the deprotonation of H8 in lipid bilayers of DLPC or DMPC, and for H16 in DOPC. When residues L8 and L16 both are substituted to yield GWALP23-H^8,16, the ²H NMR spectra show, interestingly, no titration dependence from pH 2-8, yet bilayer thickness-dependent orientation differences. The helix with H8 and H16 is found to adopt a transmembrane orientation in thin bilayers of DLPC, a combination of transmembrane and surface orientations in DMPC, and then a complete transition to a surface bound orientation in the thicker DPoPC and DOPC lipid bilayers. In the surface orientations, alanine A7 no longer fits within the core helix. These results along with previous studies with different locations of histidine residues suggest that lipid hydrophobic thickness is a first determinant and pH a second determinant for the helical orientation, along with possible side-chain snorkeling, when the His residues are incorporated into the hydrophobic region of a lipid membrane-associated helix.

Downhill (Un)Folding Coupled to Binding as a Mechanism for Engineering Broadband Protein Conformational Transducers

Canonical proteins fold and function as conformational switches that toggle between their folded (on) and unfolded (off) states, a mechanism that also provides the basis for engineering transducers for biosensor applications. One of the limitations of such transducers, however, is their relatively narrow operational range, limited to ligand concentrations 20-fold below or above their C₅₀. Previously, we discovered that certain fast-folding proteins lose/gain structure gradually (downhill folding), which led us to postulate their operation as conformational rheostats capable of processing inputs/outputs in analog fashion. Conformational rheostats could make transducers with extended sensitivity. Here we investigate this hypothesis by engineering pH transducing into the naturally pH insensitive, downhill folding protein gpW. Particularly, we engineered histidine grafts into its hydrophobic core to induce unfolding via histidine ionization. We designed and tested the effects of ionization via computational modeling and studied experimentally the four most promising single grafts and two double grafts. All tested mutants become reversible pH transducers in the 4-9 range, and their response increases proportionally to how buried the histidine graft is. Importantly, the pH-dependent reversible (un)folding occurs in rheostatic fashion, so the engineered transducers can detect up to 6 orders of magnitude in [H⁺] for single grafts, and even more for double grafts. Our results demonstrate that downhill (un)folding coupled to binding produces the gradual, analog responses to the ligand (here H⁺) that are expected of conformational rheostats, and which make them a powerful mechanism for engineering transducers with sensitivity over many orders of magnitude in ligand concentration (broadband).

A seven-residue deletion in PrP leads to generation of a spontaneous prion formed from C-terminal C1 fragment of PrP

Prions result from a drastic conformational change of the host-encoded cellular prion protein (PrP), leading to the formation of β-sheet-rich, insoluble, and protease-resistant self-replicating assemblies (PrP^Sc). The cellular and molecular mechanisms involved in spontaneous prion formation in sporadic and inherited human prion diseases or equivalent animal diseases are poorly understood, in part because cell models of spontaneously forming prions are currently lacking. Here, extending studies on the role of the H2 α-helix C terminus of PrP, we found that deletion of the highly conserved ¹⁹⁰HTVTTTT¹⁹⁶ segment of ovine PrP led to spontaneous prion formation in the RK13 rabbit kidney cell model. On long-term passage, the mutant cells stably produced proteinase K (PK)-resistant, insoluble, and aggregated assemblies that were infectious for naïve cells expressing either the mutant protein or other PrPs with slightly different deletions in the same area. The electrophoretic pattern of the PK-resistant core of the spontaneous prion (Δ^Spont) contained mainly C-terminal polypeptides akin to C1, the cell-surface anchored C-terminal moiety of PrP generated by natural cellular processing. RK13 cells expressing solely the Δ190-196 C1 PrP construct, in the absence of the full-length protein, were susceptible to Δ^Spont prions. Δ^Spont infection induced the conversion of the mutated C1 into a PK-resistant and infectious form perpetuating the biochemical characteristics of Δ^Spont prion. In conclusion, this work provides a unique cell-derived system generating spontaneous prions and provides evidence that the 113 C-terminal residues of PrP are sufficient for a self-propagating prion entity.

In Situ pH Measurements of Individual Levitated Microdroplets Using Aerosol Optical Tweezers

The pH of microscale reaction environments controls numerous physicochemical processes, requiring a real-time pH microprobe. We present highly accurate real-time pH measurements of microdroplets using aerosol optical tweezers (AOT) and analysis of the whispering gallery modes (WGMs) contained in the cavity-enhanced Raman spectra. Uncertainties ranging from ±0.03 to 0.06 in pH for picoliter droplets are obtained through averaging Raman frames acquired at 0.5 Hz over 3.3 min. The high accuracy in pH determination is achieved by combining two independent measurements uniquely provided by the AOT approach: the anion concentration ratio from the spontaneous Raman spectra, and the total solute concentration from the refractive index retrieved from WGM analysis of the stimulated cavity-enhanced Raman spectra. pH can be determined over a range of -0.36 to 0.76 using the aqueous sodium bisulfate system. This technique enables direct measurements of pH-dependent chemical and physical changes experienced by individual microparticles and exploration of the role of pH in the chemical behavior of confined microenvironments.

Structural insight into conformational change in prion protein by breakage of electrostatic network around H187 due to its protonation

A conformational change from normal prion protein(PrP^C) to abnormal prion protein(PrP^SC) induces fatal neurodegenerative diseases. Acidic pH is well-known factors involved in the conformational change. Because the protonation of H187 is strongly linked to the change in PrP stability, we examined the charged residues R156, E196, and D202 around H187. Interestingly, there have been reports on pathological mutants, such as H187R, E196A, and D202N. In this study, we focused on how an acidic pH and pathological mutants disrupt this electrostatic network and how this broken network destabilizes PrP structure. To do so, we performed a temperature-based replica-exchange molecular dynamics (T-REMD) simulation using a cumulative 252 μs simulation time. We measured the distance between amino acids comprising four salt bridges (R156–E196/D202 and H187–E196/D202). Our results showed that the spatial configuration of the electrostatic network was significantly altered by an acidic pH and mutations. The structural alteration in the electrostatic network increased the RMSF value around the first helix (H1). Thus, the structural stability of H1, which is anchored to the H2–H3 bundle, was decreased. It induces separation of R156 from the electrostatic network. Analysis of the anchoring energy also shows that two salt-bridges (R156-E196/D202) are critical for PrP stability.

A designer molecular chaperone against transmissible spongiform encephalopathy slows disease progression in mice and macaques

Transmissible spongiform encephalopathies (TSEs) are fatal neurodegenerative diseases that lack therapeutic solutions. Here, we show that the molecular chaperone (N,N'-([cyclohexylmethylene]di-4,1-phenylene)bis(2-[1-pyrrolidinyl]acetamide)), designed via docking simulations, molecular dynamics simulations and quantum chemical calculations, slows down the progress of TSEs. In vitro, the designer molecular chaperone stabilizes the normal cellular prion protein, eradicates prions in infected cells, prevents the formation of drug-resistant strains and directly inhibits the interaction between prions and abnormal aggregates, as shown via real-time quaking-induced conversion and in vitro conversion NMR. Weekly intraperitoneal injection of the chaperone in prion-infected mice prolonged their survival, and weekly intravenous administration of the compound in macaques infected with bovine TSE slowed down the development of neurological and psychological symptoms and reduced the concentration of disease-associated biomarkers in the animals' cerebrospinal fluid. The de novo rational design of chaperone compounds could lead to therapeutics that can bind to different prion protein strains to ameliorate the pathology of TSEs.

G-score: A function to solve the puzzle of modeling the protonation states of β-secretase binding pocket

The population density concept has emerged as a proposal for the analysis of molecular dynamics results, the key characteristic of population density is the evaluation of the simultaneous occurrence of a set of relevant parameters for a system. However, despite its statistical strength, selection of the tolerance level for the comparison of different models may appear as arbitrary. This work introduces the G-score, a function which summarizes and categorizes the results of population density analysis. Additionally, it incorporates parameters based on rmsd and dihedral angles, besides the protein-protein and protein-ligand interatomic distances conventionally used, which complement each other to provide a better description of the behavior of the system. These newly-proposed tools were applied to determine the most probable protonation state of the aspartic dyad of BACE1, Asp93 and Asp289, in the presence of three types of transition state inhibitors namely: reduced amides, tertiary carbinamines and hydroxyethylamines. The results show a full agreement between G-score values and population density charts, with the advantage of allowing a quick and direct comparison among all the considered models. We anticipate that the simplicity of calculating the parameters employed in this study will permit the extensive use of population density and the G-score for other molecular systems.

Anomalous Properties of Lys Residues Buried in the Hydrophobic Interior of a Protein Revealed with 15N-Detect NMR Spectroscopy

Ionizable residues buried in hydrophobic environments in proteins are essential for many fundamental biochemical processes. These residues titrate with anomalous pK_a values that are challenging to reproduce with structure-based calculations owing to the conformational reorganization coupled to their ionization. Detailed characterization of this conformational reorganization is of interest; unfortunately, the properties of buried Lys residues are difficult to study experimentally. Here we demonstrate the utility of ¹⁵N NMR spectroscopy to gain insight into the protonation state, state of hydration and conformational dynamics of the Nζ amino group of buried Lys residues. The experiments were applied to five variants of staphylococcal nuclease, with internal Lys residues that titrate with pK_a values ranging from 6.2 to 8.1. Direct detection of buried Lys residues with these NMR spectroscopy methods will enable correlation between thermodynamic and structural data as well as unprecedented examination of how conformational transitions coupled to their ionization affect their pK_a values.

In Vitro Approach To Identify Key Amino Acids in Low Susceptibility of Rabbit Prion Protein to Misfolding

Prion diseases, or transmissible spongiform encephalopathies (TSEs), are a group of rare progressive neurodegenerative disorders caused by an abnormally folded prion protein (PrP^Sc). This is capable of transforming the normal cellular prion protein (PrP^C) into new infectious PrP^Sc Interspecies prion transmissibility studies performed by experimental challenge and the outbreak of bovine spongiform encephalopathy that occurred in the late 1980s and 1990s showed that while some species (sheep, mice, and cats) are readily susceptible to TSEs, others are apparently resistant (rabbits, dogs, and horses) to the same agent. To study the mechanisms of low susceptibility to TSEs of certain species, the mouse-rabbit transmission barrier was used as a model. To identify which specific amino acid residues determine high or low susceptibility to PrP^Sc propagation, protein misfolding cyclic amplification (PMCA), which mimics PrP^C-to-PrP^Sc conversion with accelerated kinetics, was used. This allowed amino acid substitutions in rabbit PrP and accurate analysis of misfolding propensities. Wild-type rabbit recombinant PrP could not be misfolded into a protease-resistant self-propagating isoform in vitro despite seeding with at least 12 different infectious prions from diverse origins. Therefore, rabbit recombinant PrP mutants were designed to contain every single amino acid substitution that distinguishes rabbit recombinant PrP from mouse recombinant PrP. Key amino acid residue substitutions were identified that make rabbit recombinant PrP susceptible to misfolding, and using these, protease-resistant misfolded recombinant rabbit PrP was generated. Additional studies characterized the mechanisms by which these critical amino acid residue substitutions increased the misfolding susceptibility of rabbit PrP.IMPORTANCE Prion disorders are invariably fatal, untreatable diseases typically associated with long incubation periods and characteristic spongiform changes associated with neuronal loss in the brain. Development of any treatment or preventative measure is dependent upon a detailed understanding of the pathogenesis of these diseases, and understanding the mechanism by which certain species appear to be resistant to TSEs is critical. Rabbits are highly resistant to naturally acquired TSEs, and even under experimental conditions, induction of clinical disease is not easy. Using recombinant rabbit PrP as a model, this study describes critical molecular determinants that confer this high resistance to transmissible spongiform encephalopathies.