Deuteration of nonexchangeable protons on proteins affects their thermal stability, side-chain dynamics, and hydrophobicity

Unveiling the relationship between heat-resistant structure characteristics and inhibitory activity in colored highland barley proteinaceous α-amylase inhibitors

Natural α-amylase inhibitors (α-AIs) serve as food processing additives, capable of mitigating postprandial blood glucose levels, but heat resistance limits their application in high-temperature processing. This study delved into the correlation between protein structural characteristics and heat-resistance of colored highland barley (CHB) α-AIs and evaluated the inhibitory activity during chemical modification and in vitro digestion. Results demonstrated that CHB α-AIs were glycoproteins, the inhibitory activity retention rate of black highland barley α-AI salted-out with 0-60 % (NH4)2SO4 (BK1 α-AI) was 56.23 % ± 0.64 %. The protein structure analysis revealed that the preservation of three-dimensional structure was attributed to hydrogen bonds and hydrophobic interactions, and disulfide bonds played a crucial role in maintaining protein folding and activity. Succinylation increased the content of disulfide bonds after heating, and the inhibitory activity retention rate of α-AI noodles increased from 37.72 % ± 2.49 % to 42.79 % ± 0.39 %. These findings provide a theoretical foundation for the application of α-AI in thermally processed foods.

Water slowing down drives the occurrence of the low temperature dynamical transition in microgels

The protein dynamical transition marks an increase in atomic mobility and the onset of anharmonic motions at a critical temperature (T d), which is considered relevant for protein functionality. This phenomenon is ubiquitous, regardless of protein composition, structure and biological function and typically occurs at large protein content, to avoid water crystallization. Recently, a dynamical transition has also been reported in non-biological macromolecules, such as poly(N-isopropyl acrylamide) (PNIPAM) microgels, bearing many similarities to proteins. While the generality of this phenomenon is well-established, the role of water in the transition remains a subject of debate. In this study, we use atomistic molecular dynamics (MD) simulations and elastic incoherent neutron scattering (EINS) experiments with selective deuteration to investigate the microscopic origin of the dynamical transition and distinguish water and PNIPAM roles. While a standard analysis of EINS experiments would suggest that the dynamical transition occurs in PNIPAM and water at a similar temperature, simulations reveal a different perspective, also qualitatively supported by experiments. From room temperature down to about 180 K, PNIPAM exhibits only modest changes of dynamics, while water, being mainly hydration water under the probed extreme confinement, significantly slows down and undergoes a mode-coupling transition from diffusive to activated. Our findings therefore challenge the traditional view of the dynamical transition, demonstrating that it occurs in proximity of the water mode-coupling transition, shedding light on the intricate interplay between polymer and water dynamics.

Investigation of the structure of protein-polymer conjugates in solution reveals the impact of protein deuteration and the size of the polymer on its thermal stability

The conjugation of proteins with polymers offers immense biotechnological potential by creating novel macromolecules. This article presents experimental findings on the structural properties of maltose-binding protein (MBP) conjugated with linear biodegradable polyphosphoester polymers with different molecular weights. We studied isotopic effects on both proteins and polymers. Circular dichroism and fluorescence spectroscopy and small-angle neutron scattering reveal that the conjugation process destabilizes the protein, affecting the secondary more than the tertiary structure, even at room temperature, and that the presence of two domains in the MBP may contribute to its observed instability. Notably, unfolding temperatures differ between native MBP and the conjugates. In particular, this study sheds light on the complex interplay of factors such as the deuteration influencing protein stability and conformational changes in the conjugation processes. The perdeuteration influences the hydrogen bond network and hydrophobic interactions in the case of the MBP protein. The perdeuteration of the protein influences the hydrogen bond network and hydrophobic interactions. This is evident in the decreased thermal stability of deuterated MBP protein, in the conjugate, especially with high-molecular-mass polymers.

Structural characterisation of α-synuclein-membrane interactions and the resulting aggregation using small angle scattering

The presence of amyloid fibrils is a hallmark of several neurodegenerative diseases. Some amyloidogenic proteins, such as α-synuclein and amyloid β, interact with lipids, and this interaction can strongly favour the formation of amyloid fibrils. In particular the primary nucleation step, i.e. the de novo formation of amyloid fibrils, has been shown to be accelerated by lipids. However, the exact mechanism of this acceleration is still mostly unclear. Here we use a range of scattering methods, such as dynamic light scattering (DLS) and small angle X-ray and neutron scattering (SAXS and SANS) to obtain structural information on the binding of α-synuclein to model membranes formed from negatively charged lipids and their co-assembly into amyloid fibrils. We find that the model membranes take an active role in the reaction. The binding of α synuclein to the model membranes immediately induces a major structural change in the lipid assembly, which leads to a break-up into small and mostly disc- or rod-like lipid-protein particles. This transition can be reversed by temperature changes or proteolytic protein removal. Incubation of the small lipid-α-synuclein particles for several hours, however, leads to amyloid fibril formation, whereby the lipids are incorporated into the amyloid fibrils.

