Search

Chaperones TAPBPR and tapasin associate with class-I major histocompatibility complexes (MHC-I) to promote optimization (editing) of peptide cargo. Here, we use solution NMR to investigate the mechanism of peptide exchange. We identify TAPBPR-induced conformational changes on conserved MHC-I molecular surfaces, consistent with our independently determined X-ray structure of the complex. Dynamics present in the empty MHC-I are stabilized by TAPBPR, and become progressively dampened with increasing peptide occupancy. Incoming peptides are recognized according to the global stability of the final pMHC-I product, and anneal in a native-like conformation to be edited by TAPBPR. Our results demonstrate an inverse relationship between MHC-I peptide occupancy and TAPBPR binding affinity, where the lifetime and structural features of transiently bound peptides controls the regulation of a conformational switch, located near the TAPBPR binding site, which triggers TAPBPR release. These results suggest a similar mechanism for the function of tapasin in the peptide-loading complex., Graphical abstract

NMR spectroscopy of membrane proteins involved in electron transport is difficult due to the presence of both the lipids and paramagnetic centers. Here we report the solution NMR study of the NADPH-cytochrome P450 oxidoreductase (POR) in its reduced and oxidized states. We interrogate POR, first, in its truncated soluble form (70 kDa), which is followed by experiments with the full-length protein incorporated in a lipid nanodisc (240 kDa). To overcome paramagnetic relaxation in the reduced state of POR as well as the signal broadening due to its high molecular weight, we utilized the methyl-TROSY approach. Extrinsic 13C-methyl groups were introduced by modifying the engineered surface-exposed cysteines with methyl-methanethiosulfonate. Chemical shift dispersion of the resonances from different sites in POR was sufficient to monitor differential effects of the reduction–oxidation process and conformation changes in the POR structure related to its function. Despite the high molecular weight of the POR-nanodisc complex, the surface-localized 13C-methyl probes were sufficiently mobile to allow for signal detection at 600 MHz without perdeuteration. This work demonstrates a potential of the solution methyl-TROSY in analysis of structure, dynamics, and function of POR, which may also be applicable to similar paramagnetic and flexible membrane proteins.

Chitosan interaction with chitosanase was examined through analysis of spectral line shapes in the NMR HSQC titration experiments. We established that the substrate, chitosan hexamer, binds to the enzyme through the three-state induced-fit mechanism with fast formation of the encounter complex followed by slow isomerization of the bound-state into the final conformation. Mapping of the chemical shift perturbations in two sequential steps of the mechanism highlighted involvement of the substrate-binding subsites and the hinge region in the binding reaction. Equilibrium parameters of the three-state model agreed with the overall thermodynamic dissociation constant determined by ITC. This study presented the first kinetic evidence of the induced-fit mechanism in the glycoside hydrolases.

The ability of high-resolution NMR spectroscopy to readout the response of molecular interactions at multiple atomic sites presents a unique capability to define thermodynamic equilibrium constants and kinetic rate constants for complex, multiple-step biological interactions. Nonetheless, the extraction of the relevant equilibrium binding and rate constants requires the appropriate analysis of not only a readout that follows the equilibrium concentrations of typical binding titration curves, but also the lineshapes of NMR spectra. To best take advantage of NMR data for characterizing molecular interactions, we developed NmrLineGuru, a software tool with a user-friendly graphical user interface (GUI) to model two-state, three-state, and four-state binding processes. Application of NmrLineGuru is through stand-alone GUIs, with no dependency on other software and no scripted input. NMR spectra can be fitted or simulated starting with user-specified input parameters and a chosen kinetic model. The ability to both simulate and fit NMR spectra provides the user the opportunity to not only determine the binding parameters that best reproduce the measured NMR spectra for the selected kinetic model, but to also query the possibility that alternative models agree with the data. NmrLineGuru is shown to provide an accurate, quantitative analysis of complex molecular interactions.

The family of cystathionine-β-synthase (CBS)-pair domain divalent metal cation transport mediators (CNNMs) is composed of four integral membrane proteins associated with Mg

Zeolitic imidazolate frameworks (ZIFs) with intrinsic photoactivity represent an excellent scaffold for a single-site photocatalytic center. However, their application in visible-light-driven photocatalysis is hampered by their light-harvesting ability due to limited coverage in the realm of the solar spectrum and relatively low extinction coefficient. In this work, we report the first example of the enhanced light-harvesting property of ZIF-67 via energy transfer from a molecular photosensitizer, that is, the RuN3 complex, sensitized on the surface of ZIF-67. Using transient absorption spectroscopy, we show that energy transfer occurs from the excited RuN3 to ZIF-67 with ∼86.9% efficiency. More importantly, this RuN3/ZIF-67 hybrid system exhibits significantly enhanced photocatalytic activity for H2 production from water, which can be attributed to efficient energy transfer from RuN3.

