A novel neural network technique for analysis and classification of EM single-particle images.

We propose a novel self-organizing neural network for the unsupervised classification of electron microscopy (EM) images of biological macromolecules. The radical novelty of the algorithm lies in its rigorous mathematical formulation that, starting from a large set of possibly very noisy input data, finds a set of "representative" data items, organized onto an ordered output map, such that the probability density of this set of representative items resembles at its possible best the probability density of the input data. In a way, it summarizes large amounts of information into a concise description that rigorously keeps the basic pattern of the input data distribution. In this application to the field of three-dimensional EM of single particles, two different data sets have been used; one comprised 2458 rotational power spectra of individual negative stain images of the G40P helicase of Bacillus subtilis bacteriophage SPP1, and the other contained 2822 cryoelectron images of SV40 large T-antigen. Our experimental results prove that this technique is indeed very successful, providing the user with the capability of exploring complex patterns in a succinct, informative, and objective manner. The above facts, together with the consideration that the integration of this new algorithm with commonly used software packages is immediate, prompt us to propose it as a valuable new tool in the analysis of large collections of noisy data.

The DnaB.DnaC complex: a structure based on dimers assembled around an occluded channel.

Replicative helicases are motor proteins that unwind DNA at replication forks. Escherichia coli DnaB is the best characterized member of this family of enzymes. We present the 26 A resolution three-dimensional structure of the DnaB hexamer in complex with its loading partner, DnaC, obtained from cryo-electron microscopy. Analysis of the volume brings insight into the elaborate way the two proteins interact, and provides a structural basis for control of the symmetry state and inactivation of the helicase by DnaC. The complex is arranged on the basis of interactions among DnaC and DnaB dimers. DnaC monomers are observed for the first time to arrange as three dumb-bell-shaped dimers that interlock into one of the faces of the helicase. This could be responsible for the freezing of DnaB in a C(3) architecture by its loading partner. The central channel of the helicase is almost occluded near the end opposite to DnaC, such that even single-stranded DNA could not pass through. We propose that the DnaB N-terminal domain is located at this face.

pH-controlled quaternary states of hexameric DnaB helicase.

DnaB is the major helicase in the Escherichia coli replisome. It is a homohexameric enzyme that interacts with many other replisomal proteins and cofactors. It is usually loaded onto a single strand of DNA at origins of replication from its complex with its loading partner DnaC, then translocates in the 5' to 3' direction, unwinding duplex DNA in an NTP-driven process. Quaternary polymorphism has been described for the DnaB oligomer, a feature it has in common with some other hexameric helicases. In the present work, electron microscopy and in- depth rotational analysis studies of negatively stained specimens has allowed the establishment of conditions that govern the transition between the two different rotational symmetry states (C(3) and C(6)) of DnaB. It is shown: (a) that the pH value of the sample buffer, within the physiological range, dictates the quaternary organisation of the DnaB oligomer; (b) that the pH-induced transition is fully reversible; (c) that the type of adenine nucleotide complexed to DnaB, whether hydrolysable or not, does not affect its quaternary architecture; (d) that the DnaB.DnaC complex exists only as particles with C(3) symmetry; and (e) that DnaC interacts only with DnaB particles that have C(3) symmetry. Structural consequences of this quaternary polymorphism, as well as its functional implications for helicase activity, are discussed.

Large T-antigen double hexamers imaged at the simian virus 40 origin of replication.

The initial step of simian virus 40 (SV40) DNA replication is the binding of the large tumor antigen (T-Ag) to the SV40 core origin. In the presence of Mg(2+) and ATP, T-Ag forms a double-hexamer complex covering the complete core origin. By using electron microscopy and negative staining, we visualized for the first time T-Ag double hexamers bound to the SV40 origin. Image processing of side views of these nucleoprotein complexes revealed bilobed particles 24 nm long and 8 to 12 nm wide, which indicates that the two T-Ag hexamers are oriented head to head. Taking into account all of the biochemical data known on the T-Ag-DNA interactions at the replication origin, we present a model in which the DNA passes through the inner channel of both hexamers. In addition, we describe a previously undetected structural domain of the T-Ag hexamer and thereby amend the previously published dimensions of the T-Ag hexamer. This domain we have determined to be the DNA-binding domain of T-Ag.

Predicting the conformational class of short and medium size loops connecting regular secondary structures: application to comparative modelling.

