Disentangling Glycan-Protein Interactions: Nuclear Magnetic Resonance (NMR) to the Rescue.

The interactions of glycans with proteins modulate many events related to health and disease. In fact, the establishment of these recognition events and their biological consequences are intimately related to the three-dimensional structures of both partners, as well as to their dynamic features and their presentation on the corresponding cell compartments. NMR techniques are unique to disentangle these characteristics and, indeed, diverse NMR-based methodologies have been developed and applied to monitor the binding events of glycans with their associate receptors. This protocol outlines the procedures to acquire, process and analyze two of the most powerful NMR methodologies employed in the NMR-glycobiology field, 1H-Saturation transfer difference (STD) and 1H,15N-Heteronuclear single quantum coherence (HSQC) titration experiments, which complementarily offer information from the glycan and protein perspective, respectively. Indeed, when combined they offer a powerful toolkit for elucidating both the structural and dynamic aspects of molecular recognition processes. This comprehensive approach enhances our understanding of glycan-protein interactions and contributes to advancing research in the chemical glycobiology field.

Analysis of the homodimeric structure of a D-Ala-D-Ala metallopeptidase, VanX, from vancomycin-resistant bacteria.

Bacteria that have acquired resistance to most antibiotics, particularly those causing nosocomial infections, create serious problems. Among these, the emergence of vancomycin-resistant enterococci was a tremendous shock, considering that vancomycin is the last resort for controlling methicillin-resistant Staphylococcus aureus. Therefore, there is an urgent need to develop an inhibitor of VanX, a protein involved in vancomycin resistance. Although the crystal structure of VanX has been resolved, its asymmetric unit contains six molecules aligned in a row. We have developed a structural model of VanX as a stable dimer in solution, primarily utilizing nuclear magnetic resonance (NMR) residual dipolar coupling. Despite the 46 kDa molecular mass of the dimer, the analyses, which are typically not as straightforward as those of small proteins around 10 kDa, were successfully conducted. We assigned the main chain using an amino acid-selective unlabeling method. Because we found that the zinc ion-coordinating active sites in the dimer structure were situated in the opposite direction to the dimer interface, we generated an active monomer by replacing an amino acid at the dimer interface. The monomer consists of only 202 amino acids and is expected to be used in future studies to screen and improve inhibitors using NMR.

Fast collective motions of backbone in transmembrane α helices are critical to water transfer of aquaporin.

Fast collective motions are widely present in biomolecules, but their functional relevance remains unclear. Herein, we reveal that fast collective motions of backbone are critical to the water transfer of aquaporin Z (AqpZ) by using solid-state nuclear magnetic resonance (ssNMR) spectroscopy and molecular dynamics (MD) simulations. A total of 212 residue site-specific dipolar order parameters and 158 15N spin relaxation rates of the backbone are measured by combining the 13C- and 1H-detected multidimensional ssNMR spectra. Analysis of these experimental data by theoretic models suggests that the small-amplitude (~10°) collective motions of the transmembrane α helices on the nanosecond-to-microsecond timescales are dominant for the dynamics of AqpZ. The MD simulations demonstrate that these collective motions are critical to the water transfer efficiency of AqpZ by facilitating the opening of the channel and accelerating the water-residue hydrogen bonds renewing in the selectivity filter region.

Studying protein stability in crowded environments by NMR.

Most proteins perform their functions in crowded and complex cellular environments where weak interactions are ubiquitous between biomolecules. These complex environments can modulate the protein folding energy landscape and hence affect protein stability. NMR is a nondestructive and effective method to quantify the kinetics and equilibrium thermodynamic stability of proteins at an atomic level within crowded environments and living cells. Here, we review NMR methods that can be used to measure protein stability, as well as findings of studies on protein stability in crowded environments mimicked by polymer and protein crowders and in living cells. The important effects of chemical interactions on protein stability are highlighted and compared to spatial excluded volume effects.

Probing the energy barriers and stages of membrane protein unfolding using solid-state NMR spectroscopy.

Understanding how the amino acid sequence dictates protein structure and defines its stability is a fundamental problem in molecular biology. It is especially challenging for membrane proteins that reside in the complex environment of a lipid bilayer. Here, we obtain an atomic-level picture of the thermally induced unfolding of a membrane-embedded α-helical protein, human aquaporin 1, using solid-state nuclear magnetic resonance spectroscopy. Our data reveal the hierarchical two-step pathway that begins with unfolding of a structured extracellular loop and proceeds to an intermediate state with a native-like helical packing. In the second step, the transmembrane domain unravels as a single unit, resulting in a heterogeneous misfolded state with high helical content but with nonnative helical packing. Our results show the importance of loops for the kinetic stabilization of the whole membrane protein structure and support the three-stage membrane protein folding model.

