Computer assignment of the backbone resonances of labelled proteins using two-dimensional correlation experiments

Summary We present ALPS (Assignment for Labelled Protein Spectra), a flexible computer program for the automatic assignment of backbone NMR resonances of ¹⁵N/¹³C-labelled proteins. The program constructs pseudoresidues from peak-picking lists of a set of two-dimensional triple resonance experiments and uses either a systematic search or a simulated annealing-based optimization to perform the assignment. This method has been successfully tested on two-dimensional triple resonance spectra of Rhodobacter capsulatus ferrocytochrome c ₂ (116 amino acids).     

Determination of the structure of [Nle7]-endothelin by 1H NMR

[Nle7]-endothelin was synthesized and studied by 1H NMR and distance geometry calculations. The NMR study was performed first in DMSO-d6 and then in 50% acetonitrile/water since this peptide aggregates in pure water. In both cases, all spin systems were identified and assigned with the aid of two-dimensional spectroscopy (2D): COSY (for scalar couplings) and NOESY (for dipolar couplings). On the basis of the acetonitrile/water NMR parameters, and using the DISGEO program, a three-dimensional structure of [Nle7]-endothelin is proposed and discussed.     

SQUID: a program for the analysis and display of data from crystallography and molecular dynamics

SQUID is a flexible computer program that allows the analysis and display of molecular coordinates from crystallography, NMR, and molecular dynamics. The program can also display two-dimensional and three-dimensional data using many graph types, as well as perform array processing of data with numerous intrinsic functions. Graphics are based on the use of “move” and “draw” instructions, allowing easy development of new device drivers, including vector plotters.     

Two-dimensional graphic analysis of DNA sequence homologies.

We describe a computer program designed to facilitate the pattern matching analysis of homologies between DNA sequences. It takes advantage of a two-dimensional plot in order to simplify the evaluation of significant structures inherited in the sequences. The program can be divided into three parts, i) algorithm for search of homologies, ii) two-dimensional graphic display of the result, iii) further graphic treatment to enhance significant structures. The power of the graphic display is presented by the following application of the program. We conducted a search for direct repeats in the mouse immunoglobulin kappa-chain genes. Both the five J DNA sequences and other shorter repeats were found. We also found a longer stretch of homology that could indicate the presence of duplicated DNA in the J4, J5 region.     

Three-dimensional structure of acyl carrier protein determined by NMR pseudoenergy and distance geometry calculations

Distance constraints from two-dimensional NMR cross-relaxation data are used to derive a three-dimensional structure for acyl carrier protein from Escherichia coli. Several approaches to structure determination are explored. The most successful proves to be an approach that combines the early stages of a distance geometry program with energy minimization in the presence of NMR constraints represented as pseudopotentials. Approximately 450 proton to proton distance constraints including 50 long-range constraints were included in these programs. Starting structures were generated at random by the distance geometry program and energies minimized by a molecular mechanics module to give final structures. Seven of the structures were deemed acceptable on the basis of agreement with experimentally determined distances. Root-mean-square deviations from the mean of these structures for backbone atoms range from 2 to 3 A. All structures show three roughly parallel helices with hydrophobic residues facing inward and hydrophilic residues facing outward. A hydrophobic cleft is recognizable and is identified as a likely site for acyl chain binding.     

Complete (1)H and (13)C NMR chemical shift assignments of mono-, di-, and trisaccharides as basis for NMR chemical shift predictions of polysaccharides using the computer program casper

