Structural dynamics of the membrane translocation domain of colicin E9 and its interaction with TolB

In order for the 61 kDa colicin E9 protein toxin to enter the cytoplasm of susceptible cells and kill them by hydrolysing their DNA, the colicin must interact with the outer membrane BtuB receptor and Tol translocation pathway of target cells. The translocation function is located in the N-terminal domain of the colicin molecule. (1)H, (1)H-(1)H-(15)N and (1)H-(13)C-(15)N NMR studies of intact colicin E9, its DNase domain, minimal receptor-binding domain and two N-terminal constructs containing the translocation domain showed that the region of the translocation domain that governs the interaction of colicin E9 with TolB is largely unstructured and highly flexible. Of the expected 80 backbone NH resonances of the first 83 residues of intact colicin E9, 61 were identified, with 43 of them being assigned specifically. The absence of secondary structure for these was shown through chemical shift analyses and the lack of long-range NOEs in (1)H-(1)H-(15)N NOESY spectra (tau(m)=200 ms). The enhanced flexibility of the region of the translocation domain containing the TolB box compared to the overall tumbling rate of the protein was identified from the relatively large values of backbone and tryptophan indole (15)N spin-spin relaxation times, and from the negative (1)H-(15)N NOEs of the backbone NH resonances. Variable flexibility of the N-terminal region was revealed by the (15)N T(1)/T(2) ratios, which showed that the C-terminal end of the TolB box and the region immediately following it was motionally constrained compared to other parts of the N terminus. This, together with the observation of inter-residue NOEs involving Ile54, indicated that there was some structural ordering, resulting most probably from the interactions of side-chains. Conformational heterogeneity of parts of the translocation domain was evident from a multiplicity of signals for some of the residues. Im9 binding to colicin E9 had no effect on the chemical shifts or other NMR characteristics of the region of colicin E9 containing the TolB recognition sequence, though the interaction of TolB with intact colicin E9 bound to Im9 did affect resonances from this region. The flexibility of the translocation domain of colicin E9 may be connected with its need to recognise protein partners that assist it in crossing the outer membrane and in the translocation event itself.     

Taking into account the conformational dynamics of macromolecules when analyzing luminescent probe signals

A mathematical model is proposed that takes into account the specific features of electronic excitation energy transfer between the molecular probes sorbed on a polymeric chain in a liquid solution. The kinetics of the process is described from the standpoint of the stochastic changes in macromolecular conformations. The results of numerical solution of the problem are described, as well as the analytical expressions obtained in the context of perturbation theory for a low transfer rate and/or fast conformational dynamics of the macrochain. A channel of nonlinear (with respect to excitation) deactivation resulting from binary annihilation of closely spaced excited centers is considered. Expressions for the efficient rate of mutual quenching and delayed annihilation fluorescence of the probe are obtained. The time dependences of characteristic luminescent signals and the parametric curves of relative fluorescence quantum yield are given.     

NMR investigation of the dynamics of tryptophan side-chains in hemoglobins

NMR relaxation measurements of 15N spin-lattice relaxation rate (R(1)), spin-spin relaxation rate (R(2)), and heteronuclear nuclear Overhauser effect (NOE) have been carried out at 11.7T and 14.1T as a function of temperature for the side-chains of the tryptophan residues of 15N-labeled and/or (2H,15N)-labeled recombinant human normal adult hemoglobin (Hb A) and three recombinant mutant hemoglobins, rHb Kempsey (betaD99N), rHb (alphaY42D/betaD99N), and rHb (alphaV96W), in the carbonmonoxy and the deoxy forms as well as in the presence and in the absence of an allosteric effector, inositol hexaphosphate (IHP). There are three Trp residues (alpha14, beta15, and beta37) in Hb A for each alphabeta dimer. These Trp residues are located in important regions of the Hb molecule, i.e. alpha14Trp and beta15Trp are located in the alpha(1)beta(1) subunit interface and beta37Trp is located in the alpha(1)beta(2) subunit interface. The relaxation experiments show that amino acid substitutions in the alpha(1)beta(2) subunit interface can alter the dynamics of beta37Trp. The transverse relaxation rate (R(2)) for beta37Trp can serve as a marker for the dynamics of the alpha(1)beta(2) subunit interface. The relaxation parameters of deoxy-rHb Kemspey (betaD99N), which is a naturally occurring abnormal human hemoglobin with high oxygen affinity and very low cooperativity, are quite different from those of deoxy-Hb A, even in the presence of IHP. The relaxation parameters for rHb (alphaY42D/betaD99N), which is a compensatory mutant of rHb Kempsey, are more similar to those of Hb A. In addition, TROSY-CPMG experiments have been used to investigate conformational exchange in the Trp residues of Hb A and the three mutant rHbs. Experimental results indicate that the side-chain of beta37Trp is involved in a relatively slow conformational exchange on the micro- to millisecond time-scale under certain experimental conditions. The present results provide new dynamic insights into the structure-function relationship in hemoglobin.     