Loops are regions of non-repetitive conformation connecting regular secondary structures. They are both the most difficult and error prone regions of a protein to solve by X-ray crystallography and the hardest regions to model using comparative procedures. Although a loop can sometimes be modelled from a homologue, very often it must be selected from outside the family. The loop prediction procedure, SLoop, attempts to identify the conformational class of the loop rather than to select a specific loop from a set of fragments extracted from known structures or generated ab initio. Templates are constructed for each of the 161 loop conformational classes that have been identified from the clustering of the structures of some 2024 loops of one to eight residues in length. A class template describes both sequence preferences and relative disposition of bounding secondary structures. During comparative modelling, the conformation of a loop can be predicted by identifying a loop class with which its sequence and disposition of bounding secondary structures are compatible. The procedure is tested on an unrelated non-redundant set of 1785 loops under stringent and lax evaluation schemes. Optimal sequence score cut-offs are identified such that the prediction rate is equal to the percentage of loops assigned to acceptable classes. Under the stringent evaluation, at the optimal sequence score cut-off, a conformation is predicted for 50% of loops of which 47% are correct, while under the lax evaluation a conformation is predicted for 63% of loops of which 54% are correct. Sequence score is shown to be a good indicator of the probability of a prediction being correct. Loop length also has a strong affect on prediction outcomes. Considering only loops of two to five residues in length, under the stringent evaluation 62% of loops are predicted with 52% of these predictions being correct while under the lax evaluation predictions are provided for 75% of loops of which 57% are correct.

Conformational analysis and clustering of short and medium size loops connecting regular secondary structures: a database for modeling and prediction.

Loops are regions of nonrepetitive conformation connecting regular secondary structures. We identified 2,024 loops of one to eight residues in length, with acceptable main-chain bond lengths and peptide bond angles, from a database of 223 protein and protein-domain structures. Each loop is characterized by its sequence, main-chain conformation, and relative disposition of its bounding secondary structures as described by the separation between the tips of their axes and the angle between them. Loops, grouped according to their length and type of their bounding secondary structures, were superposed and clustered into 161 conformational classes, corresponding to 63% of all loops. Of these, 109 (51% of the loops) were populated by at least four nonhomologous loops or four loops sharing a low sequence identity. Another 52 classes, including 12% of the loops, were populated by at least three loops of low sequence similarity from three or fewer nonhomologous groups. Loop class suprafamilies resulting from variations in the termini of secondary structures are discussed in this article. Most previously described loop conformations were found among the classes. New classes included a 2:4 type IV hairpin, a helix-capping loop, and a loop that mediates dinucleotide-binding. The relative disposition of bounding secondary structures varies among loop classes, with some classes such as beta-hairpins being very restrictive. For each class, sequence preferences as key residues were identified; those most frequently at these conserved positions than in proteins were Gly, Asp, Pro, Phe, and Cys. Most of these residues are involved in stabilizing loop conformation, often through a positive phi conformation or secondary structure capping. Identification of helix-capping residues and beta-breakers among the highly conserved positions supported our decision to group loops according to their bounding secondary structures. Several of the identified loop classes were associated with specific functions, and all of the member loops had the same function; key residues were conserved for this purpose, as is the case for the parvalbumin-like calcium-binding loops. A significant number, but not all, of the member loops of other loop classes had the same function, as is the case for the helix-turn-helix DNA-binding loops. This article provides a systematic and coherent conformational classification of loops, covering a broad range of lengths and all four combinations of bounding secondary structure types, and supplies a useful basis for modelling of loop conformations where the bounding secondary structures are known or reliably predicted.

Analysis, clustering and prediction of the conformation of short and medium size loops connecting regular secondary structures.

Loops are regions of non-repetitive conformation connecting regular secondary structures. They are both the most difficult and error prone regions of a protein to solve by X-ray crystallography and the hardest regions to model using knowledge-based procedures. While the core of a protein can be straight forwardly modelled from the structurally conserved regions of homologues of known structure, loops must be modelled from a selected homologue or from a loop chosen from outside the family. Here we present a loop prediction procedure that attempts to identify the conformational class of the loop rather than to select a specific loop from a database of fragments. The structures of some 2083 loops of one to eight residues in length were extracted from a database of 225 protein and protein domain structures. For each loop, the relative disposition of its bounding secondary structures is described by the separation between the tips of their axes, the angle and dihedral angle between their axes. From the clustering of the loops according to the root mean square deviation of their spatial fit, a total of 162 loop conformational classes, including 79% of loops, were identified. One-hundred and eight of these, involving 66% of the loops, were populated by at least four non-homologous loops or four loops sharing a low sequence identity. Another 54 classes, including 13% of the loops, were populated by at least three loops of low sequence similarity from three or fewer non-homologous groups. Most of the previously described loop conformations were found among the populated classes. For each class a template was constructed containing both sequence preferences and the relative disposition of bounding secondary structures among member loops. During comparative modelling, the conformation of a loop can be predicted by identifying a loop class with which its sequence and disposition of bounding secondary structures are compatible.