A salt bridge of the C-terminal carboxyl group regulates PHPT1 substrate affinity and catalytic activity.

PHPT1 is a histidine phosphatase that modulates signaling in eukaryotes through its catalytic activity. Here, we present an analysis of the structure and dynamics of PHPT1 through a combination of solution NMR, molecular dynamics, and biochemical experiments. We identify a salt bridge formed between the R78 guanidinium moiety and the C-terminal carboxyl group on Y125 that is critical for ligand binding. Disruption of the salt bridge by appending a glycine residue at the C-terminus (G126) leads to a decrease in catalytic activity and binding affinity for the pseudo substrate, para-nitrophenylphosphate (pNPP), as well as the active site inhibitor, phenylphosphonic acid (PPA). We show through NMR chemical shift, 15N relaxation measurements, and analysis of molecular dynamics trajectories, that removal of this salt bridge results in an active site that is altered both structurally and dynamically thereby significantly impacting enzymatic function and confirming the importance of this electrostatic interaction.

Decoding the mechanism governing the structural stability of wheat germ agglutinin and its isolated domains: A combined calorimetric, NMR, and MD simulation study.

Wheat germ agglutinin (WGA) demonstrates potential as an oral delivery agent owing to its selective binding to carbohydrates and its capacity to traverse biological membranes. In this study, we employed differential scanning calorimetry and molecular dynamics simulations to comprehensively characterize the thermal unfolding process of both the complete lectin and its four isolated domains. Furthermore, we present the nuclear magnetic resonance structures of three domains that were previously lacking experimental structures in their isolated forms. Our results provide a collective understanding of the energetic and structural factors governing the intricate unfolding mechanism of the complete agglutinin, shedding light on the specific role played by each domain in this process. The analysis revealed negligible interdomain cooperativity, highlighting instead significant coupling between dimer dissociation and the unfolding of the more labile domains. By comparing the dominant interactions, we rationalized the stability differences among the domains. Understanding the structural stability of WGA opens avenues for enhanced drug delivery strategies, underscoring its potential as a promising carrier throughout the gastrointestinal environment.

Insights on the G protein-coupled receptor helix 8 solution structure and orientation using a neurotensin receptor 1 peptide.

G-protein coupled receptors (GPCRs) are the largest class of membrane proteins encoded in the human genome with high pharmaceutical relevance and implications to human health. These receptors share a prevalent architecture of seven transmembrane helices followed by an intracellular, amphipathic helix 8 (H8) and a disordered C-terminal tail (Ctail). Technological advancements have led to over 1000 receptor structures in the last two decades, yet frequently H8 and the Ctail are conformationally heterogeneous or altogether absent. Here we synthesize a peptide comprising the neurotensin receptor 1 (NTS1) H8 and Ctail (H8-Ctail) to investigate its structural stability, conformational dynamics, and orientation in the presence of detergent and phospholipid micelles, which mimic the membrane. Circular dichroism (CD) and nuclear magnetic resonance (NMR) measurements confirm that zwitterionic 1,2-diheptanoyl-sn-glycero-3-phosphocholine is a potent stabilizer of H8 structure, whereas the commonly-used branched detergent lauryl maltose neopentyl glycol (LMNG) is unable to completely stabilize the helix - even at amounts four orders of magnitude greater than its critical micellar concentration. We then used NMR spectroscopy to assign the backbone chemical shifts. A series of temperature and lipid titrations were used to define the H8 boundaries as F376-R392 from chemical shift perturbations, changes in resonance intensity, and chemical-shift-derived phi/psi angles. Finally, the H8 azimuthal and tilt angles, defining the helix orientation relative of the membrane normal were measured using paramagnetic relaxation enhancement NMR. Taken together, our studies reveal the H8-Ctail region is sensitive to membrane physicochemical properties and is capable of more adaptive behavior than previously suggested by static structural techniques.

Heat-induced structural and chemical changes to a computationally designed miniprotein.

