Digital dermoscopy monitoring of melanocytic lesions: Two novel calculators combining static and dynamic features to identify melanoma.

BACKGROUND: Early diagnosis is the most effective intervention to improve the prognosis of cutaneous melanoma. Even though the introduction of dermoscopy has improved the diagnostic accuracy, it can still be difficult to distinguish some melanomas from benign melanocytic lesions. Digital dermoscopy monitoring can identify dynamic changes of melanocytic lesions: To date, some algorithms were proposed, but a universally accepted one is still lacking. OBJECTIVES: To identify independent predictive variables associated with the diagnosis of cutaneous melanoma and develop a multivariable dermoscopic prediction model able to discriminate benign from malignant melanocytic lesions undergoing digital dermoscopy monitoring. METHODS: We collected dermoscopic images of melanocytic lesions excised after dermoscopy monitoring and carried out static and dynamic evaluations of dermoscopic features. We built two multivariable predictive models based on logistic regression and random forest. RESULTS: We evaluated 173 lesions (65 cutaneous melanomas and 108 nevi). Forty-two melanomas were in situ, and the median thickness of invasive melanomas was 0.35 mm. The median follow-up time was 9.8 months for melanomas and 9.1 for nevi. The logistic regression and random forest models performed with AUC values of 0.87 and 0.89, respectively, were substantially higher than those of the static evaluation models (ABCD TDS score, 0.57; 7-point checklist, 0.59). Finally, we built two risk calculators, which translate the proposed models into user-friendly applications, to assist clinicians in the decision-making process. CONCLUSIONS: The present study demonstrates that the integration of dynamic and static evaluations of melanocytic lesions is a safe approach that can significantly boost the diagnostic accuracy for cutaneous melanoma. We propose two diagnostic tools that significantly increase the accuracy in discriminating melanoma from nevi during digital dermoscopy monitoring.

The effect of backbone on the small-world properties of protein contact maps.

In the last years, small-world behavior has been extensively described for proteins, when they are represented by the undirected graph defined by the inter-residue protein contacts. By adopting this representation it was possible to compute the average clustering coefficient (C) and characteristic path length (L) of protein structures, and their values were found to be similar to those of graphs characterized by small-world topology. In this comment, we analyze a large set of non-redundant protein structures (1753) and show that by randomly mimicking the protein collapse, the covalent structure of the protein chain significantly contributes to the small-world behavior of the inter-residue contact graphs. When protein graphs are generated, imposing constraints similar to those induced by the backbone connectivity, their characteristic path lengths and clustering coefficients are indistinguishable from those computed using the real contact maps showing that L and C values cannot be used for 'protein fingerprinting'. Moreover we verified that these results are independent of the selected protein representations, residue composition and protein secondary structures.

Microbiota and environmental stress: how pollution affects microbial communities in Manila clams.

Given the crucial role of microbiota in host development, health, and environmental interactions, genomic analyses focusing on host-microbiota interactions should certainly be considered in the investigation of the adaptive mechanisms to environmental stress. Recently, several studies suggested that microbiota associated to digestive tract is a key, although still not fully understood, player that must be considered to assess the toxicity of environmental contaminants. Bacteria-dependent metabolism of xenobiotics may indeed modulate the host toxicity. Conversely, environmental variables (including pollution) may alter the microbial community and/or its metabolic activity leading to host physiological alterations that may contribute to their toxicity. Here, 16s rRNA gene amplicon sequencing has been applied to characterize the hepatopancreas microbiota composition of the Manila clam, Ruditapes philippinarum. The animals were collected in the Venice lagoon area, which is subject to different anthropogenic pressures, mainly represented by the industrial activities of Porto Marghera (PM). Seasonal and geographic differences in clam microbiotas were explored and linked to host response to chemical stress identified in a previous study at the transcriptome level, establishing potential interactions among hosts, microbes, and environmental parameters. The obtained results showed the recurrent presence of putatively detoxifying bacterial taxa in PM clams during winter and over-representation of several metabolic pathways involved in xenobiotic degradation, which suggested the potential for host-microbial synergistic detoxifying actions. Strong interaction between seasonal and chemically-induced responses was also observed, which partially obscured such potentially synergistic actions. Seasonal variables and exposure to toxicants were therefore shown to interact and substantially affect clam microbiota, which appeared to mirror host response to environmental variation. It is clear that understanding how animals respond to chemical stress cannot ignore a key component of such response, the microbiota.

