The structural landscape and diversity of Pyricularia oryzae MAX effectors revisited.

Magnaporthe AVRs and ToxB-like (MAX) effectors constitute a family of secreted virulence proteins in the fungus Pyricularia oryzae (syn. Magnaporthe oryzae), which causes blast disease on numerous cereals and grasses. In spite of high sequence divergence, MAX effectors share a common fold characterized by a ß-sandwich core stabilized by a conserved disulfide bond. In this study, we investigated the structural landscape and diversity within the MAX effector repertoire of P. oryzae. Combining experimental protein structure determination and in silico structure modeling we validated the presence of the conserved MAX effector core domain in 77 out of 94 groups of orthologs (OG) identified in a previous population genomic study. Four novel MAX effector structures determined by NMR were in remarkably good agreement with AlphaFold2 (AF2) predictions. Based on the comparison of the AF2-generated 3D models we propose a classification of the MAX effectors superfamily in 20 structural groups that vary in the canonical MAX fold, disulfide bond patterns, and additional secondary structures in N- and C-terminal extensions. About one-third of the MAX family members remain singletons, without strong structural relationship to other MAX effectors. Analysis of the surface properties of the AF2 MAX models also highlights the high variability within the MAX family at the structural level, potentially reflecting the wide diversity of their virulence functions and host targets.

Does a Similar 3D Structure Mean a Similar Folding Pathway? The Presence of a C-Terminal α-Helical Extension in the 3D Structure of MAX60 Drastically Changes the Folding Pathway Described for Other MAX-Effectors from Magnaporthe oryzae.

Does a similar 3D structure mean a similar folding pathway? This question is particularly meaningful when it concerns proteins sharing a similar 3D structure, but low sequence identity or homology. MAX effectors secreted by the phytopathogenic fungus Magnaporthe oryzae present such characteristics. They share a common 3D structure, a ß-sandwich with the same topology for all the family members, but an extremely low sequence identity/homology. In a previous study, we have investigated the folding of two MAX effectors, AVR-Pia and AVR-Pib, using High-Hydrostatic-Pressure NMR and found that they display a similar folding pathway, with a common folding intermediate. In the present work, we used a similar strategy to investigate the folding conformational landscape of another MAX effector, MAX60, and found a very different folding intermediate. Our analysis strongly supports that the presence of a C-terminal α-helical extension in the 3D structure of MAX60 could be responsible for its different folding pathway.

Combining High-Pressure NMR and Geometrical Sampling to Obtain a Full Topological Description of Protein Folding Landscapes: Application to the Folding of Two MAX Effectors from Magnaporthe oryzae.

Despite advances in experimental and computational methods, the mechanisms by which an unstructured polypeptide chain regains its unique three-dimensional structure remains one of the main puzzling questions in biology. Single-molecule techniques, ultra-fast perturbation and detection approaches and improvement in all-atom and coarse-grained simulation methods have greatly deepened our understanding of protein folding and the effects of environmental factors on folding landscape. However, a major challenge remains the detailed characterization of the protein folding landscape. Here, we used high hydrostatic pressure 2D NMR spectroscopy to obtain high-resolution experimental structural information in a site-specific manner across the polypeptide sequence and along the folding reaction coordinate. We used this residue-specific information to constrain Cyana3 calculations, in order to obtain a topological description of the entire folding landscape. This approach was used to describe the conformers populating the folding landscape of two small globular proteins, AVR-Pia and AVR-Pib, that belong to the structurally conserved but sequence-unrelated MAX effectors superfamily. Comparing the two folding landscapes, we found that, in spite of their divergent sequences, the folding pathway of these two proteins involves a similar, inescapable, folding intermediate, even if, statistically, the routes used are different.

Unravelling the Adaptation Mechanisms to High Pressure in Proteins.

Life is thought to have appeared in the depth of the sea under high hydrostatic pressure. Nowadays, it is known that the deep biosphere hosts a myriad of life forms thriving under high-pressure conditions. However, the evolutionary mechanisms leading to their adaptation are still not known. Here, we show the molecular bases of these mechanisms through a joint structural and dynamical study of two orthologous proteins. We observed that pressure adaptation involves the decoupling of protein-water dynamics and the elimination of cavities in the protein core. This is achieved by rearranging the charged residues on the protein surface and using bulkier hydrophobic residues in the core. These findings will be the starting point in the search for a complete genomic model explaining high-pressure adaptation.

