The BLAST programs are widely used tools for searching protein and DNA databases for sequence similarities. For protein comparisons, a variety of definitional, algorithmic and statistical refinements described here permits the execution time of the BLAST programs to be decreased substantially while enhancing their sensitivity to weak similarities. A new criterion for triggering the extension of word hits, combined with a new heuristic for generating gapped alignments, yields a gapped BLAST program that runs at approximately three times the speed of the original. In addition, a method is introduced for automatically combining statistically significant alignments produced by BLAST into a position-specific score matrix, and searching the database using this matrix. The resulting Position-Specific Iterated BLAST (PSIBLAST) program runs at approximately the same speed per iteration as gapped BLAST, but in many cases is much more sensitive to weak but biologically relevant sequence similarities. PSI-BLAST is used to uncover several new and interesting members of the BRCT superfamily.

/pdf/gapped-blast-and-psi-blast-a-new-generation-of-protein-2ybl67x5yx.pdf

Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

Technologies to measure whole-genome mRNA abundances1,2,3 and methods to organize and display such data4,5,6,7,8,9,10 are emerging as valuable tools for systems-level exploration of transcriptional regulatory networks. For instance, it has been shown that mRNA data from 118 genes, measured at several time points in the developing hindbrain of mice, can be hierarchically clustered into various patterns (or 'waves') whose members tend to participate in common processes5. We have previously shown that hierarchical clustering can group together genes whose cis-regulatory elements are bound by the same proteins in vivo6. Hierarchical clustering has also been used to organize genes into hierarchical dendograms on the basis of their expression across multiple growth conditions7. The application of Fourier analysis to synchronized yeast mRNA expression data has identified cell-cycle periodic genes, many of which have expected cis-regulatory elements8. Here we apply a systematic set of statistical algorithms, based on whole-genome mRNA data, partitional clustering and motif discovery, to identify transcriptional regulatory sub-networks in yeast—without any a priori knowledge of their structure or any assumptions about their dynamics. This approach uncovered new regulons (sets of co-regulated genes) and their putative cis-regulatory elements. We used statistical characterization of known regulons and motifs to derive criteria by which we infer the biological significance of newly discovered regulons and motifs. Our approach holds promise for the rapid elucidation of genetic network architecture in sequenced organisms in which little biology is known.

/pdf/systematic-determination-of-genetic-network-architecture-e5mdjear7j.pdf

Systematic determination of genetic network architecture

A wealth of protein and DNA sequence data is being generated by genome projects and other sequencing efforts. A crucial barrier to deciphering these sequences and understanding the relations among them is the difficulty of detecting subtle local residue patterns common to multiple sequences. Such patterns frequently reflect similar molecular structures and biological properties. A mathematical definition of this "local multiple alignment" problem suitable for full computer automation has been used to develop a new and sensitive algorithm, based on the statistical method of iterative sampling. This algorithm finds an optimized local alignment model for N sequences in N-linear time, requiring only seconds on current workstations, and allows the simultaneous detection and optimization of multiple patterns and pattern repeats. The method is illustrated as applied to helix-turn-helix proteins, lipocalins, and prenyltransferases.

/pdf/detecting-subtle-sequence-signals-a-gibbs-sampling-strategy-56bu9neqln.pdf

Detecting Subtle Sequence Signals: A Gibbs Sampling Strategy for Multiple Alignment

Molecular recognition is central to all biological processes. For the past 50 years, Koshland's 'induced fit' hypothesis has been the textbook explanation for molecular recognition events. However, recent experimental evidence supports an alternative mechanism. 'Conformational selection' postulates that all protein conformations pre-exist, and the ligand selects the most favored conformation. Following binding the ensemble undergoes a population shift, redistributing the conformational states. Both conformational selection and induced fit appear to play roles. Following binding by a primary conformational selection event, optimization of side chain and backbone interactions is likely to proceed by an induced fit mechanism. Conformational selection has been observed for protein-ligand, protein-protein, protein-DNA, protein-RNA and RNA-ligand interactions. These data support a new molecular recognition paradigm for processes as diverse as signaling, catalysis, gene regulation and protein aggregation in disease, which has the potential to significantly impact our views and strategies in drug design, biomolecular engineering and molecular evolution.

/pdf/the-role-of-dynamic-conformational-ensembles-in-biomolecular-33rwd16x5b.pdf

The role of dynamic conformational ensembles in biomolecular recognition

The clustered regularly interspaced short palindromic repeats (CRISPR)-associated enzyme Cas9 is an RNA-guided endonuclease that uses RNA-DNA base-pairing to target foreign DNA in bacteria. Cas9-guide RNA complexes are also effective genome engineering agents in animals and plants. Here we use single-molecule and bulk biochemical experiments to determine how Cas9-RNA interrogates DNA to find specific cleavage sites. We show that both binding and cleavage of DNA by Cas9-RNA require recognition of a short trinucleotide protospacer adjacent motif (PAM). Non-target DNA binding affinity scales with PAM density, and sequences fully complementary to the guide RNA but lacking a nearby PAM are ignored by Cas9-RNA. Competition assays provide evidence that DNA strand separation and RNA-DNA heteroduplex formation initiate at the PAM and proceed directionally towards the distal end of the target sequence. Furthermore, PAM interactions trigger Cas9 catalytic activity. These results reveal how Cas9 uses PAM recognition to quickly identify potential target sites while scanning large DNA molecules, and to regulate scission of double-stranded DNA.

/pdf/dna-interrogation-by-the-crispr-rna-guided-endonuclease-cas9-37k0aiok4x.pdf

DNA interrogation by the CRISPR RNA-guided endonuclease Cas9

http://galton.uchicago.edu/~lalley/Courses/312/Berg1989.pdf

Facilitated target location in biological systems

Selection of DNA Binding Sites by Regulatory Proteins Statistical-mechanical Theory and Application to Operators and Promoters

3. WHAT IS DIFFUSION CONTROL?" 134 3.1 The Spherical Case " 135 3.2 Other Geometries 141

Diffusion-Controlled Macromolecular Interactions

In this paper we summarize the various factors that must be considered in establishing the operational specificity of the binding of a protein regulator of gene expression to a DNA target site. We consider informational (combinatorial) aspects of binding-site specification, actual recognition mechanisms, and the thermodynamics of target-site selection against a background of competing pseudospecific and non-(sequence)-specific DNA binding sites. The results provide insight into the design, specification, and possibly the evolution of regulatory proteins and their chromosomal binding targets, as well as into practical aspects of the design of regulatory-protein isolation schemes and physicochemical regulatory considerations in vivo.

https://www.pnas.org/content/pnas/83/6/1608.full.pdf

Otto G. Berg

Papers

Facilitated target location in biological systems

Selection of DNA Binding Sites by Regulatory Proteins Statistical-mechanical Theory and Application to Operators and Promoters

Diffusion-Controlled Macromolecular Interactions

On the specificity of DNA-protein interactions.

Selection of DNA binding sites by regulatory proteins. II. The binding specificity of cyclic AMP receptor protein to recognition sites.