Characterization of 9-aminoacridine interaction with chromatophore membranes and modelling of the probe response to artificially induced transmembrane delta pH values.

We analyze the adsorption of the fluorescent monoamine 9-aminoacridine to the membrane phase of photosynthetic chromatophores, in the physiological interval of pH values ranging from 5.5 to 8.5 and at ionic strengths of 0.005 and 0.150 M. The interaction of the probe with the membrane phase is described with S-shaped isotherms of the Hill type and is modulated by electrostatic effects as modelled with the Gouy-Chapman-Boltzman theory. This description is consistent with different values of the surface change density of the chromatophore membranes decreasing from about 1.3 x 10(-3) to about 0.5 x 10(-3) e-/A2, on changing the pH from 8.5/7.5 to 6.5/5.5, respectively. Furthermore we show that, when the free concentrations of the probe in the inner and outer vesicle compartments are computed from the adsorbing isotherms at the proper pH values, the model considering the equilibrium distribution of the neutral monoamine following the onset of a delta pH is sufficient to describe the dependence of the artificially induced transmembrane delta pH values on the observed quenching of the probe fluorescence.

Topology prediction for helical transmembrane proteins at 86% accuracy.

Previously, we introduced a neural network system predicting locations of transmembrane helices (HTMs) based on evolutionary profiles (PHDhtm, Rost B, Casadio R, Fariselli P, Sander C, 1995, Protein Sci 4:521-533). Here, we describe an improvement and an extension of that system. The improvement is achieved by a dynamic programming-like algorithm that optimizes helices compatible with the neural network output. The extension is the prediction of topology (orientation of first loop region with respect to membrane) by applying to the refined prediction the observation that positively charged residues are more abundant in extra-cytoplasmic regions. Furthermore, we introduce a method to reduce the number of false positives, i.e., proteins falsely predicted with membrane helices. The evaluation of prediction accuracy is based on a cross-validation and a double-blind test set (in total 131 proteins). The final method appears to be more accurate than other methods published: (1) For almost 89% (+/-3%) of the test proteins, all HTMs are predicted correctly. (2) For more than 86% (+/-3%) of the proteins, topology is predicted correctly. (3) We define reliability indices that correlate with prediction accuracy: for one half of the proteins, segment accuracy raises to 98%; and for two-thirds, accuracy of topology prediction is 95%. (4) The rate of proteins for which HTMs are predicted falsely is below 2% (+/-1%). Finally, the method is applied to 1,616 sequences of Haemophilus influenzae. We predict 19% of the genome sequences to contain one or more HTMs. This appears to be lower than what we predicted previously for the yeast VIII chromosome (about 25%).

Transmembrane helices predicted at 95% accuracy.

We describe a neural network system that predicts the locations of transmembrane helices in integral membrane proteins. By using evolutionary information as input to the network system, the method significantly improved on a previously published neural network prediction method that had been based on single sequence information. The input data were derived from multiple alignments for each position in a window of 13 adjacent residues: amino acid frequency, conservation weights, number of insertions and deletions, and position of the window with respect to the ends of the protein chain. Additional input was the amino acid composition and length of the whole protein. A rigorous cross-validation test on 69 proteins with experimentally determined locations of transmembrane segments yielded an overall two-state per-residue accuracy of 95%. About 94% of all segments were predicted correctly. When applied to known globular proteins as a negative control, the network system incorrectly predicted fewer than 5% of globular proteins as having transmembrane helices. The method was applied to all 269 open reading frames from the complete yeast VIII chromosome. For 59 of these, at least two transmembrane helices were predicted. Thus, the prediction is that about one-fourth of all proteins from yeast VIII contain one transmembrane helix, and some 20%, more than one.

Quantum mechanical analysis of oxygenated and deoxygenated states of hemocyanin: theoretical clues for a plausible allosteric model of oxygen binding.