The de novo design of miniprotein inhibitors has recently emerged as a new technology to create proteins that bind with high affinity to specific therapeutic targets. Their size, ease of expression, and apparent high stability makes them excellent candidates for a new class of protein drugs. However, beyond circular dichroism melts and hydrogen/deuterium exchange experiments, little is known about their dynamics, especially at the elevated temperatures they seemingly tolerate quite well. To address that and gain insight for future designs, we have focused on identifying unintended and previously overlooked heat-induced structural and chemical changes in a particularly stable model miniprotein, EHEE_rd2_0005. Nuclear magnetic resonance (NMR) studies suggest the presence of dynamics on multiple time and temperature scales. Transiently elevating the temperature results in spontaneous chemical deamidation visible in the NMR spectra, which we validate using both capillary electrophoresis and mass spectrometry (MS) experiments. High temperatures also result in greatly accelerated intrinsic rates of hydrogen exchange and signal loss in NMR heteronuclear single quantum coherence spectra from local unfolding. These losses are in excellent agreement with both room temperature hydrogen exchange experiments and hydrogen bond disruption in replica exchange molecular dynamics simulations. Our analysis reveals important principles for future miniprotein designs and the potential for high stability to result in long-lived alternate conformational states.

Solid-state NMR spectroscopic analysis for structure determination of a zinc-bound catalytic amyloid fibril.

Zinc ions are commonly involved in enzyme catalysis and protein structure stabilization, but their coordination geometry of zinc-protein complex is rarely determined. Here, in this chapter, we introduce a systematic solid-state NMR approach to determine the oligomeric assembly and Zn2+ coordination geometry of a de novo designed amyloid fibrils that catalyze zinc dependent ester hydrolysis. NMR chemical shifts and intermolecular contacts confirm that the peptide forms parallel-in-register ß-sheets, with the two forms of Zn2+ bound histidines in each peptide. The amphiphilic parallel ß-sheets assemble into stacked bilayers that are stabilized by hydrophobic side chains between ß-sheets. The conformations of the histidine side chains, determined by 13C-15N distance measurements, reveal how histidines protrude from the ß-sheet. 1H-15N correlation spectra show that the single-Zn2+ coordinated histidine associated with dynamic water. The resulting structure provides insight into how metal ions contribute to stabilizing the protein structure and driving its catalytic reactivity.

Boron catalysis in a designer enzyme.

Enzymes play an increasingly important role in improving the benignity and efficiency of chemical production, yet the diversity of their applications lags heavily behind chemical catalysts as a result of the relatively narrow range of reaction mechanisms of enzymes. The creation of enzymes containing non-biological functionalities facilitates reaction mechanisms outside nature's canon and paves the way towards fully programmable biocatalysis1-3. Here we present a completely genetically encoded boronic-acid-containing designer enzyme with organocatalytic reactivity not achievable with natural or engineered biocatalysts4,5. This boron enzyme catalyses the kinetic resolution of hydroxyketones by oxime formation, in which crucial interactions with the protein scaffold assist in the catalysis. A directed evolution campaign led to a variant with natural-enzyme-like enantioselectivities for several different substrates. The unique activation mode of the boron enzyme was confirmed using X-ray crystallography, high-resolution mass spectrometry (HRMS) and 11B NMR spectroscopy. Our study demonstrates that genetic-code expansion can be used to create evolvable enantioselective enzymes that rely on xenobiotic catalytic moieties such as boronic acids and access reaction mechanisms not reachable through catalytic promiscuity of natural or engineered enzymes.

1Hß chemical shift-based phase modulated NMR methods for fast identification of amino acid types in proteins.

Here we report two phase modulated NMR experiments: PM-2D HN(CACBHB) and PM-2D HN(HB), that use 1Hß chemical shifts to rapidly identify amino acid type in proteins. The magnetization on the 1Hß spins during the experiments is allowed to evolve for a fixed evolution period that results in phase modulation (positive or negative) of the cross peaks corresponding to various amino acid residues on their 2D HN projections, resembling a typical 2D [1H-15N]-HSQC spectrum. All amino acids except glycine can be categorized into three discernible groups based on their 1Hß chemical shifts, resulting in unique phase patterns at different fixed evolution periods for 1Hß, thus facilitating their identification. Remarkably, the PM-2D HN(HB) stands out among all amino acid type identification NMR techniques for its applicability with cost-effective and most routinely employed 15N-labeled protein samples for NMR studies. Furthermore, when combined effectively with the 13Cß chemical shift-based phase modulated NMR method (PM-2D HN(CACB)), these methods resolved the identification of large groups of amino acids into relatively smaller groups. Moreover, these techniques can accelerate the sequence-specific sequential resonance assignment (SSRA) process and would help in fast tracking of assigned NMR signals exhibiting chemical shift perturbation on the 2D [1H-15N]-HSQC spectrum of proteins during various experiments (e.g., temperature change, pH change, and protein or ligand or cofactor binding) as well as in site-directed mutagenesis.

