Rotational Diffusion Anisotropy and Local Backbone Dynamics of Carbon Monoxide-Bound Rhodobacter capsulatus Cytochrome c‘

Changes in dynamic and static structures of the HIV-1 p24 capsid protein N-domain caused by amino-acid substitution are associated with its viral viability

HIV-1 capsid is comprised of over a hundred p24 protein molecules, arranged as either pentamers or hexamers. Three p24 mutants with amino acid substitutions in capsid N-terminal domain protein were examined: G60W (α3-4 loop), M68T (helix 4), and P90T (α4-5 loop), which exhibited no viability for biological activity. One common structural feature of the three p24 N-domain mutants, examined by NMR, was the long-range effect of more β-structures at the β2-strand in the N-terminal region compared with the wild-type. In addition, the presence of fewer helical structures was observed in M68T and P90T, beyond the broad area from helix 1 to the C-terminal part of helix 4. This suggests that both N-terminal beta structures and helices play important roles in the formation of p24 hexamers and pentamers. Next, compared with P90T, we examined cis-conformation or trans-conformation of wild-type adopted by isomerization at G89-P90. Since P90T mutant adopts only a trans-conformation, comparison of chemical shifts and signal intensities between each spectra revealed that the major peaks (about 85%) in the spectrum of wild-type correspond to trans-conformation. Furthermore, it was indicated that the region in cis-conformation (minor; 15%) was more stabilized than that observed in trans-conformation, based on the analyses of heteronuclear Overhauser effect as well as the order-parameter. Therefore, it was concluded that the cis-conformation is more favorable than the trans-conformation for the interaction between the p24 N-terminal domain and cyclophilin-A. This is because HIV-1 with a P90T protein, which adopts only a trans-conformation, is associated with non-viability of biological activity.

The PipX Protein, When Not Bound to Its Targets, Has Its Signaling C-Terminal Helix in a Flexed Conformation

PipX, an 89-residue protein, acts as a coactivator of the global nitrogen regulator NtcA in cyanobacteria. NtcA-PipX interactions are regulated by 2-oxoglutarate (2-OG), an inverse indicator of the ammonia abundance, and by P_II, a protein that binds to PipX at low 2-OG concentrations. The structure of PipX, when bound to NtcA or P_II, consists of an N-terminal, five-stranded β-sheet (conforming a Tudor-like domain), and two long α-helices. These helices adopt either a flexed conformation, where they are in close contact and in an antiparallel mutual orientation, also packing against the β-sheet, or an open conformation (observed only in the P_II-PipX complex) where the last α-helix moves apart from the rest of the protein. The aim of this work was to study the structure and dynamics of isolated PipX in solution by NMR. The backbone chemical shifts, the hydrogen-exchange, and the NOE patterns indicated that the isolated, monomeric PipX structure was formed by an N-terminal five-stranded β-sheet and two C-terminal α-helices. Furthermore, the observed NOEs between the two helices, and of α-helix2 with β-strand2 suggested that PipX adopted a flexed conformation. The β-strands 1 and 5 were highly flexible, as shown by the lack of interstrand backbone-backbone NOEs; in addition, the ¹⁵N-dynamics indicated that the C terminus of β-strand4 and the following β-turn (Phe42-Thr47), and the C-cap of α-helix1 (Arg70-Asn71) were particularly mobile. These two regions could act as hinges, allowing PipX to interact with its partners, including PlmA in the newly recognized P_II-PipX-PlmA ternary complex.

Effect of Glu12-His89 Interaction on Dynamic Structures in HIV-1 p17 Matrix Protein Elucidated by NMR

