Protein Allostery at Atomic Resolution

Insights into mechanisms of MALT1 allostery from NMR and AlphaFold dynamic analyses

Mucosa-associated lymphoid tissue lymphoma-translocation protein 1 (MALT1) is an attractive target for the development of modulatory compounds in the treatment of lymphoma and other cancers. While the three-dimensional structure of MALT1 has been previously determined through X-ray analysis, its dynamic behaviour in solution has remained unexplored. We present here dynamic analyses of the apo MALT1 form along with the E549A mutation. This investigation used NMR 15N relaxation and NOE measurements between side-chain methyl groups. Our findings confirm that MALT1 exists as a monomer in solution, and demonstrate that the domains display semi-independent movements in relation to each other. Our dynamic study, covering multiple time scales, along with the assessment of conformational populations by Molecular Dynamic simulations, Alpha Fold modelling and PCA analysis, put the side chain of residue W580 in an inward position, shedding light at potential mechanisms underlying the allosteric regulation of this enzyme.

NMR and Single-Molecule FRET Insights into Fast Protein Motions and Their Relation to Function

Proteins often undergo large-scale conformational transitions, in which secondary and tertiary structure elements (loops, helices, and domains) change their structures or their positions with respect to each other. Simple considerations suggest that such dynamics should be relatively fast, but the functional cycles of many proteins are often relatively slow. Sophisticated experimental methods are starting to tackle this dichotomy and shed light on the contribution of large-scale conformational dynamics to protein function. In this review, we focus on the contribution of single-molecule Förster resonance energy transfer and nuclear magnetic resonance (NMR) spectroscopies to the study of conformational dynamics. We briefly describe the state of the art in each of these techniques and then point out their similarities and differences, as well as the relative strengths and weaknesses of each. Several case studies, in which the connection between fast conformational dynamics and slower function has been demonstrated, are then introduced and discussed. These examples include both enzymes and large protein machines, some of which have been studied by both NMR and fluorescence spectroscopies.

Targeting SHP2 Cryptic Allosteric Sites for Effective Cancer Therapy

SHP2, a pivotal component downstream of both receptor and non-receptor tyrosine kinases, has been underscored in the progression of various human cancers and neurodevelopmental disorders. Allosteric inhibitors have been proposed to regulate its autoinhibition. However, oncogenic mutations, such as E76K, convert SHP2 into its open state, wherein the catalytic cleft becomes fully exposed to its ligands. This study elucidates the dynamic properties of SHP2 structures across different states, with a focus on the effects of oncogenic mutation on two known binding sites of allosteric inhibitors. Through extensive modeling and simulations, we further identified an alternative allosteric binding pocket in solution structures. Additional analysis provides insights into the dynamics and stability of the potential site. In addition, multi-tier screening was deployed to identify potential binders targeting the potential site. Our efforts to identify a new allosteric site contribute to community-wide initiatives developing therapies using multiple allosteric inhibitors to target distinct pockets on SHP2, in the hope of potentially inhibiting or slowing tumor growth associated with SHP2.

NMR tools to detect protein allostery

Allostery is a fundamental mechanism of cellular homeostasis by intra-protein communication between distinct functional sites. It is an internal process of proteins to steer interactions not only with each other but also with other biomolecules such as ligands, lipids, and nucleic acids. In addition, allosteric regulation is particularly important in enzymatic activities. A major challenge in structural and molecular biology today is unraveling allosteric sites in proteins, to elucidate the detailed mechanism of allostery and the development of allosteric drugs. Here we summarize the recently developed tools and approaches which enable the elucidation of regulatory hotspots and correlated motion in biomolecules, focusing primarily on solution-state nuclear magnetic resonance spectroscopy (NMR). These tools open an avenue towards a rational understanding of the mechanism of allostery and provide essential information for the design of allosteric drugs.

Protein Allostery Study in Cells Using NMR Spectroscopy

Protein allostery is commonly observed in vitro. But how protein allostery behaves in cells is unknown. In this work, a protein monomer-dimer equilibrium system was built with the allosteric effect on the binding characterized using NMR spectroscopy through mutations away from the dimer interface. A chemical shift linear fitting method was developed that enabled us to accurately determine the dissociation constant. A total of 28 allosteric mutations were prepared and grouped to negative allosteric, nonallosteric, and positive allosteric modulators. ∼ 50% of mutations displayed the allosteric-state changes when moving from a buffered solution into cells. For example, there were no positive allosteric modulators in the buffered solution but eight in cells. The change in protein allostery is correlated with the interactions between the protein and the cellular environment. These interactions presumably drive the surrounding macromolecules in cells to transiently bind to the monomer and dimer mutational sites and change the free energies of the two species differently which generate new allosteric effects. These surrounding macromolecules create a new protein allostery pathway that is only present in cells.