Molecular evolution and domain structure of plasminogen-related growth factors (HGF/SF and HGF1/MSP).

Plasminogen-related growth factors, a new family of polypeptide growth factors with the basic domain organization and mechanism of activation of the blood proteinase plasminogen, include hepatocyte growth factor/scatter factor (HGF/SF), a potent effector of the growth, movement, and differentiation of epithelia and endothelia, and hepatocyte growth factor-like/macrophage stimulating protein (HGF1/MSP), an effector of macrophage chemotaxis and phagocytosis. Phylogeny of the serine proteinase domains and analysis of intron-exon boundaries and kringle sequences indicate that HGF/SF, HGF1/MSP, plasminogen, and apolipoprotein (a) have evolved from a common ancestral gene that consisted of an N-terminal domain corresponding to plasminogen activation peptide (PAP), 3 copies of the kringle domain, and a serine proteinase domain. Models of the N domains of HGF/SF, HGF1/MSP, and plasminogen, characterized by the presence of 4 conserved Cys residues forming a loop in a loop, have been modeled based on disulfide-bond constraints. There is a distinct pattern of charged and hydrophobic residues in the helix-strand-helix motif proposed for the PAP domain of HGF/SF; these may be important for receptor interaction. Three-dimensional structures of the 4 kringle and the serine proteinase domains of HGF/SF were constructed by comparative modeling using the suite of programs COMPOSER and were energy minimized. Docking of a lysine analogue indicates a putative lysine-binding pocket within kringle 2 (and possibly another in kringle 4). The models suggest a mechanism for the formation of a noncovalent HGF/SF homodimer that may be responsible for the activation of the Met receptor. These data provide evidence for the divergent evolution and structural similarity of plasminogen, HGF/SF, and HGF1/MSP, and highlight a new strategy for growth factor evolution, namely the adaptation of a proteolytic enzyme to a role in receptor activation.

RNA-mediated specificity of DNA packaging into hybrid lambda/phi 29 proheads.

A small RNA (pRNA, 174 nt) is known to be essential for DNA packaging in bacteriophage phi 29. However, in an in vitro DNA packaging system based on hybrid lambda/phi 29 proheads (made up of head proteins from phage lambda and connectors from phage phi 29), the specificity of DNA packaging is lost, and different RNA molecules fulfil the requirements for DNA packaging, albeit with less efficiency than phi 29 pRNA. Competition assays with RNAs from different sources have shown that phi 29 connectors bind preferentially pRNA. An increase in the efficiency of phi 29 DNA packaging into hybrid proheads induced by phi 29 pRNA is observed because, when phi 29 pRNA is incubated with hybrid proheads, phi 29 DNA is packaged more efficiently than other DNAs of similar length. Furthermore, when hybrid proheads carrying phi 29 pRNA are incubated with a mixture of DNAs from different sources, phi 29 DNA is selectively packaged, thus indicating that phi 29 pRNA determines the specificity of DNA packaging.

Characterization of an RNA-binding domain in the bacteriophage phi 29 connector.

The connector of bacteriophage phi 29 is known to promote the viral prohead assembly, to bind DNA, and to drive DNA packaging into preformed viral shells in an RNA-dependent process. In this report, the phi 29 connector protein, p10, is shown to bind RNA in a sequence-independent fashion, and to possess an RNA recognition motif comprised approximately the region between residues 21 and 94 of the p10 sequence. Substitution mutants in specific amino acids of the RNA-binding domain obtained by site-directed mutagenesis showed that amino acids Phe23, His57, Phe59, and Tyr61 are critical for RNA binding and, subsequently, for DNA packaging into proheads. Proteolytic modified forms of the phi 29 connector have allowed us to conclude that the DNA- and RNA-binding domains are separated within the p10 sequence. It is also shown that RNA is stably associated to DNA-filled proheads during the DNA-packaging process.

DNA conformational change induced by the bacteriophage phi 29 connector.

Translocation of viral DNA inwards and outwards of the capsid of double-stranded DNA bacteriophages occurs through the connector, a key viral structure that is known to interact with DNA. It is shown here that phage phi 29 connector binds both linear and circular double-stranded DNA. However, DNA-mediated protection of phi 29 connectors against Staphylococcus aureus endoprotease V8 digestion suggests that binding to linear DNA is more stable than to circular DNA. Endoprotease V8-protection assays also suggest that the length of the linear DNA required to produce a stable phi 29 connector-DNA interaction is, at least, twice longer than the phi 29 connector channel. This result is confirmed by experiments of phi 29 connector-protection of DNA against DNase I digestion. Furthermore, DNA circularization assays indicate that phi 29 connectors restrain negative supercoiling when bound to linear DNA. This DNA conformational change is not observed upon binding to circular DNA and it could reflect the existence of some left-handed DNA coiling or DNA untwisting inside of the phi 29 connector channel.