The computer program casper uses (1)H and (13)C NMR chemical shift data of mono- to trisaccharides for the prediction of chemical shifts of oligo- and polysaccharides. In order to improve the quality of these predictions the (1)H and (13)C, as well as (31)P when applicable, NMR chemical shifts of 30 mono-, di-, and trisaccharides were assigned. The reducing sugars gave two distinct sets of NMR resonances due to the α- and β-anomeric forms. In total 35 (1)H and (13)C NMR chemical shift data sets were obtained from the oligosaccharides. One- and two-dimensional NMR experiments were used for the chemical shift assignments and special techniques were employed in some cases such as 2D (1)H,(13)C-HSQC Hadamard Transform methodology which was acquired approximately 45 times faster than a regular t(1) incremented (1)H,(13)C-HSQC experiment and a 1D (1)H,(1)H-CSSF-TOCSY experiment which was able to distinguish spin-systems in which the target protons were only 3.3Hz apart. The (1)H NMR chemical shifts were subsequently refined using total line-shape analysis with the PERCH NMR software. The acquired NMR data were then utilized in the casper program (http://www.casper.organ.su.se/casper/) for NMR chemical shift predictions of the O-antigen polysaccharides from Klebsiella O5, Shigella flexneri serotype X, and Salmonella arizonae O62. The data were compared to experimental data of the polysaccharides from the two former strains and the lipopolysaccharide of the latter strain showing excellent agreement between predicted and experimental (1)H and (13)C NMR chemical shifts.     

Solution-state 2D NMR of Ball-milled Plant Cell Wall Gels in DMSO-d 6

Although finely divided ball-milled whole cell walls do not completely dissolve in dimethylsulfoxide (DMSO), they readily swell producing a gel. Solution-state two-dimensional (2D) nuclear magnetic resonance (NMR) of this gel, produced directly in the NMR tube, provides an interpretable structural fingerprint of the polysaccharide and lignin components of the wall without actual solubilization, and without structural modification beyond that inflicted by the ball milling and ultrasonication steps. Since the cellulose is highly crystalline and difficult to swell, the component may be under-represented in the spectra. The method however provides a more rapid method for comparative structural evaluation of plant cell walls than is currently available. With the new potential for chemometric analysis using the 2D NMR fingerprint, this method may find application as a secondary screen for selecting biomass lines and for optimizing biomass processing and conversion efficiencies.     

Molecular conformational space analysis using computer graphics: Going beyond FRODO

The molecular graphics program FRODO has been modified to support analytical animation of molecular dynamics trajectories. The enhanced program, mdFRODO, supports all features available in FRODO and is interfaced to GROMOS. A variety of analytical animation modes is included. Extensive coloring and atom selection features are implemented to aid the user in distinguishing features of interest in a set of conformations. Molecular conformational space can be analyzed efficiently and comprehended. Animations may be viewed in stereo, and the animated object can be overlaid with any of the standard FRODO objects. The mdFRODO program is of wide use in molecular dynamics, X-ray crystallography and two-dimensional NMR work. Examples illustrating various aspects of collective motion in protein molecules are given and discussed.     

A general Monte Carlo/simulated annealing algorithm for resonance assignment in NMR of uniformly labeled biopolymers

We describe a general computational approach to site-specific resonance assignments in multidimensional NMR studies of uniformly ¹⁵N,¹³C-labeled biopolymers, based on a simple Monte Carlo/simulated annealing (MCSA) algorithm contained in the program MCASSIGN2. Input to MCASSIGN2 includes lists of multidimensional signals in the NMR spectra with their possible residue-type assignments (which need not be unique), the biopolymer sequence, and a table that describes the connections that relate one signal list to another. As output, MCASSIGN2 produces a high-scoring sequential assignment of the multidimensional signals, using a score function that rewards good connections (i.e., agreement between relevant sets of chemical shifts in different signal lists) and penalizes bad connections, unassigned signals, and assignment gaps. Examination of a set of high-scoring assignments from a large number of independent runs allows one to determine whether a unique assignment exists for the entire sequence or parts thereof. We demonstrate the MCSA algorithm using two-dimensional (2D) and three-dimensional (3D) solid state NMR spectra of several model protein samples (α-spectrin SH3 domain and protein G/B1 microcrystals, HET-s_218–289 fibrils), obtained with magic-angle spinning and standard polarization transfer techniques. The MCSA algorithm and MCASSIGN2 program can accommodate arbitrary combinations of NMR spectra with arbitrary dimensionality, and can therefore be applied in many areas of solid state and solution NMR.     

Conformational studies of the O-specific polysaccharide of Shigella flexneri 5a and of four related synthetic pentasaccharide fragments using NMR and molecular modeling

As part of a program for the development of synthetic vaccines against the pathogen Shigella flexneri, we used a combination of NMR and molecular modeling methods to study the conformations of the O-specific polysaccharide (O-SP) of S. flexneri 5a and of four related synthetic pentasaccharide fragments. The NMR study, based on the analysis of 1H and 13C chemical shifts, the evaluation of inter-residue distances, and the measurement of one- and three-bond heteronuclear coupling constants, showed that the conformation of one of the four pentasaccharides is similar to that of the native O-SP in solution. Interestingly, inhibition enzyme-linked immunosorbent assay demonstrated that a protective monoclonal antibody specific for S. flexneri 5a has a greater affinity for this pentasaccharide than for the others. We carried out a complete conformational search on the pentasaccharides using the CICADA algorithm interfaced with MM3 force field. We calculated Boltzmann-averaged inter-residue distances and 3JC,H coupling constants for the different conformational families and compared the results with NMR data for all pentasaccharides. Our experimental data are consistent with only one conformational family. We also used molecular modeling data to build models of the O-SP with the molecular builder program POLYS. The models that are in agreement with NMR data adopt right-handed 3-fold helical structures in which the branched glucosyl residue points outwards.     

Spectral fitting for signal assignment and structural analysis of uniformly <Superscript>13</Superscript>C-labeled solid proteins by simulated annealing based on chemical shifts and spin dynamics

We describe an approach for the signal assignment and structural analysis with a suite of two-dimensional (13)C-(13)C magic-angle-spinning solid-state NMR spectra of uniformly (13)C-labeled peptides and proteins. We directly fit the calculated spectra to experimental ones by simulated annealing in restrained molecular dynamics program CNS as a function of atomic coordinates. The spectra are calculated from the conformation dependent chemical shift obtained with SHIFTX and the cross-peak intensities computed for recoupled dipolar interactions. This method was applied to a membrane-bound 14-residue peptide, mastoparan-X. The obtained C', C(alpha) and C(beta) chemical shifts agreed with those reported previously at the precisions of 0.2, 0.7 and 0.4 ppm, respectively. This spectral fitting program also provides backbone dihedral angles with a precision of about 50 degrees from the spectra even with resonance overlaps. The restraints on the angles were improved by applying protein database program TALOS to the obtained chemical shifts. The peptide structure provided by these restraints was consistent with the reported structure at the backbone RMSD of about 1 A.     

Conformational and model-building studies of the hairpin form of the mismatched DNA octamer d(m5C-G-m5C-G-T-G-m5C-G)

The hairpin form of the mismatched octamer d(m5C-G-m5C-G-T-G-m5C-G) was studied by means of NMR spectroscopy. In a companion study it is shown that the hairpin form of this DNA fragment consists of a structure with a stem of three Watson-Crick-type base pairs and a loop consisting of only two nucleotides. The non-exchangeable proton resonances were assigned by means of two-dimensional correlation spectroscopy and two-dimensional nuclear Overhauser effect spectroscopy. Proton-proton coupling constants were used for the conformational analysis of the deoxyribose ring and for some of the backbone torsion angles. From the two-dimensional NMR spectra and the coupling-constant analysis it is concluded that: (i) the stem of the hairpin exhibits B-DNA characteristics; (ii) the sugar rings are not conformationally pure, but display a certain amount of conformational flexibility; (iii) the stacking interaction in the stem of the hairpin is elongated from the 3'-side in a more or less regular fashion with the two loop nucleotides; (iv) at the 5'-side of the stem a stacking discontinuity occurs between the stem and the loop; (v) at the 5'-side of the stem the loop is closed by means of a sharp backbone turn which involves unusual gamma' and beta+ torsion angles in residue dG(6). The NMR results led to the construction of a hairpin-loop model which was energy-minimized by means of a molecular-mechanics program. The results clearly show that a DNA hairpin-loop structure in which the loop consists of only two nucleotides bridging the minor groove in a straightforward fashion, (i) causes no undue steric strain, and (ii) involves well-known conformational principles throughout the course of the backbone. The hairpin form of the title compound is compared with the hairpin form of d(A-T-C-C-T-A-T4-T-A-G-G-A-T), in which the central -T4- part forms a loop of four nucleotides. Both models display similarities as far as stacking interactions are concerned.     

The program XEASY for computer-supported NMR spectral analysis of biological macromolecules

Summary A new program package, XEASY, was written for interactive computer support of the analysis of NMR spectra for three-dimensional structure determination of biological macromolecules. XEASY was developed for work with 2D, 3D and 4D NMR data sets. It includes all the functions performed by the precursor program EASY, which was designed for the analysis of 2D NMR spectra, i.e., peak picking and support of sequence-specific resonance assignments, cross-peak assignments, cross-peak integration and rate constant determination for dynamic processes. Since the program utilizes the X-window system and the Motif widget set, it is portable on a wide range of UNIX workstations. The design objective was to provide maximal computer support for the analysis of spectra, while providing the user with complete control over the final resonance assignments. Technically important features of XEASY are the use and flexible visual display of strips, i.e., two-dimensional spectral regions that contain the relevant parts of 3D or 4D NMR spectra, automated sorting routines to narrow down the selection of strips that need to be interactively considered in a particular assignment step, a protocol of resonance assignments that can be used for reliable bookkeeping, independent of the assignment strategy used, and capabilities for proper treatment of spectral folding and efficient transfer of resonance assignments between spectra of different types and different dimensionality, including projected, reduced-dimensionality triple-resonance experiments.Abbreviations 1D, 2D, 3D, 4D one-, two-, three-, four-dimensional - NOE nuclear Overhauser enhancement - NOESY nuclear Overhauser enhancement spectroscopy - TOCSY total correlation spectroscopy - COSY correlation spectroscopy - TPPI time-proportional phase incrementation     

NMR-pseudoenergy approach to the solution structure of acyl carrier protein

A method for protein structure determination from two-dimensional NMR cross-relaxation data is presented and explored by using short amino acid segments from acyl carrier protein as a test case. The method is based on a molecular mechanics program and incorporates NMR distance constraints in the form of a pseudoenergy term that accurately reflects the distance-dependent precision of NMR cross-relaxation data. When it is used in an indiscriminant fashion, the method has a tendency to produce structures representing local energy minima near starting structures, rather than structures representing a global energy minimum. However, stepwise inclusion of energy terms, beginning with a function heavily weighted by backbone distance constraints, appears to simplify the potential energy surface to a point where convergence to a common backbone structure from a variety of starting structures is possible. In the case of the segment from residues 3 to 15 in acyl carrier protein, a nearly perfect alpha-helix is produced starting with a linear chain, an alpha-helical chain, or a chain having residues with alternating linear and alpha-helical backbone torsional angles. In the case of the segment from residues 26 to 36 a structure having a right-handed loop is produced.     

Combined use of COSY and double quantum two-dimensional NMR spectroscopy for elucidation of spin systems in polymyxin B

Two-dimensional single quantum correlation NMR spectroscopy (COSY) and two-dimensional double quantum NMR spectroscopy (2QT) are used to study spin systems in the 1H NMR spectrum of polymyxin B. Because of different frequency relationships, the two types of two-dimensional NMR experiments are found to be highly complementary. This is demonstrated by combined use of COSY and 2QT spectroscopy to obtain a complete analysis of the complicated spectral overlap which occurs in the 1H NMR spectrum of polymyxin B.     

AURELIA,a program for computer-aided analysis of multidimensional NMR spectra

Summary AURELIA is an advanced program for the computer-aided evaluation of two-, three- and four-dimensional NMR spectra of any type of molecule. It can be used for the analysis of spectra of small molecules as well as for evaluation of complicated spectra of biological macromolecules such as proteins. AURELIA is highly interactive and offers a large number of tools, such as artefact reduction, cluster and multiplet analysis, spin system searches, resonance assignments, automated calculation of volumes in multidimensional spectra, calculation of distances with different approaches, including the full relaxation matrix approach, Bayesian analysis of peak features, correlation of molecular structures with NMR data, comparison of spectra via spectral algebra and pattern match techniques, automated sequential assignments on the basis of triple resonance spectra, and automatic strip calculation. In contrast to most other programs, many tasks are performed automatically.     

Complete H and C NMR chemical shift assignments of mono-, di-, and trisaccharides as basis for NMR chemical shift predictions of polysaccharides using the computer program casper

The computer program casper uses ¹H and ¹³C NMR chemical shift data of mono- to trisaccharides for the prediction of chemical shifts of oligo- and polysaccharides. In order to improve the quality of these predictions the ¹H and ¹³C, as well as ³¹P when applicable, NMR chemical shifts of 30 mono-, di-, and trisaccharides were assigned. The reducing sugars gave two distinct sets of NMR resonances due to the α- and β-anomeric forms. In total 35 ¹H and ¹³C NMR chemical shift data sets were obtained from the oligosaccharides. One- and two-dimensional NMR experiments were used for the chemical shift assignments and special techniques were employed in some cases such as 2D ¹H,¹³C-HSQC Hadamard Transform methodology which was acquired approximately 45 times faster than a regular t₁ incremented ¹H,¹³C-HSQC experiment and a 1D ¹H,¹H-CSSF-TOCSY experiment which was able to distinguish spin-systems in which the target protons were only 3.3 Hz apart. The ¹H NMR chemical shifts were subsequently refined using total line-shape analysis with the PERCH NMR software. The acquired NMR data were then utilized in the casper program (http://www.casper.organ.su.se/casper/) for NMR chemical shift predictions of the O-antigen polysaccharides from Klebsiella O5, Shigella flexneri serotype X, and Salmonella arizonae O62. The data were compared to experimental data of the polysaccharides from the two former strains and the lipopolysaccharide of the latter strain showing excellent agreement between predicted and experimental ¹H and ¹³C NMR chemical shifts.     

Prediction of protein secondary structure based on physical theory. Histones

Secondary structures of histones H1, H2A, H2B, H3, H4 and H5 have been calculated by the computer program ALB based on a molecular theory of protein secondary structure. The predicted secondary structures of all histones are predominantly alpha-helical. The calculated secondary structure of linker histones H1 and H5 is close to that previously obtained from two-dimensional NMR data. For each of the core histones (H2A, H2B, H3, H4) one long alpha-helix and several short ones have been predicted. These long helices can be identified with rods in the low-resolution electron density map.     

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.