Structure-based design of model proteins

A structure-based, sequence-design procedure is proposed in which one considers a set of decoy structures that compete significantly with the target structure in being low energy conformations. The decoy structures are chosen to have strong overlaps in contacts with the putative native state. The procedure allows the design of sequences with large and small stability gaps in a random-bond heteropolymer model in both two and three dimensions by an appropriate assignment of the contact energies to both the native and nonnative contacts. The design procedure is also successfully applied to the two-dimensional HP model. Proteins 31:10–20, 1998. © 1998 Wiley-Liss, Inc.     

Identifying Chromosomal Fragile Sites from a Hierarchical-Clustering Point of View

Identification of fragile sites is a way to investigate the genetic abnormalities that are the hallmark of cancer and play an important role in carcinogenesis. Manifestation of nonrandom breakage at a chromosome band is an essential criterion for determination of the fragility of the band. In this article, a new detection procedure is proposed. This new procedure takes the relationship of one site with the others into consideration and can be applied to tests of the randomness of breakpoints under either the proportional probability model (PPM) or the equiprobability model (EPM). The procedure can form a grouping structure that classifies all sites into several clusters. It is applied to identification of fragile sites for a real data set for Chinese patients with colorectal carcinoma for illustration of the proposed method.     

Some Probabilistic Results on the Nonrandomness of Simple Sequence Repeats in DNA Sequences

Some probabilistic results on simple sequence repeats (SSRs) in DNA sequences are derived and used to quantify the nonrandomness of SSRs as an index of nonrandomness. The applicability of the index of nonrandomness is illustrated using several examples from the literature on selected human diseased genes.     

Structural model for family 32 of glycosyl-hydrolase enzymes

A structural model is presented for family 32 of the glycosyl-hydrolase enzymes based on the beta-propeller fold. The model is derived from the common prediction of two different threading methods, TOPITS and THREADER. In addition, we used a correlated mutation analysis and prediction of active-site residues to corroborate the proposed model. Physical techniques (circular dichroism and differential scanning calorimetry) confirmed two aspects of the prediction, the proposed all-beta fold and the multi-domain structure. The most reliable three-dimensional model was obtained using the structure of neuraminidase (1nscA) as template. The analysis of the position of the active site residues in this model is compatible with the catalytic mechanism proposed by Reddy and Maley (J. Biol. Chem. 271:13953–13958, 1996), which includes three conserved residues, Asp, Glu, and Cys. Based on this analysis, we propose the participation of one more conserved residue (Asp 162) in the catalytic mechanism. The model will facilitate further studies of the physical and biochemical characteristics of family 32 of the glycosyl-hydrolases. Proteins 33:383–395, 1998. © 1998 Wiley-Liss, Inc.     

A new model for how O6-methylguanine-DNA methyltransferase binds DNA

Human methyltransferase (hAT) catalyzes the transfer of an alkyl group from the 6-position of guanine to an active site Cys residue. The physiological role of hAT is the repair of alkylated guanine residues in DNA. However, the repair of methylated or chloroethylated guanine bases negates the effects of certain chemotherapeutic agents. A model of how hAT binds DNA might be useful in the design of compounds that could inactivate hAT. We have used computer modeling studies to generate such a model. The model utilizes a helix-loop-wing DNA binding motif found in Mu transposase. The model incorporates a flipped out guanine base in order to bring the methylated oxygen atom close to the active site Cys residue. The model is consistent with a variety of chemical and biochemical data. Proteins 32:3–6, 1998. © 1998 Wiley-Liss, Inc.     

Detecting evolutionary trends from molecular data

The measures of compositional nonrandomness to be discussed as to their physical significance and to their power of detecting evolutionary significant variations are (see article)(pi a priori probability for amino acid i, ni its number of occurrences in a protein of length L). As a concrete example, the pi are here supposed to represent equal frequencies of all non-stop codons. For each quantity, four levels are defined: The base level, with optimal (i.e. minimal nonrandomness) composition, admitting non-integer values of ni; the integer level with optimal integer composition; the noise level, represented by a typical random cain; and the real protein level. On all these levels, S, which is the measure with the most direct physical sense, shows the smoothest behavior with the smallest relative fluctuations and thus the highest resolution.     

Models of multilocus recombination: nonrandomness in chiasma number and crossover positions.

