Protein Allostery at Atomic Resolution.

Protein allostery is a phenomenon involving the long range coupling between two distal sites in a protein. In order to elucidate allostery at atomic resoluion on the ligand-binding WW domain of the enzyme Pin1, multistate structures were calculated from exact nuclear Overhauser effect (eNOE). In its free form, the protein undergoes a microsecond exchange between two states, one of which is predisposed to interact with its parent catalytic domain. In presence of the positive allosteric ligand, the equilibrium between the two states is shifted towards domain-domain interaction, suggesting a population shift model. In contrast, the allostery-suppressing ligand decouples the side-chain arrangement at the inter-domain interface thereby reducing the inter-domain interaction. As such, this mechanism is an example of dynamic allostery. The presented distinct modes of action highlight the power of the interplay between dynamics and function in the biological activity of proteins.

Proteome-wide analysis of phospho-regulated PDZ domain interactions.

A key function of reversible protein phosphorylation is to regulate protein-protein interactions, many of which involve short linear motifs (3-12 amino acids). Motif-based interactions are difficult to capture because of their often low-to-moderate affinities. Here, we describe phosphomimetic proteomic peptide-phage display, a powerful method for simultaneously finding motif-based interaction and pinpointing phosphorylation switches. We computationally designed an oligonucleotide library encoding human C-terminal peptides containing known or predicted Ser/Thr phosphosites and phosphomimetic variants thereof. We incorporated these oligonucleotides into a phage library and screened the PDZ (PSD-95/Dlg/ZO-1) domains of Scribble and DLG1 for interactions potentially enabled or disabled by ligand phosphorylation. We identified known and novel binders and characterized selected interactions through microscale thermophoresis, isothermal titration calorimetry, and NMR We uncover site-specific phospho-regulation of PDZ domain interactions, provide a structural framework for how PDZ domains accomplish phosphopeptide binding, and discuss ligand phosphorylation as a switching mechanism of PDZ domain interactions. The approach is readily scalable and can be used to explore the potential phospho-regulation of motif-based interactions on a large scale.

The Sign of Nuclear Magnetic Resonance Chemical Shift Difference as a Determinant of the Origin of Binding Selectivity: Elucidation of the Position Dependence of Phosphorylation in Ligands Binding to Scribble PDZ1.

The use of nuclear magnetic resonance chemical shift perturbation to monitor changes taking place around the binding site of a ligand-protein interaction is a routine and widely applied methodology in the field of protein biochemistry. Shifts are often acquired by titrating various concentrations of ligand to a fixed concentration of the receptor and may serve the purpose, among others, of determining affinity constants, locating binding surfaces, or differentiating between binding mechanisms. Shifts are quantified by the so-called combined chemical shift difference. Although the directionality of shift changes is often used for detailed analysis of specific cases, the approach has not been adapted in standard chemical shift monitoring. This is surprising as it would not require additional effort. Here, we demonstrate the importance of the sign of the chemical shift difference induced by ligand-protein interaction. We analyze the sign of the 15N/1H shift changes of the PDZ1 domain of Scribble upon interaction with two pairs of phosphorylated and unphosphorylated peptides. We find that detailed differences in the molecular basis of this PDZ-ligand interaction can be obtained from our analysis to which the classical method of combined chemical shift perturbation analysis is insensitive. In addition, we find a correlation between affinity and millisecond motions. Application of the methodology to Cyclophilin a, a cis-trans isomerase, reveals molecular details of peptide recognition. We consider our directionality vector chemical shift analysis as a method of choice when distinguishing the molecular origin of binding specificities of a class of similar ligands, which is often done in drug discovery.

NOE-Derived Methyl Distances from a 360†kDa Proteasome Complex.

