New methods of structure refinement for macromolecular structure determination by NMR

Recent advances in multidimensional NMR methodology have permitted solution structures of proteins in excess of 250 residues to be solved. In this paper, we discuss several methods of structure refinement that promise to increase the accuracy of macromolecular structures determined by NMR. These methods include the use of a conformational database potential and direct refinement against three-bond coupling constants, secondary 13C shifts, 1H shifts, T1/T2 ratios, and residual dipolar couplings. The latter two measurements provide long range restraints that are not accessible by other solution NMR parameters.

Functional dynamics in the active site of the ribonuclease binase

Binase, a member of a family of microbial guanyl-specific ribonucleases, catalyzes the endonucleotic cleavage of single-stranded RNA. It shares 82% amino acid identity with the well-studied protein barnase. We used NMR spectroscopy to study the millisecond dynamics of this small enzyme, using several methods including the measurement of residual dipolar couplings in solution. Our data show that the active site of binase is flanked by loops that are flexible at the 300-μs time scale. One of the catalytic residues, His-101, is located on such a flexible loop. In contrast, the other catalytic residue, Glu-72, is located on a β-sheet, and is static. The residues Phe-55, part of the guanine base recognition site, and Tyr-102, stabilizing the base, are the most dynamic. Our findings suggest that binase possesses an active site that has a well-defined bottom, but which has sides that are flexible to facilitate substrate access/egress, and to deliver one of the catalytic residues. The motion in these loops does not change on complexation with the inhibitor d(CGAG) and compares well with the maximum kcat (1,500 s−1) of these ribonucleases. This observation indicates that the NMR-measured loop motions reflect the opening necessary for product release, which is apparently rate limiting for the overall turnover.

Excitonic couplings and electronic coherence in bridged naphthalene dimers

The electronic excitations of naphthalene and a family of bridged naphthalene dimers are calculated and analyzed by using the Collective Electronic Oscillator method combined with the oblique Lanczos algorithm. All experimentally observed trends in absorption profiles and radiative lifetimes are reproduced. Each electronic excitation is linked to the corresponding real-space transition density matrix, which represents the motions of electrons and holes created in the molecule by photon absorption. Two-dimensional plots of these matrices help visualize the degree of exciton localization and explain the dependence of the electronic interaction between chromophores on their separation.

NMR scalar couplings across Watson–Crick base pair hydrogen bonds in DNA observed by transverse relaxation-optimized spectroscopy

This paper describes the NMR observation of 15N—15N and 1H—15N scalar couplings across the hydrogen bonds in Watson–Crick base pairs in a DNA duplex, hJNN and hJHN. These couplings represent new parameters of interest for both structural studies of DNA and theoretical investigations into the nature of the hydrogen bonds. Two dimensional [15N,1H]-transverse relaxation-optimized spectroscopy (TROSY) with a 15N-labeled 14-mer DNA duplex was used to measure hJNN, which is in the range 6–7 Hz, and the two-dimensional hJNN-correlation-[15N,1H]-TROSY experiment was used to correlate the chemical shifts of pairs of hydrogen bond-related 15N spins and to observe, for the first time, hJHN scalar couplings, with values in the range 2–3.6 Hz. TROSY-based studies of scalar couplings across hydrogen bonds should be applicable for large molecular sizes, including protein-bound nucleic acids.

Residual delay maps unveil global patterns of atmospheric nonlinearity and produce improved local forecasts

We use residual-delay maps of observational field data for barometric pressure to demonstrate the structure of latitudinal gradients in nonlinearity in the atmosphere. Nonlinearity is weak and largely lacking in tropical and subtropical sites and increases rapidly into the temperate regions where the time series also appear to be much noisier. The degree of nonlinearity closely follows the meridional variation of midlatitude storm track frequency. We extract the specific functional form of this nonlinearity, a V shape in the lagged residuals that appears to be a basic feature of midlatitude synoptic weather systems associated with frontal passages. We present evidence that this form arises from the relative time scales of high-pressure versus low-pressure events. Finally, we show that this nonlinear feature is weaker in a well regarded numerical forecast model (European Centre for Medium-Range Forecasts) because small-scale temporal and spatial variation is smoothed out in the grided inputs. This is significant, in that it allows us to demonstrate how application of statistical corrections based on the residual-delay map may provide marked increases in local forecast accuracy, especially for severe weather systems.

