On the properties and sequence context of structurally ambivalent fragments in proteins

The goal of this work is to characterize structurally ambivalent fragments in proteins. We have searched the Protein Data Bank and identified all structurally ambivalent peptides (SAPs) of length five or greater that exist in two different backbone conformations. The SAPs were classified in five distinct categories based on their structure. We propose a novel index that provides a quantitative measure of conformational variability of a sequence fragment. It measures the context-dependent width of the distribution of (phi,xi) dihedral angles associated with each amino acid type. This index was used to analyze the local structural propensity of both SAPs and the sequence fragments contiguous to them. We also analyzed type-specific amino acid composition, solvent accessibility, and overall structural properties of SAPs and their sequence context. We show that each type of SAP has an unusual, type-specific amino acid composition and, as a result, simultaneous intrinsic preferences for two distinct types of backbone conformation. All types of SAPs have lower sequence complexity than average. Fragments that adopt helical conformation in one protein and sheet conformation in another have the lowest sequence complexity and are sampled from a relatively limited repertoire of possible residue combinations. A statistically significant difference between two distinct conformations of the same SAP is observed not only in the overall structural properties of proteins harboring the SAP but also in the properties of its flanking regions and in the pattern of solvent accessibility. These results have implications for protein design and structure prediction.     

Investigations into sequence and conformational dependence of backbone entropy,inter-basin dynamics and the Flory isolated-pair hypothesis for peptides

The populations and transitions between Ramachandran basins are studied for combinations of the standard 20 amino acids in monomers, dimers and trimers using an implicit solvent Langevin dynamics algorithm and employing seven commonly used force-fields. Both the basin populations and inter-conversion rates are influenced by the nearest neighbor's conformation and identity, contrary to the Flory isolated-pair hypothesis. This conclusion is robust to the choice of force-field, even though the use of different force-fields produces large variations in the populations and inter-conversion rates between the dominant helical, extended beta, and polyproline II basins. The computed variation of conformational and dynamical properties with different force-fields exceeds the difference between explicit and implicit solvent calculations using the same force-field. For all force-fields, the inter-basin transitions exhibit a directional dependence, with most transitions going through extended beta conformation, even when it is the least populated basin. The implications of these results are discussed in the context of estimates for the backbone entropy of single residues, and for the ability of all-atom simulations to reproduce experimental protein folding data.     

Effect of backbone cyclization on protein folding stability: chain entropies of both the unfolded and the folded states are restricted

Circular versions of a large number of proteins have been designed by connecting the N and C termini via peptide linkers. A motivation for these designs is the assumed enhancement in folding stability, because backbone cyclization reduces the chain entropy of the unfolded state. Here, it is recognized that backbone cyclization also reduces the chain entropy of a flexible peptide linker in the folded state. Specifically, the end-to-end distance of the linker is restricted to fluctuations around the average displacement between the N and C termini of the folded protein. The balance of the chain-entropy reductions in the folded and unfolded states is used to predict the change in the unfolding free energy, deltadeltaG(cycl), by backbone cyclization. Predicted values of deltadeltaG(cycl) are in quantitative agreement with results of a careful study on cyclizing the 34 residue PIN1 WW domain by linkers with two to seen residues. The experimental results of an optimal linker length l=4 and a maximum stabilization of 1.7 kcal/mol are reproduced. Calculations of deltadeltaG(cycl) for a broad selection of circular proteins suggest that the stabilizing effect of backbone cyclization is modest, reflecting entropy reductions in both the unfolded and the folded states.     

An information-theoretical measure of taxonomic diversity

Traditional diversity indices are computed from the abundances of species present and are insensitive to taxonomic differences between species. However, a community in which most species belong to the same genus is intuitively less diverse than another community with a similar number of species distributed more evenly between genera. In this paper, we propose an information-theoretical measure of taxonomic diversity that reflects both the abundances and taxonomic distinctness of the species. Unlike previous measures of taxonomic diversity, such as Rao's quadratic entropy, in this new measure the analyzed taxonomic properties are associated with the single species instead of species pairs.     

Information gain for genetic parameter estimation with incorporation of marker data

Genetic marker data has been increasingly incorporated into segregation analysis, as combined segregation and linkage analysis has been performed more frequently. In this article, we study the extent of information gains with incorporation of marker data in segregation analysis, a topic that has not been investigated rigorously. Specifically, the current study is to investigate the influence of marker data on genetic model parameter estimation. A variance matrix criterion (as the inverse of the Fisher information matrix) and a relative entropy criterion (a measure of flatness of expected log-likelihood surface) are used to quantify the information gains. Our results indicate that substantial information gain can be achieved with the incorporation of marker data. The amount of variance reduction increases as the heterozygosity of the linked marker increases and as the trait gets closer to the linked marker(s). Incorporation of marker data in larger pedigrees also yields greater information gains based on both criteria. The effect of pedigree structure is also studied.     

Control of biogenic H(2)S production with nitrite and molybdate

