Spectral Density Mapping Research Articles

Insight into the molecular basis of the soluble storage-state to fibrous-state structural transformation of spider pyriform silk

Pyriform (or piriform) silk is a fibrous spider silk used in a composite material to attach web silks to substrates and disparate silk types together. Argiope argentata pyriform silk is formed from a protein with a central repetitive region consisting of a 234-amino acid unit (the "Py unit") repeated 21 times. We have previously shown that the Py unit behaves modularly, with a structured core of 6 α-helices that retains its structuring upon truncation of N- and C-terminal intrinsically disordered segments. Here, we apply solution-state nuclear magnetic resonance (NMR) spectroscopy to investigate the structured core of the Py unit. Based on 15N spin relaxation measurements and reduced spectral density mapping at two field strengths, the α-helical core of the Py unit exhibits minimal backbone motion beyond the tumbling of the globular domain except for one interhelical linker that exhibits a notably higher degree of dynamics. Evaluation of rotational diffusion based on 15N spin relaxation and translational diffusion by pulsed field gradient diffusion NMR demonstrates that the core of the Py unit behaves as a compact structure, with hydrodynamic behaviour consistent with a well-packed protein. Pyriform silk protein undergoes a structural transformation upon conversion from the soluble-state to the fibrous-state, with a loss of αhelical content and gain of β-sheet content. The more dynamic segment we observe midway through the structured core of the Py core, especially in the context of the intrinsically disordered segments that flank the core, would provide an initiation point for decompaction, structural transformation, and/or intermolecular interactions.

Deuterium spin relaxation of fractionally deuterated ribonuclease H using paired 475 and 950MHz NMR spectrometers.

Deuterium (2H) spin relaxation of 13CH2D methyl groups has been widely applied to investigate picosecond-to-nanosecond conformational dynamics in proteins by solution-state NMR spectroscopy. The B0 dependence of the 2H spin relaxation rates is represented by a linear relationship between the spectral density function at three discrete frequencies J(0), J(ωD) and J(2ωD). In this study, the linear relation between 2H relaxation rates at B0 fields separated by a factor of two and the interpolation of rates at intermediate frequencies are combined for a more robust approach for spectral density mapping. The general usefulness of the approach is demonstrated on a fractionally deuterated (55%) and alternate 13C-12C labeled sample of E. coli RNase H. Deuterium relaxation rate constants (R1, R1ρ, RQ, RAP) were measured for 57 well-resolved 13CH2D moieties in RNase H at 1H frequencies of 475MHz, 500MHz, 900MHz, and 950MHz. The spectral density mapping of the 475/950MHz data combination was performed independently and jointly to validate the expected relationship between data recorded at B0 fields separated by a factor of two. The final analysis was performed by jointly analyzing 475/950MHz rates with 700MHz rates interpolated from 500/900MHz data to yield six J(ωD) values for each methyl peak. The J(ω) profile for each peak was fit to the original (τM, Sf2, τf) or extended model-free function (τM, Sf2, Ss2, τf, τs) to obtain optimized dynamic parameters.

Electroencephalograms (EEG) During Mental Arithmetic Task Performance

Abstract: The aim of the present study is to collect the electrical activity of the brain using electroencephalogram (EEG) data from subjects performing a cognitive workload task. During the study, the subjects were involved in intensecognitive activity while performing mental calculations (serial subtractions). The background EEG was recorded from each subject. The primary aim for collecting this dataset is to conduct fluctuation analysis of the EEG duringcognitive activity and to compare the result to the data provided by conventional methods, such as Fourier power spectral density mapping and coherence. It can also beutilized for studying the time-scale characteristics of the involvement of different brain areas in cognition processes and nonlinear characteristics of brain dynamics.

Model-Free or Not?

Relaxation in nuclear magnetic resonance is a powerful method for obtaining spatially resolved, timescale-specific dynamics information about molecular systems. However, dynamics in biomolecular systems are generally too complex to be fully characterized based on NMR data alone. This is a familiar problem, addressed by the Lipari-Szabo model-free analysis, a method that captures the full information content of NMR relaxation data in case all internal motion of a molecule in solution is sufficiently fast. We investigate model-free analysis, as well as several other approaches, and find that model-free, spectral density mapping, LeMaster’s approach, and our detector analysis form a class of analysis methods, for which behavior of the fitted parameters has a well-defined relationship to the distribution of correlation times of motion, independent of the specific form of that distribution. In a sense, they are all “model-free.” Of these methods, only detectors are generally applicable to solid-state NMR relaxation data. We further discuss how detectors may be used for comparison of experimental data to data extracted from molecular dynamics simulation, and how simulation may be used to extract details of the dynamics that are not accessible via NMR, where detector analysis can be used to connect those details to experiments. We expect that combined methodology can eventually provide enough insight into complex dynamics to provide highly accurate models of motion, thus lending deeper insight into the nature of biomolecular dynamics.

