Search

Motivation Protein–Protein Interactions (PPI) are essentials for most cellular processes and thus, unveiling how proteins interact is a crucial question that can be better understood by identifying which residues are responsible for the interaction. Computational approaches are orders of magnitude cheaper and faster than experimental ones, leading to proliferation of multiple methods aimed to predict which residues belong to the interface of an interaction. Results We present BIPSPI, a new machine learning-based method for the prediction of partner-specific PPI sites. Contrary to most binding site prediction methods, the proposed approach takes into account a pair of interacting proteins rather than a single one in order to predict partner-specific binding sites. BIPSPI has been trained employing sequence-based and structural features from both protein partners of each complex compiled in the Protein–Protein Docking Benchmark version 5.0 and in an additional set independently compiled. Also, a version trained only on sequences has been developed. The performance of our approach has been assessed by a leave-one-out cross-validation over different benchmarks, outperforming state-of-the-art methods. Availability and implementation BIPSPI web server is freely available at http://bipspi.cnb.csic.es. BIPSPI code is available at https://github.com/bioinsilico/BIPSPI. Docker image is available at https://hub.docker.com/r/bioinsilico/bipspi/. Supplementary information Supplementary data are available at Bioinformatics online. [ABSTRACT FROM AUTHOR]

Background: Tilt series are commonly used in electron tomography as a means of collecting three-dimensional information from two-dimensional projections. A common problem encountered is the projection alignment prior to 3D reconstruction. Current alignment techniques usually employ gold particles or image derived markers to correctly align the images. When these markers are not present, correlation between adjacent views is used to align them. However, sequential pairwise correlation is prone to bias and the resulting alignment is not always optimal. Results: In this paper we introduce an algorithm to find regions of the tilt series which can be tracked within a subseries of the tilt series. These regions act as landmarks allowing the determination of the alignment parameters. We show our results with synthetic data as well as experimental cryo electron tomography. Conclusion: Our algorithm is able to correctly align a single-tilt tomographic series without the help of fiducial markers thanks to the detection of thousands of small image patches that can be tracked over a short number of images in the series. [ABSTRACT FROM AUTHOR]

We propose a novel nonnegative matrix factorization model that aims at finding localized, part-based, representations of nonnegative multivariate data items. Unlike the classical nonnegative matrix factorization (NMF) technique, this new model, denoted "nonsmooth nonnegative matrix factorization" (nsNMF), corresponds to the optimization of an unambiguous cost function designed to explicitly represent sparseness, in the form of nonsmoothness, which is controlled by a single parameter. In general, this method produces a set of basis and encoding vectors that are not only capable of representing the original data, but they also extract highly localized patterns, which generally lend themselves to improved interpretability. The properties of this new method are illustrated with several data sets. Comparisons to previously published methods show that the new nsNMF method has some advantages in keeping faithfulness to the data in the achieving a high degree of sparseness for both the estimated basis and the encoding vectors and in better interpretability of the factors. [ABSTRACT FROM AUTHOR]

Single-particle analysis by Cryo-electron microscopy (CryoEM) provides direct access to the conformation of each macromolecule. However, the image's signal-to-noise ratio is low, and some form of classification is usually performed at the image processing level to allow structural modeling. Classical classification methods imply the existence of a discrete number of structural conformations. However, new heterogeneity algorithms introduce a novel reconstruction paradigm, where every state is represented by a lower number of particles, potentially just one, allowing the estimation of conformational landscapes representing the different structural states a biomolecule explores. In this work, we present a novel deep learning-based method called HetSIREN. HetSIREN can fully reconstruct or refine a CryoEM volume in real space based on the structural information summarized in a conformational latent space. The unique characteristics that set HetSIREN apart start with the definition of the approach as a real space-based only method, a fact that allows spatially focused analysis, but also the introduction of a novel network architecture specifically designed to make use of meta-sinusoidal activations, with proven high analytics capacities. Continuing with innovations, HetSIREN can also refine the pose parameters of the images at the same time that it conditions the network with prior information/constraints on the maps, such as Total Variation and L 1 denoising, ultimately yielding cleaner volumes with high-quality structural features. Finally, but very importantly, HetSIREN addresses one of the most confusing issues in heterogeneity analysis, as it is the fact that real structural heterogeneity estimation is entangled with pose estimation (and to a lesser extent with CTF estimation), in this way, HetSIREN introduces a novel encoding architecture able to decouple pose and CTF information from the conformational landscape, resulting in more accurate and interpretable conformational latent spaces. We present results on computer-simulated data, public data from EMPIAR, and data from experimental systems currently being studied in our laboratories. An important finding is the sensitivity of the structure and dynamics of the SARS-CoV-2 Spike protein on the storage temperature., Competing Interests: Competing Interests The authors declare that they have no competing interests.