Nuclear magnetic resonance spectroscopy is the prime tool to probe structure and dynamics of biomolecules at atomic resolution. A serious challenge for that method is the size limit imposed on molecules to be studied. Standard studies are typically restricted to ca. 30-40 kDa. More recent developments lead to spin relaxation measurements in methyl groups in single proteins or protein complexes as large as a mega-Dalton, which directly allow the extraction of angular information or experiments with paramagnetic samples. However, these probes are all of indirect nature in contrast to the most intuitive and easy-to-interpret structural/dynamics restraint, the internuclear distance, which can be measured by nuclear Overhauser enhancement (NOE). Herein, we demonstrate time-averaged distance measurements on the 360 kDa half proteasome from Thermoplasma acidophilium. The approach is based on exact quantification of the NOE (eNOE). Our findings open up an avenue for such measurements on very large molecules. These restraints will help in a detailed determination of conformational changes upon perturbation such as ligand binding.

The Exact Nuclear Overhauser Enhancement: Recent Advances.

Although often depicted as rigid structures, proteins are highly dynamic systems, whose motions are essential to their functions. Despite this, it is difficult to investigate protein dynamics due to the rapid timescale at which they sample their conformational space, leading most NMR-determined structures to represent only an averaged snapshot of the dynamic picture. While NMR relaxation measurements can help to determine local dynamics, it is difficult to detect translational or concerted motion, and only recently have significant advances been made to make it possible to acquire a more holistic representation of the dynamics and structural landscapes of proteins. Here, we briefly revisit our most recent progress in the theory and use of exact nuclear Overhauser enhancements (eNOEs) for the calculation of structural ensembles that describe their conformational space. New developments are primarily targeted at increasing the number and improving the quality of extracted eNOE distance restraints, such that the multi-state structure calculation can be applied to proteins of higher molecular weights. We then review the implications of the exact NOE to the protein dynamics and function of cyclophilin A and the WW domain of Pin1, and finally discuss our current research and future directions.

Contribution of specific residues of the ß-solenoid fold to HET-s prion function, amyloid structure and stability.

The [Het-s] prion of the fungus Podospora anserina represents a good model system for studying the structure-function relationship in amyloid proteins because a high resolution solid-state NMR structure of the amyloid prion form of the HET-s prion forming domain (PFD) is available. The HET-s PFD adopts a specific ß-solenoid fold with two rungs of ß-strands delimiting a triangular hydrophobic core. A C-terminal loop folds back onto the rigid core region and forms a more dynamic semi-hydrophobic pocket extending the hydrophobic core. Herein, an alanine scanning mutagenesis of the HET-s PFD was conducted. Different structural elements identified in the prion fold such as the triangular hydrophobic core, the salt bridges, the asparagines ladders and the C-terminal loop were altered and the effect of these mutations on prion function, fibril structure and stability was assayed. Prion activity and structure were found to be very robust; only a few key mutations were able to corrupt structure and function. While some mutations strongly destabilize the fold, many substitutions in fact increase stability of the fold. This increase in structural stability did not influence prion formation propensity in vivo. However, if an Ala replacement did alter the structure of the core or did influence the shape of the denaturation curve, the corresponding variant showed a decreased prion efficacy. It is also the finding that in addition to the structural elements of the rigid core region, the aromatic residues in the C-terminal semi-hydrophobic pocket are critical for prion propagation. Mutations in the latter region either positively or negatively affected prion formation. We thus identify a region that modulates prion formation although it is not part of the rigid cross-ß core, an observation that might be relevant to other amyloid models.

Enzyme Selectivity Fine-Tuned through Dynamic Control of a Loop.

Allostery has been revealed as an essential property of all proteins. For enzymes, shifting of the structural equilibrium distribution at one site can have substantial impacts on protein dynamics and selectivity. Promising sites of remotely shifting such a distribution by changing the dynamics would be at flexible loops because relatively large changes may be achieved with minimal modification of the protein. A ligand-selective change of binding affinity to the active site of cyclophilin is presented involving tuning of the dynamics of a highly flexible loop. Binding affinity is increased upon substitution of double Gly to Ala at the hinge regions of the loop. Quenching of the motional amplitudes of the loop slightly rearranges the active site. In particular, key residues for binding Phe60 and His126 adopt a more fixed orientation in the bound protein. Our system may serve as a model system for studying the effects of various time scales of loop motion on protein function tuned by mutations.