Structure and conformational dynamics of base excision repair DNA glycosylases

DNA glycosylases play a key role in DNA repair, which maintains the integrity of the cell genome. The structures of many DNA glycosylases have been solved to date. The review considers these structures and the dynamics of DNA glycosylase interactions with DNA. The available data suggest that lesion recognition by DNA glycosylases is a highly dynamic process that is accompanied by multiple conformational changes in the enzyme and DNA substrate.     

The use of TROSY for detection and suppression of conformational exchange NMR line broadening in biological macromolecules

The interference between conformational exchange-induced time-dependent variations of chemical shifts in a pair of scalar coupled ¹H and ¹⁵N spins is used to construct novel TROSY-type NMR experiments to suppress NMR signal loss in [¹⁵N,¹H]-correlation spectra of a 14-mer DNA duplex free in solution and complexed with the Antp homeodomain. An analysis of double- and zero-quantum relaxation rates of base ¹H–¹⁵N moieties showed that for certain residues the contribution of conformational exchange-induced transverse relaxation might represent a dominant relaxation mechanism, which, in turn, can be effectively suppressed by TROSY. The use of the new TROSY method for exchange-induced transverse relaxation optimization is illustrated with two new experiments, 2D ^h1 J _HN,^h2 J _NN-quantitative [¹⁵N,¹H]-TROSY to measure ^h1 J _HN and ^h2 J _NN scalar coupling constants across hydrogen bonds in nucleic acids, and 2D (^h2 J _NN+^h1 J _NH)-correlation-[¹⁵N,¹H]-TROSY to correlate ¹H^N chemical shifts of bases with the chemical shifts of the tertiary ¹⁵N spins across hydrogen bonds using the sum of the trans-hydrogen bond coupling constants in nucleic acids.     

Torsion angle dynamics: Reduced variable conformational sampling enhances crystallographic structure refinement

A reduced variable conformational sampling strategy for macromolecules based on molecular dynamics in torsion angle space is evaluated using crystallographic refinement as a prototypical search problem. Bae and Haug's algorithm for constrained dynamics [Bae, D. S., Haug, E. J. A recursive formulation for constrained mechanical system dynamics. Mech. Struct. Mach. 15:359–382, 1987], originally developed for robotics, was used. Their formulation solves the equations of motion exactly for arbitrary holonomic constraints, and hence differs from commonly used approximation algorithms. It uses gradients calculated in Cartesian coordinates, and thus also differs from internal coordinate formulations. Molecular dynamics can be carried out at significantly higher temperatures due to the elimination of the high frequency bond and angle vibrations. The sampling strategy presented here combines high temperature torsion angle dynamics with repeated trajectories using different initial velocities. The best solutions can be identified by the free R value, or the R value if experimental phase information is appropriately included in the refinement. Applications to crystallographic refinement show a significantly increased radius of convergence over conventional techniques. For a test system with diffraction data to 2 Å resolution, slow-cooling protocols fail to converge if the backbone atom root mean square (rms) coordinate deviation from the crystal structure is greater than 1.25 Å, but torsion angle refinement can correct backbone atom rms coordinate deviations up to approximately 1.7 Å. © 1994 Wiley-Liss, Inc.     

Alternate-site isotopic labeling of ribonucleotides for NMR studies of ribose conformational dynamics in RNA

Heteronuclear NMR spin relaxation studies of conformational dynamics are coming into increasing use to help understand the functions of ribozymes and other RNAs. Due to strong magnetic interactions within the ribose ring, however, these studies have thus far largely been limited to ¹³C and ¹⁵N resonances on the nucleotide base side chains. We report here the application of the alternate-site ¹³C isotopic labeling scheme, pioneered by LeMaster for relaxation studies of amino acid side chains, to nucleic acid systems. We have used different strains of E. coli to prepare mononucleotides containing ¹³C label in one of two patterns: Either C1′ or C2′ in addition to C4′, termed (1′/2′,4′) labeling, or nearly complete labeling at the C2′ and C4′ sites only, termed (2′,4′) labeling. These patterns provide isolated H spin systems on the labeled carbon atoms and thus allow spin relaxation studies without interference from scalar or dipolar coupling. Using relaxation studies of AMP dissolved in glycerol at varying temperature to produce systems with correlation times characteristic of different size RNAs, we demonstrate the removal of errors due to interaction in T ₁ measurements of larger nucleic acids and in T _1ρ measurements in RNA molecules. By extending the applicability of spin relaxation measurements to backbone ribose groups, this technology should greatly improve the flexibility and completeness of NMR analyses of conformational dynamics in RNA.     