Photonic-Radar Based Multiple-Target Tracking Under Complex Traffic-Environments

Recent developments in the state-of-the-art Intelligent Transportation Systems enable autonomous vehicles to offer significant safety services to take appropriate and prompt actions to avoid any probable unfortunate road-hazard. As the utmost functions of the advanced driving assistance system-equipped autonomous vehicles governed by the equipped radar, therefore, the radar system should have the ability to track multiple-targets accurately with high radar-resolutions. Unlike the microwave-radar, the photonic-radar comes out as an attractive candidate owing to provide wide-spectra to attain improved and precise radar-resolutions at low-power requirements along with extended target-range even under severe atmospheric fluctuations. Therefore, a linear frequency-modulated continuous-wave photonic-radar is developed in this work to carry out a radar cross-section-based tracking of multiple mobile-targets in the presence of fog, cloud, and rain. Besides it, some complex real-time traffic-scenarios consisting of multiple mobile-targets make the target-detection, data-association, and classification processes more complicated. Therefore, this work is tested for different multiple-mobile targets in different complicated traffic-scenarios modeled by using MATLAB™ software. The performance of the demonstrated photonic-radar is assessed through the power spectral density and range-Doppler mapping measurements. Furthermore, a comparison of the developed photonic-radar is also established with conventional microwave-radar to present a comparative analysis.

A gratuitous β-Lactamase inducer uncovers hidden active site dynamics of the Staphylococcus aureus BlaR1 sensor domain.

Increasing evidence shows that active sites of proteins have non-trivial conformational dynamics. These dynamics include active site residues sampling different local conformations that allow for multiple, and possibly novel, inhibitor binding poses. Yet, active site dynamics garner only marginal attention in most inhibitor design efforts and exert little influence on synthesis strategies. This is partly because synthesis requires a level of atomic structural detail that is frequently missing in current characterizations of conformational dynamics. In particular, while the identity of the mobile protein residues may be clear, the specific conformations they sample remain obscure. Here, we show how an appropriate choice of ligand can significantly sharpen our abilities to describe the interconverting binding poses (conformations) of protein active sites. Specifically, we show how 2-(2’-carboxyphenyl)-benzoyl-6-aminopenicillanic acid (CBAP) exposes otherwise hidden dynamics of a protein active site that binds β-lactam antibiotics. When CBAP acylates (binds) the active site serine of the β-lactam sensor domain of BlaR1 (BlaRS), it shifts the time scale of the active site dynamics to the slow exchange regime. Slow exchange enables direct characterization of inter-converting protein and bound ligand conformations using NMR methods. These methods include chemical shift analysis, 2-d exchange spectroscopy, off-resonance ROESY of the bound ligand, and reduced spectral density mapping. The active site architecture of BlaRS is shared by many β-lactamases of therapeutic interest, suggesting CBAP could expose functional motions in other β-lactam binding proteins. More broadly, CBAP highlights the utility of identifying chemical probes common to structurally homologous proteins to better expose functional motions of active sites.

Enhanced spectral density mapping through combined multiple-field deuterium 13CH2D methyl spin relaxation NMR spectroscopy.

Quadrupolar relaxation of 2H (D) nuclear spins is a powerful probe of conformational dynamics in biological macromolecules. Deuterium relaxation rate constants are determined by the spectral density function for reorientation of the C-D bond vector at zero, single-quantum, and double-quantum 2H frequencies. In the present work, 2H relaxation rate constants were measured for an E. coli ribonuclease H [U-2H, 15N] ILV-[13CH2D] sample using 400, 500, 800, and 900 MHz NMR spectrometers and analyzed by three approaches to determine spectral density values. First, data recorded at each static magnetic field were analyzed independently. Second, data recorded at 400 and 800 MHz were analyzed jointly and data recorded at other fields were analyzed independently. Third, data recorded at 400 and 500 MHz were interpolated to 450 MHz, and the resulting two pairs of data, corresponding to 400 MHz/800 MHz and 450 MHz/900 MHz, were analyzed jointly. The second and third approaches rely on the identity between the double quantum frequency at the lower field and the single quantum frequency at the higher field. Spectral density values for 32 of the 48 resolvable ILV methyl resonances were fit by the Lipari-Szabo model-free formalism and used to validate the three methods. The three spectral density mapping methods performed equally well in cross validation with data recorded at 700 MHz. However, the third method yielded approximately 10-15% more precise estimates of model-free parameters and consequently provides a general strategy for analysis of 2H spin relaxation data in biological macromolecules.

