Structural, energetic, and dynamic responses of the native state ensemble of staphylococcal nuclease to cavity-creating mutations

Diverse models of cavity engineering in enzyme modification: Creation, filling, and reshaping

Most enzyme modification strategies focus on designing the active sites or their surrounding structures. Interestingly, a large portion of the enzymes (60%) feature active sites located within spacious cavities. Despite recent discoveries, cavity-mediated enzyme engineering remains crucial for enhancing enzyme properties and unraveling folding-unfolding mechanisms. Cavity engineering influences enzyme stability, catalytic activity, specificity, substrate recognition, and docking. This article provides a comprehensive review of various cavity engineering models for enzyme modification, including cavity creation, filling, and reshaping. Additionally, it also discusses feasible tools for geometric analysis, functional assessment, and modification of cavities, and explores potential future research directions in this field. Furthermore, a promising universal modification strategy for cavity engineering that leverages state-of-the-art technologies and methodologies to tailor cavities according to the specific requirements of industrial production conditions is proposed.

Principles and practical applications of structure-based vaccine design

Viral proteins fold into a variety of structures as they perform their functions. Structure-based vaccine design aims to exploit knowledge of an antigen’s architecture to stabilize it in a vulnerable conformation. We summarize the general principles of structure-based vaccine design, with a focus on the major types of sequence modifications: proline, disulfide, cavity-filling, electrostatic and hydrogen-bond substitution, as well as domain deletion. We then review recent applications of these principles to vaccine-design efforts across five viral families: Coronaviridae, Orthomyxoviridae, Paramyxoviridae, Pneumoviridae, and Filoviridae. Outstanding challenges include continued application of proven design principles to pathogens of interest, as well as development of new strategies for those pathogens that resist traditional techniques.

A hierarchy of coupling free energies underlie the thermodynamic and functional architecture of protein structures

Protein sequences and structures evolve by satisfying varied physical and biochemical constraints. This multi-level selection is enabled not just by the patterning of amino acids on the sequence, but also via coupling between residues in the native structure. Here, we employ an energetically detailed statistical mechanical model with millions of microstates to extract such long-range structural correlations, i.e. thermodynamic coupling free energies, from a diverse family of protein structures. We find that despite the intricate and anisotropic distribution of coupling patterns, the majority of residues (>70%) are only marginally coupled contributing to functional motions and catalysis. Physical origins of ‘sectors’, determinants of native ensemble heterogeneity in extant, ancient and designed proteins, and the basis for allostery emerge naturally from coupling free energies. The statistical framework highlights how evolutionary selection and optimization occur at the level of global interaction network for a given protein fold impacting folding, function, and allosteric outputs.

•

Evolution of protein structures occurs at the level of global interaction network.

•

More than 70% of the protein residues are weakly or marginally coupled.

•

Functional ‘sector’ regions are a manifestation of marginal coupling.

•

Coupling indices vary across the entire proteins in extant-ancient and natural-designed pairs.

•

The proposed methodology can be used to understand allostery and epistasis.

Comparative Assessment of NMR Probes for the Experimental Description of Protein Folding Pathways with High-Pressure NMR

Multidimensional NMR intrinsically provides multiple probes that can be used for deciphering the folding pathways of proteins: NH amide and CαHα groups are strategically located on the backbone of the protein, while CH₃ groups, on the side-chain of methylated residues, are involved in important stabilizing interactions in the hydrophobic core. Combined with high hydrostatic pressure, these observables provide a powerful tool to explore the conformational landscapes of proteins. In the present study, we made a comparative assessment of the NH, CαHα, and CH₃ groups for analyzing the unfolding pathway of ∆+PHS Staphylococcal Nuclease. These probes yield a similar description of the folding pathway, with virtually identical thermodynamic parameters for the unfolding reaction, despite some notable differences. Thus, if partial unfolding begins at identical pressure for these observables (especially in the case of backbone probes) and concerns similar regions of the molecule, the residues involved in contact losses are not necessarily the same. In addition, an unexpected slight shift toward higher pressure was observed in the sequence of the scenario of unfolding with CαHα when compared to amide groups.

