NMR line shapes and multi-state binding equilibria

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.     

Fluorescence2D: Software for Accelerated Acquisition and Analysis of Two-Dimensional Fluorescence Spectra

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.     

Dielectric parameters of ascitic and pleural fluids in the microwave range in different nosologies

The dielectric parameters of ascitic and pleural fluids formed in the human body in oncological and nononcological diseases of different nosology have been estimated in the range between 400 MHz and 1.2 GHZ. The dependence of refractive and absorption indices of ascitic and pleural liquids on the signal frequency and mass concentration of dissolved substances was found. Common regularities and distinctions in the behavior of their dielectric properties were revealed.     

Temperature dependence of the backbone dynamics of ribonuclease A in the ground state and bound to the inhibitor 5'-phosphothymidine (3'-5')pyrophosphate adenosine 3'-phosphate

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.     

In Vitro Biosynthesis and Chemical Identification of UDP-N-acetyl-d-quinovosamine (UDP-d-QuiNAc)

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 His₆-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.     

Conservation of Flexible Residue Clusters among Structural and Functional Enzyme Homologues

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.     

Folding under inequilibrium conditions as a possible reason for partial irreversibility of heat-denatured proteins: computer simulation study

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.     

Application of methyl-TROSY to a large paramagnetic membrane protein without perdeuteration: <Superscript>13</Superscript>C-MMTS-labeled NADPH-cytochrome P450 oxidoreductase

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 ¹³C-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 ¹³C-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.     

Three-dimensional structure of the complexin/SNARE complex

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 Ca(2+)-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 alpha-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 Ca(2+)-evoked neurotransmitter release.     

Assignments of backbone 1H, 13C and 15N resonances in H-Ras (1–166) complexed with GppNHp at physiological pH

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.