Characterizing Aciniform Silk Repetitive Domain Backbone Dynamics and Hydrodynamic Modularity.

Spider aciniform (wrapping) silk is a remarkable fibrillar biomaterial with outstanding mechanical properties. It is a modular protein consisting, in Argiope trifasciata, of a core repetitive domain of 200 amino acid units (W units). In solution, the W units comprise a globular folded core, with five α-helices, and disordered tails that are linked to form a ~63-residue intrinsically disordered linker in concatemers. Herein, we present nuclear magnetic resonance (NMR) spectroscopy-based 15N spin relaxation analysis, allowing characterization of backbone dynamics as a function of residue on the ps–ns timescale in the context of the single W unit (W1) and the two unit concatemer (W2). Unambiguous mapping of backbone dynamics throughout W2 was made possible by segmental NMR active isotope-enrichment through split intein-mediated trans-splicing. Spectral density mapping for W1 and W2 reveals a striking disparity in dynamics between the folded core and the disordered linker and tail regions. These data are also consistent with rotational diffusion behaviour where each globular domain tumbles almost independently of its neighbour. At a localized level, helix 5 exhibits elevated high frequency dynamics relative to the proximal helix 4, supporting a model of fibrillogenesis where this helix unfolds as part of the transition to a mixed α-helix/β-sheet fibre.

Spectral density mapping at multiple magnetic fields suitable for [formula omitted] NMR relaxation studies

Standard spectral density mapping protocols, well suited for the analysis of 15N relaxation rates, introduce significant systematic errors when applied to 13C relaxation data, especially if the dynamics is dominated by motions with short correlation times (small molecules, dynamic residues of macromolecules). A possibility to improve the accuracy by employing cross-correlated relaxation rates and on measurements taken at several magnetic fields has been examined. A suite of protocols for analyzing such data has been developed and their performance tested. Applicability of the proposed protocols is documented in two case studies, spectral density mapping of a uniformly labeled RNA hairpin and of a selectively labeled disaccharide exhibiting highly anisotropic tumbling. Combination of auto- and cross-correlated relaxation data acquired at three magnetic fields was applied in the former case in order to separate effects of fast motions and conformational or chemical exchange. An approach using auto-correlated relaxation rates acquired at five magnetic fields, applicable to anisotropically moving molecules, was used in the latter case. The results were compared with a more advanced analysis of data obtained by interpolation of auto-correlated relaxation rates measured at seven magnetic fields, and with the spectral density mapping of cross-correlated relaxation rates. The results showed that sufficiently accurate values of auto- and cross-correlated spectral density functions at zero and 13C frequencies can be obtained from data acquired at three magnetic fields for uniformly 13C-labeled molecules with a moderate anisotropy of the rotational diffusion tensor. Analysis of auto-correlated relaxation rates at five magnetic fields represents an alternative for molecules undergoing highly anisotropic motions.

The Dynamic Landscape of the Full-Length HIV-1 Transactivator of Transcription.

The type 1 human immunodeficiency virus (HIV-1) transactivator of transcription (Tat) is a small RNA-binding protein essential for viral gene expression and replication. It has also been shown to bind to a large number of human proteins and to modulate many different cellular activities. We have used nuclear magnetic resonance (NMR) spectroscopy and hydrogen exchange chemistry to measure backbone dynamics over the millisecond to picosecond time scales. Sequential backbone assignment was facilitated by several isotope labeling schemes, including uniform labeling, site-specific labeling, and unlabeling. (15)N NMR relaxation parameters were measured and analyzed by reduced spectral density mapping and the Lipari-Szabo Model-Free approach to characterize the backbone dynamics on the picosecond to nanosecond time scale. The results indicate that the protein exists in an extended disordered conformational ensemble. NMR relaxation dispersion profiles show that on the millisecond time scale no conformational exchange is detected for any of the residues, supporting the model of a disordered backbone. NMR chemical shift differences from random coil values suggest that some segments of the protein have a modest propensity to fold; comparison to X-ray diffraction structures of Tat complexes indicates that some segments of the protein function through an induced-fit mechanism whereas other segments likely operate by conformational selection. Surprisingly, measured hydrogen exchange rates are higher than predicted for a disordered polymer, but this is explained as being caused by the high net charge on the protein that enhances base-catalyzed hydrogen exchange. The dynamics results provide a deeper understanding of the protein conformational ensemble and form a foundation for future studies of the conformational changes that accompany the formation of the superelongation complex that activates viral transcription.

