Intra-residue methyl-methyl correlations for valine and leucine residues in large proteins from a 3D-HMBC-HMQC experiment

Solution-state methyl NMR spectroscopy of large non-deuterated proteins enabled by deep neural networks

Methyl-TROSY nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for characterising large biomolecules in solution. However, preparing samples for these experiments is demanding and entails deuteration, limiting its use. Here we demonstrate that NMR spectra recorded on protonated, uniformly 13C labelled samples can be processed using deep neural networks to yield spectra that are of similar quality to typical deuterated methyl-TROSY spectra, potentially providing information for proteins that cannot be produced in bacterial systems. We validate the methodology experimentally on three proteins with molecular weights in the range 42-360 kDa. We further demonstrate the applicability of our methodology to 3D NOESY spectra of Escherichia coli Malate Synthase G (81 kDa), where observed NOE cross-peaks are in good agreement with the available structure. The method represents an advance in the field of using deep learning to analyse complex magnetic resonance data and could have an impact on the study of large biomolecules in years to come.

Structural basis of substrate recognition and allosteric activation of the proapoptotic mitochondrial HtrA2 protease

The mitochondrial serine protease HtrA2 is a human homolog of the Escherichia coli Deg-proteins exhibiting chaperone and proteolytic roles. HtrA2 is involved in both apoptotic regulation via its ability to degrade inhibitor-of-apoptosis proteins (IAPs), as well as in cellular maintenance as part of the cellular protein quality control machinery, by preventing the possible toxic accumulation of aggregated proteins. In this study, we use advanced solution NMR spectroscopy methods combined with biophysical characterization and biochemical assays to elucidate the crucial role of the substrate recognizing PDZ domain. This domain regulates the protease activity of HtrA2 by triggering an intricate allosteric network involving the regulatory loops of the protease domain. We further show that divalent metal ions can both positively and negatively modulate the activity of HtrA2, leading to a refined model of HtrA2 regulation within the apoptotic pathway.

Activation of caspase-9 on the apoptosome as studied by methyl-TROSY NMR

Mitochondrial apoptotic signaling cascades lead to the formation of the apoptosome, a 1.1-MDa heptameric protein scaffold that recruits and activates the caspase-9 protease. Once activated, caspase-9 cleaves and activates downstream effector caspases, triggering the onset of cell death through caspase-mediated proteolysis of cellular proteins. Failure to activate caspase-9 enables the evasion of programmed cell death, which occurs in various forms of cancer. Despite the critical apoptotic function of caspase-9, the structural mechanism by which it is activated on the apoptosome has remained elusive. Here, we used a combination of methyl-transverse relaxation-optimized NMR spectroscopy, protein engineering, and biochemical assays to study the activation of caspase-9 bound to the apoptosome. In the absence of peptide substrate, we observed that both caspase-9 and its isolated protease domain (PD) only very weakly dimerize with dissociation constants in the millimolar range. Methyl-NMR spectra of isotope-labeled caspase-9, within the 1.3-MDa native apoptosome complex or an engineered 480-kDa apoptosome mimic, reveal that the caspase-9 PD remains monomeric after recruitment to the scaffold. Binding to the apoptosome, therefore, organizes caspase-9 PDs so that they can rapidly and extensively dimerize only when substrate is present, providing an important layer in the regulation of caspase-9 activation. Our work highlights the unique role of NMR spectroscopy to structurally characterize protein domains that are flexibly tethered to large scaffolds, even in cases where the molecular targets are in excess of 1 MDa, as in the present example.

Intrinsic structural dynamics dictate enzymatic activity and inhibition

Enzymes are known to sample various conformations, many of which are critical for their biological function. However, structural characterizations of enzymes predominantly focus on the most populated conformation. As a result, single-point mutations often produce structures that are similar or essentially identical to those of the wild-type enzyme despite large changes in enzymatic activity. Here, we show for mutants of a histone deacetylase enzyme (HDAC8) that reduced enzymatic activities, reduced inhibitor affinities, and reduced residence times are all captured by the rate constants between intrinsically sampled conformations that, in turn, can be obtained independently by solution NMR spectroscopy. Thus, for the HDAC8 enzyme, the dynamic sampling of conformations dictates both enzymatic activity and inhibitor potency. Our analysis also dissects the functional role of the conformations sampled, where specific conformations distinct from those in available structures are responsible for substrate and inhibitor binding, catalysis, and product dissociation. Precise structures alone often do not adequately explain the effect of missense mutations on enzymatic activity and drug potency. Our findings not only assign functional roles to several conformational states of HDAC8 but they also underscore the paramount role of dynamics, which will have general implications for characterizing missense mutations and designing inhibitors.

Structure and dynamics of the mitochondrial DNA-compaction factor Abf2 from S. cerevisiae