Summary The family of cystathionine-β-synthase (CBS)-pair domain divalent metal cation transport mediators (CNNMs) is composed of four integral membrane proteins associated with Mg2+ transport. Structurally, CNNMs contain large cytosolic regions composed of a CBS-pair and a cyclic nucleotide-binding homology (CNBH) domain. How these regulate Mg2+ transport activity is unknown. Here, we determined the crystal structures of cytosolic fragments in two conformations: Mg2+-ATP-analog bound and ligand free. The structures reveal open and closed conformations with functionally important contacts not observed in structures of the individual domains. We also identified a second Mg2+-binding region in the CBS-pair domain and a different dimerization interface for the CNBH domain. Analytical ultracentrifugation and isothermal titration calorimetry experiments revealed a tight correlation between Mg2+-ATP binding and protein dimerization. Mutations that blocked either function prevented cellular Mg2+ efflux activity. The results suggest Mg2+ efflux is regulated by conformational changes associated with Mg2+-ATP binding to CNNM CBS-pair domains.

NADPH-cytochrome P450 oxidoreductase (POR) is a membrane protein in the endoplasmic reticulum of eukaryotic cells. POR is as a key reducing partner for a number of cytochrome P450 proteins involved in different metabolic degradation and signaling pathways. Preparation of the full-length recombinant POR expressed in bacteria has been reported and, typically, involved the use of Triton X-100 detergent for extraction of the overexpressed POR from bacterial membranes. However, extraction efficiency is always relatively low hindering structural studies, particularly—the NMR spectroscopy requiring isotopic enrichment. In this paper, we assessed the effect of a variety of detergents and additives on the efficiency of the membrane-extraction step in POR preparation protocol. We evaluated non-ionic detergents with the variable hydrophobicity (Triton X-100, X-114, and X-405) and structure (Triton X-100, TWEEN-20, Brij-35), a zwitterionic/non-ionic detergent combination (Triton X-100 and CHAPS), as well as a range of alkylamines and polyamines as additives to the conventional extraction buffer containing Triton X-100. None of the detergents or detergent-additive combinations yielded better extraction efficiency than the conventional protocol with the Triton X-100. Lack of variation of the extraction yield allows to hypothesize that the conventional protocol extracts all of the available natively-folded monomeric POR while the remaining fraction is possibly an unfolded aggregated POR, which did not insert in the membranes during expression. We propose that the yield of soluble POR may be increased by a careful optimization of expression conditions while monitoring the distribution of POR between soluble and insoluble fractions in the detergent extraction step.

Supported phospholipid bilayers are a convenient model of cellular membranes in studies of membrane biophysics and protein-lipid interactions. Traditionally, supported lipid bilayers are formed on a flat surface of a glass slide to be observed through fluorescence microscopes. This paper describes a method to enable fluorescence detection from the supported lipid bilayers using standard horizontal-beam spectrofluorometers instead of the microscopes. In the proposed approach, the supported lipid bilayers are formed on the inner optical surfaces of the standard fluorescence microcell. To enable observation of the bilayer absorbed on the cell wall, the microcell is placed in a standard fluorometer cell holder and specifically oriented to expose the inner cell walls to both excitation and emission channels with a help of the custom cell adaptor. The signal intensity from supported bilayers doped with 1 % (mol) of rhodamine-labeled lipid in the standard 3-mm optical microcell was equivalent to fluorescence of the 70-80 nM reference solution of rhodamine recorded in a commercial microcell adaptor. Because no modifications to the instruments are required in this method, a variety of steady-state and time-domain fluorescence measurements of the supported phospholipid bilayers may be performed with the spectral resolution using standard horizontal-beam spectrofluorometers.

Cellular membranes are heterogeneous planar lipid bilayers displaying lateral phase separation with the nanometer-scale liquid-ordered phase (also known as "lipid rafts") surrounded by the liquid-disordered phase. Many membrane-associated proteins were found to permanently integrate into the lipid rafts, which is critical for their biological function. Isoforms H and N of Ras GTPase possess a unique ability to switch their lipid domain preference depending on the type of bound guanine nucleotide (GDP or GTP). This behavior, however, has never been demonstrated in vitro in model bilayers with recombinant proteins and therefore has been attributed to the action of binding of Ras to other proteins at the membrane surface. In this paper, we report the observation of the nucleotide-dependent switch of lipid domain preferences of the semisynthetic lipidated N-Ras in lipid raft vesicles in the absence of additional proteins. To detect segregation of Ras molecules in raft and disordered lipid domains, we measured Förster resonance energy transfer between the donor fluorophore, mant, attached to the protein-bound guanine nucleotides, and the acceptor, rhodamine-conjugated lipid, localized into the liquid-disordered domains. Herein, we established that N-Ras preferentially populated raft domains when bound to mant-GDP, while losing its preference for rafts when it was associated with a GTP mimic, mant-GppNHp. At the same time, the isolated lipidated C-terminal peptide of N-Ras was found to be localized outside of the liquid-ordered rafts, most likely in the bulk-disordered lipid. Substitution of the N-terminal G domain of N-Ras with a homologous G domain of H-Ras disrupted the nucleotide-dependent lipid domain switch.