Diverse Native Ensembles Dictate the Differential Functional Responses of Nuclear Receptor Ligand-Binding Domains

Native states of folded proteins are characterized by a large ensemble of conformations whose relative populations and interconversion dynamics determine the functional output. This is more apparent in transcription factors that have evolved to be inherently sensitive to small perturbations, thus fine-tuning gene expression. To explore the extent to which such functional features are imprinted on the folding landscape of transcription factor ligand-binding domains (LBDs), we characterize paralogous LBDs of the nuclear receptor (NR) family employing an energetically detailed and ensemble-based Ising-like statistical mechanical model. We find that the native ensembles of the LBDs from glucocorticoid receptor, PPAγ, and thyroid hormone receptor display a remarkable diversity in the width of the native wells, the number and nature of partially structured states, and hence the degree of conformational order. Monte Carlo simulations employing the full state representation of the ensemble highlight that many of the functional conformations coexist in equilibrium, whose relative populations are sensitive to both temperature and the strength of ligand binding. Allosteric modulation of the degree of structure at a coregulator binding site on ligand binding is shown to arise via a redistribution of populations in the native ensembles of glucocorticoid and PPAγ LBDs. Our results illustrate how functional requirements can drive the evolution of conformationally diverse native ensembles in paralogs.

Genetically detoxified pertussis toxin displays near identical structure to its wild-type and exhibits robust immunogenicity

The mutant gdPT R9K/E129G is a genetically detoxified variant of the pertussis toxin (PTx) and represents an attractive candidate for the development of improved pertussis vaccines. The impact of the mutations on the overall protein structure and its immunogenicity has remained elusive. Here we present the crystal structure of gdPT and show that it is nearly identical to that of PTx. Hydrogen-deuterium exchange mass spectrometry revealed dynamic changes in the catalytic domain that directly impacted NAD⁺ binding which was confirmed by biolayer interferometry. Distal changes in dynamics were also detected in S2-S5 subunit interactions resulting in tighter packing of B-oligomer corresponding to increased thermal stability. Finally, antigen stimulation of human whole blood, analyzed by a previously unreported mass cytometry assay, indicated broader immunogenicity of gdPT compared to pertussis toxoid. These findings establish a direct link between the conserved structure of gdPT and its ability to generate a robust immune response.

Controlling Structure and Dimensions of a Disordered Protein via Mutations

The dimensions of intrinsically disordered proteins (IDPs) are sensitive to small energetic-entropic differences between intramolecular and protein–solvent interactions. This is commonly observed on modulating solvent composition and temperature. However, the inherently heterogeneous conformational landscape of IDPs is also expected to be influenced by mutations that can (de)stabilize pockets of local and even global structure, native and non-native, and hence the average dimensions. Here, we show experimental evidence for the remarkably tunable landscape of IDPs by employing the DNA-binding domain of CytR, a high-sequence-complexity IDP, as a model system. CytR exhibits a range of structure and compactness upon introducing specific mutations that modulate microscopic terms, including main-chain entropy, hydrophobicity, and electrostatics. The degree of secondary structure, as monitored by far-UV circular dichroism (CD), is strongly correlated to average ensemble dimensions for 14 different mutants of CytR and is consistent with the Uversky–Fink relation. Our experiments highlight how average ensemble dimensions can be controlled via mutations even in the disordered regime, the prevalence of non-native interactions and provide testable controls for molecular simulations.

Engineering Order and Cooperativity in a Disordered Protein

