SESAME: A least-squares approach to the evaluation of protein structures computed from NMR data

Summary A method is proposed for defining a probability distribution on an ensemble of protein conformations from a 2D NOE spectrum, while at the same time back-calculating the experimental spectrum from the ensemble. This enables one to assess the relative quality and significance of the conformations, and to test the consistency of the ensemble as a whole with the experimental spectrum. The method eliminates the need to integrate the cross-peak intensities and is surprisingly insensitive to random noise in the spectrum. In this communication, these advantages are demonstrated by applying the method to simulated data, for which the correct result is already known.Abbreviations NOE nuclear Overhauser effect - BPTI basic pancreatic trypsin inhibitor - rmsd root-mean-square deviation     

Recent Advances in Computational Methods for Nuclear Magnetic Resonance Data Processing

Although three-dimensional protein structure determination using nuclear magnetic resonance (NMR) spectroscopy is a computationally costly and tedious process that would benefit from advanced computational techniques, it has not garnered much research attention from specialists in bioinformatics and computational biology. In this paper, we review recent advances in computational methods for NMR protein structure determination. We summarize the advantages of and bottlenecks in the existing methods and outline some open problems in the field. We also discuss current trends in NMR technology development and suggest directions for research on future computational methods for NMR.     

H assignment and secondary structure determination of human melanoma growth stimulating activity (MGSA) by NMR spectroscopy

The solution structure of melanoma growth stimulating activity (MGSA) has been investigated using proton NMR spectroscopy. Sequential resonance assignments have been carried out, and elements of secondary structure have been identified on the basis of NOE, coupling constant, chemical shift, and amide proton exchange data. Long-range NOEs have established that MGSA is a dimer in solution. The secondary structure and dimer interface of MGSA appear to be similar to those found previously for the homologous chemokine interleukin-8 [Clore et al. (1990) Biochemistry 29, 1689-1696]. The MGSA monomer contains a three stranded anti-parallel β-sheet arranged in a ‘Greek-key’ conformation, and a C-terininal -helix (residues 58 69).     

Folding and activity of cAMP-dependent protein kinase mutants

The catalytic subunit of cAMP-dependent protein kinase (PKA) can easily be expressed in Escherichia coli and is catalytically active. Four phosphorylation sites are known in PKA (S10, S139, T197 and S338), and the isolated recombinant protein is a mixture of different phosphorylated forms. Obtaining uniformly phosphorylated protein requires separation of the protein preparation leading to significant loss in protein yield. It is found that the mutant S10A/S139D/S338D has similar properties as the wild-type protein, whereas additional replacement of T197 with either E or D reduces protein expression yield as well as folding propensity of the protein. Due to its high sequence homology to Akt/PKB, which cannot easily be expressed in E. coli, PKA has been used as a surrogate kinase for drug design. Several mutations within the ATP binding site have been described to make PKA even more similar to Akt/PKB. Two proteins with Akt/PKB-like mutations in the ATP binding site were made (PKAB6 and PKAB8), and in addition S10, S139 and S338 phosphorylation sites have been removed. These proteins can be expressed in high yields but have reduced activity compared to the wild-type. Proper folding of all proteins was analyzed by 2D 1H, 15N-TROSY NMR experiments.     

Efficient analysis of protein 2D NMR spectra using the software packageEASY

Summary The programEASY supports the spectral analysis of biomacromolecular two-dimensional (2D) nuclear magnetic resonance (NMR) data. It provides a user-friendly, window-based environment in which to view spectra for interactive interpretation. In addition, it includes a number of automated routines for peakpicking, spin-system identification, sequential resonance assignment in polypeptide chains, and cross peak integration. In this uniform environment, all resulting parameter lists can be recorded on disk, so that the paper plots and handwritten notes which normally accompany manual assignment of spectra can be largely eliminated. For example, in a protein structure determination by 2D¹H NMR,EASY accepts the frequency domain datasets as input, and after combined use of the automated and interactive routines it can yield a listing of conformational constraints in the format required as input for the calculation of the 3D structure. The program was extensively tested with current protein structure determinations in our laboratory. In this paper, its main features are illustrated with data on the protein basic pancreatic trypsin inhibitor.     

Dimethylsulfoxide-quenched hydrogen/deuterium exchange method to study amyloid fibril structure