Cellular membranes are heterogeneous planar lipid bilayers displaying lateral phase separation with the nanometer-scale liquid-ordered phase (aka “lipid rafts” or Lo) surrounded by the liquid-disordered phase (Ld). Many membrane-associated proteins were found to stably integrate in the rafts, which is critical for their biological function. Isoforms H and N of Ras GTPase possess a unique ability to switch their lipid domain preference depending on the type of bound guanine nucleotide (GDP or GTP). This behavior, however, has never been reproducedin vitroin model bilayers with recombinant proteins, and therefore has been attributed to action of other proteins binding Ras at the membrane surface. In this paper, we report the observation of the nucleotide-dependent switch of lipid domain preferences of the semisynthetic lipidated N-Ras in raft lipid vesiclesin the absence of other proteins. To detect segregation of Ras molecules in raft and disordered lipid domains, we measured Förster Resonance Energy Transfer (FRET) between the donor fluorophore, mant, attached to the protein-bound guanine nucleotides, and the acceptor, rhodamine-conjugated lipid, localized to the liquid-disordered domains. We demonstrated that N-Ras preferentially populated raft domains when bound to mant-GDP, while losing preference for rafts when it was associated with a GTP mimic, mant-GppNHp. At the same time, the isolated lipidated C-terminal peptide of N-Ras was found localized outside of the liquid-ordered rafts, most likely—in the bulk disordered lipid.

Ras GTPase is a peripheral membrane protein central to cellular signaling of growth and proliferation. Membrane attachment is critical for a range of Ras activities, therefore, ability to make faithful in-vitro samples of a mem-brane-bound Ras for detailed biophysical studies is a highly desirable goal. In this manuscript, we are describing preparation of a large-scale sample of isotopically labeled H-Ras conjugated to lipid nanodiscs. We demonstrate that the Ras-nanodisc sample is fairly stable to allow for a range of Nuclear Magnetic Resonance (NMR) and other biophysical measurements. The need to achieve a homogeneous protein-nanodisc ratio is also emphasized.

Conformational transitions between closed and open states in the NADPH-cytochrome P450 oxidoreductase (POR) play a critical role in its electron-transport function. In this study, we determined rotational diffusion coefficients of the EDANS fluorophore attached to the cytosolic POR construct lacking the N-terminal transmembrane region. We identified two dynamic modes, slow and fast, which are interpreted as the rotational diffusion of POR as a whole and the local domain motion, respectively. Timescale of the local rotational diffusion component suggests that it may correspond to the transient opening of the fully oxidized POR structure.

NMR spectroscopy of membrane proteins involved in electron transport is difficult due to the presence of both the lipids and paramagnetic centers. Here we report the solution NMR study of the NADPH-cytochrome P450 oxidoreductase (POR) in its reduced and oxidized states. We interrogate POR, first, in its truncated soluble form (70 kDa), which is followed by experiments with the full-length protein incorporated in a lipid nanodisc (240 kDa). To overcome paramagnetic relaxation in the reduced state of POR as well as the signal broadening due to its high molecular weight, we utilized the methyl-TROSY approach. Extrinsic 13C-methyl groups were introduced by modifying the engineered surface-exposed cysteines with methyl-methanethiosulfonate (MMTS). Chemical shift dispersion of the resonances from different sites in POR was sufficient to monitor differential effects of the reduction-oxidation process and conformation changes in the POR structure related to its function. Despite the high molecular weight of the POR-nanodisc complex, the surface-localized 13C-methyl probes were sufficiently mobile to allow for signal detection at 600 MHz without perdeuteration. This work demonstrates a potential of the solution methyl-TROSY in analysis of structure, dynamics, and function of POR, which may also be applicable to similar paramagnetic and flexible membrane proteins.

This manuscript describes a workflow for analysis of transient absorption (TA) spectroscopy data using Integrative Data Analysis Platforms (IDAP) software package. Time-dependent spectral series are analyzed through evaluation of the isosbestic point and kinetics of excited state and ground-state bleach decays. Model fitting and selection based on Akaike’s Information Criterion is discussed. As a practical example, we analyze excitation decays of a common protein label, Alexa Fluor 647.

N-acetyl-d-quinovosamine (2-acetamido-2,6-dideoxy-d-glucose, QuiNAc) occurs in the polysaccharide structures of many Gram-negative bacteria. In the biosynthesis of QuiNAc-containing polysaccharides, UDP-QuiNAc is the hypothetical donor of the QuiNAc residue. Biosynthesis of UDP-QuiNAc has been proposed to occur by 4,6-dehydration of UDP-N-acetyl-d-glucosamine (UDP-GlcNAc) to UDP-2-acetamido-2,6-dideoxy-d-xylo-4-hexulose followed by reduction of this 4-keto intermediate to UDP-QuiNAc. Several specific dehydratases are known to catalyze the first proposed step. A specific reductase for the last step has not been demonstrated in vitro, but previous mutant analysis suggested that Rhizobium etli gene wreQ might encode this reductase. Therefore, this gene was cloned and expressed in Escherichia coli, and the resulting His6-tagged WreQ protein was purified. It was tested for 4-reductase activity by adding it and NAD(P)H to reaction mixtures in which 4,6-dehydratase WbpM had acted on the precursor substrate UDP-GlcNAc. Thin layer chromatography of the nucleotide sugars in the mixture at various stages of the reaction showed that WbpM converted UDP-GlcNAc completely to what was shown to be its 4-keto-6-deoxy derivative by NMR and that addition of WreQ and NADH led to formation of a third compound. Combined gas chromatography-mass spectrometry analysis of acid hydrolysates of the final reaction mixture showed that a quinovosamine moiety had been synthesized after WreQ addition. The two-step reaction progress also was monitored in real time by NMR. The final UDP-sugar product after WreQ addition was purified and determined to be UDP-d-QuiNAc by one-dimensional and two-dimensional NMR experiments. These results confirmed that WreQ has UDP-2-acetamido-2,6-dideoxy-d-xylo-4-hexulose 4-reductase activity, completing a pathway for UDP-d-QuiNAc synthesis in vitro.