Structural disorder in proteins arises from a complex interplay between weak hydrophobicity and unfavorable electrostatic interactions. The extent to which the hydrophobic effect contributes to the unique and compact native state of proteins is, however, confounded by large compensation between multiple entropic and energetic terms. Here we show that protein structural order and cooperativity arise as emergent properties upon hydrophobic substitutions in a disordered system with non-intuitive effects on folding and function. Aided by sequence-structure analysis, equilibrium, and kinetic spectroscopic studies, we engineer two hydrophobic mutations in the disordered DNA-binding domain of CytR that act synergistically, but not in isolation, to promote structure, compactness, and stability. The double mutant, with properties of a fully ordered domain, exhibits weak cooperativity with a complex and rugged conformational landscape. The mutant, however, binds cognate DNA with an affinity only marginally higher than that of the wild type, though nontrivial differences are observed in the binding to noncognate DNA. Our work provides direct experimental evidence of the dominant role of non-additive hydrophobic effects in shaping the molecular evolution of order in disordered proteins and vice versa, which could be generalized to even folded proteins with implications for protein design and functional manipulation.

Lessons from pressure denaturation of proteins

Although it is now relatively well understood how sequence defines and impacts global protein stability in specific structural contexts, the question of how sequence modulates the configurational landscape of proteins remains to be defined. Protein configurational equilibria are generally characterized by using various chemical denaturants or by changing temperature or pH. Another thermodynamic parameter which is less often used in such studies is high hydrostatic pressure. This review discusses the basis for pressure effects on protein structure and stability, and describes how the unique mechanisms of pressure-induced unfolding can provide unique insights into protein conformational landscapes.

Modulation of allosteric coupling by mutations: from protein dynamics and packing to altered native ensembles and function

A large body of work has gone into understanding the effect of mutations on protein structure and function. Conventional treatments have involved quantifying the change in stability, activity and relaxation rates of the mutants with respect to the wild-type protein. However, it is now becoming increasingly apparent that mutational perturbations consistently modulate the packing and dynamics of a significant fraction of protein residues, even those that are located >10–15 Å from the mutated site. Such long-range modulation of protein features can distinctly tune protein stability and the native conformational ensemble contributing to allosteric modulation of function. In this review, I summarize a series of experimental and computational observations that highlight the incredibly pliable nature of proteins and their response to mutational perturbations manifested via the intra-protein interaction network. I highlight how an intimate understanding of mutational effects could pave the way for integrating stability, folding, cooperativity and even allostery within a single physical framework.

Practical aspects of high-pressure NMR spectroscopy and its applications in protein biophysics and structural biology

Pressure and temperature are the two fundamental variables of thermodynamics. Temperature and chemical perturbation are central experimental tools for the exploration of macromolecular structure and dynamics. Though it has long been recognized that hydrostatic pressure offers a complementary and often unique view of macromolecular structure, stability and dynamics, it has not been employed nearly as much. For solution NMR applications the limited use of high-pressure is undoubtedly traced to difficulties of employing pressure in the context of modern multinuclear and multidimensional NMR. Limitations in pressure tolerant NMR sample cells have been overcome and enable detailed studies of macromolecular energy landscapes, hydration, dynamics and function. Here we review the practical considerations for studies of biological macromolecules at elevated pressure, with a particular emphasis on applications in protein biophysics and structural biology.

The consequences of cavity creation on the folding landscape of a repeat protein depend upon context

The effect of introducing internal cavities on protein native structure and global stability has been well documented, but the consequences of these packing defects on folding free-energy landscapes have received less attention. We investigated the effects of cavity creation on the folding landscape of the leucine-rich repeat protein pp32 by high-pressure (HP) and urea-dependent NMR and high-pressure small-angle X-ray scattering (HPSAXS). Despite a modest global energetic perturbation, cavity creation in the N-terminal capping motif (N-cap) resulted in very strong deviation from two-state unfolding behavior. In contrast, introduction of a cavity in the most stable, C-terminal half of pp32 led to highly concerted unfolding, presumably because the decrease in stability by the mutations attenuated the N- to C-terminal stability gradient present in WT pp32. Interestingly, enlarging the central cavity of the protein led to the population under pressure of a distinct intermediate in which the N-cap and repeats 1-4 were nearly completely unfolded, while the fifth repeat and the C-terminal capping motif remained fully folded. Thus, despite modest effects on global stability, introducing internal cavities can have starkly distinct repercussions on the conformational landscape of a protein, depending on their structural and energetic context.