In this work with ab initio computations, we describe relevant interactions between protein active sites and ligands, using as a test case arthropod hemocyanins. A computational analysis of models corresponding to the oxygenated and deoxygenated forms of the hemocyanin active site is performed using the Density Functional Theory approach. We characterize the electron density distribution of the binding site with and without bound oxygen in relation to the geometry, which stems out of the crystals of three hemocyanin proteins, namely the oxygenated form from the horseshoe crab Limulus polyphemus, and the deoxygenated forms, respectively, from the same source and from another arthropod, the spiny lobster Panulirus interruplus. Comparison of the three available crystals indicate structural differences at the oxygen binding site, which cannot be explained only by the presence and absence of the oxygen ligand, since the geometry of the ligand site of the deoxygenated Panulirus hemocyanin is rather similar to that of the oxygenated Limulus protein. This finding was interpreted in the frame of a mechanism of allosteric regulation for oxygen binding. However, the cooperative mechanism, which is experimentally well documented, is only partially supported by crystallographic data, since no oxygenated crystal of Panulirus hemocyanin is presently available. We address the following question: is the local ligand geometry responsible for the difference of the dicopper distance observed in the two deoxygenated forms of hemocyanin or is it necessary to advocate the allosteric regulation of the active site conformations in order to reconcile the different crystal forms? We find that the difference of the dicopper distance between the two deoxygenated hemocyanins is not due to the small differences of ligand geometry found in the crystals and conclude that it must be therefore stabilized by the whole protein tertiary structure.

Prediction of the transmembrane regions of beta-barrel membrane proteins with a neural network-based predictor.

A method based on neural networks is trained and tested on a nonredundant set of beta-barrel membrane proteins known at atomic resolution with a jackknife procedure. The method predicts the topography of transmembrane beta strands with residue accuracy as high as 78% when evolutionary information is used as input to the network. Of the transmembrane beta-strands included in the training set, 93% are correctly assigned. The predictor includes an algorithm of model optimization, based on dynamic programming, that correctly models eight out of the 11 proteins present in the training/testing set. In addition, protein topology is assigned on the basis of the location of the longest loops in the models. We propose this as a general method to fill the gap of the prediction of beta-barrel membrane proteins.

A neural network approach to evaluate fold recognition results.

Fold recognition techniques assist the exploration of protein structures, and web-based servers are part of the standard set of tools used in the analysis of biochemical problems. Despite their success, current methods are only able to predict the correct fold in a relatively small number of cases. We propose an approach that improves the selection of correct folds from among the results of two methods implemented as web servers (SAMT99 and 3DPSSM). Our approach is based on the training of a system of neural networks with models generated by the servers and a set of associated characteristics such as the quality of the sequence-structure alignment, distribution of sequence features (sequence-conserved positions and apolar residues), and compactness of the resulting models. Our results show that it is possible to detect adequate folds to model 80% of the sequences with a high level of confidence. The improvements achieved by taking into account sequence characteristics open the door to future improvements by directly including such factors in the step of model generation. This approach has been implemented as an automatic system LIBELLULA, available as a public web server at http://www.pdg.cnb.uam.es/servers/libellula.html.

Can functional regions of proteins be predicted from their coding sequences? The case study of G-protein coupled receptors.

A filter based on a set of unsupervised neural networks trained with a winner-take-all strategy discloses signals along the coding sequences of G-protein coupled receptors. By comparing with the existing experimental data it appears that these signals correlate with putative functional domains of the proteins. After protein alignment within subfamilies, signals cluster in protein regions which, according to the presently available experimental results, are described as possible functional domains of the folded proteins. The mapping procedure reveals characteristic regions in the coding sequences common and/or characteristic of the receptor subtype. This is particularly noticeable for the third cytoplasmic loop, which is likely to be involved in the molecular coupling of all the subfamilies with G-proteins. The results indicate that our mapping can highlight intrinsic representative features of the coding sequences which, in the case of G-protein coupled receptors, are characteristic of protein functional regions and suggest a possible application of the filter for predicting functional determinants in proteins starting from the coding sequence.

A high diffusion coefficient for coenzyme Q10 might be related to a folded structure.