Pressure and Chemical Unfolding of an α-Helical Bundle Protein: The GH2 Domain of the Protein Adaptor GIPC1.

When combined with NMR spectroscopy, high hydrostatic pressure is an alternative perturbation method used to destabilize globular proteins that has proven to be particularly well suited for exploring the unfolding energy landscape of small single-domain proteins. To date, investigations of the unfolding landscape of all-ß or mixed-α/ß protein scaffolds are well documented, whereas such data are lacking for all-α protein domains. Here we report the NMR study of the unfolding pathways of GIPC1-GH2, a small α-helical bundle domain made of four antiparallel α-helices. High-pressure perturbation was combined with NMR spectroscopy to unravel the unfolding landscape at three different temperatures. The results were compared to those obtained from classical chemical denaturation. Whatever the perturbation used, the loss of secondary and tertiary contacts within the protein scaffold is almost simultaneous. The unfolding transition appeared very cooperative when using high pressure at high temperature, as was the case for chemical denaturation, whereas it was found more progressive at low temperature, suggesting the existence of a complex folding pathway.

Protein kinase B controls Mycobacterium tuberculosis growth via phosphorylation of the transcriptional regulator Lsr2 at threonine 112.

Mycobacterium tuberculosis (Mtb) is able to persist in the body through months of multi-drug therapy. Mycobacteria possess a wide range of regulatory proteins, including the protein kinase B (PknB) which controls peptidoglycan biosynthesis during growth. Here, we observed that depletion of PknB resulted in specific transcriptional changes that are likely caused by reduced phosphorylation of the H-NS-like regulator Lsr2 at threonine 112. The activity of PknB towards this phosphosite was confirmed with purified proteins, and this site was required for adaptation of Mtb to hypoxic conditions, and growth on solid media. Like H-NS, Lsr2 binds DNA in sequence-dependent and non-specific modes. PknB phosphorylation of Lsr2 reduced DNA binding, measured by fluorescence anisotropy and electrophoretic mobility shift assays, and our NMR structure of phosphomimetic T112D Lsr2 suggests that this may be due to increased dynamics of the DNA-binding domain. Conversely, the phosphoablative T112A Lsr2 had increased binding to certain DNA sites in ChIP-sequencing, and Mtb containing this variant showed transcriptional changes that correspond with the change in DNA binding. In summary, PknB controls Mtb growth and adaptations to the changing host environment by phosphorylating the global transcriptional regulator Lsr2.

Combining High-Pressure Perturbation with NMR Spectroscopy for a Structural and Dynamical Characterization of Protein Folding Pathways.

High-hydrostatic pressure is an alternative perturbation method that can be used to destabilize globular proteins. Generally perfectly reversible, pressure exerts local effects on regions or domains of a protein containing internal voids, contrary to heat or chemical denaturant that destabilize protein structures uniformly. When combined with NMR spectroscopy, high pressure (HP) allows one to monitor at a residue-level resolution the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. The use of HP-NMR has long been hampered by technical difficulties. Owing to the recent development of commercially available high-pressure sample cells, HP-NMR experiments can now be routinely performed. This review summarizes recent advances of HP-NMR techniques for the characterization at a quasi-atomic resolution of the protein folding energy landscape.

Monitoring Unfolding of Titin I27 Single and Bi Domain with High-Pressure NMR Spectroscopy.

A complete description of the pathways and mechanisms of protein folding requires a detailed structural and energetic characterization of the folding energy landscape. Simulations, when corroborated by experimental data yielding global information on the folding process, can provide this level of insight. Molecular dynamics (MD) has often been combined with force spectroscopy experiments to decipher the unfolding mechanism of titin immunoglobulin-like single or multidomain, the giant multimodular protein from sarcomeres, yielding information on the sequential events during titin unfolding under stretching. Here, we used high-pressure NMR to monitor the unfolding of titin I27 Ig-like single domain and tandem. Because this method brings residue-specific information on the folding process, it can provide quasiatomic details on this process without the help of MD simulations. Globally, the results of our high-pressure analysis are in agreement with previous results obtained by the combination of experimental measurements and MD simulation and/or protein engineering, although the intermediate folding state caused by the early detachment of the AB ß-sheet, often reported in previous works based on MD or force spectroscopy, cannot be detected. On the other hand, the A'G parallel ß-sheet of the ß-sandwich has been confirmed as the Achilles heel of the three-dimensional scaffold: its disruption yields complete unfolding with very similar characteristics (free energy, unfolding volume, kinetics rate constants) for the two constructs.