Molecular chaperones TAPBPR (TAP-binding protein related) and tapasin associate with class-I major histocompatibility complex (MHC-I) molecules to promote optimization (editing) of peptide cargo. Using solution NMR, we investigate the molecular mechanism of peptide exchange performed by the 90 kDa chaperone protein complex. We identify TAPBPR-induced conformational changes on conserved MHC-I surfaces, consistent with our independently determined X-ray structure of the empty complex. Conformational dynamics present in the empty MHC-I are stabilized by TAPBPR in a peptide-deficient complex, and become progressively dampened with increasing peptide occupancy. Incoming peptides are recognized by the chaperoned groove according to the global stability of the final pMHC-I product, and anneal in a native-like conformation. Our results demonstrate an inverse relationship between MHC-I occupancy by peptide and the affinity of TAPBPR for such pMHC-I molecules, where the lifetime of transiently bound peptides controls the dynamic regulation of a conformational switch, located near the TAPBPR binding site, which triggers TAPBPR release. We further discuss the role of protein dynamics in shaping chaperone specificity towards different human and murine class-I MHC allotypes, and present the high-resolution cryoEM structure of a human A*02/TAPBPR complex with novel insights into the antigen editing mechanism. These results suggest a similar mechanism for the editing function of tapasin in the peptide-loading complex, and provide a molecular blueprint for the design of novel chaperones with tailored antigen editing functions.

NADPH-cytochrome P450 oxidoreductase (CYPOR) was shown to undergo large conformational rearrangements in its functional cycle. Using a new Forster resonance energy transfer (FRET) approach based on femtosecond transient absorption spectroscopy (TA), we determined the donor–acceptor distance distribution in the reduced and oxidized states of CYPOR. The unmatched time resolution of TA allowed the quantitative assessment of the donor–acceptor FRET, indicating that CYPOR assumes a closed conformation in both reduced and oxidized states in the absence of the redox partner. The described ultrafast TA measurements of FRET with readily available red–infrared fluorescent labels open new opportunities for structural studies in chromophore-rich proteins and their complexes.

Ras is a small monomeric GTPase acting as molecular switch in multiple cellular processes. The N-terminal G domain of Ras binds GTP or GDP accompanied by a magnesium ion, which is strictly required for GTPase activity and performs a structural role. Another ion-binding site on the opposite face of the G domain has been recently observed to specifically associate with calcium acetate in the crystal [Buhrman, G., et al. (2010) Proc. Natl. Aacd. Sci. U.S.A. 107, 4931-4936]. In this article, we report thermodynamic measurements of the affinity and specificity of the remote ion-binding site in H-Ras as observed in solution. Using (15)N-(1)H nuclear magnetic resonance spectroscopy, we determined that, in contrast to the crystalline state, the remote site in solution is specific for a divalent cation, binding both calcium and magnesium with anions playing a minimal role. The affinity of the remote site for divalent cations is in the low millimolar range and remarkably different for GDP- and GppNHp-bound forms of the G domain, indicating that the GTP-binding pocket and the remote site are allosterically coupled through the distance of more than 25 Å. Considering that the remote site is oriented toward the membrane surface in vivo, we hypothesize that its cognate biological ligand might be a positively charged group extending from a lipid or an integral membrane protein.

Biological function of proteins relies on conformational transitions and binding of specific ligands. Protein-ligand interactions are thermodynamically and kinetically coupled to conformational changes in protein structures as conceptualized by the models of pre-existing equilibria and induced fit. NMR spectroscopy is particularly sensitive to complex ligand-binding modes-NMR line-shape analysis can provide for thermodynamic and kinetic constants of ligand-binding equilibria with the site-specific resolution. However, broad use of line shape analysis is hampered by complexity of NMR line shapes in multi-state systems. To facilitate interpretation of such spectral patterns, I computationally explored systems where isomerization or dimerization of a protein (receptor) molecule is coupled to binding of a ligand. Through an extensive analysis of multiple exchange regimes for a family of three-state models, I identified signature features to guide an NMR experimentalist in recognizing specific interaction mechanisms. Results show that distinct multi-state models may produce very similar spectral patterns. I also discussed aggregation of a receptor as a possible source of spurious three-state line shapes and provided specific suggestions for complementary experiments that can ensure reliable mechanistic insight.

Molecular chaperones TAPBPR (TAP-binding protein related) and tapasin associate with major histocompatibility complex class I (MHC-I) to promote the loading of antigenic peptides through a poorly understood mechanism. Here, we use solution Nuclear Magnetic Resonance (NMR) spectroscopy to probe TAPBPR-induced changes on MHC-I. Dynamic motions present in the empty MHC groove become progressively dampened with increasing peptide occupancy, while allosteric communication between the A- and F-pockets regulates a conformational switch located near the TAPBPR binding site, which is crucial for chaperone release from the complex. Our analysis of NMR data recorded for a range of TAPBPR complexes prepared with both murine H2 and human HLA alleles complements the recent X-ray structures to provide atomic-resolution mechanistic insights into the selection of optimal peptide sequences for the displayed antigen repertoire. In particular, our results show that negative allosteric coupling between the MHC groove and chaperone binding sites allows TAPBPR to proofread MHC molecules containing a range of different peptides. Since the affinity of incoming peptides for the empty groove is greatly reduced in the chaperone complex, (micromolar range, relative to nanomolar for the free MHC), these interactions can provide a mechanism for optimizing the peptide repertoire, where only the highest-affinity peptides can drive chaperone release. Finally, our results suggest that TAPBPR may promote the dissociation of tightly bound peptides from MHC molecules, thereby further scrutinizing the displayed repertoire. These findings imply a similar mechanism for the specificity and editing function of tapasin in the peptide-loading complex.

