Methods for sequential resonance assignment in solid, uniformly 13C, 15N labelled peptides: quantification and application to antamanide

The application of adiabatic polarization-transfer experiments to resonance assignment in solid, uniformly ¹³C-¹⁵N-labelled polypeptides is demonstrated for the cyclic decapeptide antamanide. A homonuclear correlation experiment employing the DREAM sequence for adiabatic dipolar transfer yields a complete assignment of the C and aliphatic side-chain ¹³C resonances to amino acid types. The same information can be obtained from a TOBSY experiment using the recently introduced P9¹ ₁₂ TOBSY sequence, which employs the J couplings as a transfer mechanism. A comparison of the two methods is presented. Except for some aromatic phenylalanine resonances, a complete sequence-specific assignment of the ¹³C and ¹⁵N resonances in antamanide is achieved by a series of selective or broadband adiabatic triple-resonance experiments. Heteronuclear transfer by adiabatic-passage Hartmann–Hahn cross polarization is combined with adiabatic homonuclear transfer by the DREAM and rotational-resonance tickling sequences into two- and three-dimensional experiments. The performance of these experiments is evaluated quantitatively.     

13C, 15N Resonance Assignment of Parts of the HET-s Prion Protein in its Amyloid Form

The partial ¹⁵N and ¹³C solid-state NMR resonance assignment of the HET-s prion protein fragment 218–289 in its amyloid form is presented. It is based on experiments measured at MAS frequencies in the range of 20–40 kHz using exclusively adiabatic polarization-transfer schemes. The resonance assignment within each residue is based on two-dimensional ¹³C––¹³C correlation spectra utilizing the DREAM mixing scheme. The sequential linking of the assigned residues used a set of two- and three-dimensional ¹⁵N––¹³C correlation experiments. Almost all cross peaks visible in the spectra are assigned, but only resonances from 43 of the 78 amino-acid residues could be detected. The missing residues are thought to be highly disordered and/or highly dynamic giving rise to broad resonance lines that escaped detection in the experiments applied. The line widths of the observed resonances are narrow and comparable to line widths observed in micro-crystalline samples. The 43 assigned residues are located in two fragments of about 20 residues. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

The Database for Reaching Experiments and Models

Reaching is one of the central experimental paradigms in the field of motor control, and many computational models of reaching have been published. While most of these models try to explain subject data (such as movement kinematics, reaching performance, forces, etc.) from only a single experiment, distinct experiments often share experimental conditions and record similar kinematics. This suggests that reaching models could be applied to (and falsified by) multiple experiments. However, using multiple datasets is difficult because experimental data formats vary widely. Standardizing data formats promises to enable scientists to test model predictions against many experiments and to compare experimental results across labs. Here we report on the development of a new resource available to scientists: a database of reaching called the Database for Reaching Experiments And Models (DREAM). DREAM collects both experimental datasets and models and facilitates their comparison by standardizing formats. The DREAM project promises to be useful for experimentalists who want to understand how their data relates to models, for modelers who want to test their theories, and for educators who want to help students better understand reaching experiments, models, and data analysis.     

Principles of quantitative analysis of two-dimensional spectra of nuclear Overhauser effect in evaluation of protein and peptide conformation

To elucidate potentialities of two-dimensional homonuclear Overhauser effect (NOESY) spectra of peptides and proteins for their spatial structure determination, impact of experimental parameters and intrinsic properties of the investigated molecule on proton cross-peak volumes in NOESY spectra was analysed. Recommendations which could increase accuracy of cross-peak volume measurements were suggested. Influence of intrinsic properties of a molecule (spin-lattice relaxation times T1, correlation time tau C and surrounding protons) on the volume of cross-peak for particular protons was analyzed using a complete relaxation matrix of the (formula; see text) helix of gramicidin A. Nonselective relaxation time T1 of the protons was found to affect only slightly the results of cross-peak volumes computer simulation, whereas correlation time tau C and surrounding protons seriously influenced cross-peak volumes. Nevertheless, cross-peak volumes between NH, C alpha H and C beta H protons of a dipeptide fragment of the entire molecule could be accurately simulated using the relaxation matrix of the individual dipeptide. Thus local conformations (torsion angles phi, psi and chi 1) of amino acid residues could be deduced independently of one another and prior to the complete analysis of a molecular structure. The result can be obtained even in the presence of spin-diffusion at mixing times providing maximal volumes of cross-peaks in NOESY spectra.     