Sensitivity-enhanced NMR 15N R1 and R1ρ relaxation experiments for the investigation of intrinsically disordered proteins at high magnetic fields

NMR relaxation experiments provide residue-specific insights into the structural dynamics of proteins. Here, we present an optimized set of sensitivity-enhanced 15N R1 and R1ρ relaxation experiments applicable to fully protonated proteins. The NMR pulse sequences are conceptually similar to the set of TROSY-based sequences and their HSQC counterpart (Lakomek et al., J. Biomol. NMR 2012). Instead of the TROSY read-out scheme, a sensitivity-enhanced HSQC read-out scheme is used, with improved and easier optimized water suppression. The presented pulse sequences are applied on the cytoplasmic domain of the SNARE protein Synpatobrevin-2 (Syb-2), which is intrinsically disordered in its monomeric pre-fusion state. A two-fold increase in the obtained signal-to-noise ratio is observed for this intrinsically disordered protein, therefore offering a four-fold reduction of measurement time compared to the TROSY-detected version. The inter-scan recovery delay can be shortened to two seconds. Pulse sequences were tested at 600 MHz and 1200 MHz 1H Larmor frequency, thus applicable over a wide magnetic field range. A comparison between protonated and deuterated protein samples reveals high agreement, indicating that reliable 15N R1 and R1ρ rate constants can be extracted for fully protonated and deuterated samples. The presented pulse sequences will benefit not only for IDPs but also for an entire range of low and medium-sized proteins.

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.

Structural characterization of an intrinsically disordered protein complex using integrated small-angle neutron scattering and computing

Characterizing structural ensembles of intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins is essential for studying structure-function relationships. Due to the different neutron scattering lengths of hydrogen and deuterium, selective labeling and contrast matching in small-angle neutron scattering (SANS) becomes an effective tool to study dynamic structures of disordered systems. However, experimental timescales typically capture measurements averaged over multiple conformations, leaving complex SANS data for disentanglement. We hereby demonstrate an integrated method to elucidate the structural ensemble of a complex formed by two IDRs. We use data from both full contrast and contrast matching with residue-specific deuterium labeling SANS experiments, microsecond all-atom molecular dynamics (MD) simulations with four molecular mechanics force fields, and an autoencoder-based deep learning (DL) algorithm. From our combined approach, we show that selective deuteration provides additional information that helps characterize structural ensembles. We find that among the four force fields, a99SB-disp and CHARMM36m show the strongest agreement with SANS and NMR experiments. In addition, our DL algorithm not only complements conventional structural analysis methods but also successfully differentiates NMR and MD structures which are indistinguishable on the free energy surface. Lastly, we present an ensemble that describes experimental SANS and NMR data better than MD ensembles generated by one single force field and reveal three clusters of distinct conformations. Our results demonstrate a new integrated approach for characterizing structural ensembles of IDPs.

Fitting Force Field Parameters to NMR Relaxation Data

We present an approach to optimize force field parameters using time-dependent data from NMR relaxation experiments. To do so, we scan parameters in the dihedral angle potential energy terms describing the rotation of the methyl groups in proteins and compare NMR relaxation rates calculated from molecular dynamics simulations with the modified force fields to deuterium relaxation measurements of T4 lysozyme. We find that a small modification of C^γ methyl groups improves the agreement with experiments both for the protein used to optimize the force field and when validating using simulations of CI2 and ubiquitin. We also show that these improvements enable a more effective a posteriori reweighting of the MD trajectories. The resulting force field thus enables more direct comparison between simulations and side-chain NMR relaxation data and makes it possible to construct ensembles that better represent the dynamics of proteins in solution.

Probing Ligand Binding Sites on Large Proteins by Nuclear Magnetic Resonance Spectroscopy of Genetically Encoded Non-Canonical Amino Acids

