Multidimensional replica-exchange method for free-energy calculations

We have developed a new simulation algorithm for free-energy calculations. The method is a multidimensional extension of the replica-exchange method. While pairs of replicas with different temperatures are exchanged during the simulation in the original replica-exchange method, pairs of replicas with different temperatures and/or different parameters of the potential energy are exchanged in the new algorithm. This greatly enhances the sampling of the conformational space and allows accurate calculations of free energy in a wide temperature range from a single simulation run, using the weighted histogram analysis method.Comment: 13 pages, (ReVTeX), 9 figures. J. Chem. Phys. 113 (2000), in pres

Highly Charged Ion Production Using an Electrode in Biased and Floating Modes

One of the most popular ways to obtain higher beam intensities in ECR ion sources is to install an electrode (usually disc) into the plasma chamber. Examined this method in detail we found that majority of the groups observed the beam intensity improvement by supplying a suitable biased voltage to the electrode and an electron current was injected into the plasma. A few groups observed the enhancement, however, when the electrode operated at floating potential - without being an electron donor. Only a few (and sometimes contradictionary) information was found on the optimised properties of the electrodes, i.e. position, dimension, shape, material. In spite of the great success of the "biased-disc" method, the mechanism is still not completely clear. In this contribution, as one step of understanding, we examine what condition we observed the above mentioned two modes. The experiments were performed at the 18 GHz RIKEN and at the 14.5 GHz ATOMKI ECR ion sources. It was found that effect of the electrode is strongly depends on the local plasma parameters and on the position of the electrode. At certain mirror ratios and electrode positions we needed to negatively bias the electrode and inject electrons into the plasma. The electrode operated as an electron source (Electron Donor ED mode). At higher mirror ratios and other axial positions the electrode works by directly changing the plasma potential dip (Potential Tuner PT mode). These two modes were checked and successfully found both in continuos and in pulsed mode operation. In both (ED and PT) modes we generated higher highly charged ion currents in the RIKEN-ECRIS than without the electrode

Amino acid "little Big Bang": Representing amino acid substitution matrices as dot products of Euclidian vectors

<p>Abstract</p> <p>Background</p> <p>Sequence comparisons make use of a one-letter representation for amino acids, the necessary quantitative information being supplied by the substitution matrices. This paper deals with the problem of finding a representation that provides a comprehensive description of amino acid intrinsic properties consistent with the substitution matrices.</p> <p>Results</p> <p>We present a Euclidian vector representation of the amino acids, obtained by the singular value decomposition of the substitution matrices. The substitution matrix entries correspond to the dot product of amino acid vectors. We apply this vector encoding to the study of the relative importance of various amino acid physicochemical properties upon the substitution matrices. We also characterize and compare the PAM and BLOSUM series substitution matrices.</p> <p>Conclusions</p> <p>This vector encoding introduces a Euclidian metric in the amino acid space, consistent with substitution matrices. Such a numerical description of the amino acid is useful when intrinsic properties of amino acids are necessary, for instance, building sequence profiles or finding consensus sequences, using machine learning algorithms such as Support Vector Machine and Neural Networks algorithms.</p

phenix.mr_rosetta: molecular replacement and model rebuilding with Phenix and Rosetta.

The combination of algorithms from the structure-modeling field with those of crystallographic structure determination can broaden the range of templates that are useful for structure determination by the method of molecular replacement. Automated tools in phenix.mr_rosetta simplify the application of these combined approaches by integrating Phenix crystallographic algorithms and Rosetta structure-modeling algorithms and by systematically generating and evaluating models with a combination of these methods. The phenix.mr_rosetta algorithms can be used to automatically determine challenging structures. The approaches used in phenix.mr_rosetta are described along with examples that show roles that structure-modeling can play in molecular replacement

Prediction of binding hot spot residues by using structural and evolutionary parameters

In this work, we present a method for predicting hot spot residues by using a set of structural and evolutionary parameters. Unlike previous studies, we use a set of parameters which do not depend on the structure of the protein in complex, so that the predictor can also be used when the interface region is unknown. Despite the fact that no information concerning proteins in complex is used for prediction, the application of the method to a compiled dataset described in the literature achieved a performance of 60.4%, as measured by F-Measure, corresponding to a recall of 78.1% and a precision of 49.5%. This result is higher than those reported by previous studies using the same data set

Aggregation and travelling wave dynamics in a two-population model of cancer cell growth and invasion