Extending the eNOE data set of large proteins by evaluation of NOEs with unresolved diagonals.

The representation of a protein's spatial sampling at atomic resolution is fundamental for understanding its function. NMR has been established as the best-suited technique toward this goal for small proteins. However, the accessible information content rapidly deteriorates with increasing protein size. We have recently demonstrated that for small proteins distance restraints with an accuracy smaller than 0.1 Å can be obtained by replacing traditional semi-quantitative Nuclear Overhauser Effects (NOEs) with exact NOEs (eNOE). The high quality of the data allowed us to calculate structural ensembles of the small model protein GB3 consisting of multiple rather than a single state. The analysis has been limited to small proteins because NOEs of spins with unresolved diagonal peaks cannot be used. Here we propose a simple approach to translate such NOEs into correct upper distance restraints, which opens access to larger biomolecules. We demonstrate that for 16 kDa cyclophilin A the collection of such restraints extends the original 1254 eNOEs to 3471.

A high-affinity, dimeric inhibitor of PSD-95 bivalently interacts with PDZ1-2 and protects against ischemic brain damage.

Inhibition of the ternary protein complex of the synaptic scaffolding protein postsynaptic density protein-95 (PSD-95), neuronal nitric oxide synthase (nNOS), and the N-methyl-D-aspartate (NMDA) receptor is a potential strategy for treating ischemic brain damage, but high-affinity inhibitors are lacking. Here we report the design and synthesis of a novel dimeric inhibitor, Tat-NPEG4(IETDV)(2) (Tat-N-dimer), which binds the tandem PDZ1-2 domain of PSD-95 with an unprecedented high affinity of 4.6 nM, and displays extensive protease-resistance as evaluated in vitro by stability-measurements in human blood plasma. X-ray crystallography, NMR, and small-angle X-ray scattering (SAXS) deduced a true bivalent interaction between dimeric inhibitor and PDZ1-2, and also provided a dynamic model of the conformational changes of PDZ1-2 induced by the dimeric inhibitor. A single intravenous injection of Tat-N-dimer (3 nmol/g) to mice subjected to focal cerebral ischemia reduces infarct volume with 40% and restores motor functions. Thus, Tat-N-dimer is a highly efficacious neuroprotective agent with therapeutic potential in stroke.

A Structural Ensemble for the Enzyme Cyclophilin Reveals an Orchestrated Mode of Action at Atomic Resolution.

For enzyme activity, an exact structural and motional orchestration of the active site and its surroundings is believed to be key. In order to reveal such possible phenomena at atomic resolution on the basis of experimental evidence, an experimental restraint driven two-state ensemble of the prototypical enzyme cyclophilin was determined by using a recently introduced exact NOE approach. The ensemble description reveals the presence of an open and a closed state of cyclophilin, which is indicative of large-scale correlated motion. In the open state, the catalytic site is preorganized for catalysis, thus suggesting the mechanism of action to be conformational sampling, while the ligand-binding loop appears to act through an induced fit mechanism. This finding is supported by affinity measurements of a cyclophilin designed to be more open. Overall, more than 60-70 % of the side-chain conformations of cyclophilin appear to be correlated.

Folded Alpha Helical Putative New Proteins from Apilactobacillus kunkeei.