To test the existence of the salt bridge and stability of the HIV-1 p17 matrix protein, an E12A (mutated at helix 1) was established to abolish possible electrostatic interactions. The chemical shift perturbation from the comparison between wild type and E12A suggested the existence of an electrostatic interaction in wild type between E12 and H89 (located in helix 4). Unexpectedly, the studies using urea denaturation indicated that the E12A substitution slightly stabilized the protein. The dynamic structure of E12A was examined under physiological conditions by both amide proton exchange and relaxation studies. The quick exchange method of amide protons revealed that the residues with faster exchange were located at the mutated region, around A12, compared to those of the wild-type protein. In addition, some residues at the region of helix 4, including H89, exhibited faster exchange in the mutant. In contrast, the average values of the kinetic rate constants for amide proton exchange for residues located in all loop regions were slightly lower in E12A than in wild type. Furthermore, the analyses of the order parameter revealed that less flexible structures existed at each loop region in E12A. Interestingly, the structures of the regions including the alpha1-2 loop and helix 5 of E12A exhibited more significant conformational exchanges with the NMR time-scale than those of wild type. Under lower pH conditions, for further destabilization, the helix 1 and alpha2-3 loop in E12A became more fluctuating than at physiological pH. Because the E12A mutant lacks the activities for trimer formation on the basis of the analytical ultra-centrifuge studies on the sedimentation distribution of p17 (Fledderman et al. Biochemistry 49, 9551–9562, 2010), it is possible that the changes in the dynamic structures induced by the absence of the E12-H89 interaction in the p17 matrix protein contributes to a loss of virus assembly.

A Multivalent Marine Lectin from Crenomytilus grayanus Possesses Anti-cancer Activity through Recognizing Globotriose Gb3

In this study, we report the structure and function of a lectin from the sea mollusk Crenomytilus grayanus collected from the sublittoral zone of Peter the Great Bay of the Sea of Japan. The crystal structure of C. grayanus lectin (CGL) was solved to a resolution of 1.08 Å, revealing a β-trefoil fold that dimerizes into a dumbbell-shaped quaternary structure. Analysis of the crystal CGL structures bound to galactose, galactosamine, and globotriose Gb3 indicated that each CGL can bind three ligands through a carbohydrate-binding motif involving an extensive histidine- and water-mediated hydrogen bond network. CGL binding to Gb3 is further enhanced by additional side-chain-mediated hydrogen bonds in each of the three ligand-binding sites. NMR titrations revealed that the three binding sites have distinct microscopic affinities toward galactose and galactosamine. Cell viability assays showed that CGL recognizes Gb3 on the surface of breast cancer cells, leading to cell death. Our findings suggest the use of this lectin in cancer diagnosis and treatment.

Dynamics of nitric oxide synthase-calmodulin interactions at physiological calcium concentrations

The intracellular Ca²⁺ concentration is an important regulator of many cellular functions. The small acidic protein calmodulin (CaM) serves as a Ca²⁺ sensor and control element for many enzymes. Nitric oxide synthase (NOS) is one of the proteins that is activated by CaM and plays a major role in a number of key physiological and pathological processes. Previous studies have shown CaM to act like a switch that causes a conformational change in NOS to allow for the electron transfer between the reductase and oxygenase domains through a process that is thought to be highly dynamic. We have analyzed the structure and dynamics of complexes formed by peptides based on inducible NOS (iNOS) and endothelial NOS (eNOS) with CaM at Ca²⁺ concentrations that mimic the physiological basal (17 and 100 nM) and elevated levels (225 nM) found in mammalian cells using fluorescence techniques and nuclear magnetic resonance spectroscopy. The results show the CaM-NOS complexes have similar structures at physiological and fully saturated Ca²⁺ levels; however, their dynamics are remarkably different. At 225 nM Ca²⁺, the CaM-NOS complexes show overall an increase in backbone dynamics, when compared to the dynamics of the complexes at saturating Ca²⁺ concentrations. Specifically, the N-lobe of CaM in the CaM-iNOS complex displays a lower internal mobility (higher S²) and higher exchange protection compared to those of the CaM-eNOS complex. In contrast, the C-lobe of CaM in the CaM-eNOS complex is less dynamic. These results illustrate that structures of CaM-NOS complexes determined at saturated Ca²⁺ concentrations cannot provide a complete picture because the differences in intramolecular dynamics become visible only at physiological Ca²⁺ levels.

Cytochromes c': Structure, Reactivity and Relevance to Haem-Based Gas Sensing

Cytochromes c' are a group of class IIa cytochromes with pentacoordinate haem centres and are found in photosynthetic, denitrifying and methanotrophic bacteria. Their function remains unclear, although roles in nitric oxide (NO) trafficking during denitrification or in cellular defence against nitrosoative stress have been proposed. Cytochromes c' are typically dimeric with each c-type haem-containing monomer folding as a four-α-helix bundle. Their hydrophobic and crowded distal sites impose severe restrictions on the binding of distal ligands, including diatomic gases. By contrast, NO binds to the proximal haem face in a similar manner to that of the eukaryotic NO sensor, soluble guanylate cyclase and bacterial analogues. In this review, we focus on how structural features of cytochromes c' influence haem spectroscopy and reactivity with NO, CO and O2. We also discuss the relevance of cytochrome c' to understanding the mechanisms of gas binding to haem-based sensor proteins.

