Design and implementation of a qualitative simulation model of lambda phage infection

MOTIVATION: Molecular biology databases hold a large number of empirical facts about many different aspects of biological entities. That data is static in the sense that one cannot ask a database 'What effect has protein A on gene B?' or 'Do gene A and gene B interact, and if so, how?'. Those questions require an explicit model of the target organism. Traditionally, biochemical systems are modelled using kinetics and differential equations in a quantitative simulator. For many biological processes however, detailed quantitative information is not available, only qualitative or fuzzy statements about the nature of interactions. RESULTS: We designed and implemented a qualitative simulation model of lambda phage growth control in Escherichia coli based on the existing simulation environment QSim. Qualitative reasoning can serve as the basis for automatic transformation of contents of genomic databases into interactive modelling systems that can reason about the relations and interactions of biological entities.      

A sample preparation protocol for H nuclear magnetic resonance studies of water-soluble metabolites in blood and urine

We describe a general protocol for preparing protein-containing biofluids for ¹H nuclear magnetic resonance (NMR) metabolomic studies. In this protocol, untreated samples are diluted in deuterated solvents to precipitate proteins and recover metabolites quantitated relative to standard reference compounds such as 3-trimethylsilylpropionic acid (TSP) and 2,2-dimethyl-2-silapentane-5-sulfonic acid (DSS). The efficacy of this protocol was tested using a bovine serum albumin/metabolite mix and human serum samples. This sample preparation method can be readily applied to any protein-containing biofluid for ¹H NMR studies.     

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.     

Covariance NMR in higher dimensions: application to 4D NOESY spectroscopy of proteins

Elucidation of high-resolution protein structures by NMR spectroscopy requires a large number of distance constraints that are derived from nuclear Overhauser effects between protons (NOEs). Due to the high level of spectral overlap encountered in 2D NMR spectra of proteins, the measurement of high quality distance constraints requires higher dimensional NMR experiments. Although four-dimensional Fourier transform (FT) NMR experiments can provide the necessary kind of spectral information, the associated measurement times are often prohibitively long. Covariance NMR spectroscopy yields 2D spectra that exhibit along the indirect frequency dimension the same high resolution as along the direct dimension using minimal measurement time. The generalization of covariance NMR to 4D NMR spectroscopy presented here exploits the inherent symmetry of certain 4D NMR experiments and utilizes the trace metric between donor planes for the construction of a high-resolution spectral covariance matrix. The approach is demonstrated for a 4D (13)C-edited NOESY experiment of ubiquitin. The 4D covariance spectrum narrows the line-widths of peaks strongly broadened in the FT spectrum due to the necessarily short number of increments collected, and it resolves otherwise overlapped cross peaks allowing for an increase in the number of NOE assignments to be made from a given dataset. At the same time there is no significant decrease in the positive predictive value of observing a peak as compared to the corresponding 4D Fourier transform spectrum. These properties make the 4D covariance method a potentially valuable tool for the structure determination of larger proteins and for high-throughput applications in structural biology.     

Determining molecular binding sites on human serum albumin by displacement of oleic acid

An NMR method was developed for determining binding sites of small molecules on human serum albumin (HSA) by competitive displacement of (13)C-labeled oleic acid. This method is based on the observation that in the crystal structure of HSA complexed with oleic acid, two principal drug-binding sites, Sudlow's sites I (warfarin) and II (ibuprofen), are also occupied by fatty acids. In two-dimensional [(1)H,(13)C]heteronuclear single quantum coherence NMR spectra, seven distinct resonances were observed for the (13)C-methyl-labeled oleic acid as a result of its binding to HSA. Resonances corresponding to the major drug-binding sites were identified through competitive displacement of molecules that bind specifically to each site. Thus, binding of molecules to these sites can be followed by their displacement of oleic acids. Furthermore, the amount of bound ligand at each site can be determined from changes in resonance intensities. For molecules containing fluorine, binding results were further validated by direct observations of the bound ligands using (19)F NMR. Identifying the binding sites for drug molecules on HSA can aid in determining the structure-activity relationship of albumin binding and assist in the design of molecules with altered albumin binding.     

Refolding out of guanidine hydrochloride is an effective approach for high-throughput structural studies of small proteins