The emergence of new proteins is a central question in biology. Most tertiary protein folds known to date appear to have an ancient origin, but it is clear from bioinformatic analyses that new proteins continuously emerge in all organismal groups. However, there is a paucity of experimental data on new proteins regarding their structure and biophysical properties. We performed a detailed phylogenetic analysis and identified 48 putative open reading frames in the honeybee-associated bacterium Apilactobacillus kunkeei for which no or few homologs could be identified in closely-related species, suggesting that they could be relatively new on an evolutionary time scale and represent recently evolved proteins. Using circular dichroism-, fluorescence- and nuclear magnetic resonance (NMR) spectroscopy we investigated six of these proteins and show that they are not intrinsically disordered, but populate alpha-helical dominated folded states with relatively low thermodynamic stability (0-3 kcal/mol). The NMR and biophysical data demonstrate that small new proteins readily adopt simple folded conformations suggesting that more complex tertiary structures can be continuously re-invented during evolution by fusion of such simple secondary structure elements. These findings have implications for the general view on protein evolution, where de novo emergence of folded proteins may be a common event.

Shotgun Metagenomic Analyses of Microbial Assemblages in the Aquatic Ecosystem of Winam Gulf of Lake Victoria, Kenya Reveals Multiclass Pollution.

Lake Victoria, the second-largest freshwater lake in the world, provides an important source of food and income, particularly fish for both domestic consumption and for export market. In recent years, Lake Victoria has suffered massive pollution from both industrial and wastewater discharge. Microplastic biomes, pharmaceutical residues, drugs of abuse, heavy metals, agrochemicals, and personal care products are ubiquitous in the aquatic ecosystem of Winam Gulf. These pollutants are known to alter microbial assemblages in aquatic ecosystems with far-reaching ramification including a calamitous consequence to human health. Indeed, some of these pollutants have been associated with human cancers and antimicrobial resistance. There is a paucity of data on the microbial profiles of this important but heavily polluted aquatic ecosystem. The current study sought to investigate the metagenomic profiles of microbial assemblages in the Winam Gulf ecosystem. Water and sediment samples were collected from several locations within the study sites. Total genomic DNA pooled from all sampling sites was extracted and analyzed by whole-genome shotgun sequencing. Analyses revealed three major kingdoms: bacteria, archaea and eukaryotes belonging to 3 phyla, 13 classes, 14 families, 9 orders, 14 genera, and 10 species. Proteobacteria, Betaproteobacteria, Comamonadaceae, Burkholdariales, and Arcobacter were the dominated phyla, class, family, order, genera, and species, respectively. The Kyoto Encyclopedia of Genes and Genomes indicated the highest number of genes involved in metabolism. The presence of carbohydrate metabolism genes and enzymes was used to infer organic pollutions from sewage and agricultural runoffs. Similarly, the presence of xylene and nutrotoluene degradation genes and enzyme was used to infer industrial pollution into the lake. Drug metabolism genes lend credence to the possibility of pharmaceutical pollutants in water. Taken together, there is a clear indication of massive pollution. In addition, carbohydrate-active enzymes were the most abundant and included genes in glycoside hydrolases. Shotgun metagenomic analyses conveyed an understanding of the microbial communities of the massively polluted aquatic ecosystem of Winam Gulf, Lake Vicoria, Kenya. The current study documents the presence of multiclass pollutants in Lake Victoria and reveals information that might be useful for a potential bioremediation strategy using the native microbial communities.

Identification of motif-based interactions between SARS-CoV-2 protein domains and human peptide ligands pinpoint antiviral targets.

The virus life cycle depends on host-virus protein-protein interactions, which often involve a disordered protein region binding to a folded protein domain. Here, we used proteomic peptide phage display (ProP-PD) to identify peptides from the intrinsically disordered regions of the human proteome that bind to folded protein domains encoded by the SARS-CoV-2 genome. Eleven folded domains of SARS-CoV-2 proteins were found to bind 281 peptides from human proteins, and affinities of 31 interactions involving eight SARS-CoV-2 protein domains were determined (KD ∼ 7-300 µM). Key specificity residues of the peptides were established for six of the interactions. Two of the peptides, binding Nsp9 and Nsp16, respectively, inhibited viral replication. Our findings demonstrate how high-throughput peptide binding screens simultaneously identify potential host-virus interactions and peptides with antiviral properties. Furthermore, the high number of low-affinity interactions suggest that overexpression of viral proteins during infection may perturb multiple cellular pathways.

