Structural and functional innovations in the real-time evolution of new (ßα)8 barrel enzymes.

New genes can arise by duplication and divergence, but there is a fundamental gap in our understanding of the relationship between these genes, the evolving proteins they encode, and the fitness of the organism. Here we used crystallography, NMR dynamics, kinetics, and mass spectrometry to explain the molecular innovations that arose during a previous real-time evolution experiment. In that experiment, the (ßα)8 barrel enzyme HisA was under selection for two functions (HisA and TrpF), resulting in duplication and divergence of the hisA gene to encode TrpF specialists, HisA specialists, and bifunctional generalists. We found that selection affects enzyme structure and dynamics, and thus substrate preference, simultaneously and sequentially. Bifunctionality is associated with two distinct sets of loop conformations, each essential for one function. We observed two mechanisms for functional specialization: structural stabilization of each loop conformation and substrate-specific adaptation of the active site. Intracellular enzyme performance, calculated as the product of catalytic efficiency and relative expression level, was not linearly related to fitness. Instead, we observed thresholds for each activity above which further improvements in catalytic efficiency had little if any effect on growth rate. Overall, we have shown how beneficial substitutions selected during real-time evolution can lead to manifold changes in enzyme function and bacterial fitness. This work emphasizes the speed at which adaptive evolution can yield enzymes with sufficiently high activities such that they no longer limit the growth of their host organism, and confirms the (ßα)8 barrel as an inherently evolvable protein scaffold.

Measurement of protein backbone 13CO and 15N relaxation dispersion at high resolution.

Peak overlap in crowded regions of two-dimensional spectra prevents characterization of dynamics for many sites of interest in globular and intrinsically disordered proteins. We present new three-dimensional pulse sequences for measurement of Carr-Purcell-Meiboom-Gill relaxation dispersions at backbone nitrogen and carbonyl positions. To alleviate increase in the measurement time associated with the additional spectral dimension, we use non-uniform sampling in combination with two distinct methods of spectrum reconstruction: compressed sensing and co-processing with multi-dimensional decomposition. The new methodology was validated using disordered protein CD79A from B-cell receptor and an SH3 domain from Abp1p in exchange between its free form and bound to a peptide from the protein Ark1p. We show that, while providing much better resolution, the 3D NUS experiments give the similar accuracy and precision of the dynamic parameters to ones obtained using traditional 2D experiments. Furthermore, we show that jackknife resampling of the spectra yields robust estimates of peak intensities errors, eliminating the need for recording duplicate data points.

Comprehensive analysis of NMR data using advanced line shape fitting.