Intramolecular dynamics of low molecular weight protein tyrosine phosphatase in monomer-dimer equilibrium studied by NMR: a model for changes in dynamics upon target binding

Low molecular weight protein tyrosine phosphatase (LMW-PTP) dimerizes in the phosphate-bound state in solution with a dissociation constant of K(d)=1.5(+/-0.1)mM and an off-rate on the order of 10(4)s(-1). 1H and 15N NMR chemical shifts identify the dimer interface, which is in excellent agreement with that observed in the crystal structure of the dimeric S19A mutant. Two tyrosine residues of each molecule interact with the active site of the other molecule, implying that the dimer may be taken as a model for a complex between LMW-PTP and a target protein. 15N relaxation rates for the monomeric and dimeric states were extrapolated from relaxation data acquired at four different protein concentrations. Relaxation data of satisfactory precision were extracted for the monomer, enabling model-free analyses of backbone fluctuations on pico- to nanosecond time scales. The dimer relaxation data are of lower quality due to extrapolation errors and the possible presence of higher-order oligomers at higher concentrations. A qualitative comparison of order parameters in the monomeric and apparent dimeric states shows that loops forming the dimer interface become rigidified upon dimerization. Qualitative information on monomer-dimer exchange and intramolecular conformational exchange was obtained from the concentration dependence of auto- and cross-correlated relaxation rates. The loop containing the catalytically important Asp129 fluctuates between different conformations in both the monomeric and dimeric (target bound) states. The exchange rate compares rather well with that of the catalyzed reaction step, supporting existing hypotheses that catalysis and enzyme dynamics may be coupled. The side-chain of Trp49, which is important for substrate specificity, exhibits conformational dynamics in the monomer that are largely quenched upon formation of the dimer, suggesting that binding is associated with the selection of a single side-chain conformer.     

Nudged elastic band calculation of minimal energy paths for the conformational change of a GG non-canonical pair

The nudged elastic band (NEB) technique has been implemented in AMBER to calculate low-energy paths for conformational changes. A novel simulated annealing protocol that does not require an initial hypothesis for the path is used to sample low-energy paths. This was used to study the conformational change of an RNA cis Watson-Crick/Hoogsteen GG non-canonical pair, with one G syn around the glycosidic bond and the other anti. A previous solution structure, determined by NMR-constrained modeling, demonstrated that the GG pairs change from (syn)G-(anti)G to (anti)G-(syn)G in the context of duplex r(GCAGGCGUGC) on the millisecond timescale. The set of low-energy paths found by NEB show that each G flips independently around the glycosidic bond, with the anti G flipping to syn first. Guanine bases flip without opening adjacent base-pairs by protruding into the major groove, accommodated by a transient change by the ribose to C2'-exo sugar pucker. Hydrogen bonds between bases and the backbone, which lower the energetic barrier to flipping, are observed along the path. The results show the plasticity of RNA base-pairs in helices, which is important for biological processes, including mismatch repair, protein recognition, and translation. The modeling of the GG conformational change also demonstrates that NEB can be used to discover non-trivial paths for macromolecules and therefore NEB can be used as an exploratory method for predicting putative conformational change paths.     

Effects of calcium binding on the side-chain methyl dynamics of calbindin D9k: a 2H NMR relaxation study

The effects of Ca(2+) binding on the side-chain methyl dynamics of calbindin D(9k) have been characterized by (2)H NMR relaxation rate measurements. Longitudinal, transverse in-phase, quadrupolar order, transverse anti-phase and double quantum relaxation rates are reported for both the apo and Ca(2+)-loaded states of the protein at two magnetic field strengths. The relatively large size of the data set allows for a detailed analysis of the underlying conformational dynamics by spectral density mapping and model-free fitting procedures. The results reveal a correlation between a methyl group's distance from the Ca(2+) binding sites and its conformational dynamics. Several methyl groups segregate into two limiting classes, one proximal and the other distal to the binding sites. Methyl groups in these two classes respond differently to Ca(2+) binding, both in terms of the timescale and amplitude of their fluctuations. Ca(2+) binding elicits a partial immobilization among methyl groups in the proximal class, which is consistent with previous studies of calbindin's backbone dynamics. The distal class, however, exhibits a trend that could not be inferred from the backbone data in that its mobility actually increases with Ca(2+) binding. We have introduced the term polar dynamics to describe this type of organization across the molecule. The trend may represent an important mechanism by which calbindin D(9k) achieves high affinity binding while minimizing the corresponding loss of conformational entropy.     