A general method to analyze the structure of a supramolecular complex of amyloid fibrils at amino acid residue resolution has been developed. This method combines the NMR-detected hydrogen/deuterium (H/D) exchange technique to detect hydrogen-bonded amide groups and the ability of the aprotic organic solvent dimethylsulfoxide (DMSO) to dissolve amyloid fibrils into NMR-observable, monomeric components while suppressing the undesired H/D exchange reaction. Moreover, this method can be generally applied to amyloid fibrils to elucidate the distribution of hydrogen-bonded amino acid residues in the three-dimensional molecular organization in the amyloid fibrils. In this study, we describe theoretical considerations in the H/D exchange method to obtain the structural information of proteins, and the DMSO-quenched H/D exchange method to study a supramolecular complex of amyloid fibrils. A possible application of this method to study the interaction of a protein/peptide with phospholipid membrane is also discussed.     

Challenges at the frontiers of structural biology

Knowledge of the three-dimensional structures of proteins is the key to unlocking the full potential of genomic information. There are two distinct directions along which cutting-edge research in structural biology is currently moving towards this goal. On the one hand, tightly focused long-term research in individual laboratories is leading to the determination of the structures of macromolecular assemblies of ever-increasing size and complexity. On the other hand, large consortia of structural biologists, inspired by the pace of genome sequencing, are developing strategies to determine new protein structures rapidly, so that it will soon be possible to predict reasonably accurate structures for most protein domains. We anticipate that a small number of complex systems, studied in depth, will provide insights across the field of biology with the aid of genome-based comparative structural analysis.     

The role of dynamics in modulating ligand exchange in intracellular lipid binding proteins

Lipids are essential for many biological processes and crucial in the pathogenesis of several diseases. Intracellular lipid-binding proteins (iLBPs) provide mobile hydrophobic binding sites that allow hydrophobic or amphipathic lipid molecules to penetrate into and across aqueous layers. Thus iLBPs mediate the lipid transport within the cell and participate to a spectrum of tissue-specific pathways involved in lipid homeostasis. Structural studies have shown that iLBPs' binding sites are inaccessible from the bulk, implying that substrate binding should involve a conformational change able to produce a ligand entry portal. Many studies have been reported in the last two decades on iLBPs indicating that their dynamics play a pivotal role in regulating ligand binding and targeted release. The ensemble of reported data has not been reviewed until today. This review is thus intended to summarize and possibly generalize the results up to now described, providing a picture which could help to identify the missing notions necessary to improve our understanding of the role of dynamics in iLBPs' molecular recognition. Such notions would clarify the chemistry of lipid binding to iLBPs and set the basis for the development of new drugs.     

A computer-based protocol for semiautomated assignments and 3D structure determination of proteins

Summary A strategy is presented for the semiautomated assignment and 3D structure determination of proteins from heteronuclear multidimensional Nuclear Magnetic Resonance (NMR) data. This approach involves the computer-based assignment of the NMR signals, identification of distance restraints from nuclear Overhauser effects, and generation of 3D structures by using the NMR-derived restraints. The protocol is described in detail and illustrated on a resonance assignment and structure determination of the FK506 binding protein (FKBP, 107 amino acids) complexed to the immunosuppressant, ascomycin. The 3D structures produced from this automated protocol attained backbone and heavy atom rmsd of 1.17 and 1.69 Å, respectively. Although more highly resolved structures of the complex have been obtained by standard interpretation of NMR data (Meadows et al. (1993) Biochemistry, 32, 754–765), the structures generated with this automated protocol required minimal manual intervention during the spectral assignment and 3D structure calculations stages. Thus, the protocol may yield an approximate order of magnitude reduction in the time required for the generation of 3D structures of proteins from NMR data.     

Determination of melarsoprol in biological fluids by high-performance liquid chromatography and characterisation of two stereoisomers by nuclear magnetic resonance spectroscopy

The analysis of melarsoprol in whole blood, plasma, urine and cerebrospinal fluid is described. Extraction was made with a mixture of chloroform and acetonitrile followed by back-extraction into phosphoric acid. A reversed-phase liquid chromatography system with ultraviolet detection was used. The relative standard deviation was 1% at concentrations around 10 μmol/l and 3–6% at the lower limit of determination (9 nmol/l in plasma, 93 nmol/l in whole blood, 45 nmol/l in urine and 10 nmol/l in cerebrospinal fluid). Melarsoprol is not a stable compound and samples to be stored for longer periods of time should be kept at −70°C. Plasma samples can be stored at −20°C for upt to 2 months. Chromatography showed that melarsoprol contains two components. Using nuclear magnetic resonance spectroscopy the two components were shown to be diastereomers which slowly equilibrate by inversion of the configuration at the As atom.     