Flexible and rigid structures in HIV-1 p17 matrix protein monitored by relaxation and amide proton exchange with NMR

The HIV-1 p17 matrix protein is a multifunctional protein that interacts with other molecules including proteins and membranes. The dynamic structure between its folded and partially unfolded states can be critical for the recognition of interacting molecules. One of the most important roles of the p17 matrix protein is its localization to the plasma membrane with the Gag polyprotein. The myristyl group attached to the N-terminus on the p17 matrix protein functions as an anchor for binding to the plasma membrane. Biochemical studies revealed that two regions are important for its function: D14-L31 and V84-V88. Here, the dynamic structures of the p17 matrix protein were studied using NMR for relaxation and amide proton exchange experiments at the physiological pH of 7.0. The results revealed that the α12-loop, which includes the 14-31 region, was relatively flexible, and that helix 4, including the 84-88 region, was the most protected helix in this protein. However, the residues in the α34-loop near helix 4 had a low order parameter and high exchange rate of amide protons, indicating high flexibility. This region is probably flexible because this loop functions as a hinge for optimizing the interactions between helices 3 and 4. The C-terminal long region of K113-Y132 adopted a disordered structure. Furthermore, the C-terminal helix 5 appeared to be slightly destabilized due to the flexible C-terminal tail based on the order parameters. Thus, the dynamic structure of the p17 matrix protein may be related to its multiple functions.

Structural characterization of CHCHD5 and CHCHD7: two atypical human twin CX9C proteins

Twin CX(9)C proteins constitute a large protein family among all eukaryotes; are putative substrates of the mitochondrial Mia40-dependent import machinery; contain a coiled coil-helix-coiled coil-helix (CHCH) fold stabilized by two disulfide bonds as exemplified by three structures available for this family. However, they considerably differ at the primary sequence level and this prevents an accurate prediction of their structural models. With the aim of expanding structural information on CHCH proteins, here we structurally characterized human CHCHD5 and CHCHD7. While CHCHD5 has two weakly interacting CHCH domains which sample a range of limited conformations as a consequence of hydrophobic interactions, CHCHD7 has a third helix hydrophobically interacting with an extension of helix α2, which is part of the CHCH domain. Upon reduction of the disulfide bonds both proteins become unstructured exposing hydrophobic patches, with the result of protein aggregation/precipitation. These results suggest a model where the molecular interactions guiding the protein recognition between Mia40 and the disulfide-reduced CHCHD5 and CHCHD7 substrates occurs in vivo when the latter proteins are partially embedded in the protein import pore of the outer membrane of mitochondria.

Promiscuous binding at the crossroads of numerous cancer pathways: insight from the binding of glutaminase interacting protein with glutaminase L

The glutaminase interacting protein (GIP) is composed of a single PDZ domain that interacts with a growing list of partner proteins, including glutaminase L, that are involved in a number of cell signaling and cancer pathways. Therefore, GIP makes a good target for structure-based drug design. Here, we report the solution structures of both free GIP and GIP bound to the C-terminal peptide analogue of glutaminase L. This is the first reported nuclear magnetic resonance structure of GIP in a complex with one of its binding partners. Our analysis of both free GIP and GIP in a complex with the glutaminase L peptide provides important insights into how a promiscuous binding domain can have affinity for multiple binding partners. Through a detailed chemical shift perturbation analysis and backbone dynamics studies, we demonstrate here that the binding of the glutaminase L peptide to GIP is an allosteric event. Taken together, the insights reported here lay the groundwork for the future development of a specific inhibitor for GIP.

Nuclear magnetic resonance provides a quantitative description of protein conformational flexibility on physiologically important time scales