Cytological evidence indicates that genetic interference can be partitioned into two empirical components: nonrandomness in the number of chiasmata that occur and nonrandomness in the locations of multiple chiasmata. Previous studies have incorporated the first effect into genetic models for analyzing multipoint data. An extension to this approach is presented which allows for the second component of interference by modeling the probability density function of the locations of multiple crossovers. Results of reanalyses of multilocus data for the Drosophila X chromosome show that models that incorporate only the first effect give a better fit to these data than do standard mapping functions and that the extended model significantly improved the fit by decreasing the predicted frequency of multiple crossovers in nearby regions. Our results demonstrate that chiasma-based models of multilocus recombination, which are unique in incorporating direct estimates of the frequency of multiple crossovers for a chromosome region, can provide a powerful and realistic means of accounting for genetic interference when applied to the problems of gene localization, locus ordering, and exclusion mapping.     

Point perturbation analysis of experimental data

A new method of data analysis is proposed. The method is based on discrete perturbation of experimental data points, which is used to probe the metric of the parameter hyperspace. Perturbation-induced fluctuations in the residual values are analysed by discrete Fourier transform to yield the autocorrelation function and a relaxation length for each experimental point. This parameter provides a quantitative measure of correlation and hence nonrandomness of residuals. The method is applied to the analysis of measurements of the shear viscosity of a 2,6-lutidine/water mixture near the critical point, and to the oxygen and carbon monoxide binding reactions to human hemoglobin. Relaxation profiles are constructed for several experimental data sets. Departure from random behavior in the residuals is discussed in connection with the theoretical interpretations of the phenomenon under consideration.     

Influence of Exon Duplication on Intron and Exon Phase Distribution

Nonrandomness in the intron and exon phase distributions in a sample of 305 human genes has been found and analyzed. It was shown that exon duplications had a significant effect on the exon phase nonrandomness. All of the nonrandomness is probably due to both the processes of exon duplication and shuffling. A quantitative estimation of exon duplications in the human genome and their influence on the intron and exon phase distributions has been analyzed. According to our estimation, the proportion of duplicated exons in the human genome constitutes at least 6% of the total. Generalizing the particular case of exon duplication to the more common event of exon shuffling, we modeled and analyzed the influence of exon shuffling on intron phase distribution. Received: 28 March 1997 / Accepted: 9 July 1997     

pH-induced co-operative effects in hysteretic enzymes. 1. A theoretical model of a new type of co-operative behaviour controlled by pH

A new model which provides an explanation for pH-induced co-operativity of hysteretic enzymes is proposed. The essence of the model is that a region, or a domain, of the enzyme undergoes a spontaneous 'slow' conformational change which does not affect the geometry of the active site. The region which undergoes this spontaneous conformational transition bears an ionizable group. Kinetic co-operativity occurs if the pK of this ionizable group changes upon this conformational transition. Thus co-operativity does not arise from a distortion of the active site. An interesting prediction of the model is that at 'extreme' pH values co-operativity must be suppressed. Although the kinetic equation pertaining to the model is of the 2:2 type, co-operativity is not maximum or minimum at half-saturation of the enzyme by the substrate, as occurs with 2:2 binding isotherms. A new index of maximum or minimum kinetic co-operativity, whether this extreme occurs at half-saturation or not, has been proposed which allows the change of kinetic co-operativity to be followed as a function of pH. It is believed that this model will be useful in explaining the behaviour of enzymes attached to biological polyelectrolytes, such as membranes or cell envelopes.     

Accessibility to internal cavities and ligand binding sites monitored by protein crystallographic thermal factors

Protein structures are flexible both in solution and in the solid state. X-ray crystallographically determined thermal factors monitor the flexibility of protein atoms. A method utilizing such factors is proposed to delineate protein regions through which a ligand can exchange between binding site and bulk solvent. It is based on the assumption that thermally excited protein regions are excellent candidates for opening a ligand channel. Computationally simple and inexpensive, the method analyzes directions from which thermal factors can propagate within the protein, resulting in thermal motion paths (TMPs). Applications to engineered T4 lysozymes, where an artificial internal cavity can host hydrophobic molecules, and to sperm whale myoglobins, where the active site is completely buried, yielded results in agreement with other independent structural observations and with previous hypotheses. Further new features could also be suggested. The proposed TMP analysis could aid molecular dynamics simulation studies as well as time-resolved and site-directed mutagenesis experimental studies, especially given its modest computational expense and its direct roots in experimental results based on thermal factors determined in high-resolution crystallographic studies. Proteins 31:201–213,1998. © 1998 Wiley-Liss, Inc.     