Solution 1H NMR study of the accommodation of the side chain of n-butyl-etiohemin-I incorporated into the active site of cyano-metmyoglobin

In order to identify the most readily deformable portion of the heme pocket in myoglobin, equine myoglobin was reconstituted with a meso-n-butyl substituent on centrosymmetric etiohemin-I. Solution ¹H NMR data for the low-spin iron(III) cyanide complex of oxidized myoglobin that include 2D nuclear Overhauser enhancement spectroscopy contacts, paramagnetic relaxation, and dipolar shifts resulting from magnetic anisotropy show that the heme binds uniquely to the iron in a manner that arranges the methyl and ethyl substituents on a given pyrrole in a clockwise manner when viewed from the proximal side, and with the n-butyl group seated at the canonical -meso position of native protohemin-IX. The butyl group is oriented sharply toward the proximal side and its protein contacts demonstrate that it is oriented largely into the xenon hole in myoglobin. The location of the n-butyl group on the proximal side near the vacancies places it within the region found to be most flexible in molecular dynamics simulation. A small, counterclockwise rotation of the pyrrole N–Fe–N vector of n-butyl-etiohemin-I relative to that for native protohemin, indicated by both the prosthetic group methyl contact shift pattern and the prosthetic group contacts to heme pocket residues, is proposed to allow the xenon hole to accommodate better the n-butyl group. In contrast to previous work, which showed that a bulky polar substituent on etiohemin-I preferentially seats at the canonical -meso position, the nonpolar n-butyl group selects the -meso position.Electronic Supplementary Material Supplementary material is available for this article at .     

Sampling and efficiency of metric matrix distance geometry: A novel partial metrization algorithm

Summary In this paper, we present a reassessment of the sampling properties of the metric matrix distance geometry algorithm, which is in wide-spread use in the determination of three-dimensional structures from nuclear magnetic resonance (NMR) data. To this end, we compare the conformational space sampled by structures generated with a variety of metric matrix distance geometry protocols. As test systems we use an unconstrained polypeptide, and a small protein (rabbit neutrophil defensin peptide 5) for which only few tertiary distances had been derived from the NMR data, allowing several possible folds of the polypeptide chain. A process called metrization in the preparation of a trial distance matrix has a very large effect on the sampling properties of the algorithm. It is shown that, depending on the metrization protocol used, metric matrix distance geometry can have very good sampling properties'indeed, both for the unconstrained model system and the NMR-structure case. We show that the sampling properties are to a great degree determined by the way in which the first few distances are chosen within their bounds. Further, we present a new protocol (partial metrization) that is computationally more efficient but has the same excellent sampling properties. This novel protocol has been implemented in an expanded new release of the program X-PLOR with distance geometry capabilities.     

Solution structure of betacellulin,a new member of EGF-family ligands

The solution structure of the EGF-like domain of betacellulin (BTCe), a newly discovered member of the epidermal growth factor (EGF) family, has been determined using two-dimensional nuclear magnetic resonance spectroscopy. This is the first report to identify the solution structure of the EGF-family ligand monomers that interact with both ErbB-1 and ErbB-4. The solution structure of BTCe was calculated using 538 NMR-derived restraints. The overall structure of BTCe was stabilized by three disulfide bonds, a hydrophobic core, and 23 hydrogen bonds. It appears that BTCe is comprised of five beta-strands and one short 3(10) helical turn. The secondary structural elements of BTCe are basically similar to those of the other EGF-family proteins, except that several significant variations of the structural properties were found. It is suggested that the structural variations between BTCe and the other EGF-family ligands may affect the specific receptor-recognition properties of EGF-family ligands.     

The MUMO (minimal under-restraining minimal over-restraining) method for the determination of native state ensembles of proteins

While reliable procedures for determining the conformations of proteins are available, methods for generating ensembles of structures that also reflect their flexibility are much less well established. Here we present a systematic assessment of the ability of ensemble-averaged molecular dynamics simulations with ensemble-averaged NMR restraints to simultaneously reproduce the average structure of proteins and their associated dynamics. We discuss the effects that under-restraining (overfitting) and over-restraining (underfitting) have on the structures generated in ensemble-averaged molecular simulations. We then introduce the MUMO (minimal under-restraining minimal over-restraining) method, a procedure in which different observables are averaged over a different number of molecules. As both over-restraining and under-restraining are significantly reduced in the MUMO method, it is possible to generate ensembles of conformations that accurately characterize both the structure and the dynamics of native states of proteins. The application of the MUMO method to the protein ubiquitin yields a high-resolution structural ensemble with an RDC Q-factor of 0.19.