Interactions outside the boundaries of the canonical binding groove of a PDZ domain influence ligand binding.

The postsynaptic density protein-95/discs large/zonula occludens-1 (PDZ) domain is a protein-protein interaction module with a shallow binding groove where protein ligands bind. However, interactions that are not part of this canonical binding groove are likely to modulate peptide binding. We have investigated such interactions beyond the binding groove for PDZ3 from PSD-95 and a peptide derived from the C-terminus of the natural ligand CRIPT. We found via nuclear magnetic resonance experiments that up to eight residues of the peptide ligand interact with the PDZ domain, showing that the interaction surface extends far outside of the binding groove as defined by the crystal structure. PDZ3 contains an extra structural element, a C-terminal helix (α3), which is known to affect affinity. Deletion of this helix resulted in the loss of several intermolecular nuclear Overhauser enhancements from peptide residues outside of the binding pocket, suggesting that α3 forms part of the extra binding surface in wild-type PDZ3. Site-directed mutagenesis, isothermal titration calorimetry, and fluorescence intensity experiments confirmed the importance of both α3 and the N-terminal part of the peptide for the affinity. Our data suggest a general mechanism in which different binding surfaces outside of the PDZ binding groove could provide sites for specific interactions.

Biophysical characterization of the complex between human papillomavirus E6 protein and synapse-associated protein 97.

The E6 protein of human papillomavirus (HPV) exhibits complex interaction patterns with several host proteins, and their roles in HPV-mediated oncogenesis have proved challenging to study. Here we use several biophysical techniques to explore the binding of E6 to the three PDZ domains of the tumor suppressor protein synapse-associated protein 97 (SAP97). All of the potential binding sites in SAP97 bind E6 with micromolar affinity. The dissociation rate constants govern the different affinities of HPV16 and HPV18 E6 for SAP97. Unexpectedly, binding is not mutually exclusive, and all three PDZ domains can simultaneously bind E6. Intriguingly, this quaternary complex has the same apparent hydrodynamic volume as the unliganded PDZ region, suggesting that a conformational change occurs in the PDZ region upon binding, a conclusion supported by kinetic experiments. Using NMR, we discovered a new mode of interaction between E6 and PDZ: a subset of residues distal to the canonical binding pocket in the PDZ(2) domain exhibited noncanonical interactions with the E6 protein. This is consistent with a larger proportion of the protein surface defining binding specificity, as compared with that reported previously.

Sequence-specific long range networks in PSD-95/discs large/ZO-1 (PDZ) domains tune their binding selectivity.

Protein-protein interactions mediated by modular protein domains are critical for cell scaffolding, differentiation, signaling, and ultimately, evolution. Given the vast number of ligands competing for binding to a limited number of domain families, it is often puzzling how specificity can be achieved. Selectivity may be modulated by intradomain allostery, whereby a remote residue is energetically connected to the functional binding site via side chain or backbone interactions. Whereas several energetic pathways, which could mediate intradomain allostery, have been predicted in modular protein domains, there is a paucity of experimental data to validate their existence and roles. Here, we have identified such functional energetic networks in one of the most common protein-protein interaction modules, the PDZ domain. We used double mutant cycles involving site-directed mutagenesis of both the PDZ domain and the peptide ligand, in conjunction with kinetics to capture the fine energetic details of the networks involved in peptide recognition. We performed the analysis on two homologous PDZ-ligand complexes and found that the energetically coupled residues differ for these two complexes. This result demonstrates that amino acid sequence rather than topology dictates the allosteric pathways. Furthermore, our data support a mechanism whereby the whole domain and not only the binding pocket is optimized for a specific ligand. Such cross-talk between binding sites and remote residues may be used to fine tune target selectivity.

Side-chain interactions form late and cooperatively in the binding reaction between disordered peptides and PDZ domains.