On the Validation of Protein Force Fields Based on Structural Criteria.

Molecular dynamics simulations depend critically on the quality of the force field used to describe the interatomic interactions and the extent to which it has been validated for use in a specific application. Using a curated test set of 52 high-resolution structures, 39 derived from X-ray diffraction and 13 solved using NMR, we consider the extent to which different parameter sets of the GROMOS protein force field can be distinguished based on comparing a range of structural criteria, including the number of backbone hydrogen bonds, the number of native hydrogen bonds, polar and nonpolar solvent-accessible surface area, radius of gyration, the prevalence of secondary structure elements, J-coupling constants, nuclear Overhauser effect (NOE) intensities, positional root-mean-square deviations (RMSD), and the distribution of backbone ϕ and ψ dihedral angles. It is shown that while statistically significant differences between the average values of individual metrics could be detected, these were in general small. Furthermore, improvements in agreement in one metric were often offset by loss of agreement in another. The work establishes a framework and test set against which protein force fields can be validated. It also highlights the danger of inferring the relative quality of a given force field based on a small range of structural properties or small number of proteins.

High-resolution heteronuclear correlations between spin-1/2 and half-integer quadrupolar nuclei under fast MAS solid-state NMR.

High isotropic resolution is essential for the structural elucidation of samples with multiple sites. In this study, utilizing the benefits of TRAPDOR-based heteronuclear multiple quantum coherence (T-HMQC) and a pair of one rotor period long cosine amplitude modulated low-power (cos-lp) pulse-based symmetric-split-t1 multiple-quantum magic angle spinning (MQMAS) methods, we have developed a proton-detected 2D 35Cl/1H T-HMQC-MQMAS pulse sequence under fast MAS (70 kHz) to achieve high-resolution in the indirect dimension of the spin-3/2 (35Cl) nuclei connected via protons. As T-HMQC polarizes not only single-quantum central transition (SQCT) but also triple-quantum (TQ) coherences, the proposed 2D pulse sequence is implemented via selection of two coherence pathways (SQCT→TQ →SQCT and TQ → SQCT→TQ) resulting in the 35Cl isotropic dimension and is superior to the existing double-quantum satellite-transition (DQST) T-HMQC in terms of resolution.

Improved analysis of NMR chemical shift perturbations through an error estimation method.

In solution NMR, chemical shift perturbation (CSP) experiments are widely employed to study intermolecular interactions. However, excluding the nonsignificant peak shift is difficult because little is known about errors in CSP. Here, to address this issue, we introduce a method for estimating errors in CSP based on the noise level. First, we developed a technique that involves line shape fitting to estimate errors in peak position via Monte Carlo simulations. Second, this technique was applied to estimate errors in CSP. In intermolecular interaction analysis of VAP-A with SNX2, error estimation of CSP enabled the evaluation of small but significant changes in peak position and yielded detailed insights that are unattainable with conventional CSP analysis. Third, this technique was successfully applied to estimate errors in residual dipolar couplings. In conclusion, our error estimation method improves CSP analysis by excluding the nonsignificant peak shift.

1H, 13C, and 15N resonance assignments of the La Motif of the human La-related protein 1.

Human La-related protein 1 (HsLARP1) is involved in post-transcriptional regulation of certain 5' terminal oligopyrimidine (5'TOP) mRNAs as well as other mRNAs and binds to both the 5'TOP motif and the 3'-poly(A) tail of certain mRNAs. HsLARP1 is heavily involved in cell proliferation, cell cycle defects, and cancer, where HsLARP1 is significantly upregulated in malignant cells and tissues. Like all LARPs, HsLARP1 contains a folded RNA binding domain, the La motif (LaM). Our current understanding of post-transcriptional regulation that emanates from the intricate molecular framework of HsLARP1 is currently limited to small snapshots, obfuscating our understanding of the full picture on HsLARP1 functionality in post-transcriptional events. Here, we present the nearly complete resonance assignment of the LaM of HsLARP1, providing a significant platform for future NMR spectroscopic studies.