Mitochondria are essential organelles that produce most of the energy via the oxidative phosphorylation (OXPHOS) system in all eukaryotic cells. Several essential subunits of the OXPHOS system are encoded by the mitochondrial genome (mtDNA) despite its small size. Defects in mtDNA maintenance and expression can lead to severe OXPHOS deficiencies and are amongst the most common causes of mitochondrial disease. The mtDNA is packaged as nucleoprotein structures, referred to as nucleoids. The conserved mitochondrial proteins, ARS-binding factor 2 (Abf2) in the budding yeast Saccharomyces cerevisiae (S. cerevisiae) and mitochondrial transcription factor A (TFAM) in mammals, are nucleoid associated proteins (NAPs) acting as condensing factors needed for packaging and maintenance of the mtDNA. Interestingly, gene knockout of Abf2 leads, in yeast, to the loss of mtDNA and respiratory growth, providing evidence for its crucial role. On a structural level, the condensing factors usually contain two DNA binding domains called high-mobility group boxes (HMG boxes). The co-operating mechanical activities of these domains are crucial in establishing the nucleoid architecture by bending the DNA strands. Here we used advanced solution NMR spectroscopy methods to characterize the dynamical properties of Abf2 together with its interaction with DNA. We find that the two HMG-domains react notably different to external cues like temperature and salt, indicating distinct functional properties. Biophysical characterizations show the pronounced difference of these domains upon DNA-binding, suggesting a refined picture of the Abf2 functional cycle.

Disulfide-Bond-Induced Structural Frustration and Dynamic Disorder in a Peroxiredoxin from MAS NMR

Disulfide bond formation is fundamentally important for protein structure and constitutes a key mechanism by which cells regulate the intracellular oxidation state. Peroxiredoxins (PRDXs) eliminate reactive oxygen species such as hydrogen peroxide through a catalytic cycle of Cys oxidation and reduction. Additionally, upon Cys oxidation PRDXs undergo extensive conformational rearrangements that may underlie their presently structurally poorly defined functions as molecular chaperones. Rearrangements include high molecular-weight oligomerization, the dynamics of which are, however, poorly understood, as is the impact of disulfide bond formation on these properties. Here we show that formation of disulfide bonds along the catalytic cycle induces extensive μs time scale dynamics, as monitored by magic-angle spinning NMR of the 216 kDa-large Tsa1 decameric assembly and solution-NMR of a designed dimeric mutant. We ascribe the conformational dynamics to structural frustration, resulting from conflicts between the disulfide-constrained reduction of mobility and the desire to fulfill other favorable contacts.

Progress toward automated methyl assignments for methyl-TROSY applications

Methyl-TROSY spectroscopy has extended the reach of solution-state NMR to supra-molecular machineries over 100 kDa in size. Methyl groups are ideal probes for studying structure, dynamics, and protein-protein interactions in quasi-physiological conditions with atomic resolution. Successful implementation of the methodology requires accurate methyl chemical shift assignment, and the task still poses a significant challenge in the field. In this work, we outline the current state of technology for methyl labeling, data collection, data analysis, and nuclear Overhauser effect (NOE)-based automated methyl assignment approaches. We present MAGIC-Act and MAGIC-View, two Python extensions developed as part of the popular NMRFAM-Sparky package, and MAGIC-Net a standalone structure-based network analysis program. MAGIC-Act conducts statistically driven amino acid typing, Leu/Val pairing guided by 3D HMBC-HMQC, and NOESY cross-peak symmetry checking. MAGIC-Net provides model-based NOE statistics to aid in selection of a methyl labeling scheme. The programs provide a versatile, semi-automated framework for rapid methyl assignment.

Nearest-neighbor NMR spectroscopy: categorizing spectral peaks by their adjacent nuclei

Methyl-NMR enables atomic-resolution studies of structure and dynamics of large proteins in solution. However, resonance assignment remains challenging. The problem is to combine existing structural informational with sparse distance restraints and search for the most compatible assignment among the permutations. Prior classification of peaks as either from isoleucine, leucine, or valine reduces the search space by many orders of magnitude. However, this is hindered by overlapped leucine and valine frequencies. In contrast, the nearest-neighbor nuclei, coupled to the methyl carbons, resonate in distinct frequency bands. Here, we develop a framework to imprint additional information about passively coupled resonances onto the observed peaks. This depends on simultaneously orchestrating closely spaced bands of resonances along different magnetization trajectories, using principles from control theory. For methyl-NMR, the method is implemented as a modification to the standard fingerprint spectrum (the 2D-HMQC). The amino acid type is immediately apparent in the fingerprint spectrum. There is no additional relaxation loss or an increase in experimental time. The method is validated on biologically relevant proteins. The idea of generating new spectral information using passive, adjacent resonances is applicable to other contexts in NMR spectroscopy.

Protein-ligand structure determination with the NMR molecular replacement tool, NMR2

We recently reported on a new method called NMR Molecular Replacement that efficiently derives the structure of a protein-ligand complex at the interaction site. The method was successfully applied to high and low affinity complexes covering ligands from peptides to small molecules. The algorithm used in the NMR Molecular Replacement program has until now not been described in detail. Here, we present a complete description of the NMR Molecular Replacement implementation as well as several new features that further reduce the time required for structure elucidation.

Automated assignment of methyl NMR spectra from large proteins

As structural biology trends towards larger and more complex biomolecular targets, a detailed understanding of their interactions and underlying structures and dynamics is required. The development of methyl-TROSY has enabled NMR spectroscopy to provide atomic-resolution insight into the mechanisms of large molecular assemblies in solution. However, the applicability of methyl-TROSY has been hindered by the laborious and time-consuming resonance assignment process, typically performed with domain fragmentation, site-directed mutagenesis, and analysis of NOE data in the context of a crystal structure. In response, several structure-based automatic methyl assignment strategies have been developed over the past decade. Here, we present a comprehensive analysis of all available methods and compare their input data requirements, algorithmic strategies, and reported performance. In general, the methods fall into two categories: those that primarily rely on inter-methyl NOEs, and those that utilize methyl PRE- and PCS-based restraints. We discuss their advantages and limitations, and highlight the potential benefits from standardizing and combining different methods.