3d Trosy-HncaCodedcb and Trosy-HncaCodedco Experiments: Triple Resonance nmr Experiments With two Sequential Connectivity Pathways and High Sensitivity

The concept of chemical shift-coding monitors chemical shifts in multi-dimensional NMR experiments without additional polarization transfer elements and evolution periods. The chemical shifts are coded in the line-shape of the cross-peak through an apparent scalar coupling dependent upon the chemical shift. This concept is applied to the three-dimensional triple-resonance experiment HNCA adding the information of (13)C(beta) or (13)C' chemical shifts. On average, the proposed TROSY-HNCA(coded)CB experiment is a factor of 2 less sensitive than the HNCA experiment. However, it contains correlations via the chemical shifts of both (13)C(alpha) and (13)C(beta), and provides up to three times higher resolution along the (13)C(alpha) chemical shift axis. Thus, it dramatically reduces ambiguities in linking the spin systems of adjacent residues in the protein sequence during the sequential assignment. The TROSY-HNCA(coded)CO experiment which is conceptually similar contains correlations via the chemical shifts of (13)C(alpha) and (13)C' without major signal losses. The proposed triple resonance experiments are applied to a approximately 70% (2)H, approximately 85% (13)C,(15)N labeled protein with a molecular weight of 25 kDa.     

Coherence transfer in deoxyribose sugars produced by isotropic mixing: an improved intraresidue assignment strategy for the two-dimensional NMR spectra of DNA

Two-dimensional isotropic mixing spectroscopy has been used to confirm assignments of the deoxyribose sugar protons in the H NMR spectrum of the DNA oligonucleotides d(CGCGTTTTCGCG) and [d(GCCGTGGCCACGGC)]2. The broad-band decoupling sequence MLEV-16 was applied during the mixing period to induce isotropic coupling within the spin systems, resulting in net transfer of coherence throughout the coupled spin networks. Nearly all 1', 2', 2', 3', and 4' protons of a given nucleotide could be identified on the basis of through-bond scalar connectivities. In addition, in the hairpin, a number of connectivities to 5'/5' protons were found. The dependence of cross-peak intensity on the length of radio-frequency irradiation for several different coherence transfer orders is presented, and implications for optimization are discussed.     

Protein–nucleotide contacts in motor proteins detected by DNP-enhanced solid-state NMR

DNP (dynamic nuclear polarization)-enhanced solid-state NMR is employed to directly detect protein–DNA and protein–ATP interactions and identify the residue type establishing the intermolecular contacts. While conventional solid-state NMR can detect protein–DNA interactions in large oligomeric protein assemblies in favorable cases, it typically suffers from low signal-to-noise ratios. We show here, for the oligomeric DnaB helicase from Helicobacter pylori complexed with ADP and single-stranded DNA, that this limitation can be overcome by using DNP-enhanced spectroscopy. Interactions are established by DNP-enhanced ³¹P–¹³C polarization-transfer experiments followed by the recording of a 2D ¹³C–¹³C correlation experiment. The NMR spectra were obtained in less than 2 days and allowed the identification of residues of the motor protein involved in nucleotide binding.     

High pressure NMR reveals a variety of fluctuating conformers in beta-lactoglobulin