A Hybrid Nanotube Stamp System in Intracellular Protein Delivery for Cancer Treatment and NMR Analytical Techniques.

In contrast to intracellular gene transfer, the direct delivery of expressed proteins is a significantly challenging yet essential technique for elucidating cellular functions, including protein complex structure, liquid-liquid phase separation, therapeutic applications, and reprogramming. In this study, we developed a hybrid nanotube (HyNT) stamp system that physically inserts the HyNTs into adhesive cells, enabling the injection of target molecules through HyNT ducts. This system demonstrates the capability to deliver multiple proteins, such as lactate oxidase (LOx) and ubiquitin (UQ), to approximately 1.8 × 107 adhesive cells with a delivery efficiency of 89.9% and a viability of 97.1%. The delivery of LOx enzyme into HeLa cancer cells induced cell death, while enzyme-delivered healthy cells remained viable. Furthermore, our stamp system can deliver an isotope-labeled UQ into adhesive cells for detection by nuclear magnetic resonance (NMR).

Structural characterization of E22G Aß1-42 fibrils via1H detected MAS NMR.

Amyloid fibrils have been implicated in the pathogenesis of several neurodegenerative diseases, the most prevalent example being Alzheimer's disease (AD). Despite the prevalence of AD, relatively little is known about the structure of the associated amyloid fibrils. This has motivated our studies of fibril structures, extended here to the familial Arctic mutant of Aß1-42, E22G-Aß1-42. We found E22G-AßM0,1-42 is toxic to Escherichia coli, thus we expressed E22G-Aß1-42 fused to the self-cleavable tag NPro in the form of its EDDIE mutant. Since the high surface activity of E22G-Aß1-42 makes it difficult to obtain more than sparse quantities of fibrils, we employed 1H detected magic angle spinning (MAS) nuclear magnetic resonance (NMR) experiments to characterize the protein. The 1H detected 13C-13C methods were first validated by application to fully protonated amyloidogenic nanocrystals of GNNQQNY, and then applied to fibrils of the Arctic mutant of Aß, E22G-Aß1-42. The MAS NMR spectra indicate that the biosynthetic samples of E22G-Aß1-42 fibrils comprise a single conformation with 13C chemical shifts extracted from hCH, hNH, and hCCH spectra that are very similar to those of wild type Aß1-42 fibrils. These results suggest that E22G-Aß1-42 fibrils have a structure similar to that of wild type Aß1-42.

One touch is all it takes: the supramolecular interaction between ubiquitin and lanthanide complexes revisited by paramagnetic NMR and molecular dynamics.

The supramolecular interaction between lanthanide complexes and proteins is at the heart of numerous chemical and biological studies. Some of these complexes have demonstrated remarkable interaction properties with proteins or peptides in solution and in the crystalline state. Here we have used the paramagnetism of lanthanide ions to characterize the affinity of two lanthanide complexes for ubiquitin. As the interaction process is dynamic, the acquired NMR data only reflect the time average of the different steps. We have used molecular dynamics (MD) simulations to get a deeper insight into the detailed interaction scenario at the microsecond scale. This NMR/MD approach enabled us to establish that the tris-dipicolinate complex interacts specifically with arginines and lysines, while the crystallophore explores the protein surface through weak interactions with carboxylates. These observations shed new light on the dynamic interaction properties of these complexes, which will ultimately enable us to propose a crystallization mechanism.

Assessing the Performance of Peptide Force Fields for Modeling the Solution Structural Ensembles of Cyclic Peptides.

Molecular dynamics simulation is a powerful tool for characterizing the solution structural ensembles of cyclic peptides. However, the ability of simulation to recapitulate experimental results and make accurate predictions largely depends on the force fields used. In our work here, we evaluate the performance of seven state-of-the-art force fields in recapitulating the experimental NMR results in water of 12 benchmark cyclic peptides, consisting of 6 cyclic pentapeptides, 4 cyclic hexapeptides, and 2 cyclic heptapeptides. The results show that RSFF2+TIP3P, RSFF2C+TIP3P, and Amber14SB+TIP3P exhibit similar and the best performance, all recapitulating the NMR-derived structure information on 10 cyclic peptides. Amber19SB+OPC successfully recapitulates the NMR-derived structure information on 8 cyclic peptides. In contrast, OPLS-AA/M+TIP4P, Amber03+TIP3P, and Amber14SBonlysc+GB-neck2 could only recapitulate the NMR-derived structure information on 5 cyclic peptides, the majority of which are not well-structured.