N6-(((trimethylsilyl)-methoxy)carbonyl)-l-lysine (TMSK) and N6-trifluoroacetyl-l-lysine (TFAK) are non-canonical amino acids, which can be installed in proteins by genetic encoding. In addition, we describe a new aminoacyl-tRNA synthetase specific for N6-(((trimethylsilyl)methyl)-carbamoyl)-l-lysine (TMSNK), which is chemically more stable than TMSK. Using the dimeric SARS-CoV-2 main protease (Mpro) as a model system with three different ligands, we show that the 1H and 19F nuclei of the solvent-exposed trimethylsilyl and CF3 groups produce intense signals in the nuclear magnetic resonance (NMR) spectrum. Their response to active-site ligands differed significantly when positioned near rather than far from the active site. Conversely, the NMR probes failed to confirm the previously reported binding site of the ligand pelitinib, which was found to enhance the activity of Mpro by promoting the formation of the enzymatically active dimer. In summary, the amino acids TMSK, TMSNK, and TFAK open an attractive path for site-specific NMR analysis of ligand binding to large proteins of limited stability and at low concentrations.

Slow methyl axes motions in perdeuterated villin headpiece subdomain probed by cross-correlated NMR relaxation measurements

Protein methyl groups can participate in multiple motional modes on different time scales. Sub-nanosecond to nano-second time scale motions of methyl axes are particularly challenging to detect for small proteins in solutions. In this work we employ NMR relaxation interference between the methyl H-H/H-C dipole-dipole interactions [Sun&Tugarinov, J. Magn. Reason. 2012] to characterize methyl axes motions as a function of temperature in a small model protein villin headpiece subdomain (HP36), in which all non-exchangeable protons are deuterated with the exception of methyl groups of leucine and valine residues. The data points to the existence of slow motional modes of methyl axes on sub-nanosecond to nanosecond time scales. Further, at high temperatures for which the overall tumbling of the protein is on the order of 2 ns, we observe a coupling between the slow internal motion and the overall molecular tumbling, based on the anomalous order parameters and their temperature-dependent trends. The addition of 28%(w/w) glycerol-d8 increases the viscosity of the solvent and separates the timescales of internal and overall tumbling, thus permitting for another view of the necessity of the coupling assumption for these sites at high temperatures.

Deuteration for High-Resolution Detection of Protons in Protein Magic Angle Spinning (MAS) Solid-State NMR

Proton detection developed in the last 20 years as the method of choice to study biomolecules in the solid state. In perdeuterated proteins, proton dipolar interactions are strongly attenuated, which allows yielding of high-resolution proton spectra. Perdeuteration and backsubstitution of exchangeable protons is essential if samples are rotated with MAS rotation frequencies below 60 kHz. Protonated samples can be investigated directly without spin dilution using proton detection methods in case the MAS frequency exceeds 110 kHz. This review summarizes labeling strategies and the spectroscopic methods to perform experiments that yield assignments, quantitative information on structure, and dynamics using perdeuterated samples. Techniques for solvent suppression, H/D exchange, and deuterium spectroscopy are discussed. Finally, experimental and theoretical results that allow estimation of the sensitivity of proton detected experiments as a function of the MAS frequency and the external B₀ field in a perdeuterated environment are compiled.

The impact of folding modes and deuteration on the atomic resolution structure of hen egg-white lysozyme

The biological function of a protein is intimately related to its structure and dynamics, which in turn are determined by the way in which it has been folded. In vitro refolding is commonly used for the recovery of recombinant proteins that are expressed in the form of inclusion bodies and is of central interest in terms of the folding pathways that occur in vivo. Here, biophysical data are reported for in vitro-refolded hydrogenated hen egg-white lysozyme, in combination with atomic resolution X-ray diffraction analyses, which allowed detailed comparisons with native hydrogenated and refolded perdeuterated lysozyme. Distinct folding modes are observed for the hydrogenated and perdeuterated refolded variants, which are determined by conformational changes to the backbone structure of the Lys97-Gly104 flexible loop. Surprisingly, the structure of the refolded perdeuterated protein is closer to that of native lysozyme than that of the refolded hydrogenated protein. These structural differences suggest that the observed decreases in thermal stability and enzymatic activity in the refolded perdeuterated and hydrogenated proteins are consequences of the macromolecular deuteration effect and of distinct folding dynamics, respectively. These results are discussed in the context of both in vitro and in vivo folding, as well as of lysozyme amyloidogenesis.

Inspecting the structural characteristics of chiral drug penicillamine under different pH conditions using Raman optical activity spectroscopy and DFT calculations