Adaptive behavior with protozoa-inspired dynamics

A protozoa-inspired system for controlling autonomous robots is explained. The proposed system uses the nonlinear dynamics of the molecule concentrations in the cytoplasm of a protozoan. The model includes many parameters which are optimized using the genetic algorithm. A mobile robot obstacle avoidance example is shown with simulation results. The proposed protozoa-inspired system opens interesting possibilities to create new species by combining individual models by symbiotic or colony relationships. Received: 13 March 1996 / Accepted in revised form: 17 July 1998     

A hierarchical method for generating low-energy conformers of a protein-ligand complex

A novel dynamical protocol for finding the low-energy conformations of a protein-ligand complex is described. The energy functions examined consist of an empirical force field with four different dielectric screening models; the generalized Born/surface area model also is examined. Application of the method to three complexes of known crystal structure provides insights into the energy functions used for selecting low-energy docked conformations and into the structure of the binding-energy surface. Evidence is presented that the local energy minima of a ligand in a binding site are arranged in a hierarchical fashion. This observation motivates the construction of a hierarchical docking algorithm that substantially enriches the population of ligand conformations close to the crystal conformation. The algorithm is also adapted to permit docking into a flexible binding site and preliminary tests of this method are presented. Proteins 33:475–495, 1998. © 1998 Wiley-Liss, Inc.     

Nonrandom location of IS1 elements in the genomes of natural isolates of Escherichia coli

We have studied the spatial distribution of IS1 elements in the genomes of natural isolates comprising the ECOR reference collection of Escherichia coli. We find evidence for nonrandomness at three levels. Many pairs of IS1 elements are in much closer proximity (< 10 kb) than can be accounted for by chance. IS1 elements in close proximity were identified by long-range PCR amplification of the genomic sequence between them. Each amplified region was sequenced and its map location determined by database screening of DNA hybridization. Among the ECOR strains with at least two IS1 elements, 54% had one or more pairs of elements separated by < 10 kb. We propose that this type of clustering is a result of "local hopping," in which we assume that a significant proportion of tranposition events leads to the insertion of a daughter IS element in the vicinity of the parental element. A second level of nonrandomness is found in strains with a modest number of IS1 elements that are mapped through the use of inverse PCR to amplify flanking genomic sequences: in these strains, the insertion sites tend to be clustered over a smaller region of chromosome than would be expected by chance. A third level of nonrandomness is observed in the composite distribution of IS elements across strains: among 20 mapped IS1 elements, none were found in the region of 48-77 minutes, a significant gap. One region of the E. coli chromosome, at 98 min, had a cluster of IS1 elements in seven ECOR strains of diverse phylogenetic origin. We deduce from sequence analysis that this pattern of distribution is a result of initial insertion in the most recent common ancestor of these strains and therefore not a hot spot of insertion. Analysis using long- range PCR with primers for IS2 and IS3 also yielded pairs of elements in close proximity, suggesting that these elements may also occasionally transpose by local hopping.      

A liquid-medium-based protocol for rapid regeneration from embryogenic soybean cultures

Soybean [Glycine max (L.) Merrill] somatic embryos of the cultivar Jack underwent histodifferentiation in liquid Murashige and Skoog (MS) medium with 3% maltose, or according to the standard published procedure employing solidified MS media, permitting the recovery of an average of 8.1 and 3.9 embryos/mg of embryogenic tissue, respectively. Cotyledon-stage embryos that developed in liquid medium were ready for desiccation within 4 weeks, while the embryos from the standard procedure required a maturation step for an additional 4 weeks. Comparison of embryo development in MS medium with maltose or FN Lite-based medium without growth regulators and supplemented with maltose or an equimolar amount of sucrose revealed that sucrose promotes faster embryo histodifferentiation and maturation, and allows the recovery of up to 50% more mature, cotyledon-stage embryos within 3 weeks. The use of this liquid-medium-based protocol relative to the standard procedure led to a fourfold increase in the number of cotyledon-stage embryos recovered from other genotypes tested. In many cases, however, the percent germination was lower. Application of this new procedure also made it possible to harvest transgenic seed 9 months following biolistic bombardment, as compared to the 13 months required when the standard solid-medium-based protocol was used. Received: 1 December 1997 / Revision received: 27 April 1998 / Accepted: 20 May 1998     