The ability to use conformational flexibility is a hallmark of enzyme function. Here we show that protein motions and catalytic activity in a RNase are coupled and display identical solvent isotope effects. Solution NMR relaxation experiments identify a cluster of residues, some distant from the active site, that are integral to this motion. These studies implicate a single residue, histidine-48, as the key modulator in coupling protein motion with enzyme function. Mutation of H48 to alanine results in loss of protein motion in the isotope-sensitive region of the enzyme. In addition, k cat decreases for this mutant and the kinetic solvent isotope effect on k cat , which was 2.0 in WT, is near unity in H48A. Despite being located 18 Å from the enzyme active site, H48 is essential in coordinating the motions involved in the rate-limiting enzymatic step. These studies have identified, of ≈160 potential exchangeable protons, a single site that is integral in the rate-limiting step in RNase A enzyme function.

Two indicators of conformational variability of Escherichia coli acyl carrier protein (ACP) have been investigated, namely backbone dynamics and chemical shift variations of ACP. Hydrophobic interactions between the 4'-PP prosthetic group and the hydrophobic pocket enclosed by the amphipathic helices resulted in chemical shift perturbations in the residues near the prosthetic group binding sites and contact sites in the hydrophobic pockets upon conversion from apo- to holo-forms. At pH 7.9, destabilization of ACP due to negative charge repulsions and the deprotonated state of His 75 resulted in observed chemical shift changes in the C-terminal region. Model-free analysis showed that the alpha(1)alpha(2) loop region near the prosthetic group binding site in ACP shows the greatest flexibility (lowest S(2) values) and this result may suggest these flexibilities are required for structural rearrangements when the acyl chain binds to the prosthetic group of ACP. Flexibility of ACP shown in this study is essential for its ability to interact with functionally different enzyme partners specifically and weakly in the rapid delivery of acyl chain from one partner to another.

During neurotransmitter release, the neuronal SNARE proteins synaptobrevin/VAMP, syntaxin, and SNAP-25 form a four-helix bundle, the SNARE complex, that pulls the synaptic vesicle and plasma membranes together possibly causing membrane fusion. Complexin binds tightly to the SNARE complex and is essential for efficient Ca2+-evoked neurotransmitter release. A combined X-ray and TROSY-based NMR study now reveals the atomic structure of the complexin/SNARE complex. Complexin binds in an antiparallel α-helical conformation to the groove between the synaptobrevin and syntaxin helices. This interaction stabilizes the interface between these two helices, which bears the repulsive forces between the apposed membranes. These results suggest that complexin stabilizes the fully assembled SNARE complex as a key step that enables the exquisitely high speed of Ca2+-evoked neurotransmitter release.

The small GTPase Ras is an important signaling molecule acting as a molecular switch in eukaryotic cells. Recent findings of global conformational exchange and a putative allosteric binding site in the G domain of Ras opened an avenue to understanding novel aspects of Ras function. To facilitate detailed NMR studies of Ras in physiological solution conditions, we performed backbone resonance assignments of Ras bound to slowly hydrolysable GTP mimic, guanosine 5'-[s, γ-imido]triphosphate at pH 7.2. Out of 163 non-proline residues of the G domain, signals from backbone amide proton, nitrogen and carbon spins of 127 residues were confidently assigned with the remaining unassigned residues mostly located at the exchange-broadened effectors interface.

Using computer simulations we have studied possible effects of heating and cooling at different scan rates on unfolding and refolding of macromolecules. We have shown that even the simplest two-state reversible transition can behave irreversibly when an unfavorable combination of cooling rate, relaxation time and activation energy of refolding occurs. On the basis of this finding we suppose that apparent irreversibility of some proteins denatured by heat may result from slow relaxation on cooling rather than thermodynamic instability and/or irreversible alterations of the polypeptide chain. Using this kinetic reversible two-state model, we estimated the effects of the scan rate and kinetic parameters of the macromolecule on its unfolding–refolding process. A few recommendations are suggested on how to reach maximal possible recovery after denaturation if refolding appears to be under kinetic control.

Proteins exist in a delicate balance between the native and unfolded states, where thermodynamic stability may be sacrificed to attain the flexibility required for efficient catalysis, binding, or allosteric control. Partition-defective 6 (Par-6) regulates the Par polarity complex by transmitting a GTPase signal through the Cdc42/Rac interaction binding PSD-95/Dlg/ZO-1 (CRIB-PDZ) module that alters PDZ ligand binding. Allosteric activation of the PDZ is achieved by local rearrangement of the L164 and K165 side chains to stabilize the interdomain CRIB:PDZ interface and reposition a conserved element of the ligand binding pocket. However, microsecond to millisecond dynamics measurements revealed that L164/K165 exchange requires a larger rearrangement than expected. The margin of thermodynamic stability for the PDZ domain is modest (∼3 kcal/mol) and further reduced by transient interactions with the disordered CRIB domain. Measurements of local structural stability revealed that tertiary contacts within the PDZ are disrupted by a partial unfolding transition that enables interconversion of the L/K switch. The unexpected participation of partial PDZ unfolding in the allosteric mechanism of Par-6 suggests that native-state unfolding may be essential for the function of other marginally stable proteins.