High pressure 1H/15N two-dimensional NMR spectroscopy has been used to study conformational fluctuation in bovine beta-lactoglobulin at pH 2.0 and 36 degrees C. Pressure dependencies of 1H and 15N chemical shifts and cross-peak intensities were analyzed at more than 80 independent atom sites between 30 and 2000 bar. Unusually large and non-linear chemical shift pressure dependencies are found for residues centering in the hydrophobic core region, suggesting the existence of low-lying excited native states (N') of the protein. Measurement of 1H/15N cross-peak intensities at individual amide sites as a function of pressure suggests that unfolding events occur independently in two sides of the beta-barrel, i.e. the hydrophobic core side (betaF-H) (producing I2) and the non-core side (betaB-E) (producing I1). At 1 bar the stability is higher for the core region (DeltaG0 = 6.5(+/-2.0) kcal/mol) than for the non-core region (4.6(+/-1.3) kcal/mol), but at high pressure the stability is reversed due to a larger DeltaV value of unfolding for the core region (90.0(+/-35.2) ml/mol) than that for the non-core region (57.4(+/-14.4) ml/mol), possibly due to an uneven distribution of cavities. The DeltaG0 profile along the amino acid sequence obtained from the pressure experiment is found to coincide well with that estimated from hydrogen exchange experiments. Altogether, the high pressure NMR experiment has revealed a variety of fluctuating conformers of beta-lactoglobulin, notably N, N', I1, I2 and the totally unfolded conformer U. Fluctuation of N to I1 and I2 conformers with open barrel structures could be a common design of lipocalin family proteins which bind various hydrophobic compounds in its barrel structure.     

针对DREAM 基因序列外显子的siRNA 设计

目的:针对DREAM基因中外显子序列,设计并筛选出起作用的siRNA,为进一步研究以DREAM为靶标的基因治疗提供依据。方法:通过利用生物信息学的方法设计出66条潜在的针对DREAM基因中外显子序列的siRNA序列。潜在的序列根据G+C含量分析及Genbank BLAST检测,我们筛选了三条较理想的DREAM基因siRNA靶标,并将其构建到pENTR/H1/TO载体中,转染细胞后提取总蛋白并用Western blot方法检测DREAM蛋白的表达水平。结果:示外显子9中5’末端位置为1253的序列能够抑制80%的DREAM的蛋白表达。结论:这一研究结果为进一步siRNA类药物的实验研究提供了理论基础。     

Improved TROSY-HNCA experiment with suppression of conformational exchange induced relaxation

A general method for improving of the sensitivity of the TROSY-type triple resonance experiments in the presence of conformational exchange-induced (CSX) relaxation is proposed based on the use of CPMG-INEPT (Müller et al., J. Am. Chem. Soc., 1995, 117, 11043–11048) during the N–C polarization transfer periods. Significantly improved sensitivity is demonstrated for the majority of cross-peaks in the new [¹⁵N,¹H]-TROSY-XY-HNCA experiment, measured with partially folded RNase AS-Protein, with negligible loss of sensitivity for resonances unaffected by CSX relaxation. In addition, a comparison of cross-peak amplitudes in [¹⁵N,¹N]-TROSY-XY-HNCA and conventional [¹⁵N,¹H]-TROSY-HNCA spectra provides a quick and sensitive estimation of the CSX relaxation contribution.     

Refinement of the spatial structure of the gramicidin A ion channel]

The spatial structure of the gramicidin A (GA) transmembrane ion-channel was refined on the base of cross-peak volumes measured in NOESY spectra (mixing time tau m = 100 and 200 ms). The refinement methods included the comparison of experimental cross-peak volumes with those calculated for low-energy GA conformations, dynamic averaging of the low-energy conformation set and restrained energy minimization. Accuracy of the spatial structure determination was estimated by the penalty function Fr defined as a root mean square deviation of interproton distances corresponding to the calculated and experimental cross-peak volumes. As the initial conformation we used the right-handed pi 6,3 LD pi 6,3 LD helix established on the base of NMR data regardless of the cross-peak volumes. The conformation is in a good agreement with NOE cross-peak volumes (Fr 0.2 to 0.5 A depending on NOESY spectrum). For a number of NOEs formed by the side chain protons, distances errors were found as much as 0.5-2.0 A. Restrained energy minimization procedure had little further success. However some of these errors were eliminated by the change in torsional angle chi 2 of D-Leu12 and dynamic averaging of the Val7 side chain conformations. Apparently, majority of deviations of the calculated and experimental cross-peak volumes are due to the intramolecular mobility of GA and cannot be eliminated within the framework of rigid globule model. In summary the spatial structure of GA ion-channel can be thought as a set of low-energy conformations, differing by the side chain torsion angles chi 1 Val7 and chi 2 D-Leu4 and D-Leu10 and the orientation of the C-terminal ethanolamine group. Root mean square differences between the atomic coordinates of conformations are in the range of 0.3-0.8 A.     