Elimination of residual metastatic prostate cancer after surgery and adjunctive cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) blockade immunotherapy

Cancer relapse after surgery is a common occurrence, most frequently resulting from the outgrowth of minimal residual disease in the form of metastases. We examined the effectiveness of cytotoxic T lymphocyte-associated antigen 4 (CTLA-4) blockade as an adjunctive immunotherapy to reduce metastatic relapse after primary prostate tumor resection. For these studies, we developed a murine model in which overt metastatic outgrowth of TRAMP-C2 (C2) prostate cancer ensues after complete primary tumor resection. Metastatic relapse in this model occurs reliably and principally within the draining lymph nodes in close proximity to the primary tumor, arising from established metastases present at the time of surgery. Using this model, we demonstrate that adjunctive CTLA-4 blockade administered immediately after primary tumor resection reduces metastatic relapse from 97.4 to 44%. Consistent with this, lymph nodes obtained 2 weeks after treatment reveal marked destruction or complete elimination of C2 metastases in 60% of mice receiving adjunctive anti-CTLA-4 whereas 100% of control antibody-treated mice demonstrate progressive C2 lymph node replacement. Our study demonstrates the potential of adjunctive CTLA-4 blockade immunotherapy to reduce cancer relapse emanating from minimal residual metastatic disease and may have broader implications for improving the capability of immunotherapy by combining such forms of therapy with other cytoreductive measures including surgery.

Antigen-driven CD4+ T cell and HIV-1 dynamics: Residual viral replication under highly active antiretroviral therapy

Antigen-induced stimulation of the immune system can generate heterogeneity in CD4+ T cell division rates capable of explaining the temporal patterns seen in the decay of HIV-1 plasma RNA levels during highly active antiretroviral therapy. Posttreatment increases in peripheral CD4+ T cell counts are consistent with a mathematical model in which host cell redistribution between lymph nodes and peripheral blood is a function of viral burden. Model fits to patient data suggest that, although therapy reduces HIV replication below replacement levels, substantial residual replication continues. This residual replication has important consequences for long-term therapy and the evolution of drug resistance and represents a challenge for future treatment strategies.

Changes in brain cell shape create residual extracellular space volume and explain tortuosity behavior during osmotic challenge

Diffusion of molecules in brain extracellular space is constrained by two macroscopic parameters, tortuosity factor λ and volume fraction α. Recent studies in brain slices show that when osmolarity is reduced, λ increases while α decreases. In contrast, with increased osmolarity, α increases, but λ attains a plateau. Using homogenization theory and a variety of lattice models, we found that the plateau behavior of λ can be explained if the shape of brain cells changes nonuniformly during the shrinking or swelling induced by osmotic challenge. The nonuniform cellular shrinkage creates residual extracellular space that temporarily traps diffusing molecules, thus impeding the macroscopic diffusion. The paper also discusses the definition of tortuosity and its independence of the measurement frame of reference.

NMR investigation of ferricytochrome c unfolding: Detection of an equilibrium unfolding intermediate and residual structure in the denatured state

Horse ferricytochrome c (cyt c) undergoes exchange of one of its axial heme ligands (Met-80) for one or more non-native ligands under denaturing conditions. We have used 1H NMR spectroscopy to detect two conformations of paramagnetic cyt c with non-native heme ligation through a range of urea concentrations. One non-native form is an equilibrium unfolding intermediate observed under partially denaturing conditions and is attributed to replacement of Met-80 with one or more Lys side chains. The second non-native form, in which the native Met ligand is replaced by a His, is observed under strongly denaturing conditions. Thermodynamic analysis of these data indicates a relatively small ΔG (17 kJ/mol) for the transition from native to the Lys-ligated intermediate and a significantly larger ΔG (47 kJ/mol) for the transition from native to the His-ligated species. Although CD and fluorescence data indicate that the equilibrium unfolding of cyt c is a two-state process, these NMR results implicate an intermediate with His-Lys ligation.

Characterization of residual structure in the thermally denatured state of barnase by simulation and experiment: Description of the folding pathway