High-Resolution Mapping of a Repeat Protein Folding Free Energy Landscape.

A complete description of the pathways and mechanisms of protein folding requires a detailed structural and energetic characterization of the conformational ensemble along the entire folding reaction coordinate. Simulations can provide this level of insight for small proteins. In contrast, with the exception of hydrogen exchange, which does not monitor folding directly, experimental studies of protein folding have not yielded such structural and energetic detail. NMR can provide residue specific atomic level structural information, but its implementation in protein folding studies using chemical or temperature perturbation is problematic. Here we present a highly detailed structural and energetic map of the entire folding landscape of the leucine-rich repeat protein, pp32 (Anp32), obtained by combining pressure-dependent site-specific 1H-15N HSQC data with coarse-grained molecular dynamics simulations. The results obtained using this equilibrium approach demonstrate that the main barrier to folding of pp32 is quite broad and lies near the unfolded state, with structure apparent only in the C-terminal region. Significant deviation from two-state unfolding under pressure reveals an intermediate on the folded side of the main barrier in which the N-terminal region is disordered. A nonlinear temperature dependence of the population of this intermediate suggests a large heat capacity change associated with its formation. The combination of pressure, which favors the population of folding intermediates relative to chemical denaturants; NMR, which allows their observation; and constrained structure-based simulations yield unparalleled insight into protein folding mechanisms.

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.

Prioritizing targets for structural biology through the lens of proteomics: the archaeal protein TGAM_1934 from Thermococcus gammatolerans.