The effects of the metabolic inhibitors, sodium nitrite and ammonium molybdate, on production of H₂S by a pure culture of the sulfate-reducing bacterium (SRB) Desulfovibrio sp. strain Lac6 and a consortium of SRB, enriched from produced water of a Canadian oil field, were investigated. Addition of 0.1 mM nitrite or 0.024 mM molybdate at the start of growth prevented the production of H₂S by strain Lac6. With exponentially growing cultures, higher levels of inhibitors, 0.25 mM nitrite or 0.095 mM molybdate, were required to suppress the production of H₂S. Simultaneous addition of nitrite and molybdate had a synergistic effect: at time 0, 0.05 mM nitrite and 0.01 mM molybdate, whereas during the exponential phase, 0.1 mM nitrite and 0.047 mM molybdate were sufficient to stop H₂S production. With an exponentially growing consortium of SRB, enriched from produced water of the Coleville oil field, much higher levels of inhibitors, 4 mM nitrite or 0.47 mM molybdate, were needed to stop the production of H₂S. The addition of these inhibitors had no effect on the composition of the microbial community, as shown by reverse sample genome probing. The results indicate that the efficiency of inhibitors in containment of SRB depends on the composition and metabolic state of the microbial community. Journal of Industrial Microbiology & Biotechnology (2001) 26, 350–355. Received 02 August 2000/ Accepted in revised form 17 April 2001     

Solute complexation degree with human serum albumin: biochromatographic approach

A mathematical model was developed for the study of the D,L-dansylamino acid retention mechanism in reversed-phase liquid chromatography using a C18 column as a stationary phase and human serum albumin (HSA) as an eluent modifier. The solute retention factor is dependent on the HSA concentration in the eluent as well as the binding constant of the guest-HSA complex. A determination of the degree of complexation n(c) (the percent of the complexed guest) could be carried out. Different Van 't Hoff plot shapes of the degree of complexation were observed with different eluent pH, confirming a change in the solute complexation mechanism for physiological pH (between 7-7.5). Enthalpy-entropy compensation was also analysed in relation to this mathematical model to confirm the solute complexation behavior with HSA. These results finally confirmed that at physiological pH and temperature (approximately 35 degrees C) values the HSA was in a favorable structural conformation for its binding with a great majority of drugs.     

Automated large-volume sample stacking procedure to detect labeled peptides at picomolar concentration using capillary electrophoresis and laser-induced fluorescence detection

We have developed an automated large-volume sample stacking (LVSS) procedure to detect fluorescein isothiocyanate-labeled peptides in the picomolar range. The injection duration is 10 min at 50 mbar to fill 62% of the capillary volume to the detection cell. The calculated limit of detection (S/N=3), filling 1% of the capillary volume, is 74 pM for bradykinin and 45 pM for L-enkephalin with samples diluted in water and analyzed in a 50 mM borate buffer, pH 9.2. With the automated LVSS system, the limits of detection are 7 pM for bradykinin, 3 pM for L-enkephalin and 2 pM for substance P. LVSS is shown to be quantitative from 500 to 10 pM.     

Optimizing resolution in multidimensional NMR by three-way decomposition

Resolution depends on the number of points sampled in a FID; in indirectly detected dimensions it is an important determinant of the total experiment time. Based on the high redundancy present in NMR data, we propose the following timesaving scheme for three-dimensional spectra. An extensive grid of discrete t1- and t2-values is used, which increases resolution while preserving the spectral width. Total experiment time is reduced by avoiding the recording of t3-FIDs for selected pairs of t1 and t2; typically the recording is omitted for about 75% of the (t1,t2) combinations. These data sets are referred to as sparse, and post-experimental processing making optimal use of spectral redundancy provides the missing, non-recorded data. We have previously shown that three-way decomposition (TWD) within the MUNIN approach provides a practical way to process dense NMR data sets. Here, a novel TWD algorithm [Ibraghimov, (2002) Numer. Linear Algebra Appl. 9, 551–565] is used to complement a sparselyrecorded time-domain data set by providing the missing FIDs for all (t1,t2) combinations omitted in the experiment. A necessary condition is that for each t1-value at least a few FIDs are recorded, and similar for each t2-value. The method is demonstrated on non-uniformly sampled ¹⁵N-NOESY-HSQC data sets recorded for the 14 kD protein azurin. The spectra obtained by TWD, reconstruction and ordinary transform to frequency-domain are, in spite of the large number of signals and the high dynamic range typical for NOESYs, highly similar to a corresponding reference spectrum, for which all (t1,t2) combinations were recorded.     

A new approach to the design of uniquely folded thermally stable proteins

A new computer program (CORE) is described that predicts core hydrophobic sequences of predetermined target protein structures. A novel scoring function is employed, which for the first time incorporates parameters directly correlated to free energies of unfolding (deltaGu), melting temperatures (Tm), and cooperativity. Metropolis-driven simulated annealing and low-temperature Monte Carlo sampling are used to optimize this score, generating sequences predicted to yield uniquely folded, stable proteins with cooperative unfolding transitions. The hydrophobic core residues of four natural proteins were predicted using CORE with the backbone structure and solvent exposed residues as input. In the two smaller proteins tested (Gbeta1, 11 core amino acids; 434 cro, 10 core amino acids), the native sequence was regenerated as well as the sequence of known thermally stable variants that exhibit cooperative denaturation transitions. Previously designed sequences of variants with lower thermal stability and weaker cooperativity were not predicted. In the two larger proteins tested (myoglobin, 32 core amino acids; methionine aminopeptidase, 63 core amino acids), sequences with corresponding side-chain conformations remarkably similar to that of native were predicted.