Spectral density mapping protocols for analysis of molecular motions in disordered proteins

Spectral density mapping represents the method of choice for investigations of molecular motions of intrinsically disordered proteins (IDPs). However, the current methodology has been developed for well-folded proteins. In order to find conditions for a reliable analysis of relaxation of IDPs, accuracy of the current reduced spectral density mapping protocols applied to IDPs was examined and new spectral density mapping methods employing cross-correlated relaxation rates have been designed. Various sources of possible systematic errors were analyzed theoretically and the presented approaches were tested on a partially disordered protein, delta subunit of bacterial RNA polymerase. Results showed that the proposed protocols provide unbiased description of molecular motions of IDPs and allow to separate slow exchange from fast dynamics.

Model-Free Spectral Density Mapping Applied to Dynamics of Rhodopsin in Solid-State NMR Spectroscopy

Crystal structures of rhodopsin are available, yet details of the activation mechanism remain unknown [1,2]. We applied solid-state 2H NMR to investigate structural and dynamical changes occurring in the process of rhodopsin activation. From the 2H NMR spectra, molecular mobility can be obtained by calculating the segmental order parameters from the residual quadrupolar couplings (RQCs). Moreover 2H nuclear spin relaxation rates related to the dynamics can be measured together with the RQCs [2]. Site-specific 2H labels were introduced into different methyl groups of retinal, and relaxation rate measurements were performed as a function of temperature (−30 to −150°C) [3]. Model-free analysis employed an irreducible representation of the combined 2H NMR line shape and relaxation data. Fluctuations of the irreducible components with respect to the average values are characterized by the individual spectral densities of motion evaluated at characteristic frequencies: J0(0), J1(ω0), and J2(2ω0), where ω0 is the nuclear resonance frequency [4]. Differences in the spectral densities manifest details of the methyl group motions within the retinal binding pocket at low temperature. At the high temperature limit, J1(ω0) and J2(2ω0) are insensitive to details of motion and collapse to a universal curve, thus substantiating the validity of the model-free analysis. Further analysis of J1(ω0) and J2(2ω0) involved simultaneous temperature-dependent fitting in terms of both diffusion and jump models for the methyl group dynamics. We conclude that spectral density analysis in terms of fluctuations of nuclear spin Hamiltonian will help us understand the activation mechanism and molecular dynamics of rhodopsin and related G protein-coupled receptors.[1] A.V. Struts et al. (2011) PNAS 108, 8263-8268. [2] M.F. Brown et al. (2010) BBA 1798, 177-193. [3] A.V. Struts et al. (2011) NSMB 18, 392-394. [4] M.F. Brown (1982) JCP 77, 1576-1599.

Influence of Heme Post-Translational Modification and Distal Ligation on the Backbone Dynamics of a Monomeric Hemoglobin

The cyanobacterium Synechococcus sp. PCC 7002 uses a hemoglobin of the truncated lineage (GlbN) in the detoxification of reactive species generated in the assimilation of nitrate. In view of a sensing or enzymatic role, several states of GlbN are of interest with respect to its structure-activity relationship. Nuclear magnetic resonance spectroscopy was applied to compare the structure and backbone dynamics of six GlbN forms differing in their oxidation state [Fe(II) or Fe(III)], distal ligand to the iron (histidine, carbon monoxide, or cyanide), or heme post-translational modification (b heme or covalently attached heme). Structural properties were assessed with pseudocontact shift calculations. (15)N relaxation data were analyzed by reduced spectral density mapping (picosecond to nanosecond motions) and by inspection of elevated R(2) values (microsecond to millisecond motions). On the picosecond to nanosecond time scale, GlbN exhibited little flexibility and was unresponsive to the differences among the various forms. Regions of slightly higher mobility were the CE turn, the EF loop, and the H-H' kink. In contrast, fluctuations on the microsecond to millisecond time scale depended on the form. Cyanide binding to the ferric state did not enhance motions, whereas reduction to the ferrous bis-histidine state resulted in elevated R(2) values for several amides. This response was attributed, at least in part, to a weakening of the distal histidine coordination. Carbon monoxide binding quenched some of these fluctuations. The results emphasized the role of the distal ligand in dictating backbone flexibility and illustrated the multiple ways in which motions are controlled by the hemoglobin fold.