ORFans are hypothetical proteins lacking any significant sequence similarity with other proteins. Here, we highlighted by quantitative proteomics the TGAM_1934 ORFan from the hyperradioresistant Thermococcus gammatolerans archaeon as one of the most abundant hypothetical proteins. This protein has been selected as a priority target for structure determination on the basis of its abundance in three cellular conditions. Its solution structure has been determined using multidimensional heteronuclear NMR spectroscopy. TGAM_1934 displays an original fold, although sharing some similarities with the 3D structure of the bacterial ortholog of frataxin, CyaY, a protein conserved in bacteria and eukaryotes and involved in iron-sulfur cluster biogenesis. These results highlight the potential of structural proteomics in prioritizing ORFan targets for structure determination based on quantitative proteomics data. The proteomic data and structure coordinates have been deposited to the ProteomeXchange with identifier PXD000402 (http://proteomecentral.proteomexchange.org/dataset/PXD000402) and Protein Data Bank under the accession number 2mcf, respectively.

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.

Cavities determine the pressure unfolding of proteins.

It has been known for nearly 100 years that pressure unfolds proteins, yet the physical basis of this effect is not understood. Unfolding by pressure implies that the molar volume of the unfolded state of a protein is smaller than that of the folded state. This decrease in volume has been proposed to arise from differences between the density of bulk water and water associated with the protein, from pressure-dependent changes in the structure of bulk water, from the loss of internal cavities in the folded states of proteins, or from some combination of these three factors. Here, using 10 cavity-containing variants of staphylococcal nuclease, we demonstrate that pressure unfolds proteins primarily as a result of cavities that are present in the folded state and absent in the unfolded one. High-pressure NMR spectroscopy and simulations constrained by the NMR data were used to describe structural and energetic details of the folding landscape of staphylococcal nuclease that are usually inaccessible with existing experimental approaches using harsher denaturants. Besides solving a 100-year-old conundrum concerning the detailed structural origins of pressure unfolding of proteins, these studies illustrate the promise of pressure perturbation as a unique tool for examining the roles of packing, conformational fluctuations, and water penetration as determinants of solution properties of proteins, and for detecting folding intermediates and other structural details of protein-folding landscapes that are invisible to standard experimental approaches.

Volumetric properties underlying ligand binding in a monomeric hemoglobin: a high-pressure NMR study.

The 2/2 hemoglobin of the cyanobacterium Synechococcus sp. PCC 7002, GlbN, coordinates the heme iron with two histidines and exists either with a b heme or with a covalently attached heme. The binding of exogenous ligands displaces the distal histidine and induces a conformational rearrangement involving the reorganization of internal void volumes. The formation of passageways within the resulting conformation is thought to facilitate ligand exchange and play a functional role. Here we monitored the perturbation induced by pressure on the ferric bis-histidine and cyanide-bound states of GlbN using (1)H-(15)N HSQC NMR spectroscopy. We inspected the outcome with a statistical analysis of 170 homologous 2/2 hemoglobin sequences. We found that the compression landscape of GlbN, as represented by the variation of an average chemical shift parameter, was highly sensitive to ligand swapping and heme covalent attachment. Stabilization of rare conformers was observed at high pressures and consistent with cavity redistribution upon ligand binding. In all states, the EF loop was found to be exceptionally labile to pressure, suggesting a functional role as a semi-flexible hinge between the adjacent helices. Finally, coevolved clusters presented a common pattern of compensating pressure responses. The high-pressure dissection combined with protein sequence analysis established locations with volumetric signatures relevant to residual communication of 2/2 hemoglobins. This article is part of a Special Issue entitled: Oxygen Binding and Sensing Proteins.

Measuring residual dipolar couplings at high hydrostatic pressure: robustness of alignment media to high pressure.

Among other perturbations, high hydrostatic pressure has proven to be a mild yet efficient way to unfold proteins. Combining pressure perturbation with NMR spectroscopy allows for a residue-per-residue description of folding reactions. Accessing the full power of NMR spectroscopy under pressure involves the investigation of conformational sampling using orientational restraints such as residual dipolar couplings (RDCs) under conditions of partial alignment. The aim of this study was to identify and characterize stable and pressure resistant alignment media for measurement of RDCs at high pressure. Four alignment media were tested. A C12E5/n-hexanol alcohol mixture remains stable from 1 to 2,500 bar, whereas Pf1 phage and DNA nanotubes undergo a reversible transition between 300 and 900 bar. Phospholipid bicelles are stable only until 300 bar at ambient temperature. Hence, RDCs can be measured at high pressure, and their interpretation will provide atomic details of the structural and dynamic perturbations on unfolded or partially folded states of proteins under pressure.

1H, 13C and 15N backbone and side-chain resonance assignments of the human oncogenic protein NCYM.

NCYM is a cis-antisense gene of MYCN oncogene and encodes an oncogenic protein that stabilizes MYCN via inhibition of GSK3b. High NCYM expression levels are associated with poor clinical outcomes in human neuroblastomas, and NCYM overexpression promotes distant metastasis in animal models of neuroblastoma. Using vacuum-ultraviolet circular dichroism and small-angle X-ray scattering, we previously showed that NCYM has high flexibility with partially folded structures; however, further structural characterization is required for the design of anti-cancer agents targeting NCYM. Here we report the 1H, 15N and 13C nuclear magnetic resonance assignments of NCYM. Secondary structure prediction using Secondary Chemical Shifts and TALOS-N analysis demonstrates that the structure of NCYM is essentially disordered, even though residues in the central region of the peptide clearly present a propensity to adopt a dynamic helical structure. This preliminary study provides foundations for further analysis of interaction between NCYM and potential partners.

Structural and functional characterization of Bc28.1, major erythrocyte-binding protein from Babesia canis merozoite surface.

Babesiosis (formerly known as piroplasmosis) is a tick-borne disease caused by the intraerythrocytic development of protozoa parasites from the genus Babesia. Like Plasmodium falciparum, the agent of malaria, or Toxoplasma gondii, responsible for human toxoplasmosis, Babesia belongs to the Apicomplexa family. Babesia canis is the agent of the canine babesiosis in Europe. Clinical manifestations of this disease range from mild to severe and possibly lead to death by multiple organ failure. The identification and characterization of parasite surface proteins represent major goals, both for the understanding of the Apicomplexa invasion process and for the vaccine potential of such antigens. Indeed, we have already shown that Bd37, the major antigenic adhesion protein from Babesia divergens, the agent of bovine babesiosis, was able to induce complete protection against various parasite strains. The major merozoite surface antigens of Babesia canis have been described as a 28-kDa membrane protein family, anchored at the surface of the merozoite. Here, we demonstrate that Bc28.1, a major member of this multigenic family, is expressed at high levels at the surface of the merozoite. This protein is also found in the parasite in vitro culture supernatants, which are the basis of effective vaccines against canine babesiosis. We defined the erythrocyte binding function of Bc28.1 and determined its high resolution solution structure using NMR spectroscopy. Surprisingly, although these proteins are thought to play a similar role in the adhesion process, the structure of Bc28.1 from B. canis appears unrelated to the previously published structure of Bd37 from B. divergens. Site-directed mutagenesis experiments also suggest that the mechanism of the interaction with the erythrocyte membrane could be different for the two proteins. The resolution of the structure of Bc28 represents a milestone for the characterization of the parasite erythrocyte binding and its interaction with the host immune system.

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

The effects of cavity-creating mutations on the structural flexibility, local and global stability, and dynamics of the folded state of staphylococcal nuclease (SNase) were examined with NMR spectroscopy, MD simulations, H/D exchange, and pressure perturbation. Effects on global thermodynamic stability correlated well with the number of heavy atoms in the vicinity of the mutated residue. Variants with substitutions in the C-terminal domain and the interface between α and ß subdomains showed large amide chemical shift variations relative to the parent protein, moderate, widespread, and compensatory perturbations of the H/D protection factors and increased local dynamics on a nanosecond time scale. The pressure sensitivity of the folded states of these variants was similar to that of the parent protein. Such observations point to the capacity of the folded proteins to adjust to packing defects in these regions. In contrast, cavity creation in the ß-barrel subdomain led to minimal perturbation of the structure of the folded state, However, significant pressure dependence of the native state amide resonances, along with strong effects on native state H/D exchange are consistent with increased probability of population of excited state(s) for these variants. Such contrasted responses to the creation of cavities could not be anticipated from global thermodynamic stability or crystal structures; they depend on the local structural and energetic context of the substitutions.

Effect of internal cavities on folding rates and routes revealed by real-time pressure-jump NMR spectroscopy.

The time required to fold proteins usually increases significantly under conditions of high pressure. Taking advantage of this general property of proteins, we combined P-jump experiments with NMR spectroscopy to examine in detail the folding reaction of staphylococcal nuclease (SNase) and of some of its cavity-containing variants. The nearly 100 observables that could be measured simultaneously collectively describe the kinetics of folding as a function of pressure and denaturant concentration with exquisite site-specific resolution. SNase variants with cavities in the central core of the protein exhibit a highly heterogeneous transition-state ensemble (TSE) with a smaller solvent-excluded void volume than the TSE of the parent SNase. This heterogeneous TSE experiences Hammond behavior, becoming more native-like (higher molar volume) with increasing denaturant concentration. In contrast, the TSE of the L125A variant, which has a cavity at the secondary core, is only slightly different from that of the parent SNase. Because pressure acts mainly to eliminate solvent-excluded voids, which are heterogeneously distributed throughout structures, it perturbs the protein more selectively than chemical denaturants, thereby facilitating the characterization of intermediates and the consequences of packing on folding mechanisms. Besides demonstrating how internal cavities can affect the routes and rates of folding of a protein, this study illustrates how the combination of P-jump and NMR spectroscopy can yield detailed mechanistic insight into protein folding reactions with exquisite site-specific temporal information.

Equivalence between Euler angle conventions for the description of tensorial interactions in liquid NMR: application to different software programs.

Long-range orientational restraints derived from alignment or rotational diffusion tensors have greatly contributed to the expansion of applications in biomolecular NMR. The orientation of the principal axis system of these tensors is usually described by the so-called Euler angles. However, no clear consensus has emerged concerning the convention of the associated orthogonal rotations. As a result, the different programs that derive or predict them have adopted different conventions, which make comparison between their results difficult. Moreover, the rotation schemes are seldom completely described. Here, we summarize the different conventions, determine which ones are adopted by commonly used software packages, and establish the formal equivalencies between the different calculated Euler angles.