International audience; Conformational flexibility between structural ensembles is an essential component of enzyme function. Although the broad dynamical landscape of proteins is known to promote a number of functional events on multiple time scales, it is yet unknown whether structural and functional enzyme homologues rely on the same concerted residue motions to perform their catalytic function. It is hypothesized that networks of contiguous and flexible residue motions occurring on the biologically relevant millisecond time scale evolved to promote and/or preserve optimal enzyme catalysis. In this study, we use a combination of NMR relaxation dispersion, model-free analysis, and ligand titration experiments to successfully capture and compare the role of conformational flexibility between two structural homologues of the pancreatic ribonuclease family: RNase A and eosinophil cationic protein (or RNase 3). In addition to conserving the same catalytic residues and structural fold, both homologues show similar yet functionally distinct clusters of millisecond dynamics, suggesting that conformational flexibility can be conserved among analogous protein folds displaying low sequence identity. Our work shows that the reduced conformational flexibility of eosinophil cationic protein can be dynamically and functionally reproduced in the RNase A scaffold upon creation of a chimeric hybrid between the two proteins. These results support the hypothesis that conformational flexibility is partly required for catalytic function in homologous enzyme folds, further highlighting the importance of dynamic residue sectors in the structural organization of proteins.

Peptidyl prolyl cis-trans isomerization acts as an effective molecular timer that plays significant roles in biological and pathological processes. Enzymes such as Pin1 catalyze cis-trans isomerization, accelerating the otherwise slow isomerization rate into time scales relevant for cellular signaling. Here we have combined NMR line shape analysis, fluorescence spectroscopy, and isothermal titration calorimetry to determine the kinetic and thermodynamic parameters describing the trans-specific interaction between the binding domain of Pin1 (WW domain) and a key cis-trans molecular switch in the amyloid precursor protein cytoplasmic tail. A three-state model, in which the cis-trans isomerization equilibrium is coupled to the binding equilibrium through the trans isomer, was found to fit the data well. The trans isomer binds the WW domain with ∼22 μM affinity via very fast association (approaching the diffusion limit) and dissociation rates. The common structural and electrostatic characteristics of Pin1 substrates, which contain a phosphorylated serine/threonine-proline motif, suggest that very rapid binding kinetics are a general feature of Pin1 interactions with other substrates. The fast binding kinetics of the WW domain allows rapid response of Pin1 to the dynamic events of phosphorylation and dephosphorylation in the cell that alter the relative populations of diverse Pin1 substrates. Furthermore, our results also highlight the vastly different rates at which slow uncatalyzed cis-trans isomerization and fast isomer-specific binding events occur. These results, along with the experimental methods presented herein, should guide future experiments aimed at the thermodynamic and kinetic characterization of cis-trans molecular switches and isomer-specific interactions involved in various biological processes.

We have recently discovered an allosteric switch in Ras, bringing an additional level of complexity to this GTPase whose mutants are involved in nearly 30% of cancers. Upon activation of the allosteric switch, there is a shift in helix 3/loop 7 associated with a disorder to order transition in the active site. Here, we use a combination of multiple solvent crystal structures and computational solvent mapping (FTMap) to determine binding site hot spots in the “off” and “on” allosteric states of the GTP-bound form of H-Ras. Thirteen sites are revealed, expanding possible target sites for ligand binding well beyond the active site. Comparison of FTMaps for the H and K isoforms reveals essentially identical hot spots. Furthermore, using NMR measurements of spin relaxation, we determined that K-Ras exhibits global conformational dynamics very similar to those we previously reported for H-Ras. We thus hypothesize that the global conformational rearrangement serves as a mechanism for allosteric coupling between the effector interface and remote hot spots in all Ras isoforms. At least with respect to the binding sites involving the G domain, H-Ras is an excellent model for K-Ras and probably N-Ras as well. Ras has so far been elusive as a target for drug design. The present work identifies various unexplored hot spots throughout the entire surface of Ras, extending the focus from the disordered active site to well-ordered locations that should be easier to target.

The phosphorylation-specific peptidyl-prolyl isomerase Pin1 catalyzes the isomerization of the peptide bond preceding a proline residue between cis and trans isomers. To best understand the mechanisms of Pin1 regulation, rigorous enzymatic assays of isomerization are required. However, most measures of isomerase activity require significant constraints on substrate sequence and only yield rate constants for the cis isomer, $$ k_{cat}^{cis} $$ and apparent Michaelis constants, $$ K_{M}^{App} $$ . By contrast, NMR lineshape analysis is a powerful tool for determining microscopic rates and populations of each state in a complex binding scheme. The isolated catalytic domain of Pin1 was employed as a first step towards elucidating the reaction scheme of the full-length enzyme. A 24-residue phosphopeptide derived from the amyloid precurser protein intracellular domain (AICD) phosphorylated at Thr668 served as a biologically-relevant Pin1 substrate. Specific 13C labeling at the Pin1-targeted proline residue provided multiple reporters sensitive to individual isomer binding and on-enzyme catalysis. We have performed titration experiments and employed lineshape analysis of phosphopeptide 13C–1H constant time HSQC spectra to determine $$ k_{cat}^{cis} $$ , $$ k_{cat}^{trans} $$ , $$ K_{D}^{cis} $$ , and $$ K_{D}^{trans} $$ for the catalytic domain of Pin1 acting on this AICD substrate. The on-enzyme equilibrium value of [E·trans]/[E·cis] = 3.9 suggests that the catalytic domain of Pin1 is optimized to operate on this substrate near equilibrium in the cellular context. This highlights the power of lineshape analysis for determining the microscopic parameters of enzyme catalysis, and demonstrates the feasibility of future studies of Pin1-PPIase mutants to gain insights on the catalytic mechanism of this important enzyme.