Residual structure in the denatured state of a protein may contain clues about the early events in folding. We have simulated by molecular dynamics the denatured state of barnase, which has been studied by NMR spectroscopy. An ensemble of 104 structures was generated after 2 ns of unfolding and following for a further 2 ns. The ensemble was heterogeneous, but there was nonrandom, residual structure with persistent interactions. Helical structure in the C-terminal portion of helix α1 (residues 13–17) and in helix α2 as well as a turn and nonnative hydrophobic clustering between β3 and β4 were observed, consistent with NMR data. In addition, there were tertiary contacts between residues in α1 and the C-terminal portion of the β-sheet. The simulated structures allow the rudimentary NMR data to be fleshed out. The consistency between simulation and experiment inspires confidence in the methods. A description of the folding pathway of barnase from the denatured to the native state can be constructed by combining the simulation with experimental data from φ value analysis and NMR.

A nicked duplex decamer DNA with a PEG6 tether

A dumbbell double-stranded DNA decamer tethered with a hexaethylene glycol linker moiety (DDSDPEG), with a nick in the centre of one strand, has been synthesised. The standard NMR methods, E.COSY, TOCSY, NOESY and HMQC, were used to measure 1H, 31P and T1 spectral parameters. Molecular modelling using rMD-simulated annealing was used to compute the structure. Scalar couplings and dipolar contacts show that the molecule adopts a right-handed B-DNA helix in 38 mM phosphate buffer at pH 7. Its high melting temperature confirms the good base stacking and stability of the duplex. This is partly attributed to the presence of the PEG6 linker at both ends of the duplex that restricts the dynamics of the stem pentamers and thus stabilises the oligonucleotide. The inspection of the global parameters shows that the linker does not distort the B-DNA geometry. The computed structure suggests that the presence of the nick is not disturbing the overall tertiary structure, base pair geometry or duplex base pairing to a substantial extent. The nick has, however, a noticeable impact on the local geometry at the nick site, indicated clearly by NMR analysis and reflected in the conformational parameters of the computed structure. The 1H spectra also show much sharper resonances in the presence of K+ indicating that conformational heterogeneity of DDSDPEG is reduced in the presence of potassium as compared to sodium or caesium ions. At the same time the 1H resonances have longer T1 times. This parameter is suggested as a sensitive gauge of stabilisation.

Molecular dynamics of a grid-mounted molecular dipolar rotor in a rotating electric field

Classical molecular dynamics is applied to the rotation of a dipolar molecular rotor mounted on a square grid and driven by rotating electric field E(ν) at T ≃ 150 K. The rotor is a complex of Re with two substituted o-phenanthrolines, one positively and one negatively charged, attached to an axial position of Rh\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{2}^{4+}}}\end{equation*}\end{document} in a [2]staffanedicarboxylate grid through 2-(3-cyanobicyclo[1.1.1]pent-1-yl)malonic dialdehyde. Four regimes are characterized by a, the average lag per turn: (i) synchronous (a < 1/e) at E(ν) = |E(ν)| > Ec(ν) [Ec(ν) is the critical field strength], (ii) asynchronous (1/e < a < 1) at Ec(ν) > E(ν) > Ebo(ν) > kT/μ, [Ebo(ν) is the break-off field strength], (iii) random driven (a ≃ 1) at Ebo(ν) > E(ν) > kT/μ, and (iv) random thermal (a ≃ 1) at kT/μ > E(ν). A fifth regime, (v) strongly hindered, W > kT, Eμ, (W is the rotational barrier), has not been examined. We find Ebo(ν)/kVcm−1 ≃ (kT/μ)/kVcm−1 + 0.13(ν/GHz)1.9 and Ec(ν)/kVcm−1 ≃ (2.3kT/μ)/kVcm−1 + 0.87(ν/GHz)1.6. For ν > 40 GHz, the rotor behaves as a macroscopic body with a friction constant proportional to frequency, η/eVps ≃ 1.14 ν/THz, and for ν < 20 GHz, it exhibits a uniquely molecular behavior.

Nuclear magnetic dipole interactions in field-oriented proteins: information for structure determination in solution.

The measurement of dipolar contributions to the splitting of 15N resonances of 1H-15N amide pairs in multidimensional high-field NMR spectra of field-oriented cyanometmyoglobin is reported. The splittings appear as small field-dependent perturbations of normal scalar couplings. Assignment of more than 90 resonances to specific sequential sites in the protein allows correlation of the dipolar contributions with predictions based on the known susceptibility and known structure of the protein. Implications as an additional source of information for protein structure determination in solution are discussed.