Electron tomography is a technique to obtain three-dimensional structural information of samples. However, the technique is limited by shifts occurring during acquisition that need to be corrected before the reconstruction process. In 2009, we proposed an approach for post-acquisition alignment of tilt series images. This approach was marker-free, based on patch tracking and integrated in free software. Here, we present improvements to the method to make it more reliable, stable and accurate. In addition, we modified the image formation model underlying the alignment procedure to include different deformations occurring during acquisition. We propose a new way to correct these computed deformations to obtain reconstructions with reduced artifacts. The new approach has demonstrated to improve the quality of the final 3D reconstruction, giving access to better defined structures for different transmission electron tomography methods: resin embedded STEM-tomography and cryo-TEM tomography. The method is freely available in TomoJ software., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2020 The Author(s).)

Resolution (global and local) is one of the most reported metrics of quality measurement in Single Particle Analysis (SPA). However, in electron tomography, the situation is different and its computation is not straightforward. Typically, resolution estimation is global and, therefore, reduces the assessment of a whole tomogram to a single number. However, it is known that tomogram quality is spatially variant. Still, up to our knowledge, a method to estimate local quality metrics in tomography is lacking. This work introduces MonoTomo , a method developed to estimate locally in a tomogram the highest reliable frequency component, expressed as a form of local resolution. The fundamentals lie in a local analysis of the density map via monogenic signals, which, in analogy to MonoRes , allows for local estimations. Results with experimental data show that the local resolution range that MonoTomo casts agrees with reported resolution values for experimental data sets, with the advantage of providing a local estimation. A range of applications of MonoTomo are suggested for further exploration., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2019 Published by Elsevier Inc.)

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.

RP4 TrbB, an essential component of the conjugative transfer apparatus of the broad-host-range plasmid RP4, is a member of the PulE protein superfamily involved in multicomponent machineries transporting macromolecules across the bacterial envelope. PulE-like proteins share several well conserved motifs, most notable a nucleoside triphosphate binding motif (P-loop). Helicobacter pylori HP0525 also belongs to the PulE superfamily and is encoded by the pathogenicity island cag, involved in the inflammatory response of infected gastric epithelial cells in mammals. The native molecular masses of TrbB and HP0525 as determined by gel filtration and glycerol gradient centrifugation suggested a homohexameric structure in the presence of ATP and Mg(2+). In the absence of nucleotides and bivalent cations, TrbB behaved as a tetramer whereas the hexameric state of HP0525 remained unaffected. Electron microscopy and image processing demonstrated that TrbB and HP0525 form ring-shaped complexes (diameter: 12 nm) with a central region (diameter: 3 nm) of low electron density when incubated in the presence of ATP and Mg(2+). However, the TrbB average image appeared to be more elliptical with strong twofold rotational symmetry whereas HP0525 complexes are regular hexagons. Six well defined triangle-shaped areas of high electron density were distinguishable in both cases. Covalent crosslinking of TrbB suggests that the hexameric ring is composed from a trimer of dimers, because only dimeric, tetrameric, and hexameric species were detectable. The toroidal structure of TrbB and HP0525 suggests that both proteins catalyze a repetitive process, most probably translocating a cognate substrate across the inner membrane.

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.