The 100-protein NMR spectra dataset: A resource for biomolecular NMR data analysis

Multidimensional NMR spectra are the basis for studying proteins by NMR spectroscopy and crucial for the development and evaluation of methods for biomolecular NMR data analysis. Nevertheless, in contrast to derived data such as chemical shift assignments in the BMRB and protein structures in the PDB databases, this primary data is in general not publicly archived. To change this unsatisfactory situation, we present a standardized set of solution NMR data comprising 1329 2-4-dimensional NMR spectra and associated reference (chemical shift assignments, structures) and derived (peak lists, restraints for structure calculation, etc.) annotations. With the 100-protein NMR spectra dataset that was originally compiled for the development of the ARTINA deep learning-based spectra analysis method, 100 protein structures can be reproduced from their original experimental data. The 100-protein NMR spectra dataset is expected to help the development of computational methods for NMR spectroscopy, in particular machine learning approaches, and enable consistent and objective comparisons of these methods.

Integrated approaches for the recognition of small molecule inhibitors for Toll-like receptor 4

Toll-like receptors (TLRs) are pattern recognition receptors present on the surface of cells playing a crucial role in innate immunity. One of the TLRs, TLR4, recognizes LPS (Lipopolysaccharide) as its ligand leading to the release of anti-inflammatory mediators as well as pro-inflammatory cytokines through signal transduction and domain recruitment. TLR4 homodimerizes at its intracellular TIR (Toll/interleukin-1 receptor) domain that helps in the recruitment of the TRAM/TICAM2 (TIR domain-containing adaptor molecule 2) molecule. TRAM also contains TIR domain which in turn, dimerizes and functions as an adapter protein to further recruit TRIF/TICAM1 (TIR domain-containing adaptor molecule 1) protein for mediating downstream signaling. Apart from LPS, TLR4 also recognizes endogenous ligands like fibrinogen, HMGB1, and hyaluronan in autoimmune conditions and sepsis. We employed computational approaches to target TRAM and recognize small molecule inhibitors from small molecules of natural origin, as contained in the Super Natural II database. Finally, cell reporter assays and NMR studies enabled the identification of promising lead compounds. Hence, this study aims to attenuate the signaling of the TLR4-TRAM-TRIF cascade in these auto-inflammatory conditions.

Water irradiation devoid pulses enhance the sensitivity of 1H,1H nuclear Overhauser effects

The nuclear Overhauser effect (NOE) is one of NMR spectroscopy's most important and versatile parameters. NOE is routinely utilized to determine the structures of medium-to-large size biomolecules and characterize protein-protein, protein-RNA, protein-DNA, and protein-ligand interactions in aqueous solutions. Typical [¹H,¹H] NOESY pulse sequences incorporate water suppression schemes to reduce the water signal that dominates ¹H-detected spectra and minimize NOE intensity losses due to unwanted polarization exchange between water and labile protons. However, at high- and ultra-high magnetic fields, the excitation of the water signal during the execution of the NOESY pulse sequences may cause significant attenuation of NOE cross-peak intensities. Using an evolutionary algorithm coupled with artificial intelligence, we recently designed high-fidelity pulses [Water irrAdiation DEvoid (WADE) pulses] that elude water excitation and irradiate broader bandwidths relative to commonly used pulses. Here, we demonstrate that WADE pulses, implemented into the 2D [¹H,¹H] NOESY experiments, increase the intensity of the NOE cross-peaks for labile and, to a lesser extent, non-exchangeable protons. We applied the new 2D [¹H,¹H] WADE-NOESY pulse sequence to two well-folded, medium-size proteins, i.e., the K48C mutant of ubiquitin and the Raf kinase inhibitor protein. We observed a net increase of the NOE intensities varying from 30 to 170% compared to the commonly used NOESY experiments. The new WADE pulses can be easily engineered into 2D and 3D homo- and hetero-nuclear NOESY pulse sequences to boost their sensitivity.

Light-triggered Modulation of Supramolecular Chirality