Ligand-induced changes in the conformational dynamics of a bacterial cytotoxic endonuclease

Knowledge about the conformational dynamics of a protein is key to understanding its biochemical and biophysical properties. In the present work we investigated the dynamic properties of the enzymatic domain of DNase colicins via time-resolved fluorescence and anisotropy decay analysis in combination with steady-state acrylamide quenching experiments. The dynamic properties of the apoenzyme were compared to those of the E9 DNase ligated to the transition metal ion Zn(2+) and the natural inhibitor Im9. We further investigated the contributions of each of the two tryptophans within the E9 DNase (Trp22 and Trp58) using two single-tryptophan mutants (E9 W22F and E9 W58F). Wild-type E9 DNase, E9 W22F, and E9 W58F, as well as Im9, showed multiple lifetime decays. The time-resolved and steady-state fluorescence results indicated that complexation of E9 DNase with Zn(2+) induces compaction of the E9 DNase structure, accompanied by immobilization of Trp22 along with a reduced solvent accessibility for both tryptophans. Im9 binding resulted in immobilization of Trp22 along with a decrease in the longest lifetime component. In contrast, Trp58 experienced less restriction on complexation of E9 DNase with Im9 and showed an increase in the longest lifetime component. Furthermore, the results point out that the Im9-induced changes in the conformational dynamics of E9 DNase are predominant and occur independently of the Zn(2+)-induced conformational effects.     

Light-induced conformational changes in full-length Arabidopsis thaliana cryptochrome

Cryptochromes (CRYs) are widespread flavoproteins with homology to photolyases (PHRs), a class of blue-light-activated DNA repair enzymes. Unlike PHRs, both plant and animal CRYs have a C-terminal domain. This cryptochrome C-terminal (CCT) domain mediates interactions with other proteins, while the PHR-like domain converts light energy into a signal via reduction and radical formation of the flavin adenine dinucleotide cofactor. However, the mechanism by which the PHR-like domain regulates the CCT domain is not known. Here, we applied the pulsed-laser-induced transient grating method to detect conformational changes induced by blue-light excitation of full-length Arabidopsis thaliana cryptochrome 1 (AtCRY1). A significant reduction in the diffusion coefficient of AtCRY1 was observed upon photoexcitation, indicating that a large conformational change occurs in this monomeric protein. AtCRY1 containing a single mutation (W324F) that abolishes an intra-protein electron transfer cascade did not exhibit this conformational change. Moreover, the conformational change was much reduced in protein lacking the CCT domain. Thus, we conclude that the observed large conformational changes triggered by light excitation of the PHR-like domain result from C-terminal domain rearrangement. This inter-domain modulation would be critical for CRYs' ability to transduce a blue-light signal into altered protein-protein interactions for biological activity. Lastly, we demonstrate that the transient grating technique provides a powerful method for the direct observation and understanding of photoreceptor dynamics.     

Investigations into sequence and conformational dependence of backbone entropy,inter-basin dynamics and the Flory isolated-pair hypothesis for peptides

The populations and transitions between Ramachandran basins are studied for combinations of the standard 20 amino acids in monomers, dimers and trimers using an implicit solvent Langevin dynamics algorithm and employing seven commonly used force-fields. Both the basin populations and inter-conversion rates are influenced by the nearest neighbor's conformation and identity, contrary to the Flory isolated-pair hypothesis. This conclusion is robust to the choice of force-field, even though the use of different force-fields produces large variations in the populations and inter-conversion rates between the dominant helical, extended beta, and polyproline II basins. The computed variation of conformational and dynamical properties with different force-fields exceeds the difference between explicit and implicit solvent calculations using the same force-field. For all force-fields, the inter-basin transitions exhibit a directional dependence, with most transitions going through extended beta conformation, even when it is the least populated basin. The implications of these results are discussed in the context of estimates for the backbone entropy of single residues, and for the ability of all-atom simulations to reproduce experimental protein folding data.     

Characterization of slow conformational dynamics in solids: dipolar CODEX