We measured the lateral diffusion of different coenzyme Q homologues and analogues in model lipid vesicles using the fluorescence collisional quenching technique with pyrene derivatives and found diffusion coefficients in the range of 10(-6) cm2/s. Theoretical diffusion coefficients for these highly hydrophobic components were calculated according to the free volume theory. An important parameter in the free volume theory is the relative dimension between diffusant and solvent: a molecular dynamics computer simulation of the coenzymes yielded their most probable geometries and volumes and revealed surprisingly similar sizes of the short and long homologues, due to a folded structure of the isoprenoid chain in the latter, with a length for coenzyme Q10 of 21 A. Using this information we were able to calculate diffusion coefficients in the range of 10(-6) cm2/s, in good agreement with those found experimentally.

Neural networks predict protein folding and structure: artificial intelligence faces biomolecular complexity.

In the genomic era DNA sequencing is increasing our knowledge of the molecular structure of genetic codes from bacteria to man at a hyperbolic rate. Billions of nucleotides and millions of aminoacids are already filling the electronic files of the data bases presently available, which contain a tremendous amount of information on the most biologically relevant macromolecules, such as DNA, RNA and proteins. The most urgent problem originates from the need to single out the relevant information amidst a wealth of general features. Intelligent tools are therefore needed to optimise the search. Data mining for sequence analysis in biotechnology has been substantially aided by the development of new powerful methods borrowed from the machine learning approach. In this paper we discuss the application of artificial feedforward neural networks to deal with some fundamental problems tied with the folding process and the structure-function relationship in proteins.

Protein structure prediction and biomolecular recognition: from protein sequence to peptidomimetic design with the human beta3 integrin.

Computational tools can bridge the gap between sequence and protein 3D structure based on the notion that information is to be retrieved from the databases and that knowledge-based methods can help in approaching a solution of the protein-folding problem. To this aim our group has implemented neural network-based predictors capable of performing with some success in different tasks, including predictions of the secondary structure of globular and membrane proteins, the topology of membrane proteins and porins and stable alpha-helical segments suited for protein design. Moreover we have developed methods for predicting contact maps in proteins and the probability of finding a cysteine in a disulfide bridge, tools which can contribute to the goal of predicting the 3D structure starting from the sequence (the so called ab initio prediction). All our predictors take advantage of evolution information derived from the structural alignments of homologous (evolutionary related) proteins and taken from the sequence and structure databases. When it is necessary to build models for proteins of unknown spatial structure, which have very little homology with other proteins of known structure, non-standard techniques need to be developed and the tools for protein structure predictions may help in protein modeling. The results of a recent simulation performed in our lab highlights the role of high performing computing technology and of tools of computational biology in protein modeling and peptidomimetic design.

The prediction of protein-protein interacting sites in genome-wide protein interaction networks: the test case of the human cell cycle.

In this paper we aim at investigating possible correlations between the number of putative interaction patches of a given protein, as inferred by an algorithm that we have developed, and its degree (number of edges of the protein node in a protein interaction network). We focus on the human cell cycle that, as compared with other biological processes, comprises the largest number of proteins whose structure is known at atomic resolution both as monomers and as interacting complexes. For predicting interaction patches we specifically develop a HM-SVM based method reaching 71% overall accuracy with a correlation coefficient value equal to 0.43 on a non redundant set of protein complexes. To test the biological meaning of our predictions, we also explore whether interacting patches contain energetically important residues and/or disease related mutations and find that predicted patches are endowed with both features. Based on this, we propose that mapping the protein with all the predicted interaction patches bridges the molecule to the interactome at the cell level. To test our hypothesis we downloaded interaction data from interaction data bases and find that the number of predicted interaction patches significantly correlates (Pearson correlation value >0.3) with the number of the known interactions (edges) per protein in the human interactome, as contained in MINT and IntAct. We also show that the correlation increases (Pearson correlation value >0.5) when the subcellular co-localization and the co-expression levels of the interacting partners are taken into account.

Prediction of the number of residue contacts in proteins.

Knowing the number of residue contacts in a protein is crucial for deriving constraints useful in modeling protein folding and/or scoring remote homology search. Here we focus on the prediction of residue contacts and show that this figure can be predicted with a neural network based method. The accuracy of the prediction is 12 percentage points higher than that of a simple statistical method. The neural network is used to discriminate between two different states of residue contacts, characterized by a contact number higher or lower than the average value of the residue distribution. When evolutionary information is taken into account, our method correctly predicts 69% of the residue states in the data base and it adds to the prediction of residue solvent accessibility. The predictor is available at htpp://www.biocomp.unibo.it