Dynamically controlling the supramolecular chirality is of great significance in development of functional chiral materials, which is thus essential for the specific function implementation. As an external energy input, light is remote and accurate for modulating chiral assemblies. In non-polarized light control, some photochemically reactive units (e. g., azobenzene, ɑ-cyanostilbene, spiropyran, anthracene) or photo-induced directionally rotating molecular motors were designed to drive chiral transfer or amplification. Besides, photoexcitation induced assembly based physical approach was also explored recently to regulate supramolecular chirality beyond photochemical reactions. In addition, circularly polarized light was applied to induce asymmetric arrangement of organic molecules and asymmetric photochemical synthesis of inorganic metallic nanostructures, in which both wavelength and handedness of circularly polarized light have effects on the induced supramolecular chirality. Although light-triggered chiral assemblies have been widely applied in photoelectric materials, biomedical fields, soft actuator, chiral catalysis and chiral sensing, there is a lack of systematic review on this topic. In this review, we summarized the recent studies and perspectives in the constructions and applications of light-responsive chiral assembled systems, aiming to provide better knowledge for the development of multifunctional chiral nanomaterials.

Atomic resolution protein allostery from the multi-state structure of a PDZ domain

Recent methodological advances in solution NMR allow the determination of multi-state protein structures and provide insights into structurally and dynamically correlated protein sites at atomic resolution. This is demonstrated in the present work for the well-studied PDZ2 domain of protein human tyrosine phosphatase 1E for which protein allostery had been predicted. Two-state protein structures were calculated for both the free form and in complex with the RA-GEF2 peptide using the exact nuclear Overhauser effect (eNOE) method. In the apo protein, an allosteric conformational selection step comprising almost 60% of the domain was detected with an "open" ligand welcoming state and a "closed" state that obstructs the binding site by changing the distance between the β-sheet 2, α-helix 2, and sidechains of residues Lys38 and Lys72. The observed induced fit-type apo-holo structural rearrangements are in line with the previously published evolution-based analysis covering ~25% of the domain with only a partial overlap with the protein allostery of the open form. These presented structural studies highlight the presence of a dedicated highly optimized and complex dynamic interplay of the PDZ2 domain owed by the structure-dynamics landscape.

Advances in the exact nuclear Overhauser effect 2018-2022

The introduction of the exact nuclear Overhauser enhancement (eNOE) methodology to solution-state nuclear magnetic resonance (NMR) spectroscopy results in tighter distance restraints from NOEs than in convention analysis. These improved restraints allow for higher resolution in structure calculation and even the disentanglement of different conformations of macromolecules. While initial work primarily focused on technical development of the eNOE, structural studies aimed at the elucidation of spatial sampling in proteins and nucleic acids were published in parallel prior to 2018. The period of 2018-2022 saw a continued series of technical innovation, but also major applications addressing biological questions. Here, we review both aspects, covering topics from the implementation of non-uniform sampling of NOESY buildups, novel pulse sequences, adaption of the eNOE to solid-state NMR, advances in eNOE data analysis, and innovations in structural ensemble calculation, to applications to protein, RNA, and DNA structure elucidation.

A Specific pSer/Thr-Pro Motif Generates Interdomain Communication Bifurcations of Two Modes of Pin1 in Solution Nuclear Magnetic Resonance

Peptides mediate the interdomain communication of Pin1 (peptidyl-prolyl cis-trans isomerase) and can regulate its conformation and biochemical functions, providing an idea for drug design using Pin1. Two template peptide sequences have been widely used in the extended or compact state of Pin1 (Cdc25C, E-Q-P-L-pT-P-V-T-D-L; Pintide, W-F-Y-pS-P-R). The way in which specific pSer/Thr-Pro peptides regulate interdomain communication to achieve the opposite state is not clear. In this study, we subdivided the sequence composition of eight types of modified peptides and investigated the interaction with Pin1 by solution nuclear magnetic resonance and molecular dynamics. Demonstrating sequence dependence on the pSer-Pro or pThr-Pro motif and different residues in anchoring the WW domain, the Pin peptide (Pintide, PintideT, Pin25C, and Pin25CT) transmits this concentration accumulation to the PPIase domain, thus exhibiting two anchoring tendencies. However, the Cdc peptide (Cdc25C, Cdc25CS, Cdctide, and CdctideS) has a low binding energy that makes it difficult for the conformation to reach a steady state. In addition, Pin1 is influenced by both compact and extended states, regulated precisely by the sequence as well as by threonine or serine. These results provide new insight into the interdomain communication of Pin1 via pSer/Thr-Pro peptide binding.

On the use of residual dipolar couplings in multi-state structure calculation of two-domain proteins

Residual dipolar couplings (RDCs) are powerful nuclear magnetic resonance (NMR) probes for the structure calculation of biomacromolecules. Typically, an alignment tensor that defines the orientation of the entire molecule relative to the magnetic field is determined either before refinement of individual bond vectors or simultaneously with this refinement. For single-domain proteins this approach works well since all bond vectors can be described within the same coordinate frame, which is given by the alignment tensor. However, novel approaches are sought after for systems where no universal alignment tensor can be used. Here, we present an approach that can be applied to two-domain proteins that enables the calculation of multiple states within each domain as well as with respect to the relative positions of the two domains.