Water Model Tuning for Improved Reproduction of Rotational Diffusion and NMR Spectral Density

A water model for molecular simulation was optimized to improve the reproduction of translational and rotational diffusion of pure water and proteins. The SPC/E(b) model was developed from the original SPC/E model with a slight increase of the O-H bond length of 1%. This tuning has significantly improved the translational and rotational diffusion when compared to the experimental values, whereas only small changes were observed in the other thermodynamic properties examined. The overall tumbling correlation times (τ(p)) from ubiquitin, protein G, bovine pancreatic trypsin inhibitor, and barstar C42/80A were successfully reproduced using the SPC/E(b) model. Calculated site-specific spectral densities of the main chain amide bond rotation in ubiquitin and protein G were in good agreement with those derived from nuclear magnetic resonance reduced spectral density mapping. The SPC/E(bT) model was also developed with temperature-dependent bond-length tuning to facilitate reproduction of the experimental τ(p) around room temperature.

Insights into internal dynamics of 6-phosphogluconolactonase from Trypanosoma brucei studied by nuclear magnetic resonance and molecular dynamics

Nuclear magnetic resonance is used to investigate the backbone dynamics in 6-phosphogluconolactonase from Trypanosoma brucei (Tb6PGL) with (holo-) and without (apo-) 6-phosphogluconic acid as ligand. Relaxation data were analyzed using the model-free approach and reduced spectral density mapping. Comparison of predictions, based on 77 ns molecular dynamics simulations, with the observed relaxation rates gives insight into dynamical properties of the protein and their alteration on ligand binding. Data indicate dynamics changes in the vicinity of the binding site. More interesting is the presence of perturbations located in remote regions of this well-structured globular protein in which no large-amplitude motions are involved. This suggests that delocalized changes in dynamics that occur upon binding could be a general feature of protein-target interactions.

NMR backbone dynamics studies of human PED/PEA-15 outline protein functional sites

PED/PEA-15 (phosphoprotein enriched in diabetes/phosphoprotein enriched in astrocytes) is a ubiquitously expressed protein and a key regulator of cell growth and glucose metabolism. PED/PEA-15 mediates both homotypic and heterotypic interactions and is constituted by an N-terminal canonical death effector domain and a C-terminal tail. In the present study, the backbone dynamics of PED/PEA-15 via (15)N R(1) and R(2) and steady-state [(1)H]-(15)N NOE measurements is reported. The dynamic parameters were analyzed using both Lipari-Szabo model-free formalism and a reduced spectral density mapping approach. The results obtained define a polar and charged surface of the death effector domain characterized by internal motions in the micro- to millisecond timescale, which is crucial for the multiple heterotypic functional protein-protein interactions in which PED/PEA-15 is involved. The present study contributes to a better understanding of the molecular basis of the PED/PEA-15 functional interactions and provides a more detailed surface for the design and development of PED/PEA-15 binders.

Interpretation of biomolecular NMR spin relaxation parametersThis paper is one of a selection of papers published in this special issue entitled “Canadian Society of Biochemistry, Molecular & Cellular Biology 52nd Annual Meeting — Protein Folding: Principles and Diseases” and has undergone the Journal's usual peer review process.

Biomolecular nuclear magnetic resonance (NMR) spin relaxation experiments provide exquisite information on the picosecond to nanosecond timescale motions of bond vectors. Spin-lattice (T1) and spin-spin (T2) relaxation times and the steady-state nuclear Overhauser effect (NOE) are the first set of parameters extracted from typical 15N or 13C NMR relaxation experiments. Therefore, verifying that T1, T2, and NOE are consistent with theoretical predictions is an important step before carrying out the more detailed model-free and reduced spectral density mapping analyses commonly employed. In this mini-review, we discuss the essential motional parameters used to describe biomolecular dynamics in the context of a variety of examples of folded and intrinsically disordered proteins and peptides in aqueous and membrane mimetic environments. Estimates of these parameters can be used as input for an online interface, introduced herein, allowing plotting of trends of T1, T2, and NOE with magnetic field strength. The plots may serve as a first-check to the spectroscopist preparing to embark on a detailed NMR relaxation analysis.