The second ADF/cofilin actin-binding site exists in F-actin, the cofilin-G-actin complex, but not in G-actin.

ADF/cofilins are actin binding proteins that bind actin close to both the N- and C-termini (site 1), and we have found a second cofilin binding site (site 2) centered around helix 112-125 [Renoult, C., Ternent, D., Maciver, S.K., Fattoum, A., Astier, C., Benyamin, Y. & Roustan, C. (1999) J. Biol. Chem. 274, 28893-28899]. We proposed a model in which ADF/cofilin intercalated between subdomains 1 and 2 of two longitudinally associated actin monomers within the actin:cofilin cofilament, explaining the change in twist that ADF/cofilins induce in the filament [McGough, A. Pope, B., Chiu, W. & Weeds, A. (1998) J. Cell Biol. 138, 771-781]. Here, we have determined the fuller extent of the cofilin footprint on site 1 of actin. Site 1 is primarily the G-actin binding site. Experiments with both peptide mimetics and fluorescently labeled cofilin suggest that site 2 only becomes available for cofilin binding within the filament, possibly due to motion between subdomains 1 and 2 within an actin monomer. We have detected motion between subdomains 1 and 2 of G-actin by FRET induced by cofilin, to reveal the second cofilin-binding site. This motion may also explain how cofilins inhibit the nucleotide exchange of actin, and why the actin:cofilin complex is polymerizable without dissociation.     

Determination of drug levels in human serum by circular dichroism

A study of the applicability of circular dichroism (CD) for the determination of drug levels in human serum is described and a new method for the quantitative determination of optically active absorbing drugs having Cotton effects at wavelengths above 250 nm in human serum and/or plasma is proposed. The principal advantages of this method are speed, economy, and simplicity, no derivatization or chromatographic separation steps being needed. The validity of the CD determination was confirmed by analysis of variance, β-lactam antibiotics being chosen as model drugs. In addition, the validation studies performed confirm the accuracy and precision of the proposed method. For β-lactam antibiotics lacking Cotton effects above 250 nm, an alternative method based on the extraction of the drug from serum is considered. Chirality 10:507–512, 1998. © 1998 Wiley-Liss, Inc.     

Heterochromatin protein 1 binds transgene arrays

Heterochromatin protein 1 (HP1) of Drosophila and its homologs in vertebrates are key components of constitutive heterochromatin. Here we provide cytological evidence for the presence of heterochromatin within a euchromatic chromosome arm by immunolocalization of HP1 to the site of a silenced transgene repeat array. The amount of HP1 associated with arrays in polytene chromosomes is correlated with the array size. Inverted transposons within an array or increased proximity of an array to blocks of naturally occurring heterochromatin may increase transgene silencing without increasing HP1 labeling. Less dense anti-HP1 labeling is found at transposon arrays in which there is no transgene silencing. The results indicate that HP1 targets the chromatin of transposon insertions and binds more densely at a site with repeated sequences susceptible to heterochromatin formation. Received: 26 June 1998; in revised form: 6 July 1998 / Accepted: 12 July 1998     

On the statistical analysis of multiple-choice feeding preference experiments

A stopping rule for an experiment defines when (under what conditions) the experiment is terminated. I investigated the stopping rules used in numerous multiple–choice feeding-preference experiments and also examined a recently proposed method for analyzing the data arising from such experiments. All of the surveyed experiments imposed stopping rules which result in a random total food consumption. If an acceptable quantification of preference is relative consumption of different food types, then the proposed analysis will likely misstate the information about preference conveyed by the data. This is due to the fact that the method may confound differences in preferences among food types with differences in the total consumption across trials. I discuss this issue in detail and present an alternative procedure which is appropriate under all stopping regimes when preference is quantified through relative consumption. The procedure I suggest uses an index which is a multivariate generalization of the preference index suggested by Kogan and Goeden (Ann Entomol Soc 1970; 63: 1175–1180) and Kogan (Ann Entomol Soc 1972; 65: 675–683) and which is analogous to a selection index for discrete food units proposed by Manly (Biometrics 1974; 30: 281–294). Received: 29 November 1997 / Accepted: 20 April 1998