The investigation of the structural characteristics of chiral drugs in physiological environments is a challenging research topic, which may lead to a better understanding of how the drugs work. Raman optical activity (ROA) spectroscopy in combination with density functional theory (DFT) calculations was exploited to inspect the structural changes in penicillamine under different acid-base states in aqueous solutions. The B3LYP/aug-cc-PVDZ method was employed and the implicit solvation model density (SMD) was considered for describing the solvation effect in H₂O. The conformations of penicillamine varied with pH, but penicillamine was liable to stabilize in the form of the P_C conformation (the sulfur atom is in a trans orientation with respect to carboxylate) in most cases for both D- and L-isomers. The relationship between the conformations of penicillamine and the ROA peaks, as well as peak assignments, were comprehensively studied and elucidated. In the fingerprint region, two ROA couplets and one ROA triplet with different patterns were recognized. The intensity, sign and frequency of the corresponding peaks also changed with varying pH. Deuteration was carried out to identify the vibrational modes, and the ROA peaks of the deuterated amino group in particular are sensitive to change in the ambient environment. The results are expected not only to serve as a reference for the interpretation of the ROA spectra of penicillamine and other chiral drugs with analogous structures but also to evaluate the structural changes of chiral molecules in physiological environments, which will form the basis of further exploration of the effects of structural characteristics on the pharmacological and toxicological properties of chiral drugs.

Fitting Side-Chain NMR Relaxation Data Using Molecular Simulations

Proteins display a wealth of dynamical motions that can be probed using both experiments and simulations. We present an approach to integrate side-chain NMR relaxation measurements with molecular dynamics simulations to study the structure and dynamics of these motions. The approach, which we term ABSURDer (average block selection using relaxation data with entropy restraints), can be used to find a set of trajectories that are in agreement with relaxation measurements. We apply the method to deuterium relaxation measurements in T4 lysozyme and show how it can be used to integrate the accuracy of the NMR measurements with the molecular models of protein dynamics afforded by the simulations. We show how fitting of dynamic quantities leads to improved agreement with static properties and highlight areas needed for further improvements of the approach.

Structural insights into protein folding, stability and activity using in vivo perdeuteration of hen egg-white lysozyme

This structural and biophysical study exploited a method of perdeuterating hen egg-white lysozyme based on the expression of insoluble protein in Escherichia coli followed by in-column chemical refolding. This allowed detailed comparisons with perdeuterated lysozyme produced in the yeast Pichia pastoris, as well as with unlabelled lysozyme. Both perdeuterated variants exhibit reduced thermal stability and enzymatic activity in comparison with hydrogenated lysozyme. The thermal stability of refolded perdeuterated lysozyme is 4.9°C lower than that of the perdeuterated variant expressed and secreted in yeast and 6.8°C lower than that of the hydrogenated Gallus gallus protein. However, both perdeuterated variants exhibit a comparable activity. Atomic resolution X-ray crystallographic analyses show that the differences in thermal stability and enzymatic function are correlated with refolding and deuteration effects. The hydrogen/deuterium isotope effect causes a decrease in the stability and activity of the perdeuterated analogues; this is believed to occur through a combination of changes to hydrophobicity and protein dynamics. The lower level of thermal stability of the refolded perdeuterated lysozyme is caused by the unrestrained Asn103 peptide-plane flip during the unfolded state, leading to a significant increase in disorder of the Lys97-Gly104 region following subsequent refolding. An ancillary outcome of this study has been the development of an efficient and financially viable protocol that allows stable and active perdeuterated lysozyme to be more easily available for scientific applications.

Deuteration of nonexchangeable protons on proteins affects their thermal stability, side-chain dynamics, and hydrophobicity

We have investigated the effect of deuteration of non-exchangeable protons on protein global thermal stability, hydrophobicity, and local flexibility using well-known thermostable model systems such as the villin headpiece subdomain (HP36) and the third immunoglobulin G-binding domain of protein G (GB3). Reversed-phase high-performance liquid chromatography (RP-HPLC) measurements as a function of temperature probe global thermal stability in the presence of acetonitrile, while differential scanning calorimetry determines thermal stability in solution. Both indicate small but measurable changes in the order of several degrees. RP-HPLC also permitted quantification of the effect of deuteration of just three core phenylalanine side chains of HP36. NMR dynamics investigation has focused on methyl axes motions using cross-correlated relaxation measurements. The analysis of order parameters provided a complex picture indicating that deuteration generally increases motional amplitudes of sub-nanosecond motion in GB3 but decreases those in HP36. Combined with earlier dynamics measurements at C_α -C_β sites and backbone sites of GB3, which probed slower time scales, the results point to the need to probe multiple atoms in the protein and variety of time scales to the discern the full complexity of the effects of deuteration on dynamics.