A complete description of biomolecular activity requires an understanding of the nature and the role of protein conformational dynamics. In recent years, novel nuclear magnetic resonance-based techniques that provide hitherto inaccessible detail concerning biomolecular motions occurring on physiologically important time scales have emerged. Residual dipolar couplings (RDCs) provide precise information about time- and ensemble-averaged structural and dynamic processes with correlation times up to the millisecond and thereby encode key information for understanding biological activity. In this review, we present the application of two very different approaches to the quantitative description of protein motion using RDCs. The first is purely analytical, describing backbone dynamics in terms of diffusive motions of each peptide plane, using extensive statistical analysis to validate the proposed dynamic modes. The second is based on restraint-free accelerated molecular dynamics simulation, providing statistically sampled free energy-weighted ensembles that describe conformational fluctuations occurring on time scales from pico- to milliseconds, at atomic resolution. Remarkably, the results from these two approaches converge closely in terms of distribution and absolute amplitude of motions, suggesting that this kind of combination of analytical and numerical models is now capable of providing a unified description of protein conformational dynamics in solution.

Optimisation of NMR dynamic models I. Minimisation algorithms and their performance within the model-free and Brownian rotational diffusion spaces

The key to obtaining the model-free description of the dynamics of a macromolecule is the optimisation of the model-free and Brownian rotational diffusion parameters using the collected R (1), R (2) and steady-state NOE relaxation data. The problem of optimising the chi-squared value is often assumed to be trivial, however, the long chain of dependencies required for its calculation complicates the model-free chi-squared space. Convolutions are induced by the Lorentzian form of the spectral density functions, the linear recombinations of certain spectral density values to obtain the relaxation rates, the calculation of the NOE using the ratio of two of these rates, and finally the quadratic form of the chi-squared equation itself. Two major topological features of the model-free space complicate optimisation. The first is a long, shallow valley which commences at infinite correlation times and gradually approaches the minimum. The most severe convolution occurs for motions on two timescales in which the minimum is often located at the end of a long, deep, curved tunnel or multidimensional valley through the space. A large number of optimisation algorithms will be investigated and their performance compared to determine which techniques are suitable for use in model-free analysis. Local optimisation algorithms will be shown to be sufficient for minimisation not only within the model-free space but also for the minimisation of the Brownian rotational diffusion tensor. In addition the performance of the programs Modelfree and Dasha are investigated. A number of model-free optimisation failures were identified: the inability to slide along the limits, the singular matrix failure of the Levenberg-Marquardt minimisation algorithm, the low precision of both programs, and a bug in Modelfree. Significantly, the singular matrix failure of the Levenberg-Marquardt algorithm occurs when internal correlation times are undefined and is greatly amplified in model-free analysis by both the grid search and constraint algorithms. The program relax ( http://www.nmr-relax.com ) is also presented as a new software package designed for the analysis of macromolecular dynamics through the use of NMR relaxation data and which alleviates all of the problems inherent within model-free analysis.

Set theory formulation of the model-free problem and the diffusion seeded model-free paradigm

Model-free analysis of NMR relaxation data, which describes the motion of individual atoms, is a problem intricately linked to the Brownian rotational diffusion of the macromolecule. The diffusion tensor parameters strongly influence the optimisation of the various model-free models and the subsequent model selection between them. Finding the optimal model of the dynamics of the system among the numerous diffusion and model-free models is hence quite complex. Using set theory, the entirety of this global problem has been encapsulated by the universal set Ll, and its resolution mathematically formulated as the universal solution Ll. Ever since the original Lipari and Szabo papers the model-free dynamics of a molecule has most often been solved by initially estimating the diffusion tensor. The model-free models which depend on the diffusion parameter values are then optimised and the best model is chosen to represent the dynamics of the residue. Finally, the global model of all diffusion and model-free parameters is optimised. These steps are repeated until convergence. For simplicity this approach to Ll will be labelled the diffusion seeded model-free paradigm. Although this technique suffers from a number of problems many have been solved. All aspects of the diffusion seeded paradigm and its consequences, together with a few alternatives to the paradigm, will be reviewed through the use of set notation.

Heterologous Overexpression and Purification of Cytochrome c‘ from Rhodobacter capsulatus and a Mutant (K42E) in the Dimerization Region. Mutation Does Not Alter Oligomerization but Impacts the Heme Iron Spin State and Nitric Oxide Binding Properties

Rhodobacter capsulatus cytochrome c' (RCCP) has been overexpressed in Escherichia coli, and its spectroscopic and ligand-binding properties have been investigated. It is concluded that the heterologously expressed protein is assembled correctly, as judged by UV-vis absorption, EPR, and resonance Raman (RR) spectroscopy of the unligated protein as well as forms in which the heme is ligated by CO or NO. To probe the oligomerization state of RCCP and its potential influence on heme reactivity, we have compared the properties of wild-type RCCP with a mutant (K42E) that lacks a salt bridge at the subunit interface. Analytical ultracentrifugation indicates that wild-type and K42E proteins are both monomeric in solution, contrary to the homodimeric structure of the crystalline state. Surprisingly, the K42E mutation produces a number of changes at the heme center (nearly 20 A distant), including perturbation of the ferric spin-state equilibrium and a change in the ferrous heme-nitrosyl complex from a six-coordinate/five-coordinate mixture to a predominantly five-coordinate heme-NO species. RR spectra indicate that ferrous K42E and wild-type RCCP both have relatively high Fe-His stretching frequencies, suggesting that the more favored five-coordinate heme-nitrosyl formation in K42E is not caused by a weaker Fe2+-His bond. Nevertheless, the altered reactivity of ferrous K42E with NO, together with its modified ferric spin state, shows that structural changes originating at the dimer interface can affect the properties of the heme center, raising the exciting possibility that intermolecular encounters at the protein surface might modulate the reactivity of cytochrome c' in vivo.

A new spin probe of protein dynamics: nitrogen relaxation in 15N-2H amide groups

(15)N spin relaxation data have provided a wealth of information on protein dynamics in solution. Standard R(1), R(1)(rho), and NOE experiments aimed at (15)N[(1)H] amide moieties are complemented in this work by HA(CACO)N-type experiments allowing the measurement of nitrogen R(1) and R(1)(rho) rates at deuterated (15)N[(2)D] sites. Difference rates obtained using this approach, R(1)((15)N[(1)H]) - R(1)((15)N[(2)D]) and R(2)((15)N[(1)H]) - R(2)((15)N[(2)D]), depend exclusively on dipolar interactions and are insensitive to (15)N CSA and R(ex) relaxation mechanisms. The methodology has been tested on a sample of peptostreptococcal protein L (63 residues) prepared in 50% H(2)O-50% D(2)O solvent. The results from the new and conventional experiments are found to be consistent, with respect to both local backbone dynamics and overall protein tumbling. Combining several data sets permits evaluation of the spectral density J(omega(D) + omega(N)) for each amide site. This spectral density samples a uniquely low frequency (26 MHz at a 500 MHz field) and, therefore, is expected to be highly useful for characterizing nanosecond time scale local motions. The spectral density mapping demonstrates that, in the case of protein L, J(omega(D) + omega(N)) values are compatible with the Lipari-Szabo interpretation of backbone dynamics based on the conventional (15)N relaxation data.

Efficient and accurate determination of the overall rotational diffusion tensor of a molecule from (15)N relaxation data using computer program ROTDIF

We present a computer program ROTDIF for efficient determination of a complete rotational diffusion tensor of a molecule from NMR relaxation data. The derivation of the rotational diffusion tensor in the case of a fully anisotropic model is based on a six-dimensional search, which could be very time consuming, particularly if a grid search in the Euler angle space is involved. Here, we use an efficient Levenberg-Marquardt algorithm combined with Monte Carlo generation of initial guesses. The result is a dramatic, up to 50-fold improvement in the computational efficiency over the previous approaches. This method is demonstrated on a computer-generated and real protein systems. We also address the issue of sensitivity of the diffusion tensor determination from (15)N relaxation measurements to experimental errors in the relaxation rates and discuss possible artifacts from applying higher-symmetry tensor model and how to recognize them.

Solution Structures of a Cyanobacterial Metallochaperone

The Atx1 copper metallochaperone from Synechocystis PCC 6803, ScAtx1, interacts with two P(1)-type copper ATPases to supply copper proteins within intracellular compartments, avoiding ATPases for other metals en route. Here we report NMR-derived solution structures for ScAtx1. The monomeric apo form has a betaalphabetabetaalpha fold with backbone motions largely restricted to loop 1 containing Cys-12 and Cys-15. The tumbling rate of Cu(I)ScAtx1 (0.1-0.8 mm) implies dimers. Experimental restraints are satisfied by symmetrical dimers with Cys-12 or His-61, but not Cys-15, invading the copper site of the opposing subunit. A full sequence of copper ligands from the cell surface to thylakoid compartments is proposed, considering in vitro homodimer liganding to mimic in vivo liganding in ScAtx1-ATPase heterodimers. A monomeric high resolution structure for Cu(I)ScAtx1, with Cys-12, Cys-15, and His-61 as ligands, is calculated without violations despite the rotational correlation time. (2)J(NH) couplings in the imidazole ring of His-61 establish coordination of N(epsilon2) to copper. His-61 is analogous to Lys-65 in eukaryotic metallochaperones, stabilizing Cu(I)S(2) complexes but by binding Cu(I) rather than compensating charge. Cys-Cys-His ligand sets are an emergent theme in some copper metallochaperones, although not in related Atx1, CopZ, or Hah1. Surface charge (Glu-13) close to the metal-binding site of ScAtx1 is likely to support interaction with complementary surfaces of copper-transporting ATPases (PacS-Arg-11 and CtaA-Lys-14) but to discourage interaction with zinc ATPase ZiaA and so inhibit aberrant formation of copper-ZiaA complexes.

Direct structure determination using residual dipolar couplings: reaction-site conformation of methionine sulfoxide reductase in solution

Residual dipolar couplings (RDC) from partially aligned molecules provide long-range structural data and are thus particularly well adapted to rapid structure validation or protein fold recognition. Extensive measurements in two alignment media can also provide precise de novo structure from RDC alone. We have applied a novel combination of these approaches to the study of methionine sulfoxide reductase (MsrA) from Erwinia chrysanthemi, a 27 kDa enzyme essential for repairing oxidative stress damage. The tertiary fold was initially validated by comparing backbone RDC to expected values from the crystal structure of the homologous MsrA from Escherichia coli. Good agreement was found throughout the chain, verifying the overall topology of the molecule, with the exception of the catalytically important peptide P196-L202, where strong and systematic RDC violation was observed. No evidence for local differential mobility in this region was detected, implying that the structure of the strand differs in the two molecules. We have therefore applied the de novo approach meccano to determine the conformation of this peptide using only RDC. A single conformation is found that is in agreement with all measured data. The aligned peptide can be docked onto the expected covalence of the rest of the template molecule while respecting its strictly defined relative orientation. In contrast to the structure of MsrA from E. coli, the reactive side chain of Cys200 is oriented toward the interior of the molecule and therefore closer to the catalytic Cys53, obviating the need for previously proposed conformational reorganization prior to formation of this disulfide intermediate. This analysis requires only backbone assignment and uses unambiguously assigned and readily measurable structural data, thereby greatly economizing investigation time compared to established nuclear Overhauser effect- (nOe-) based structure calculation methods.

Temperature-dependent dynamics of the villin headpiece helical subdomain, an unusually small thermostable protein

(15)N spin relaxation experiments were used to measure the temperature-dependence of protein backbone conformational fluctuations in the thermostable helical subdomain, HP36, of the F-actin-binding headpiece domain of chicken villin. HP36 is the smallest domain of a naturally occurring protein that folds cooperatively to a compact native state. Spin-lattice, spin-spin, and heteronuclear nuclear Overhauser effect relaxation data for backbone amide (15)N spins were collected at five temperatures in the range of 275-305 K. The data were analyzed using a model-free formalism to determine generalized order parameters, S, that describe the distribution of N-H bond vector orientations in a molecular reference frame. A novel parameter, Lambda=dln(1-S)/dln T is introduced to characterize the temperature-dependence of S. An average value of Lambda=4.5 is obtained for residues in helical conformations in HP36. This value of Lambda is not reproduced by model potential energy functions commonly used to parameterize S. The maximum entropy principle was used to derive a new model potential function that reproduces both S and Lambda. Contributions to the entropy, S(r), and heat capacity, C(r)(p), from reorientational conformational fluctuations were analyzed using this potential energy function. Values of S(r) show a qualitative dependence on S similar to that obtained for the diffusion-in-a-cone model; however, quantitative differences of up to 0.5k, in which k is the Boltzmann constant, are observed. Values of C(r)(p) approach zero for small values of S and approach k for large values of S; the largest values of C(r)(p) are predicted to occur for intermediate values of S. The results suggest that backbone dynamics, as probed by relaxation measurements, make very little contribution to the heat capacity difference between folded and unfolded states for HP36.

Base flexibility in HIV-2 TAR RNA mapped by solution (15)N, (13)C NMR relaxation

Binding of the HIV tat protein to the TAR (transactivating response region) RNA element activates transcription of the HIV viral genome. The complex of TAR with argininamide serves as a model for the RNA conformation in the tat-TAR complex. The dynamics of the HIV-2 TAR-argininamide complex was investigated by measurements of the relaxation rates of protonated base carbon and nitrogen nuclei. Six auto-correlation rates as well as cross-correlation rates were measured to map the frequencies of base motion in the HIV-2 TAR-argininamide complex. These measurements reveal a broad range of dynamic heterogeneity exhibited by hexanucleotide loop, the dinucleotide bulge, and the A-form helical regions. U23 in the bulge undergoes the largest dynamic change on binding argininamide, while U25 remains flexible, reflecting the large conformational change that is triggered upon ligand binding.

Macromolecular import into Escherichia coli: the TolA C-terminal domain changes conformation when interacting with the colicin A toxin

Various macromolecules such as bacteriotoxins and phage DNA parasitize some envelope proteins of Escherichia coli to infect the bacteria. A two-step import mechanism involves the primary interaction with an outer membrane receptor or with a pilus followed by the translocation across the outer membrane. However, this second step is poorly understood. It was shown that the TolA, TolQ, and TolR proteins play a critical role in the translocation of group A colicins and filamentous bacteriophage minor coat proteins (g3p). Translocation of these proteins requires the interaction of their N-terminal domain with the C-terminal domain of TolA (TolAIII). In this work, short soluble TolAIII domains were overproduced in the cytoplasm and in the periplasm of E. coli. In TolAIII, the two cysteine residues were found to be reduced in the cytoplasmic form and oxidized in the periplasmic form. The interaction of TolAIII with the N-terminal domain of colicin A (ATh) is observed in the presence and in the absence of the disulfide bridge. The complex formation of TolAIII and ATh was found to be independent of the ionic strength. An NMR study of TolAIII, both free and bound, shows a significant structural change when interacting with ATh, in the presence or absence of the disulfide bridge. In contrast, such a structural modification was not observed when TolAIII interacts with g3p N1. These results suggest that bacteriotoxins and Ff bacteriophages parasitize E. coli using different interactions between TolA and the translocation domain of the colicin and g3p protein, respectively.

(15)N backbone dynamics of ferricytochrome b(562): comparison with the reduced protein and the R98C variant

The backbone dynamics of ferricytochrome b(562), a four-helix bundle protein from Escherichia coli, have been studied by NMR spectroscopy. The consequences of the introduction of a c-type thioether linkage between the heme and protein and the reduction to the ferrous cytochrome have also been analyzed. (15)N relaxation rates R(1) and R(2) and (1)H-(15)N NOEs were measured at proton Larmor frequencies of 500 and 600 MHz for the oxidized and reduced protein as well as for the oxidized R98C variant. In the latter protein, an "artificial" thioether covalent bond has been introduced between the heme group and the protein frame [Arnesano, F., Banci, L., Bertini, I., Ciofi-Baffoni, S., de Lumley Woodyear, T., Johnson, C. M., and Barker, P. D. (2000) Biochemistry 39, 1499-1514]. The (15)N relaxation data were analyzed with the ModelFree protocol, and the mobility parameters on the picosecond to nanosecond time scale were compared for the three species. The three forms are rather rigid as a whole, with average generalized order parameters values of 0.87 +/- 0.08 (oxidized cytochrome b(562)), 0.84 +/- 0.07 (reduced cytochrome b(562)), and 0.85 +/- 0.07 (oxidized R98C cytochrome b(562)), indicating similar mobility for each system. Lower order parameters (S(2)) are found for residues belonging to loops 1 and 2. Higher mobility, as indicated by lower order parameters, is found for heme binding helices alpha 1 and alpha 4 in the R98C variant with respect to the wild-type protein. The analysis requires a relatively long rotational correlation time (tau(m) = 9.6 ns) whose value is accounted for on the basis of the anisotropy of the molecular shape and the high phosphate concentration needed to ensure the occurrence of monomer species. A parallel study of motions in the millisecond to microsecond time scale has also been performed on oxidized wild-type and R98C cytochrome b(562). In a CPMG experiment, decay rates were analyzed in the presence of spin-echo pulse trains of variable spacing. The dynamic behavior on this time scale is similar to that observed on the sub-nanosecond time scale, showing an increased mobility in the residues connected to the heme ligands in the R98C variant. It appears that the increased protein stability of the variant, established previously, is not correlated with an increase in rigidity.

NMR structure and backbone dynamics of a concatemer of epidermal growth factor homology modules of the human low-density lipoprotein receptor

The ligand-binding region of the low-density lipoprotein (LDL) receptor is formed by seven N-terminal, imperfect, cysteine-rich (LB) modules. This segment is followed by an epidermal growth factor precursor homology domain with two N-terminal, tandem, EGF-like modules that are thought to participate in LDL binding and recycling of the endocytosed receptor to the cell surface. EGF-A and the concatemer, EGF-AB, of these modules were expressed in Escherichia coli. Correct protein folding of EGF-A and the concatemer EGF-AB was achieved in the presence or absence of calcium ions, in contrast to the LB modules, which require them for correct folding. Homonuclear and heteronuclear 1H-15N NMR spectroscopy at 17.6 T was used to determine the three-dimensional structure of the concatemer. Both modules are formed by two pairs of short, anti-parallel beta-strands. In the concatemer, these modules have a fixed relative orientation, stabilized by calcium ion-binding and hydrophobic interactions at the interface. 15N longitudinal and transverse relaxation rates, and [1H]-15N heteronuclear NOEs were used to derive a model-free description of the backbone dynamics of the molecule. The concatemer appears relatively rigid, particularly near the calcium ion-binding site at the module interface, with an average generalized order parameter of 0.85+/-0.11. Some mutations causing familial hypercholesterolemia may now be rationalized. Mutations of D41, D43 and E44 in the EGF-B calcium ion-binding region may affect the stability of the linker and thus the orientation of the tandem modules. The diminutive core also provides little structural stabilization, necessitating the presence of disulfide bonds. The structure and dynamics of EGF-AB contrast with the N-terminal LB modules, which require calcium ions both for folding to form the correct disulfide connectivities and for maintenance of the folded structure, and are connected by highly mobile linking peptides.

Magnetic susceptibility tensor and heme contact shifts determinations in the Rhodobacter capsulatus ferricytochrome c': NMR and magnetic susceptibility studies

The 1H and 15N resonances of the carbon monoxide complex of ferrocytochrome c' of Rhodobacter capsulatus, a ferrous diamagnetic heme protein, have been extensively assigned by TOCSY-HSQC, NOESY-HSQC, and HSQC-NOESY-HSQC 3D heteronuclear experiments performed on a 7 mM sample labeled with 15N. Based on short-range and medium-range NOEs and H(N) exchange rates, the secondary structure consists of four helices: helix 1 (3-29), helix 2 (33-48), helix 3 (78-101), and helix 4 (103-125). The 15N, 1HN, and 1H(alpha) chemical shifts of the CO complex form are compared to those of the previously assigned oxidized (or ferric) state. From the chemical shift differences between these redox states, the orientation and the anisotropy of the paramagnetic susceptibility tensor have been determined using the crystallographic coordinates of the ferric state. The chi-tensor is axial, and the orientation of the z-axis is approximately perpendicular to the heme plane. The paramagnetic chemical shifts of the protons of the heme ligand have been determined and decomposed into the Fermi shift and dipolar shift contributions. Magnetic susceptibility studies in frozen solutions have been performed. Fits of the susceptibility data using the model of Maltempo (Maltempo, M. M. J. Chem. Phys. 1974, 61, 2540-2547) are consistent with a rather low contribution of the S = 3/2 spin state over the range of temperatures and confirm the value of the axial anisotropy. Values in the range 10.4-12.5 cm(-1) have been inferred for the axial zero-field splitting parameter (D). Analysis of the contact shift and the susceptibility data suggests that cytochrome c' of Rb. capsulatus exhibits a predominant high-spin character of the iron in the oxidized state at room temperature.