Funding: Engineering and Physical Sciences Research Council (UK) grant numbers EP/L504932/1 (VB), EP/K033689/1 (RE).Cells adhere to each other and to the extracellular matrix (ECM) through protein molecules on the surface of the cells. The breaking and forming of adhesive bonds, a process critical in cancer invasion and metas- tasis, can be influenced by the mutation of cancer cells. In this paper, we develop a nonlocal mathematical model describing cancer cell invasion and movement as a result of integrin-controlled cell-cell adhesion and cell-matrix adhesion, for two cancer cell populations with different levels of mutation. The partial differential equations for cell dynamics are coupled with ordinary differential equations describing the extracellular matrix (ECM) degradation and the production and decay of integrins. We use this model to investigate the role of cancer mutation on the possibility of cancer clonal competition with alternating dominance, or even competitive exclusion (phenomena observed experimentally). We discuss different possible cell aggregation patterns, as well as travelling wave patterns. In regard to the travelling waves, we investigate the effect of cancer mutation rate on the speed of cancer invasion.Publisher PDFPeer reviewe

The interaction between caveolin-1 and Rho-GTPases promotes metastasis by controlling the expression of alpha5-integrin and the activation of Src, Ras and Erk

Proteins containing a caveolin-binding domain (CBD), such as the Rho-GTPases, can interact with caveolin-1 (Cav1) through its caveolin scaffold domain. Rho-GTPases are important regulators of p130Cas, which is crucial for both normal cell migration and Src kinase-mediated metastasis of cancer cells. However, although Rho-GTPases (particularly RhoC) and Cav1 have been linked to cancer progression and metastasis, the underlying molecular mechanisms are largely unknown. To investigate the function of Cav1–Rho-GTPase interaction in metastasis, we disrupted Cav1–Rho-GTPase binding in melanoma and mammary epithelial tumor cells by overexpressing CBD, and examined the loss-of-function of RhoC in metastatic cancer cells. Cancer cells overexpressing CBD or lacking RhoC had reduced p130Cas phosphorylation and Rac1 activation, resulting in an inhibition of migration and invasion in vitro. The activity of Src and the activation of its downstream targets FAK, Pyk2, Ras and extracellular signal-regulated kinase (Erk)1/2 were also impaired. A reduction in α5-integrin expression, which is required for binding to fibronectin and thus cell migration and survival, was observed in CBD-expressing cells and cells lacking RhoC. As a result of these defects, CBD-expressing melanoma cells had a reduced ability to metastasize in recipient mice, and impaired extravasation and survival in secondary sites in chicken embryos. Our data indicate that interaction between Cav1 and Rho-GTPases (most likely RhoC but not RhoA) promotes metastasis by stimulating α5-integrin expression and regulating the Src-dependent activation of p130Cas/Rac1, FAK/Pyk2 and Ras/Erk1/2 signaling cascades

Nature of protein family signatures: Insights from singular value analysis of position-specific scoring matrices

Position-specific scoring matrices (PSSMs) are useful for detecting weak homology in protein sequence analysis, and they are thought to contain some essential signatures of the protein families. In order to elucidate what kind of ingredients constitute such family-specific signatures, we apply singular value decomposition to a set of PSSMs and examine the properties of dominant right and left singular vectors. The first right singular vectors were correlated with various amino acid indices including relative mutability, amino acid composition in protein interior, hydropathy, or turn propensity, depending on proteins. A significant correlation between the first left singular vector and a measure of site conservation was observed. It is shown that the contribution of the first singular component to the PSSMs act to disfavor potentially but falsely functionally important residues at conserved sites. The second right singular vectors were highly correlated with hydrophobicity scales, and the corresponding left singular vectors with contact numbers of protein structures. It is suggested that sequence alignment with a PSSM is essentially equivalent to threading supplemented with functional information. The presented method may be used to separate functionally important sites from structurally important ones, and thus it may be a useful tool for predicting protein functions.Comment: 22 pages, 7 figures, 4 table

Enhanced and effective conformational sampling of protein molecular systems for their free energy landscapes

Protein folding and protein–ligand docking have long persisted as important subjects in biophysics. Using multicanonical molecular dynamics (McMD) simulations with realistic expressions, i.e., all-atom protein models and an explicit solvent, free-energy landscapes have been computed for several systems, such as the folding of peptides/proteins composed of a few amino acids up to nearly 60 amino-acid residues, protein–ligand interactions, and coupled folding and binding of intrinsically disordered proteins. Recent progress in conformational sampling and its applications to biophysical systems are reviewed in this report, including descriptions of several outstanding studies. In addition, an algorithm and detailed procedures used for multicanonical sampling are presented along with the methodology of adaptive umbrella sampling. Both methods control the simulation so that low-probability regions along a reaction coordinate are sampled frequently. The reaction coordinate is the potential energy for multicanonical sampling and is a structural identifier for adaptive umbrella sampling. One might imagine that this probability control invariably enhances conformational transitions among distinct stable states, but this study examines the enhanced conformational sampling of a simple system and shows that reasonably well-controlled sampling slows the transitions. This slowing is induced by a rapid change of entropy along the reaction coordinate. We then provide a recipe to speed up the sampling by loosening the rapid change of entropy. Finally, we report all-atom McMD simulation results of various biophysical systems in an explicit solvent