Role of the amino-terminal domain of bacteriophage phi 29 connector in DNA binding and packaging.

The connector of bacteriophage phi 29 is required for prohead assembly, binds DNA, and drives DNA packaging into viral proheads. Limited proteolysis of the connector protein with endoproteinase Glu-C from Staphylococcus aureus V8 and chymotrypsin showed that a domain of the NH2-terminal region is involved in DNA binding and in the subsequent packaging into preformed proheads, but not in prohead assembly. Mutants in specific amino acids of the NH2-terminal domain, obtained by directed mutagenesis techniques, showed that the Ala1-Arg2-Lys3-Arg4 region of the connector is absolutely necessary for DNA packaging into the proheads as well as for efficient DNA binding.

Characterization of a versatile in vitro DNA-packaging system based on hybrid lambda/phi 29 proheads.

We have studied the assembly of bacteriophage lambda head proteins on the phage phi 29 connector to produce in vitro chimeric proheads, whose ability to package different types of DNA depends on the physical integrity of the phi 29 connector. Terminal protein-free phi 29 as well as nonviral DNAs have been shown to be efficiently packaged by this hybrid system. An RNA, that can be provided by any of the extracts used in the complementation mixture, was required for DNA packaging, both by the hybrid system as well as by the homologous lambda system. The DNA-packaging activity of RNase-treated proheads can be restored by adding a mixture of ribosomal RNAs. There is also a requirement for a minimal length of DNA to be stably packaged. The packaging protein p16 of phi 29 can replace the lambda terminase complex in the in vitro packaging system, both with the chimeric as well as genuine lambda proheads.

Production of lambda-phi 29 phage chimeras.

Proheads of bacteriophage lambda which carry the connector of phage phi 29 instead of that of lambda have been produced in vitro. These hybrid proheads have a structure similar to that of normal lambda proheads. Furthermore, the chimeric proheads can package both lambda and phi 29 DNA. These data show that the connector domains involved in both head assembly and DNA packaging are functionally similar. The DNA-containing lambda-phi 29 proheads can be complemented in vitro with phi 29 tails to yield infective particles capable of DNA transfer.

Bacteriophage T3 connector: three-dimensional structure and comparison with other viral head-tail connecting regions.

The bacteriophage T3 connector, which consists of 12 copies of protein gp8, has been studied by image processing of electron micrographs from negatively stained ordered aggregates. A three-dimensional reconstruction of T3 connectors was obtained by collection of tilted views and using the direct Fourier method, up to 2.3 nm resolution. The reconstructed unit cell contains two connectors whose main structural features are essentially identical, but facing in opposite directions. The T3 connector has a height of about 10.9 nm, with two clearly defined domains: a wider one 14.4 nm in diameter, with 12 morphological units in the periphery, and a narrower one, 9.7 nm in diameter. There is a channel clearly defined in the narrower domain that almost closes along the wider domain. Comparison of the three-dimensional structure obtained for the connector of phages T3 and phi 29, and that of the neck extracted from phage phi 29 particles, reveals striking similarities and significant differences. A model for a general connector to account for the common functions carried out by these viral assemblies is discussed together with the possible role of the channel for DNA translocation.

Three-dimensional reconstruction of the connector of bacteriophage phi 29 at 1.8 nm resolution.

The three-dimensional reconstruction of the connector of bacteriophage phi 29 has been obtained from tilt series of negatively stained tetragonal ordered aggregates under low-dose conditions and up to a resolution of (1/1.8) nm-1. These connectors are built up as dodecamers of only one structural polypeptide (p10). Two connectors form the crystal unit cell, each one facing in the opposite direction with respect to the plane of the crystal and partially overlapping. The main features of the two connectors that build the unit cell were essentially the same, although they were negatively stained in slightly different ways, probably due to their situations with respect to the carbon-coated support grid. The main features of the phi 29 connector structure revealed by this three-dimensional reconstruction are: the existence of two clearly defined domains, one with a diameter of around 14 nm and the other narrower (diameter approximately equal to 7.5 nm); an inner hole running all along the structure (around 7 to 8 nm in height) with a cylindrical profile and an average diameter of 4 nm; a general 6-fold symmetry along the whole structure and a 12-fold one in the wider domain; a clockwise twist of the more contrasted regions of both domains from the narrower towards the wider domain (the direction of DNA encapsidation). These features are compatible with an active role for the connector in the process of DNA packaging.