NMR spectroscopy is uniquely suited for atomic resolution studies of biomolecules such as proteins, nucleic acids and metabolites, since detailed information on structure and dynamics are encoded in positions and line shapes of peaks in NMR spectra. Unfortunately, accurate determination of these parameters is often complicated and time consuming, in part due to the need for different software at the various analysis steps and for validating the results. Here, we present an integrated, cross-platform and open-source software that is significantly more versatile than the typical line shape fitting application. The software is a completely redesigned version of PINT ( https://pint-nmr.github.io/PINT/ ). It features a graphical user interface and includes functionality for peak picking, editing of peak lists and line shape fitting. In addition, the obtained peak intensities can be used directly to extract, for instance, relaxation rates, heteronuclear NOE values and exchange parameters. In contrast to most available software the entire process from spectral visualization to preparation of publication-ready figures is done solely using PINT and often within minutes, thereby, increasing productivity for users of all experience levels. Unique to the software are also the outstanding tools for evaluating the quality of the fitting results and extensive, but easy-to-use, customization of the fitting protocol and graphical output. In this communication, we describe the features of the new version of PINT and benchmark its performance.

Structural Insights into the Calcium-Mediated Allosteric Transition in the C-Terminal Domain of Calmodulin from Nuclear Magnetic Resonance Measurements.

Calmodulin is a two-domain signaling protein that becomes activated upon binding cooperatively two pairs of calcium ions, leading to large-scale conformational changes that expose its binding site. Despite significant advances in understanding the structural biology of calmodulin functions, the mechanistic details of the conformational transition between closed and open states have remained unclear. To investigate this transition, we used a combination of molecular dynamics simulations and nuclear magnetic resonance (NMR) experiments on the Ca(2+)-saturated E140Q C-terminal domain variant. Using chemical shift restraints in replica-averaged metadynamics simulations, we obtained a high-resolution structural ensemble consisting of two conformational states and validated such an ensemble against three independent experimental data sets, namely, interproton nuclear Overhauser enhancements, (15)N order parameters, and chemical shift differences between the exchanging states. Through a detailed analysis of this structural ensemble and of the corresponding statistical weights, we characterized a calcium-mediated conformational transition whereby the coordination of Ca(2+) by just one oxygen of the bidentate ligand E140 triggers a concerted movement of the two EF-hands that exposes the target binding site. This analysis provides atomistic insights into a possible Ca(2+)-mediated activation mechanism of calmodulin that cannot be achieved from static structures alone or from ensemble NMR measurements of the transition between conformations.

Fast and accurate resonance assignment of small-to-large proteins by combining automated and manual approaches.

The process of resonance assignment is fundamental to most NMR studies of protein structure and dynamics. Unfortunately, the manual assignment of residues is tedious and time-consuming, and can represent a significant bottleneck for further characterization. Furthermore, while automated approaches have been developed, they are often limited in their accuracy, particularly for larger proteins. Here, we address this by introducing the software COMPASS, which, by combining automated resonance assignment with manual intervention, is able to achieve accuracy approaching that from manual assignments at greatly accelerated speeds. Moreover, by including the option to compensate for isotope shift effects in deuterated proteins, COMPASS is far more accurate for larger proteins than existing automated methods. COMPASS is an open-source project licensed under GNU General Public License and is available for download from http://www.liu.se/forskning/foass/tidigare-foass/patrik-lundstrom/software?l=en. Source code and binaries for Linux, Mac OS X and Microsoft Windows are available.

Capturing Excited States in the Fast-Intermediate Exchange Limit in Biological Systems Using 1 H†NMR Spectroscopy.

Changes in molecular structure are essential for the function of biomolecules. Characterization of these structural fluctuations can illuminate alternative states and help in correlating structure to function. NMR relaxation dispersion (RD) is currently the only method for detecting these alternative, high-energy states. In this study, we present a versatile 1 H R1ρ RD experiment that not only extends the exchange timescales at least three times beyond the rate limits of 13 C/15 N R1ρ and ten times for CPMG experiments, but also makes use of easily accessible probes, thus allowing a general description of biologically important excited states. This technique can be used to extract chemical shifts for the structural characterization of excited states and to elucidate complex excited states.

Local destabilization of the metal-binding region in human copper-zinc superoxide dismutase by remote mutations is a possible determinant for progression of ALS.

More than 100 distinct mutations in the gene CuZnSOD encoding human copper-zinc superoxide dismutase (CuZnSOD) have been associated with familial amyotrophic lateral sclerosis (fALS), a fatal neuronal disease. Many studies of different mutant proteins have found effects on protein stability, catalytic activity, and metal binding, but without a common pattern. Notably, these studies were often performed under conditions far from physiological. Here, we have used experimental conditions of pH 7 and 37 °C and at an ionic strength of 0.2 M to mimic physiological conditions as close as possible in a sample of pure protein. Thus, by using NMR spectroscopy, we have analyzed amide hydrogen exchange of the fALS-associated I113T CuZnSOD variant in its fully metalated state, both at 25 and 37 °C, where (15)N relaxation data, as expected, reveals that CuZnSOD I113T exists as a dimer under these conditions. The local dynamics at 82% of all residues have been analyzed in detail. When compared to the wild-type protein, it was found that I113T CuZnSOD is particularly destabilized locally at the ion binding sites of loop 4, the zinc binding loop, which results in frequent exposure of the aggregation prone outer ß-strands I and VI of the ß-barrel, possibly enabling fibril or aggregate formation. A similar study (Museth, A. K., et al. (2009) Biochemistry, 48, 8817-8829) of amide hydrogen exchange at pH 7 and 25 °C on the G93A variant also revealed a selective destabilization of the zinc binding loop. Thus, a possible scenario in ALS is that elevated local dynamics at the metal binding region can result in toxic species from formation of new interactions at local ß-strands.

Fractional enrichment of proteins using [2-(13)C]-glycerol as the carbon source facilitates measurement of excited state 13Cα chemical shifts with improved sensitivity.

A selective isotope labeling scheme based on the utilization of [2-(13)C]-glycerol as the carbon source during protein overexpression has been evaluated for the measurement of excited state (13)Cα chemical shifts using Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion (RD) experiments. As expected, the fractional incorporation of label at the Cα positions is increased two-fold relative to labeling schemes based on [2-(13)C]-glucose, effectively doubling the sensitivity of NMR experiments. Applications to a binding reaction involving an SH3 domain from the protein Abp1p and a peptide from the protein Ark1p establish that accurate excited state (13)Cα chemical shifts can be obtained from RD experiments, with errors on the order of 0.06 ppm for exchange rates ranging from 100 to 1000 s(-1), despite the small fraction of (13)Cα-(13)Cß spin-pairs that are present for many residue types. The labeling approach described here should thus be attractive for studies of exchanging systems using (13)Cα spin probes.

A high-affinity, dimeric inhibitor of PSD-95 bivalently interacts with PDZ1-2 and protects against ischemic brain damage.

Inhibition of the ternary protein complex of the synaptic scaffolding protein postsynaptic density protein-95 (PSD-95), neuronal nitric oxide synthase (nNOS), and the N-methyl-D-aspartate (NMDA) receptor is a potential strategy for treating ischemic brain damage, but high-affinity inhibitors are lacking. Here we report the design and synthesis of a novel dimeric inhibitor, Tat-NPEG4(IETDV)(2) (Tat-N-dimer), which binds the tandem PDZ1-2 domain of PSD-95 with an unprecedented high affinity of 4.6 nM, and displays extensive protease-resistance as evaluated in vitro by stability-measurements in human blood plasma. X-ray crystallography, NMR, and small-angle X-ray scattering (SAXS) deduced a true bivalent interaction between dimeric inhibitor and PDZ1-2, and also provided a dynamic model of the conformational changes of PDZ1-2 induced by the dimeric inhibitor. A single intravenous injection of Tat-N-dimer (3 nmol/g) to mice subjected to focal cerebral ischemia reduces infarct volume with 40% and restores motor functions. Thus, Tat-N-dimer is a highly efficacious neuroprotective agent with therapeutic potential in stroke.

Gastrointestinal quality-of-life after cholecystectomy: indication predicts gastrointestinal symptoms and abdominal pain.

BACKGROUND: Despite the fact that cholecystectomy is a common surgical procedure, the impact on long-term gastrointestinal quality of life is not fully known. METHODS: All surgical procedures for gallstone disease performed at Mora County Hospital, Sweden, between 2 January 2002 and 2 January 2005, were registered on a standard database form. In 2007, all patients under the age of 80 years at follow-up were requested to fill in a form containing the Gastrointestinal Quality-of-Life Index (GIQLI) questionnaire and a number of additional questions. The outcome was analysed with respect to age, gender, smoking, surgical technique, and original indication for cholecystectomy. RESULTS: A total of 627 patients (447 women, 180 men) underwent cholecystectomy, including laparoscopic cholecystectomy (N = 524), laparoscopic cholecystectomy converted to open cholecystectomy (N = 43), and open cholecystectomy (N = 60). The mean time between cholecystectomy and follow-up with the questionnaire was 49 months. The participation rate was 79 %. Using multivariate analysis in the form of generalised linear modelling, the original indication for cholecystectomy in combination with gender (p = 0.0042) was found to predict the GIQLI score. Female gender in combination with biliary colic as indication for cholecystectomy correlated with low GIQLI scores. Female gender also correlated with a higher risk for pain in the right upper abdominal quadrant after cholecystectomy (p = 0.028). CONCLUSIONS: We found the original indication for cholecystectomy, together with gender, to predict gastrointestinal symptoms and abdominal pain after cholecystectomy. Careful evaluation of symptoms is important before planning elective cholecystectomy.

Transient structure and dynamics in the disordered c-Myc transactivation domain affect Bin1 binding.

The crucial role of Myc as an oncoprotein and as a key regulator of cell growth makes it essential to understand the molecular basis of Myc function. The N-terminal region of c-Myc coordinates a wealth of protein interactions involved in transformation, differentiation and apoptosis. We have characterized in detail the intrinsically disordered properties of Myc-1-88, where hierarchical phosphorylation of S62 and T58 regulates activation and destruction of the Myc protein. By nuclear magnetic resonance (NMR) chemical shift analysis, relaxation measurements and NOE analysis, we show that although Myc occupies a very heterogeneous conformational space, we find transiently structured regions in residues 22-33 and in the Myc homology box I (MBI; residues 45-65); both these regions are conserved in other members of the Myc family. Binding of Bin1 to Myc-1-88 as assayed by NMR and surface plasmon resonance (SPR) revealed primary binding to the S62 region in a dynamically disordered and multivalent complex, accompanied by population shifts leading to altered intramolecular conformational dynamics. These findings expand the increasingly recognized concept of intrinsically disordered regions mediating transient interactions to Myc, a key transcriptional regulator of major medical importance, and have important implications for further understanding its multifaceted role in gene regulation.

PINT: a software for integration of peak volumes and extraction of relaxation rates.

We present the software Peak INTegration (PINT), designed to perform integration of peaks in NMR spectra. The program is very simple to run, yet powerful enough to handle complicated spectra. Peaks are integrated by fitting predefined line shapes to experimental data and the fitting can be customized to deal with, for instance, heavily overlapped peaks. The results can be inspected visually, which facilitates systematic optimization of the line shape fitting. Finally, integrated peak volumes can be used to extract parameters such as relaxation rates and information about low populated states. The utility of PINT is demonstrated by applications to the 59 residue SH3 domain of the yeast protein Abp1p and the 289 residue kinase domain of murine EphB2.

Protein conformational exchange measured by 1H R1ρ relaxation dispersion of methyl groups.

Activated dynamics plays a central role in protein function, where transitions between distinct conformations often underlie the switching between active and inactive states. The characteristic time scales of these transitions typically fall in the microsecond to millisecond range, which is amenable to investigations by NMR relaxation dispersion experiments. Processes at the faster end of this range are more challenging to study, because higher RF field strengths are required to achieve refocusing of the exchanging magnetization. Here we describe a rotating-frame relaxation dispersion experiment for (1)H spins in methyl (13)CHD2 groups, which improves the characterization of fast exchange processes. The influence of (1)H-(1)H rotating-frame nuclear Overhauser effects (ROE) is shown to be negligible, based on a comparison of R 1ρ relaxation data acquired with tilt angles of 90° and 35°, in which the ROE is maximal and minimal, respectively, and on samples containing different (1)H densities surrounding the monitored methyl groups. The method was applied to ubiquitin and the apo form of calmodulin. We find that ubiquitin does not exhibit any (1)H relaxation dispersion of its methyl groups at 10 or 25 °C. By contrast, calmodulin shows significant conformational exchange of the methionine methyl groups in its C-terminal domain, as previously demonstrated by (1)H and (13)C CPMG experiments. The present R 1ρ experiment extends the relaxation dispersion profile towards higher refocusing frequencies, which improves the definition of the exchange correlation time, compared to previous results.

Appropriation of the MinD protein-interaction motif by the dimeric interface of the bacterial cell division regulator MinE.

MinE is required for the dynamic oscillation of Min proteins that restricts formation of the cytokinetic septum to the midpoint of the cell in gram negative bacteria. Critical for this oscillation is MinD-binding by MinE to stimulate MinD ATP hydrolysis, a function that had been assigned to the first ∼30 residues in MinE. Previous models based on the structure of an autonomously folded dimeric C-terminal fragment suggested that the N-terminal domain is freely accessible for interactions with MinD. We report here the solution NMR structure of the full-length MinE dimer from Neisseria gonorrhoeae, with two parts of the N-terminal domain forming an integral part of the dimerization interface. Unexpectedly, solvent accessibility is highly restricted for residues that were previously hypothesized to directly interact with MinD. To delineate the true MinD-binding region, in vitro assays for MinE-stimulated MinD activity were performed. The relative MinD-binding affinities obtained for full-length and N-terminal peptides from MinE demonstrated that residues that are buried in the dimeric interface nonetheless participate in direct interactions with MinD. According to results from NMR spin relaxation experiments, access to these buried residues may be facilitated by the presence of conformational exchange. We suggest that this concealment of MinD-binding residues by the MinE dimeric interface provides a mechanism for prevention of nonspecific interactions, particularly with the lipid membrane, to allow the free diffusion of MinE that is critical for Min protein oscillation.

Interactions outside the boundaries of the canonical binding groove of a PDZ domain influence ligand binding.

The postsynaptic density protein-95/discs large/zonula occludens-1 (PDZ) domain is a protein-protein interaction module with a shallow binding groove where protein ligands bind. However, interactions that are not part of this canonical binding groove are likely to modulate peptide binding. We have investigated such interactions beyond the binding groove for PDZ3 from PSD-95 and a peptide derived from the C-terminus of the natural ligand CRIPT. We found via nuclear magnetic resonance experiments that up to eight residues of the peptide ligand interact with the PDZ domain, showing that the interaction surface extends far outside of the binding groove as defined by the crystal structure. PDZ3 contains an extra structural element, a C-terminal helix (α3), which is known to affect affinity. Deletion of this helix resulted in the loss of several intermolecular nuclear Overhauser enhancements from peptide residues outside of the binding pocket, suggesting that α3 forms part of the extra binding surface in wild-type PDZ3. Site-directed mutagenesis, isothermal titration calorimetry, and fluorescence intensity experiments confirmed the importance of both α3 and the N-terminal part of the peptide for the affinity. Our data suggest a general mechanism in which different binding surfaces outside of the PDZ binding groove could provide sites for specific interactions.

Biophysical characterization of the complex between human papillomavirus E6 protein and synapse-associated protein 97.

The E6 protein of human papillomavirus (HPV) exhibits complex interaction patterns with several host proteins, and their roles in HPV-mediated oncogenesis have proved challenging to study. Here we use several biophysical techniques to explore the binding of E6 to the three PDZ domains of the tumor suppressor protein synapse-associated protein 97 (SAP97). All of the potential binding sites in SAP97 bind E6 with micromolar affinity. The dissociation rate constants govern the different affinities of HPV16 and HPV18 E6 for SAP97. Unexpectedly, binding is not mutually exclusive, and all three PDZ domains can simultaneously bind E6. Intriguingly, this quaternary complex has the same apparent hydrodynamic volume as the unliganded PDZ region, suggesting that a conformational change occurs in the PDZ region upon binding, a conclusion supported by kinetic experiments. Using NMR, we discovered a new mode of interaction between E6 and PDZ: a subset of residues distal to the canonical binding pocket in the PDZ(2) domain exhibited noncanonical interactions with the E6 protein. This is consistent with a larger proportion of the protein surface defining binding specificity, as compared with that reported previously.

Quantifying millisecond exchange dynamics in proteins by CPMG relaxation dispersion NMR using side-chain 1H probes.

A Carr-Purcell-Meiboom-Gill relaxation dispersion experiment is presented for quantifying millisecond time-scale chemical exchange at side-chain (1)H positions in proteins. Such experiments are not possible in a fully protonated molecule because of magnetization evolution from homonuclear scalar couplings that interferes with the extraction of accurate transverse relaxation rates. It is shown, however, that by using a labeling strategy whereby proteins are produced using {(13)C,(1)H}-glucose and D(2)O a significant number of 'isolated' side-chain (1)H spins are generated, eliminating such effects. It thus becomes possible to record (1)H dispersion profiles at the ß positions of Asx, Cys, Ser, His, Phe, Tyr, and Trp as well as the γ positions of Glx, in addition to the methyl side-chain moieties. This brings the total of amino acid side-chain positions that can be simultaneously probed using a single (1)H dispersion experiment to 16. The utility of the approach is demonstrated with an application to the four-helix bundle colicin E7 immunity protein, Im7, which folds via a partially structured low populated intermediate that interconverts with the folded, ground state on the millisecond time-scale. The extracted (1)H chemical shift differences at side-chain positions provide valuable restraints in structural studies of invisible, excited states, complementing backbone chemical shifts that are available from existing relaxation dispersion experiments.

Side-chain interactions form late and cooperatively in the binding reaction between disordered peptides and PDZ domains.

Intrinsically disordered proteins are very common and mediate numerous protein-protein and protein-DNA interactions. While it is clear that these interactions are instrumental for the life of the mammalian cell, there is a paucity of data regarding their molecular binding mechanisms. Here we have used short peptides as a model system for intrinsically disordered proteins. Linear free energy relationships based on rate and equilibrium constants for the binding of these peptides to ordered target proteins, PDZ domains, demonstrate that native side-chain interactions form mainly after the rate-limiting barrier for binding and in a cooperative fashion. This finding suggests that these disordered peptides first form a weak encounter complex with non-native interactions. The data do not support the recent notion that the affinities of intrinsically disordered proteins toward their targets are generally governed by their association rate constants. Instead, we observed the opposite for peptide-PDZ interactions, namely, that changes in K(d) correlate with changes in k(off).

Isotope labeling methods for large systems.

A major drawback of nuclear magnetic resonance (NMR) spectroscopy compared to other methods is that the technique has been limited to relatively small molecules. However, in the last two decades the size limit has been pushed upwards considerably and it is now possible to use NMR spectroscopy for structure calculations of proteins of molecular weights approaching 100 kDa and to probe dynamics for supramolecular complexes of molecular weights in excess of 500 kDa. Instrumental for this progress has been development in instrumentation and pulse sequence design but also improved isotopic labeling schemes that lead to increased sensitivity as well as improved spectral resolution and simplification. These are described and discussed in this chapter, focusing on labeling schemes for amide proton and methyl proton detected experiments. We also discuss labeling methods for other potentially useful positions in proteins.

Isotope labeling methods for relaxation measurements.

Nuclear magnetic spin relaxation has emerged as a powerful technique for probing molecular dynamics. Not only is it possible to use it for determination of time constant(s) for molecular reorientation but it can also be used to characterize internal motions on time scales from picoseconds to seconds. Traditionally, uniformly (15)N labeled samples have been used for these experiments but it is clear that this limits the applications. For instance, sensitivity for large systems is dramatically increased if dynamics is probed at methyl groups and structural characterization of low-populated states requires measurements on (13)Cα, (13)Cß or (13)CO or (1)Hα. Unfortunately, homonuclear scalar couplings may lead to artifacts in the latter types of experiments and selective isotopic labeling schemes that only label the desired position are necessary. Both selective and uniform labeling schemes for measurements of relaxation rates for a large number of positions in proteins are discussed in this chapter.