A solid state NMR experiment is introduced for probing relatively slow conformational exchange, based on dephasing and refocusing dipolar couplings. The method is closely related to the previously described Centerband-Only Detection of Exchange or CODEX experiment. The use of dipolar couplings for this application is advantageous because their values are known a priori from molecular structures, and their orientations and reorientations relate in a simple way to molecular geometry and motion. Furthermore the use of dipolar couplings in conjunction with selective isotopic enrichment schemes is consistent with selection for unique sites in complex biopolymers. We used this experiment to probe the correlation time for the motion of ¹³C, ¹⁵N enriched urea molecules within their crystalline lattice.     

Conformational analysis of an EP24.15 inhibitor by NMR and molecular modelling

Summary The enzyme thimet oligopeptidase (EC3.4.24.15, EP24.15) is responsible for the hydrolysis of a number of neuropeptides. Despite much research examining its substrate specificity, little is known about the conformational requirements of its active site. We have used 1D¹H and 2D TOCSY NMR experiments to assign the proton resonances of the EP24.15 inhibitor,N-[1-(R, S)-carboxy-3-phenylpropyl]-Ala-Ala-Tyr-p-aminobenzoate (cFP), and 2D ROESY NMR to investigate whether cFP exhibits any conformational preferences in CD₃OD and in aqueous CD₃OD. Molecular modelling of charged cFP in the gaseous phase generated a number of conformations that were consistent with the NMR data obtained in CD₃OD. Analogous modelling on the uncharged cFP did not result in conformations consistent with any of the NMR data, but did suggest that, under non-polar conditions, cFP could adopt a hairpin conformation which would allow simultaneous coordination of the two carboxyl groups of cFP to the zinc ion in the active site of EP24.15.     

Characterisation of a mobile protein-binding epitope in the translocation domain of colicin E9

The 61 kDa colicin E9 protein toxin enters the cytoplasm of susceptible cells by interacting with outer membrane and periplasmic helper proteins, and kills them by hydrolysing their DNA. The membrane translocation function is located in the N-terminal domain of the colicin, with a key signal sequence being a pentapeptide region that governs the interaction with the helper protein TolB (the TolB box). Previous NMR studies (Collins et al., 2002 J. Mol. Biol. 318, 787-804) have shown that the N-terminal 83 residues of colicin E9, which includes the TolB box, is largely unstructured and highly flexible. In order to further define the properties of this region we have studied a fusion protein containing residues 1-61 of colicin E9 connected to the N-terminus of the E9 DNase by an eight-residue linking sequence. 53 of the expected 58 backbone NH resonances for the first 61 residues and all of the expected 7 backbone NH resonances of the linking sequence were assigned with 3D (1)H-(13)C-(15)N NMR experiments, and the backbone dynamics of these regions investigated through measurement of (1)H-(15)N relaxation properties. Reduced spectral density mapping, extended Lipari-Szabo modelling, and fitting backbone R(2) relaxation rates to a polymer dynamics model identifies three clusters of interacting residues, each containing a tryptophan. Each of these clusters is perturbed by TolB binding to the intact colicin, showing that the significant region for TolB binding extends beyond the recognized five amino acids of the TolB box and demonstrating that the binding epitope for TolB involves a considerable degree of order within an otherwise disordered and flexible domain. Abbreviations : Im9, the immunity protein for colicin E9; E9 DNase, the endonuclease domain of colicin E9; HSQC, heteronuclear single quantum coherence; ppm, parts per million; DSS, 2,2-(dimethylsilyl)propanesulfonic acid; TSP, sodium 3-trimethylsilypropionate; T(1 - 61)-DNase fusion protein, residues 1-61 of colicin E9 connected to the N-terminus of the E9 DNase by an eight residue thrombin cleavage sequence.     

A conformational study in solution of pro-somatostatin fragments by NMR and computational methods

The results of a conformational study by nuclear magnetic spectroscopy and computational methods on a series of point-mutated synthetic peptides, containing 14 amino acid residues and mimicking the region containing the Arg-Lys dibasic cleavage site of pro-somatostatin, have confirmed the possible role of a well defined secondary structure in the recognition phenomenon by processing enzymes. The importance of the residues located near the Arg-Lys dibasic site in the C-terminal region of the pro-hormone for the cleavage of the precursor into somatostatin-14 has been confirmed. The present structural analysis indicates the occurrence of two β-turns in the 4–7 and 11–14 regions, flanking the cleavage site, for all the peptides recognized as substrates by the processing enzyme. Interestingly, in the point-mutated analogue not processed by the enzyme and containing the replacement of proline by alanine in position 5 the first β-turn is displaced by one residue and involves the Ala⁵-Arg⁸ segment. This observation may explain the lack of recognition by the maturation enzyme. © 1998 European Peptide Society and John Wiley & Sons, Ltd.