Computationally efficient gradients for relaxation matrix-based structure refinement including the accommodation of internal motions

Summary A general method for deriving analytical gradients based on NOESY cross-peak intensities is presented. This method allows for very rapid calculation of exact gradients of cross-peak intensities with respect to parameters directly related to the pairwise dipole-dipole interaction giving rise to the cross peak. In the simplest case, gradients with respect to internuclear separation can be calculated, which allows for the rational modification of distance constraints used in structural refinement. In the general case, any arbitrary level of knowledge of the internal dynamics of the molecule can be introduced, thereby providing a pathway for the experimental determination of motional parameters. The motional characteristics of the internuclear vectors defining dipole-dipole interactions are cast in model-free terms. The form of the gradient circumvents many of the limitations of gradients expressed in terms of Cartesian or dihedral variables. The gradients presented are simple, direct, and exact and require little computational effort to calculate.     

Omega 1-decoupled 1H homonuclear shift-correlated nuclear magnetic resonance spectroscopy (COSYDEC) applied to steroids

Problems of cross-peak overlap in two-dimensional 1H homonuclear shift-correlated (COSY) spectra of steroids can often be avoided by use of the omega 1-decoupled COSY (COSYDEC) method. The selection of experimental parameters is discussed, and COSYDEC spectra are illustrated for 17a-oxa-D-homoandrost-4-ene-3,17-dione (testololactone), testosterone, and 17 alpha-hydroxyprogesterone. In a good case, a COSYDEC spectrum obtained at 250 MHz allows cross-peak recognition and assignment with facility comparable to that available only at 500 MHz for normal COSY spectra.     

A relative variation-based method to unraveling gene regulatory networks

Gene regulatory network (GRN) reconstruction is essential in understanding the functioning and pathology of a biological system. Extensive models and algorithms have been developed to unravel a GRN. The DREAM project aims to clarify both advantages and disadvantages of these methods from an application viewpoint. An interesting yet surprising observation is that compared with complicated methods like those based on nonlinear differential equations, etc., methods based on a simple statistics, such as the so-called Z-score, usually perform better. A fundamental problem with the Z-score, however, is that direct and indirect regulations can not be easily distinguished. To overcome this drawback, a relative expression level variation (RELV) based GRN inference algorithm is suggested in this paper, which consists of three major steps. Firstly, on the basis of wild type and single gene knockout/knockdown experimental data, the magnitude of RELV of a gene is estimated. Secondly, probability for the existence of a direct regulation from a perturbed gene to a measured gene is estimated, which is further utilized to estimate whether a gene can be regulated by other genes. Finally, the normalized RELVs are modified to make genes with an estimated zero in-degree have smaller RELVs in magnitude than the other genes, which is used afterwards in queuing possibilities of the existence of direct regulations among genes and therefore leads to an estimate on the GRN topology. This method can in principle avoid the so-called cascade errors under certain situations. Computational results with the Size 100 sub-challenges of DREAM3 and DREAM4 show that, compared with the Z-score based method, prediction performances can be substantially improved, especially the AUPR specification. Moreover, it can even outperform the best team of both DREAM3 and DREAM4. Furthermore, the high precision of the obtained most reliable predictions shows that the suggested algorithm may be very helpful in guiding biological experiment designs.     

Experiments for correlating quaternary carbons in RNA bases