A neural network based predictor of residue contacts in proteins.

We describe a method based on neural networks for predicting contact maps of proteins using as input chemicophysical and evolutionary information. Neural networks are trained on a data set comprising the contact maps of 200 non-homologous proteins of well resolved three-dimensional structures. The systems learn the association rules between the covalent structure of each protein and its correspondent contact map by means of a standard back propagation algorithm. Validation of the predictor on the training set and on 408 proteins of known structure which are not homologous to those contained in the training set indicate that this method scores higher than statistical approaches previously described and based on correlated mutations and sequence information.

HTP: a neural network-based method for predicting the topology of helical transmembrane domains in proteins.

In this paper we describe a microcomputer program (HTP) for predicting the location and orientation of alpha-helical transmembrane segments in integral membrane proteins. HTP is a neural network-based tool which gives as output the protein membrane topology based on the statistical propensity of residues to be located in external and internal loops. This method, which uses single protein sequences as input to the network system, correctly predicts the topology of 71 out of 92 membrane proteins of putative membrane orientation, independently of the protein source.

RCNPRED: prediction of the residue co-ordination numbers in proteins.

UNLABELLED: The RCNPRED server implements a neural network-based method to predict the co-ordination numbers of residues starting from the protein sequence. Using evolutionary information as input, RCNPRED predicts the residue states of the proteins in the database with 69% accuracy and scores 12 percentage points higher than a simple statistical method. Moreover the server implements a neural network to predict the relative solvent accessibility of each residue. A protein sequence can be directly submitted to RCNPRED: residue co-ordination numbers and solvent accessibility for each chain are returned via e-mail. AVAILABILITY: Freely available to non-commercial users at http://prion.biocomp.unibo.it/rcnpred.html.

Prediction of disulfide connectivity in proteins.

MOTIVATION: A major problem in protein structure prediction is the correct location of disulfide bridges in cysteine-rich proteins. In protein-folding prediction, the location of disulfide bridges can strongly reduce the search in the conformational space. Therefore the correct prediction of the disulfide connectivity starting from the protein residue sequence may also help in predicting its 3D structure. RESULTS: In this paper we equate the problem of predicting the disulfide connectivity in proteins to a problem of finding the graph matching with the maximum weight. The graph vertices are the residues of cysteine-forming disulfide bridges, and the weight edges are contact potentials. In order to solve this problem we develop and test different residue contact potentials. The best performing one, based on the Edmonds-Gabow algorithm and Monte-Carlo simulated annealing reaches an accuracy significantly higher than that obtained with a general mean force contact potential. Significantly, in the case of proteins with four disulfide bonds in the structure, the accuracy is 17 times higher than that of a random predictor. The method presented here can be used to locate putative disulfide bridges in protein-folding. AVAILABILITY: The program is available upon request from the authors. CONTACT: Casadio@alma.unibo.it; Piero@biocomp.unibo.it.

Predicting secondary structures of membrane proteins with neural networks.

Back-propagation, feed-forward neural networks are used to predict the secondary structures of membrane proteins whose structures are known to atomic resolution. These networks are trained on globular proteins and can predict globular protein structures having no homology to those of the training set with correlation coefficients (Ci) of 0.45, 0.32 and 0.43 for alpha-helix, beta-strand and random coil structures, respectively. When tested on membrane proteins, neural networks trained on globular proteins do, on average, correctly predict (Qi) 62%, 38% and 69% of the residues in the alpha-helix, beta-strand and random coil structures. These scores rank higher than those obtained with the currently used statistical methods and are comparable to those obtained with the joint approaches tested so far on membrane proteins. The lower success score for beta-strand as compared to the other structures suggests that the sample of beta-strand patterns contained in the training set is less representative than those of alpha-helix and random coil. Our analysis, which includes the effects of the network parameters and of the structural composition of the training set on the prediction, shows that regular patterns of secondary structures can be successfully extrapolated from globular to membrane proteins.