NMR studies of a channel protein without membranes: Structure and dynamics of water-solubilized KcsA

Structural studies of polytopic membrane proteins are often hampered by the vagaries of these proteins in membrane mimetic environments and by the difficulties in handling them with conventional techniques. Designing and creating water-soluble analogues with preserved native structures offer an attractive alternative. We report here solution NMR studies of WSK3, a water-soluble analogue of the potassium channel KcsA. The WSK3 NMR structure (PDB ID code 2K1E) resembles the KcsA crystal structures, validating the approach. By more stringent comparison criteria, however, the introduction of several charged residues aimed at improving water solubility seems to have led to the possible formations of a few salt bridges and hydrogen bonds not present in the native structure, resulting in slight differences in the structure of WSK3 relative to KcsA. NMR dynamics measurements show that WSK3 is highly flexible in the absence of a lipid environment. Reduced spectral density mapping and model-free analyses reveal dynamic characteristics consistent with an isotropically tumbling tetramer experiencing slow (nanosecond) motions with unusually low local ordering. An altered hydrogen-bond network near the selectivity filter and the pore helix, and the intrinsically dynamic nature of the selectivity filter, support the notion that this region is crucial for slow inactivation. Our results have implications not only for the design of water-soluble analogues of membrane proteins but also for our understanding of the basic determinants of intrinsic protein structure and dynamics.

Entropy-driven cAMP-dependent Allosteric Control of Inhibitory Interactions in Exchange Proteins Directly Activated by cAMP

Exchange proteins directly activated by cAMP (EPACs) are guanine nucleotide-exchange factors for the small GTPases Rap1 and Rap2 and represent a key receptor for the ubiquitous cAMP second messenger in eukaryotes. The cAMP-dependent activation of apoEPAC is typically rationalized in terms of a preexisting equilibrium between inactive and active states. Structural and mutagenesis analyses have shown that one of the critical determinants of the EPAC activation equilibrium is a cluster of salt bridges formed between the catalytic core and helices alpha1 and alpha2 at the N terminus of the cAMP binding domain and commonly referred to as ionic latch (IL). The IL stabilizes the inactive states in a closed topology in which access to the catalytic domain is sterically occluded by the regulatory moiety. However, it is currently not fully understood how the IL is allosterically controlled by cAMP. Chemical shift mapping studies consistently indicate that cAMP does not significantly perturb the structure of the IL spanning sites within the regulatory region, pointing to cAMP-dependent dynamic modulations as a key allosteric carrier of the cAMP-signal to the IL sites. Here, we have therefore investigated the dynamic profiles of the EPAC1 cAMP binding domain in its apo, cAMP-bound, and Rp-cAMPS phosphorothioate antagonist-bound forms using several 15N relaxation experiments. Based on the comparative analysis of dynamics in these three states, we have proposed a model of EPAC activation that incorporates the dynamic features allosterically modulated by cAMP and shows that cAMP binding weakens the IL by increasing its entropic penalty due to dynamic enhancements.

Conformational change upon ligand binding and dynamics of the PDZ domain from leukemia‐associated Rho guanine nucleotide exchange factor

Leukemia-associated Rho guanine nucleotide exchange factor (LARG) is a RhoA-specific guanine nucleotide exchange factor (GEF) that can activate RhoA. The PDZ (PSD-95/Disc-large/ZO-1 homology) domain of LARG interacts with membrane receptors, which can relay extracellular signals to RhoA signal transduction pathways. Until now there is no structural and dynamic information about these interactions. Here we report the NMR structures of the LARG PDZ in the apo form and in complex with the plexin-B1 C-terminal octapeptide. Unobservable resonances of the residues in betaB/betaC and betaE/alphaB loops in apo state were observed in the complex state. A distinct region of the binding groove in the LARG PDZ was found to undergo conformational change compared with other PDZs. Analysis of the (15)N relaxation data using reduced spectral density mapping shows that the apo LARG PDZ (especially its ligand-binding groove) is flexible and exhibits internal motions on both picosecond to nanosecond and microsecond to millisecond timescales. Mutagenesis and thermodynamic studies indicate that the conformation of the betaB/betaC and betaE/alphaB loops affects the PDZ-peptide interaction. It is suggested that the conformational flexibility could facilitate the change of structures upon ligand binding.