Low in vivo solubility of recombinant proteins expressed in Escherichia coli can seriously hinder the purification of structural samples for large-scale proteomic NMR and X-ray crystallography studies. Previous results from our laboratory have shown that up to one half of all bacterial and archaeal proteins are insoluble when overexpressed in E. coli. Although a number of strategies may be used to increase in vivo protein solubility, there are no generally applicable methods, and the expression of each insoluble recombinant protein must be individually optimized. For this reason, we have tested a generic denaturation/refolding protein purification procedure to assess the number of structural samples that could be generated by using this methodology. Our results show that a denaturation/refolding protocol is appropriate for many small proteins (相似文献   

SOACS index: an easy NMR-based query for glycan retrieval

¹H NMR is now a standard method to determine de novo primary sequence of all sorts of glycans. These last 30 years, tens of thousands of oligosaccharide sequences have been elucidated by NMR spectroscopy in conjunction with other physico-chemical methods including mass spectrometry and gas chromatography. Most of these sequences are now compiled and available in several web databases recently unified in publicly available GlycomeDB, along with sets of experimental data. However, because the search for an exact sequence exclusively based on proton chemical shifts is sometimes delicate for NMR non-specialists, we worked out a new type of query, named SOACS, which allows the easy retrieval of existing sequences. This query is based on the readily distinguished ¹H chemical shifts from any ¹H NMR spectrum, and was designed to be usable to the widest scientist community.     

Complex glycosaminoglycans: profiling substitution patterns by two-dimensional nuclear magnetic resonance spectroscopy

Biological and pharmacological interactions of heparin and structurally related glycosaminoglycans (GAGs) such as heparan sulfate (HS) involve complex sequences of variously sulfated uronic acid and aminosugar residues. Due to their structural microheterogeneity, these sequences are usually characterized in statistical terms, by high-performance liquid chromatographic analysis of fragments obtained by enzymatic or chemical degradation. Nuclear magnetic resonance (NMR) spectroscopy is also currently used for structural characterization of GAGs. However, the use of monodimensional NMR analysis of complex GAGs is often limited by severe signal overlap that does not allow reliable quantitative measurements. Using magnetically equivalent signals, the higher resolution achieved by two-dimensional NMR methods could be also exploited for quantitative applications. In this work, heteronuclear single quantum coherence (HSQC) spectroscopy has been evaluated to determine variously substituted monosaccharide components of HS and HS mimics obtained by chemical modification of the Escherichia coli K5 polysaccharide (K5-PS) structurally related to the common biosynthetic precursor of heparin and HS. Heparin was used as a model for assessing the influence of 1H-13C spin-spin couplings on "volumes" of the corresponding signals. For major signals, the HSQC approach permitted quantification of additional structural features both in heparins and in a typical HS. The method was applied to profile the substitution patterns of K5-PS derivatives involving different degrees of N,O-sulfation and N-acetylation, including O-sulfated heparosans bearing free amino groups.     

PR-CALC: A Program for the Reconstruction of NMR Spectra from Projections

Projection-reconstruction NMR (PR-NMR) has attracted growing attention as a method for collecting multidimensional NMR data rapidly. The PR-NMR procedure involves measuring lower-dimensional projections of a higher-dimensional spectrum, which are then used for the mathematical reconstruction of the full spectrum. We describe here the program PR-CALC, for the reconstruction of NMR spectra from projection data. This program implements a number of reconstruction algorithms, highly optimized to achieve maximal performance, and manages the reconstruction process automatically, producing either full spectra or subsets, such as regions or slices, as requested. The ability to obtain subsets allows large spectra to be analyzed by reconstructing and examining only those subsets containing peaks, offering considerable savings in processing time and storage space. PR-CALC is straightforward to use, and integrates directly into the conventional pipeline for data processing and analysis. It was written in standard C+ + and should run on any platform. The organization is flexible, and permits easy extension of capabilities, as well as reuse in new software. PR-CALC should facilitate the widespread utilization of PR-NMR in biomedical research. Electronic supplementary material Electronic supplementary material is available for this article at and accessible for authorised users.     

One- and two-dimensional NMR spectral analysis of the consequences of single amino acid replacements in proteins

The traditional approach of using homologous sequences to elucidate the role of specific amino acid residues in protein structure and function becomes more meaningful as the number of differences is minimized, with the limit being alteration of a single residue. For small proteins in solution, NMR spectroscopy offers a means of obtaining detailed information about each residue and its response to a given change in the protein sequence. Extraction of this information has been aided by recent progress in spectrometer technology (higher magnetic fields, more sensitive signal detection, more sophisticated computers) and experimental strategies (new NMR pulse sequences including multiple-quantum and two-dimensional NMR methods). The set of avian ovomucoid third domains, which consists of the third domain proper plus a short leader (connecting peptide) and has a maximum of 56 amino acid residues, offers an attractive system for developing experimental methods for investigating sequence-structure and structure-function relationships in proteins. Our NMR results provide examples of sequence effects on pKa' values, average conformation, and internal motion of amino acid side chains.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

High resolution 4-D spectroscopy with sparse concentric shell sampling and FFT-CLEAN

Recent efforts to reduce the measurement time for multidimensional NMR experiments have fostered the development of a variety of new procedures for sampling and data processing. We recently described concentric ring sampling for 3-D NMR experiments, which is superior to radial sampling as input for processing by a multidimensional discrete Fourier transform. Here, we report the extension of this approach to 4-D spectroscopy as Randomized Concentric Shell Sampling (RCSS), where sampling points for the indirect dimensions are positioned on concentric shells, and where random rotations in the angular space are used to avoid coherent artifacts. With simulations, we show that RCSS produces a very low level of artifacts, even with a very limited number of sampling points. The RCSS sampling patterns can be adapted to fine rectangular grids to permit use of the Fast Fourier Transform in data processing, without an apparent increase in the artifact level. These artifacts can be further reduced to the noise level using the iterative CLEAN algorithm developed in radioastronomy. We demonstrate these methods on the high resolution 4-D HCCH-TOCSY spectrum of protein G's B1 domain, using only 1.2% of the sampling that would be needed conventionally for this resolution. The use of a multidimensional FFT instead of the slow DFT for initial data processing and for subsequent CLEAN significantly reduces the calculation time, yielding an artifact level that is on par with the level of the true spectral noise.     

Online NMR for monitoring biocatalysed reactions

Monitoring biocatalysed reactions and metabolic pathways using NMR spectroscopy is of growing interest. As a non-invasive analytical method providing simultaneous information about intracellular and extracellular constituents, it is superior to other analytical techniques and has a wide range of applications: kinetics and stoichiometrics of metabolic events, metabolic fluxes and enzyme activities can be detected in situ or after taking a sample from the biotransformation mixture. New NMR pulse sequences provide even more valuable experiments in these fields. Research topics range from the monitoring of polymer formation to fermentations producing beverages or antibiotics. Routine monitoring of industrial fermentations by NMR seems to be imminent.     

Double quantum filtering homonuclear MAS NMR correlation spectra: a tool for membrane protein studies

13C homonuclear correlation spectra based on proton driven spin diffusion (PDSD) are becoming increasingly important for obtaining distance constraints from multiply labeled biomolecules by MAS NMR. One particular challenging situation arises when such constraints are to be obtained from spectra with a large natural abundance signal background which causes detrimental diagonal peak intensities. They obscure cross peaks, and furthermore impede the calculation of a buildup rates matrix which may be used to derive distance constraints, as carried out in "NMR crystallography". Here, we combine double quantum (DQ) filtering with 13C-13C dipolar assisted rotational resonance (DARR) experiments to yield correlation spectra free of natural abundance contributions. Two experimental schemes, using DQ filtering prior to evolution (DOPE), and after mixing (DOAM), have been evaluated. Diagonal peak intensities along the spectrum diagonal are removed completely, and crosspeaks close to the diagonal are easily identifiable. For DOAM spectra with negligible mixing times, it is possible to carry out 'assignment walks' which simplify peak identification substantially. The method is demonstrated on 13C-cys labeled proteorhodopsin, a 27 kDa membrane protein. The magnetization transfer characteristics were studied using buildup curves obtained on uniformly 13C labelled crystalline tripeptide MLF. Our data show that DQ filtered DARR experiments pave the way for obtaining through space constraints for structural studies on ligands, bound to membrane receptors, or on small fragments within large proteins.     

Modification of carboplatin by Jurkat cells

Using [(1)H,(15)N] heteronuclear single quantum coherance (HSQC) NMR and (15)N-labeled carboplatin, 1, we show that Jurkat cells affect the rate of disappearance of the HSQC NMR peak in culture medium for this Pt(2+) anticancer drug. The decay or disappearance rate constant for 1 in culture medium containing cells is k(1)=k(c)[CO(3)(2-)]+k(m)+k(u)N, where k(c) is the rate constant for reaction of 1 with carbonate in the medium, k(m) is the rate constant for reaction of 1 with all other components of the medium, and k(u) is the rate constant for reaction of 1 with cells having a number density N in the medium. Since Jurkat cells only take up a small amount of the platinum present in the medium (<1%), the observed disappearance of the HSQC NMR peak for 1 cannot be due to uptake of carboplatin by the cells.     

NMR spectroscopic filtration of polypeptides and proteins in complex mixtures

Due to the inherent complexity of the natural biological environment, most studies on polypeptides, proteins and nucleic acids have so far been performed in vitro, away from physiologically relevant conditions. Nuclear magnetic resonance is an ideal technique to extend the in vitro analysis of simple model systems to the more complex biological context. This work shows how diffusion-based spectroscopic selection can be combined with isotopic labeling to tackle and optimize the NMR analysis of specific macromolecules in multicomponent mixtures. Typical media include cell-free systems containing overexpressed proteins, lysates and proteolytic mixtures. We present a few variants of diffusion-edited HSQC pulse sequences for the selective spectroscopic detection of protein and polypeptide resonances within complex mixtures containing undesired species of smaller molecular weight. Due to diffusion-based filtering, peak intensities of fast diffusing small molecules are attenuated more than peaks due to large molecules. The basic sequence, denoted as PFGSTE-HSQC, combines translational diffusion-ordering with two dimensional heteronuclear single quantum correlation spectroscopy. The GCSTE-HSQC and BPPSTE-HSQC sequences include bipolar gradients and are therefore suitable for both diffusion-based filtering and determination of diffusion coefficients of individual mixture components. Practical applications range from protein stability/folding investigations in physiologically relevant contexts to prescreening of tertiary fold and resonance assignments in structural genomics studies. A few applications of diffusion-edited HSQC to an E. coli cell lysate containing the (15)N-labeled B domain of streptococcal protein G (GB1), and to a (15)N-labeled N-acetylglycine/apomyoglobin mixture are presented. In addition, we provide specific guidelines for experimental setup and parameter optimization.     

Analysis of peptide rotational diffusion by homonuclear NMR

The analysis of the rotational diffusion of a molecule using homonuclear NMR is investigated. The homonuclear longitudinal and transverse cross-relaxation rates, which can be quantitatively measured using off-Resonance Rotating frame nuclear Overhauser Effect Spectroscopy (ROESY), are used to build a distribution, which exhibits a solid-state-like pattern characteristic of the diffusion tensor. The distributions of the antimicrobial peptide ranalexin in water and in 30% of trifluoracetic acid (TFE) are compared, and the peptide rotational diffusion is shown to be more isotropic in water than in 30% TFE. This difference is further supported by the analysis of NMR ranalexin conformers in 30% TFE, and by the analysis of a molecular dynamics simulation of peptide in water.     

Sensitivity enhancement in NMR of macromolecules by application of optimal control theory

NMR of macromolecules is limited by large transverse relaxation rates. In practice, this results in low efficiency of coherence transfer steps in multidimensional NMR experiments, leading to poor sensitivity and long acquisition times. The efficiency of coherence transfer can be maximized by design of relaxation optimized pulse sequences using tools from optimal control theory. In this paper, we demonstrate that this approach can be adopted for studies of large biological systems, such as the 800 kDa chaperone GroEL. For this system, the ¹H–¹⁵N coherence transfer module presented here yields an average sensitivity enhancement of 20–25% for cross-correlated relaxation induced polarization transfer (CRIPT) experiments.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1007/s10858-005-3592-0