Background: Head-tail connectors are viral substructures that are very important in the viral morphogenetic cycle, having roles in the formation of the precursor capsid (prohead), DNA packaging, tail binding to the mature head and in the infection process. Structural information on the connector would, therefore, help us to understand how this structure is related to a multiplicity of functions., Results: Recombinant bacteriophage phi29 connectors have been crystallized in two-dimensional aggregates. An average projection image and a three-dimensional map have been obtained at 8 A and 10 A resolution, respectively, from untilted and tilted images of vitrified specimens of the two-dimensional crystals. The average projection image reveals a central mass surrounding a channel with 12 appendages protruding from the central mass. The three-dimensional map reveals a wide domain surrounded by 12 appendages that interact with the prohead vertex, and a narrow domain that interacts with the bacteriophage tail. At the junction of the two domains, 12 smaller appendages are visualized. A channel runs along the axis of the connector structure and is sufficiently wide to allow a double-stranded DNA molecule to pass through., Conclusions: The propeller-like structure of the phi29 connector strengthens the notion of the connector rotating during DNA packaging. The groove formed by the two lanes of large and small appendages may act as a rail to prevent the liberation of the connector from the prohead vertex during rotation.

The self-organizing tree algorithm (SOTA) was recently introduced to construct phylogenetic trees from biological sequences, based on the principles of Kohonen's self-organizing maps and on Fritzke's growing cell structures. SOTA is designed in such a way that the generation of new nodes can be stopped when the sequences assigned to a node are already above a certain similarity threshold. In this way a phylogenetic tree resolved at a high taxonomic level can be obtained. This capability is especially useful to classify sets of diversified sequences. SOTA was originally designed to analyze pre-aligned sequences. It is now adapted to be able to analyze patterns associated to the frequency of residues along a sequence, such as protein dipeptide composition and other n-gram compositions. In this work we show that the algorithm applied to these data is able to not only successfully construct phylogenetic trees of protein families, such as cytochrome c, triosephophate isomerase, and hemoglobin alpha chains, but also classify very diversified sequence data sets, such as a mixture of interleukins and their receptors.

Background: DNA helicases play a fundamental role in all aspects of nucleic acid metabolism and defects in these enzymes have been implicated in a number of inherited human disorders. DnaB is the major replicative DNA helicase in Escherichia coli and has been used as a model system for studying the structure and function of hexameric helicases. The native protein is a hexamer of identical subunits, which in solution forms a complex with six molecules of the loading protein DnaC. DnaB is delivered from this complex onto the DNA template, with the subsequent release of DnaC. We report here the structures of the DnaB helicase hexamer and its complex with DnaC under a defined set of experimental conditions, as determined by three-dimensional cryoelectron microscopy. It was hoped that the structures would provide insight into the mechanisms of helicase activity., Results: The DnaB structure reveals that six DnaB monomers assemble as three asymmetric dimers to form a polar, ring-like hexamer. The hexamer has two faces, one displaying threefold and the other sixfold symmetry. The six DnaC protomers bind tightly to the sixfold face of the DnaB hexamer. This is the first report of a three-dimensional structure of a helicase obtained using cryoelectron microscopy, and the first report of the structure of a helicase in complex with a loading protein., Conclusions: The structures of the DnaB helicase and its complex with DnaC reveal some interesting structural features relevant to helicase function and to the assembly of the two-protein complex. The results presented here provide a basis for a more complete understanding of the structure and function of these important proteins.

We propose a new type of unsupervised, growing, self-organizing neural network that expands itself by following the taxonomic relationships that exist among the sequences being classified. The binary tree topology of this neutral network, contrary to other more classical neural network topologies, permits an efficient classification of sequences. The growing nature of this procedure allows to stop it at the desired taxonomic level without the necessity of waiting until a complete phylogenetic tree is produced. This novel approach presents a number of other interesting properties, such as a time for convergence which is, approximately, a lineal function of the number of sequences. Computer simulation and a real example show that the algorithm accurately finds the phylogenetic tree that relates the data. All this makes the neural network presented here an excellent tool for phylogenetic analysis of a large number of sequences.

We present an approach to the study of the relationships between biological sequences and structures applying image analysis methods to dotplots. We introduce a set of analytical tools based on different types of digital image-processing filters that are new within the context of dotplots. We have reformulated some of the usual approaches in dotplot analysis as mathematical operations on images within the framework of mathematical morphology. An X-Window-based implementation of this new approach has been developed and is available by anonymous FTP.