Optimization and validation of multi-state NMR protein structures using structural correlations

Recent advances in the field of protein structure determination using liquid-state NMR enable the elucidation of multi-state protein conformations that can provide insight into correlated and non-correlated protein dynamics at atomic resolution. So far, NMR-derived multi-state structures were typically evaluated by means of visual inspection of structure superpositions, target function values that quantify the violation of experimented restraints and root-mean-square deviations that quantify similarity between conformers. As an alternative or complementary approach, we present here the use of a recently introduced structural correlation measure, PDBcor, that quantifies the clustering of protein states as an additional measure for multi-state protein structure analysis. It can be used for various assays including the validation of experimental distance restraints, optimization of the number of protein states, estimation of protein state populations, identification of key distance restraints, NOE network analysis and semiquantitative analysis of the protein correlation network. We present applications for the final quality analysis stages of typical multi-state protein structure calculations.

High Resolution 31P Field Cycling NMR Reveals Unsuspected Features of Enzyme-Substrate-Cofactor Dynamics

The dynamic interactions of enzymes and substrates underpins catalysis, yet few techniques can interrogate the dynamics of protein-bound ligands. Here we describe the use of field cycling NMR relaxometry to measure the dynamics of enzyme-bound substrates and cofactors in catalytically competent complexes of GMP reductase. These studies reveal new binding modes unanticipated by x-ray crystal structures and reaction-specific dynamic networks. Importantly, this work demonstrates that distal interactions not usually considered part of the reaction coordinate can play an active role in catalysis. The commercialization of shuttling apparatus will make field cycling relaxometry more accessible and expand its use to additional nuclei, promising more intriguing findings to come.

PDBcor: An automated correlation extraction calculator for multi-state protein structures

Allostery and correlated motion are key elements linking protein dynamics with the mechanisms of action of proteins. Here, we present PDBCor, an automated and unbiased method for the detection and analysis of correlated motions from experimental multi-state protein structures. It uses torsion angle and distance statistics and does not require any structure superposition. Clustering of protein conformers allows us to extract correlations in the form of mutual information based on information theory. With PDBcor, we elucidated correlated motion in the WW domain of PIN1, the protein GB3, and the enzyme cyclophilin, in line with reported findings. Correlations extracted with PDBcor can be utilized in subsequent assays including nuclear magnetic resonance (NMR) multi-state structure optimization and validation. As a guide for the interpretation of PDBcor results, we provide a series of protein structure ensembles that exhibit different levels of correlation, including non-correlated, locally correlated, and globally correlated ensembles.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

Reducing the measurement time of exact NOEs by non-uniform sampling

We have previously reported on the measurement of exact NOEs (eNOEs), which yield a wealth of additional information in comparison to conventional NOEs. We have used these eNOEs in a variety of applications, including calculating high-resolution structures of proteins and RNA molecules. The collection of eNOEs is challenging, however, due to the need to measure a NOESY buildup series consisting of typically four NOESY spectra with varying mixing times in a single measurement session. While the 2D version can be completed in a few days, a fully sampled 3D-NOESY buildup series can take 10 days or more to acquire. This can be both expensive as well as problematic in the case of samples that are not stable over such a long period of time. One potential method to significantly decrease the required measurement time of eNOEs is to use non-uniform sampling (NUS) to decrease the number of points measured in the indirect dimensions. The effect of NUS on the extremely tight distance restraints extracted from eNOEs may be very pronounced. Therefore, we investigated the fidelity of eNOEs measured from three test cases at decreasing NUS densities: the 18.4 kDa protein human Pin1, the 4.1 kDa WW domain of Pin1 (both in 3D), and a 4.6 kDa 14mer RNA UUCG tetraloop (2D). Our results show that NUS imparted negligible error on the eNOE distances derived from good quality data down to 10% sampling for all three cases, but there is a noticeable decrease in the eNOE yield that is dependent upon the underlying sparsity, and thus complexity, of the sample. For Pin1, this transition occurred at roughly 40% while for the WW domain and the UUCG tetraloop it occurred at lower NUS densities of 20% and 10%, respectively. We rationalized these numbers through reconstruction simulations under various conditions. The extent of this loss depends upon the number of scans taken as well as the number of peaks to be reconstructed. Based on these findings, we have created guidelines for choosing an optimal NUS density depending on the number of peaks needed to be reconstructed in the densest region of a 2D or 3D NOESY spectrum.