Intrinsically disordered proteins are very common and mediate numerous protein-protein and protein-DNA interactions. While it is clear that these interactions are instrumental for the life of the mammalian cell, there is a paucity of data regarding their molecular binding mechanisms. Here we have used short peptides as a model system for intrinsically disordered proteins. Linear free energy relationships based on rate and equilibrium constants for the binding of these peptides to ordered target proteins, PDZ domains, demonstrate that native side-chain interactions form mainly after the rate-limiting barrier for binding and in a cooperative fashion. This finding suggests that these disordered peptides first form a weak encounter complex with non-native interactions. The data do not support the recent notion that the affinities of intrinsically disordered proteins toward their targets are generally governed by their association rate constants. Instead, we observed the opposite for peptide-PDZ interactions, namely, that changes in K(d) correlate with changes in k(off).

An expanded view of the protein folding landscape of PDZ domains.

Most protein domains fold in an apparently co-operative and two-state manner with only the native and denatured states significantly populated at any experimental condition. However, the protein folding energy landscape is often rugged and different transition states may be rate limiting for the folding reaction under different conditions, as seen for the PDZ protein domain family. We have here analyzed the folding kinetics of two PDZ domains and found that a previously undetected third transition state is rate limiting under conditions that stabilize the native state relative to the denatured state. In light of these results, we have re-analyzed previous folding data on PDZ domains and present a unified folding mechanism with three distinct transition states separated by two high-energy intermediates. Our data show that sequence composition tunes the relative stabilities of folding transition states within the PDZ family, while the overall mechanism is determined by topology. This model captures the kinetic folding mechanism of all PDZ domains studied to date.

Solution nuclear magnetic resonance spectroscopy of bacterial outer membrane proteins in natively excreted vesicles using engineered Escherichia coli.

Gaining structural information on membrane proteins in their native lipid environment is a long-standing challenge in molecular biology. Instead, it is common to employ membrane mimetics, which has been shown to affect protein structure, dynamics, and function severely. Here, we describe the incorporation of a bacterial outer membrane protein (OmpW) into natively excreted membrane vesicles for solution nuclear magnetic resonance (NMR) spectroscopy using a mutant Escherichia coli strain with a high outer membrane vesicle (OMV) production rate. We collected NMR spectra from both vesicles containing overexpressed OmpW and vesicles from a control strain to account for the presence of physiologically relevant outer membrane proteins in vesicles and observed distinct resonance signals from OmpW. Due to the increased production of OMVs and the use of non-uniform sampling techniques we were able to obtain high-resolution 2D (HSQC) and 3D (HNCO) NMR spectra of our target protein inside its native lipid environment. While this workflow is not yet sufficient to achieve in situ structure determination, our results pave the way for further research on vesicle-based solution NMR spectroscopy.

Structural Characterization of a Thorarchaeota Profilin Indicates Eukaryotic-Like Features but with an Extended N-Terminus.

The emergence of the first eukaryotic cell is preceded by evolutionary events, which are still highly debatable. Clues of the exact sequence of events are beginning to emerge. Recent metagenomics analyses has uncovered the Asgard super-phylum as the closest yet known archaea host of eukaryotes. Some of these have been tested and confirmed experimentally. However, the bulk of eukaryotic signature proteins predicted to be encoded by the Asgard super-phylum have not been studied, and their true functions, at least in the context of a eukaryotic cell, are still elusive. For example, there are several different variants of the profilin within each Asgardian Achaea, and there are some conflicting results of their actual roles. Here, the 3D structure of profilin from Thorarchaeota is determined by nuclear magnetic resonance spectroscopy and shows that this profilin has a eukaryotic-like profilin with a rigid core and an extended N-terminus previously implicated in polyproline binding. In addition, it is also shown that Thorarchaeota Profilin co-localizes with eukaryotic actin in cultured HeLa cells. This finding reaffirms the notion that Asgardian encoded proteins possess eukaryotic-like characteristics and strengthen the likely existence of a complex cytoskeleton already in a last eukaryotic common ancestor.