The small GTPase Ras is an important signaling molecule acting as a molecular switch in eukaryotic cells. Recent findings of global conformational exchange and a putative allosteric binding site in the G domain of Ras opened an avenue to understanding novel aspects of Ras function. To facilitate detailed NMR studies of Ras in physiological solution conditions, we performed backbone resonance assignments of Ras bound to slowly hydrolysable GTP mimic, guanosine 5'-[ß, γ-imido]triphosphate at pH 7.2. Out of 163 non-proline residues of the G domain, signals from backbone amide proton, nitrogen and carbon spins of 127 residues were confidently assigned with the remaining unassigned residues mostly located at the exchange-broadened effectors interface.

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.

Ras and its homologues are central to regulation of a multitude of cellular processes. Ras in complex with GTP binds and activates its downstream signaling partners. (31)P NMR studies indicated that the Ras-GTP conformation is heterogeneous on a millisecond time scale, but details of its conformational dynamics remain unknown. Here we present evidence that the conformational exchange process in human H-Ras complexed with GTP mimic GppNHp is global, encompassing most of the GTPase catalytic domain. The correlated character of conformational dynamics in Ras opens opportunities for understanding allosteric effects in Ras function.

Conformational flexibility of the enzyme architecture is essential for biological function. These structural transitions often encompass significant portions of the enzyme molecule. Here, we present a detailed study of functionally relevant RNase A dynamics in the wild type and a D121A mutant form by NMR spin-relaxation techniques. In the wild-type enzyme, the dynamic properties are largely conserved in the apo, enzyme−substrate, and enzyme−product complexes. In comparison, mutation of aspartic acid 121 to alanine disrupts the timing of active-site dynamics, the product-release step, and global conformational changes, indicating that D121 plays a significant role in coordinating the dynamic events in RNase A. In addition, this mutation results in 90% loss of catalytic activity despite the absence of direct participation of D121 in the chemical reaction or in interactions with the substrate. These data suggest that one role of this conserved residue is to facilitate important millisecond protein dynamics.

This work examines the robustness of fitting of parameters describing conformational exchange (k(ex), p(a/b), and Deltaomega) processes from CPMG relaxation dispersion data. We have analyzed the equations describing conformational exchange processes for the intrinsic inter-dependence of their parameters that leads to the existence of multiple equivalent solutions, which equally satisfy the experimental data. We have used Monte-Carlo simulations and fitting to the synthetic data sets as well as the direct 3-D mapping of the parameter space of k(ex), p(a/b), and Deltaomega to quantitatively assess the degree of the parameter inter-dependence. The demonstrated high correlation between parameters can preclude accurate dynamics parameter estimation from NMR spin-relaxation data obtained at a single static magnetic field. The strong parameter inter-dependence can readily be overcome through acquisition of spin-relaxation data at more than one static magnetic field thereby allowing accurate assessment of conformational exchange properties.

The interaction of the dinucleotide inhibitor 5'-phosphothymidine(3',5')pyrophosphate adenosine 3'-phosphate (pTppAp) with bovine pancreatic ribonuclease A (RNase A) was characterized by calorimetry and solution NMR spectroscopy. Calorimetric data show that binding of pTppAp to RNase A is exothermic (DeltaH = -60.1 +/- 4.1 kJ/mol) with a dissociation constant of 16 nM at 298 K. At this temperature, the binding results in an entropy loss (TDeltaS = -16.8 +/- 7.3 kJ/mol) that is more favorable than that with the product analogue, 2'-CMP (TDeltaS = -31.3 +/- 0.9 kJ/mol). Temperature-dependent calorimetric experiments give a DeltaC(p) for ligand binding of -230 +/- 100 J/mol K. Binding of pTppAp results in noticeable effects on the backbone amide chemical shifts and dynamics. Amide backbone (15)N NMR spin-relaxation studies were performed on both apo RNase A and RNase A/pTppAp as a function of temperature. At each temperature, the model-free-determined order parameters, S(2), were significantly higher for RNase A/pTppAp than for the apo enzyme indicating a decrease in the conformational entropy of the protein upon ligand binding. Furthermore, the magnitude of this difference varies along the amino acid sequence specifically locating the entropic changes. The temperature dependence of S(2) at each residue enabled assessment of the local heat capacity changes (DeltaC(p)) from ligand binding. In an overall, average sense, DeltaC(p) for the protein backbone, determined from the NMR dynamics measurements, did not differ between apo RNase A and RNase A/pTppAp indicating that backbone dynamics contribute little to DeltaC(p) for protein-ligand interactions in this system. However, residue-by-residue comparison of the temperature-dependent change in entropy (DeltaS(B)) between free and bound forms reveals nonzero contributions to DeltaC(p) at individual sites. The balance of positive and negative changes reveals a redistribution of energetics upon binding. Furthermore, experiment and semiempirical estimates suggest that a large negative DeltaC(p) should accompany binding of pTppAp, and we conclude that this contribution must arise from factors other than amide backbone dynamics.