Monitoring protein folding through high pressure NMR spectroscopy

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.

A Universal Pattern in the Percolation and Dissipation of Protein Structural Perturbations

Understanding the extent to which information is transmitted through the intramolecular interaction network of proteins upon a perturbation, that is, an allosteric effect, has long remained an unsolved problem. Through an analysis of high-resolution NMR data from the literature on 28 different proteins and 49 structural perturbations, we show that the extent of induced structural changes through mutations and molecular events including protein-protein, protein-peptide, protein-ligand binding, and post-translational modifications exhibit a near-universal exponential functional form. The extent of percolation into the protein structures can be up to 20-25 Å despite no apparent change in the 3D structures. These observations are also consistent with theoretical expectations, elementary graph theoretic analysis of protein structures, detailed molecular dynamics simulations, and experimental double-mutant cycles. Our analysis highlights that most molecular events would contribute to allosteric effects independent of protein structure, topology, or identity and provides a simple avenue to test and potentially model their effects.

V67L Mutation Fills an Internal Cavity To Stabilize RecA Mtu Intein

Inteins mediate protein splicing, which has found extensive applications in protein science and biotechnology. In the Mycobacterium tuberculosis RecA mini-mini intein (ΔΔIhh), a single valine to leucine substitution at position 67 (V67L) dramatically increases intein stability and activity. However, crystal structures show that the V67L mutation causes minimal structural rearrangements, with a root-mean-square deviation of 0.2 Å between ΔΔIhh-V67 and ΔΔIhh-L67. Thus, the structural mechanisms for V67L stabilization and activation remain poorly understood. In this study, we used intrinsic tryptophan fluorescence, high-pressure nuclear magnetic resonance (NMR), and molecular dynamics (MD) simulations to probe the structural basis of V67L stabilization of the intein fold. Guanidine hydrochloride denaturation monitored by fluorescence yielded free energy changes (ΔG_f°) of -4.4 and -6.9 kcal mol^-1 for ΔΔIhh-V67 and ΔΔIhh-L67, respectively. High-pressure NMR showed that ΔΔIhh-L67 is more resistant to pressure-induced unfolding than ΔΔIhh-V67 is. The change in the volume of folding (ΔV_f) was significantly larger for V67 (71 ± 2 mL mol^-1) than for L67 (58 ± 3 mL mol^-1) inteins. The measured difference in ΔV_f (13 ± 3 mL mol^-1) roughly corresponds to the volume of the additional methylene group for Leu, supporting the notion that the V67L mutation fills a nearby cavity to enhance intein stability. In addition, we performed MD simulations to show that V67L decreases side chain dynamics and conformational entropy at the active site. It is plausible that changes in cavities in V67L can also mediate allosteric effects to change active site dynamics and enhance intein activity.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Evolutionarily Conserved Pattern of Interactions in a Protein Revealed by Local Thermal Expansion Properties

The way in which the network of intramolecular interactions determines the cooperative folding and conformational dynamics of a protein remains poorly understood. High-pressure NMR spectroscopy is uniquely suited to examine this problem because it combines the site-specific resolution of the NMR experiments with the local character of pressure perturbations. Here we report on the temperature dependence of the site-specific volumetric properties of various forms of staphylococcal nuclease (SNase), including three variants with engineered internal cavities, as measured with high-pressure NMR spectroscopy. The strong temperature dependence of pressure-induced unfolding arises from poorly understood differences in thermal expansion between the folded and unfolded states. A significant inverse correlation was observed between the global thermal expansion of the folded proteins and the number of strong intramolecular hydrogen bonds, as determined by the temperature coefficient of the backbone amide chemical shifts. Comparison of the identity of these strong H-bonds with the co-evolution of pairs of residues in the SNase protein family suggests that the architecture of the interactions detected in the NMR experiments could be linked to a functional aspect of the protein. Moreover, the temperature dependence of the residue-specific volume changes of unfolding yielded residue-specific differences in expansivity and revealed how mutations impact intramolecular interaction patterns. These results show that intramolecular interactions in the folded states of proteins impose constraints against thermal expansion and that, hence, knowledge of site-specific thermal expansivity offers insight into the patterns of strong intramolecular interactions and other local determinants of protein stability, cooperativity, and potentially also of function.