The paper presents a set of triple-resonance two-dimensional experiments for correlating all quaternary carbons in RNA bases to one or more of the base protons. The experiments make use of either three-bond proton-carbon couplings and one selective INEPT step (the long-range selective HSQC experiment) to transfer the magnetization between a proton and the carbon of interest and back, or they rely on one- and/or two-bond heteronuclear (the H(CN)C and H(N)C experiments) or carbon-carbon (the H(C)C experiment) couplings and multiple INEPT transfer steps. The effect of the large one-bond carbon-carbon coupling in t(1) is removed by a constant time evolution or by a selective refocusing. The performance of the proposed approach is demonstrated on a 0.5 mM 25-mer RNA. The results show that the experiments are applicable to samples containing agents for weak molecular alignment. The design of the correlation experiments has been supported by ab initio calculations of scalar spin-spin couplings in the free bases and the AU and GC base pairs. The ab initio data reveal surprisingly high values of guanine (2)J(N1C5) and uracil (2)J(N3C5) couplings that are in a qualitative agreement with the experimental data. The sensitivity of the spin-spin couplings to base pairing as well as the agreement with the experiment depend strongly on the type of nuclei involved and the number of bonds separating them.     

Gene Regulatory Network Inference from Multifactorial Perturbation Data Using both Regression and Correlation Analyses

An important problem in systems biology is to reconstruct gene regulatory networks (GRNs) from experimental data and other a priori information. The DREAM project offers some types of experimental data, such as knockout data, knockdown data, time series data, etc. Among them, multifactorial perturbation data are easier and less expensive to obtain than other types of experimental data and are thus more common in practice. In this article, a new algorithm is presented for the inference of GRNs using the DREAM4 multifactorial perturbation data. The GRN inference problem among genes is decomposed into different regression problems. In each of the regression problems, the expression level of a target gene is predicted solely from the expression level of a potential regulation gene. For different potential regulation genes, different weights for a specific target gene are constructed by using the sum of squared residuals and the Pearson correlation coefficient. Then these weights are normalized to reflect effort differences of regulating distinct genes. By appropriately choosing the parameters of the power law, we constructe a 0–1 integer programming problem. By solving this problem, direct regulation genes for an arbitrary gene can be estimated. And, the normalized weight of a gene is modified, on the basis of the estimation results about the existence of direct regulations to it. These normalized and modified weights are used in queuing the possibility of the existence of a corresponding direct regulation. Computation results with the DREAM4 In Silico Size 100 Multifactorial subchallenge show that estimation performances of the suggested algorithm can even outperform the best team. Using the real data provided by the DREAM5 Network Inference Challenge, estimation performances can be ranked third. Furthermore, the high precision of the obtained most reliable predictions shows the suggested algorithm may be helpful in guiding biological experiment designs.     

Application of ANS fluorescent probes to identify hydrophobic sites on the surface of DREAM

DREAM (calsenilin or KChIP-3) is a calcium sensor involved in regulation of diverse physiological processes by interactions with multiple intracellular partners including DNA, K_v4 channels, and presenilin, however the detailed mechanism of the recognition of the intracellular partners remains unclear. To identify the surface hydrophobic surfaces on apo and Ca^2 +DREAM as a possible interaction sites for target proteins and/or specific regulators of DREAM function the binding interactions of 1,8-ANS and 2,6-ANS with DREAM were characterized by fluorescence and docking studies. Emission intensity of ANS–DREAM complexes increases upon Ca^2 + association which is consistent with an overall decrease in surface polarity. The dissociation constants for ANS binding to apoDREAM and Ca^2 +DREAM were determined to be 195 ± 20 μM and 62 ± 4 μM, respectively. Fluorescence lifetime measurements indicate that two ANS molecules bind in two independent binding sites on DREAM monomer. One site is near the exiting helix of EF-4 and the second site is located in the hydrophobic crevice between EF-3 and EF-4. 1,8-ANS displacement studies using arachidonic acid demonstrate that the hydrophobic crevice between EF-3 and EF-4 serves as a binding site for fatty acids that modulate functional properties of K_v4 channel:KChIP complexes. Thus, the C-terminal hydrophobic crevice may be involved in DREAM interactions with small hydrophobic ligands as well as other intracellular proteins.