The Fluorescence2D is free software that allows analysis of two-dimensional fluorescence spectra obtained using the accelerated “triangular” acquisition schemes. The software is a combination of Python and MATLAB-based programs that perform conversion of the triangular data, display of the two-dimensional spectra, extraction of 1D slices at different wavelengths, and output in various graphic formats.

Stability of hen lysozyme in the presence of acetonitrile (MeCN) at different pH values of the medium was studied by scanning microcalorimetry with a special emphasis on determination of reliable values of the denaturational heat capacity change. It was found that the temperature of denaturation decreases on addition of MeCN. However, the free energy extrapolation showed that below room temperature the thermodynamic stability increases at low concentrations of MeCN in spite of the general destabilizing effect at higher concentrations and temperatures. Charge-induced contribution to this stabilization was shown to be negligible (no pH-dependence was found); therefore, the most probable cause for the phenomenon is an increase of hydrophobic interactions at low temperatures in aqueous solutions containing small amounts of the organic additive. The difference in preferential solvation of native and denatured states of lysozyme was calculated from the stabilization free energy data. It was found that the change in preferential solvation strongly depends on the temperature in the water-rich region. At the higher MeCN content this dependence decreases until, at 0.06 mole fractions of MeCN, the difference in the preferential solvation between native and denatured lysozyme becomes independent of the temperature over a range of 60 K. The importance of taking into account non-ideality of a mixed solution, when analyzing preferential solvation phenomena was emphasized.

Through characterization of the solvent isotope effect on protein dynamics, we have examined determinants of the rate limitation to enzyme catalysis. A global conformational change in Ribonuclease A limits the overall rate of catalytic turnover. Here we show that this motion is sensitive to solvent deuterium content; the isotope effect is 2.2, a value equivalent to the isotope effect on the catalytic rate constant. We further demonstrate that the protein motion possesses a linear proton inventory plot, indicating that a single proton is transferred in the transition state. These results provide compelling evidence for close coupling between enzyme dynamics and function and demonstrate that characterization of the transition state for protein motion in atomic detail is experimentally accessible.

On the basis of scanning microcalorimetry data from literature and our own measurements, we have calculated the changes in preferential solvation of lysozyme upon heat denaturation in six solvent systems: water + methanol, ethanol, propanol [data from Velicelebi, G.,Sturtevant, J. M. (1979) Biochemistry 18, 1180], acetone, p-dioxane [data from Fujita, Y.,Noda, Y. (1983) Bull. Chem. Soc.Jpn. 56, 233], and dimethylsulfoxide [our data Kovrigin, E. L., Kirkitadze, M. D.,Potekhin, S. A. (1996) Biofizika 41, 549-553; Kovrigin, E. L.,Potekhin, S. A. (1996) Biofizika 41, 1201-1206]. These preferential solvation changes are (in effect) the numbers of cosolvent molecules entering or leaving the solvation shell of the protein upon denaturation. It has been shown that for a group of five substances in the initial activity range (approximately up to 0. 3) the denaturational changes of preferential solvation of lysozyme do not depend on the nature of the solvent and depend only on its activity. This suggests that lysozyme does not distinguish these substances in the initial activity range and preferential solvation has a nonspecific character. It has been shown also that preferential solvation DeltaGamma23 does not depend on the pH value at least for dimethylsulfoxide-water solutions. This indicates that the chargeable groups exposed on denaturation do not contribute significantly to preferential interaction of the protein surface with the solution components.

Ras GTPase is a molecular switch controlling a number of cellular pathways including growth, proliferation, differentiation, and apoptosis. Recent reports indicated that Ras undergoes dimerization at the membrane surface through protein-protein interactions. If firmly established this property of Ras would require profound reassessment of a large amount of published data and modification of the Ras signaling paradigm. One proposed mechanism of dimerization involves formation of salt bridges between the two GTPase domains (G domains) leading to formation of a compact dimer as observed in Ras crystal structures. In this work, we interrogated the intrinsic ability of Ras to self-associate in solution by creating conditions of high local concentration through irreversibly tethering the two G domains together at their unstructured C-terminal tails. We evaluated possible self-association in this inverted tandem conjugate via analysis of the time-domain fluorescence anisotropy and NMR chemical shift perturbations. We did not observe the increased rotational correlation time expected for the G domain dimer. Variation of the ionic strength (to modulate stability of the salt bridges) did not affect the rotational correlation time in the tandem further supporting independent rotational diffusion of two G domains. In a parallel line of experiments to detect and map weak self-association of the G domains, we analyzed NMR chemical shifts perturbations at a number of sites near the crystallographic dimer interface. The nearly complete lack of chemical shift perturbations in the tandem construct supported a simple model with the independent G domains repelled from each other by their overall negative charge. These results lead us to the conclusion that self-association of the G domains cannot be responsible for homodimerization of Ras reported in the literature.