Exploring the Protein Folding Pathway with High-Pressure NMR: Steady-State and Kinetics Studies

Defining the physical-chemical determinants of protein folding and stability, under normal and pathological conditions has constituted a major subfield in biophysical chemistry for over 50 years. Although a great deal of progress has been made in recent years towards this goal, a number of important questions remain. These include characterizing the structural, thermodynamic and dynamic properties of the barriers between conformational states on the protein energy landscape, understanding the sequence dependence of folding cooperativity, defining more clearly the role of solvation in controlling protein stability and dynamics and probing the high energy thermodynamic states in the native state basin and their role in misfolding and aggregation. Fundamental to the elucidation of these questions is a complete thermodynamic parameterization of protein folding determinants. In this chapter, we describe the use of high-pressure coupled to Nuclear Magnetic Resonance (NMR) spectroscopy to reveal unprecedented details on the folding energy landscape of proteins.

Computational prediction and experimental characterization of a "size switch type repacking" during the evolution of dengue envelope protein domain III (ED3)

Dengue viruses (DEN) are classified into four serotypes (DEN1-DEN4) exhibiting high sequence and structural similarities, and infections by multiple serotypes can lead to the deadly dengue hemorrhagic fever. Here, we aim at characterizing the thermodynamic stability of DEN envelope protein domain III (ED3) during its evolution, and we report a structural analysis of DEN4wt ED3 combined with a systematic mutational analysis of residues 310 and 387. Molecular modeling based on our DEN3 and DEN4 ED3 structures indicated that the side-chains of residues 310/387, which are Val(310)/Ile(387) and Met(310)/Leu(387) in DEN3wt and DEN4wt, respectively, could be structurally compensated, and that a "size switch type repacking" might have occurred at these sites during the evolution of DEN into its four serotypes. This was experimentally confirmed by a 10°C and 5°C decrease in the thermal stability of, respectively, DEN3 ED3 variants with Met(310)/Ile(387) and Val(310)/Leu(387), whereas the variant with Met(310)/Leu(387), which contains a double mutation, had the same stability as the wild type DEN3. Namely, the Met310Val mutation should have preceded the Leu387Ile mutation in order to maintain the tight internal packing of ED3 and thus its thermodynamic stability. This view was confirmed by a phylogenetic reconstruction indicating that a common DEN ancestor would have Met(310)/Leu(387), and the intermediate node protein, Val(310)/Leu(387), which then mutated to the Val(310)/Ile(387) pair found in the present DEN3. The hypothesis was further confirmed by the observation that all of the present DEN viruses exhibit only stabilizing amino acid pairs at the 310/387 sites.

Conformational selection and adaptation to ligand binding in T4 lysozyme cavity mutants

The studies presented here explore the relationship between protein packing and molecular flexibility using ligand-binding cavity mutants of T4 lysozyme. Although previously reported crystal structures of the mutants investigated show single conformations that are similar to the WT protein, site-directed spin labeling in solution reveals additional conformational substates in equilibrium exchange with a WT-like population. Remarkably, binding of ligands, including the general anesthetic halothane shifts the population to the WT-like state, consistent with a conformational selection model of ligand binding, but structural adaptation to the ligand is also apparent in one mutant. Distance mapping with double electron-electron resonance spectroscopy and the absence of ligand binding suggest that the new substates induced by the cavity-creating mutations represent alternate packing modes in which the protein fills or partially fills the cavity with side chains, including the spin label in one case; external ligands compete with the side chains for the cavity space, stabilizing the WT conformation. The results have implications for mechanisms of anesthesia, the response of proteins to hydrostatic pressure, and protein engineering.