Structure optimization in a three-dimensional off-lattice protein model

We studied a three-dimensional off-lattice AB model with two species of monomers, hydrophobic (A) and hydrophilic (B), and present two optimization algorithms: face-centered-cubic (FCC)-lattice pruned-enriched-Rosenbluth method (PERM) and subsequent conjugate gradient (PERM++) minimization and heuristic conjugate gradient (HCG) simulation based on "off-trap" strategy. In PERM++, we apply the PERM to the FCC-lattice to produce the initial conformation, and conjugate gradient minimization is then used to reach the minimum energy state. Both algorithms have been tested in the three-dimensional AB model for all sequences with lengths 13 < or = n < or = 55. The numerical results show that the proposed methods are very promising for finding the ground states of proteins. In several cases, we renew the putative ground states energy values.     

A branch and bound algorithm for the protein folding problem in the HP lattice model

A branch and bound algorithm is proposed for the two-dimensional protein folding problem in the HP lattice model. In this algorithm, the benefit of each possible location of hydrophobic monomers is evaluated and only promising nodes are kept for further branching at each level. The proposed algorithm is compared with other well-known methods for 10 benchmark sequences with lengths ranging from 20 to 100 monomers. The results indicate that our method is a very efficient and promising tool for the protein folding problem.     

Chaotic Artificial Bee Colony algorithm: A new approach to the problem of minimization of energy of the 3D protein structure

Prediction of the three-dimensional structure of a protein from its amino acid sequence can be considered as a global optimization problem. In this paper, the Chaotic Artificial Bee Colony (CABC) algorithm was introduced and applied to 3D protein structure prediction. Based on the 3D off-lattice AB model, the CABC algorithm combines global search and local search of the Artificial Bee Colony (ABC) algorithm with the chaotic search algorithm to avoid the problem of premature convergence and easily trapping the local optimum solution. The experiments carried out with the popular Fibonacci sequences demonstrate that the proposed algorithm provides an effective and high-performance method for protein structure prediction.     

Biophysical characterization of anticoagulant hemextin AB complex from the venom of snake Hemachatus haemachatus

Hemextin AB complex from the venom of Hemachatus haemachatus is the first known natural anticoagulant that specifically inhibits the enzymatic activity of blood coagulation factor VIIa in the absence of factor Xa. It is also the only known heterotetrameric complex of two three-finger toxins. Individually only hemextin A has mild anticoagulant activity, whereas hemextin B is inactive. However, hemextin B synergistically enhances the anticoagulant activity of hemextin A and their complex exhibits potent anticoagulant activity. In this study we characterized the nature of molecular interactions leading to the complex formation. Circular dichroism studies indicate the stabilization of β-sheet in the complex. Hemextin AB complex has an increased apparent molecular diameter in both gas and liquid phase techniques. The complex formation is enthalpically favorable and entropically unfavorable with a negative change in the heat capacity. Thus, the anticoagulant complex shows less structural flexibility than individual subunits. Both electrostatic and hydrophobic interactions are important for the complexation; the former driving the process and the latter helping in the stabilization of the tetramer. The tetramer dissociates into dimers and monomers with the increase in the ionic strength of the solution and also with increase in the glycerol concentration in the buffer. The two dimers formed under each of these conditions display distinct differences in their apparent molecular diameters and anticoagulant properties. Based on these results, we have proposed a model for this unique anticoagulant complex.     

Molecular modelling of a template substitute and monomers used in molecular imprinting for aflatoxin B1 micro-HPLC analysis

The preparation of molecularly imprinted polymers (MIPs) involves the polymerisation of functional monomers in the presence of template molecules. 5,7-Dimethoxycoumarin (DMC) was found to be a structural analogue for aflatoxin B1 (AB1) and serves as its substitute in a grafting solution for the MIP synthesis. It was found that both methacrylic acid and allylamine are functional monomers which could provide a similar binding towards AB1 and DMC.     

Conformational dynamics of a lipid-interacting protein: MD simulations of saposin B

Saposin B is a water soluble alpha-helical protein which can bind to membranes and extract selected lipids, especially cerebroside sulfates. The X-ray structure of saposin B is homodimeric. There are two conformations of the dimer in the crystal-one with a closed central cavity (the AB dimer) and one (the CD dimer) with a more open cavity. We have conducted a series of short (5 ns) molecular dynamics simulations of saposin B, starting from both the AB and CD conformations and with/without bound lipid and/or water molecules within the central hydrophobic cavity. The more open (CD) dimer showed greater conformational drift than the AB dimer. The conformational drift was also somewhat higher in the absence of bound lipid. Two more extended (30 ns) simulations of AB and CD dimers were performed and analyzed in terms of changes in intersubunit packing within the dimers. The AB dimer remained largely unchanged in conformation over the duration of the extended simulation. In contrast, the CD dimer underwent a substantial conformational change corresponding to a 'scissor' motion of the two monomers so as to compress the central cavity to a more closed conformation than that seen in the AB dimer structure. A H-bond between the Q53 and Y54 side chains of the alpha3 helices of the two opposing monomers seems to hold the dimer in this 'scissor-closed' conformation. We suggest that a cycle of conformational changes, expanding and compressing the central cavity of the saposin B dimer, may play a key role in facilitating lipid extraction from bilayers.     

Structure optimization of the two-dimensional off-lattice hydrophobic–hydrophilic model

A two-dimensional off-lattice protein model with two species of monomers, hydrophobic and hydrophilic, was studied. Low-energy configurations in the model were optimized using the improved energy landscape paving (ELP+) method. In ELP+, the energy landscape paving (ELP) was first applied to search for the low-energy states. After the ELP led to the basins of the local energy minima, the additional degree-of-freedom of bond length was introduced, and the gradient descent method was then used to search for lower energy states near the local minima. Numerical results show that the proposed methods are quite effective for finding the ground states of proteins. A comparison between ELP+ and other methods is made.     

Potential natural exposure of Mississippi sandhill cranes to aflatoxin B1.

A survey was conducted to determine if carcinogenic mycotoxins were present in foods consumed by Mississippi sandhill cranes (Grus canadensis pulla). Samples of field corn (Zea mays) (n = 111) and chufa (Cyperus esculentus) (n = 20), obtained in 1987, 1988 and 1989 on the Mississippi Sandhill Crane National Wildlife Refuge (MSCNWR) and nearby private lands were analyzed for aflatoxin B1(AB1), ochratoxin A and sterigmatocystin using thin layer chromatography. Chufa samples were negative for all three mycotoxins. Aflatoxin B1 was found in corn at concentrations from 5 to 5,000 ppb; the other mycotoxins were not found in corn. Contaminated corn was found in 72% of all corn fields, but the proportion of contaminated fields was 57 to 100% for the 3-yr period. Contamination with AB1 was greatest in corn obtained from the ground post-harvest. Overall, 32% of corn samples from the ground had levels greater than or equal to 200 ppb with a mean of 427 ppb (range = 5 to 5,000 ppb) in contaminated fields. In 1989, mean AB1 concentration in corn on the ground was 5 to 1138 ppb for individual fields. The concentration of AB1 was less than or equal to 200 ppb in all corn samples from upright stalks. The study demonstrated that AB1 is available to sandhill cranes and at levels that may pose a serious health threat.     

The pathway by which the tetrameric protein transthyretin dissociates

The homotetrameric protein transthyretin (TTR) must undergo rate-limiting dissociation to its constituent monomers in order to enable partial denaturation that allows the process of amyloidogenesis associated with human pathology to ensue. The TTR quaternary structure contains two distinct dimer interfaces, one of which creates the two binding sites for the natural ligand thyroxine. Tetramer dissociation could proceed through three distinct pathways; scission into dimers along either of the two unique quaternary interfaces followed by dimer dissociation represents two possibilities. Alternatively, the tetramer could lose monomers sequentially. To elucidate the TTR dissociation pathway, we employed two different TTR constructs, each featuring covalent attachment of proximal subunits. We demonstrate that tethering the A and B subunits of TTR with a disulfide bond (as well as the symmetrically disposed C and D subunits) allows urea-mediated dissociation of the resulting (TTR-S-S-TTR)(2) construct, affording (TTR-S-S-TTR)(1) retaining a stable 16-stranded beta-sheet structure that is equivalent to the dimer not possessing a thyroid binding site. In contrast, linking the A and C subunits employing a peptide tether (TTR-L-TTR)(2) affords a kinetically stable quaternary structure that does not dissociate or denature in urea. Both tethered constructs and wild-type TTR exhibit analogous stability based on guanidine hydrochloride denaturation curves. The latter denaturant can denature the tetramer, unlike urea, which can only denature monomeric TTR; hence urea requires dissociation to monomers to function. Under native conditions, the (TTR-S-S-TTR)(2) construct is able to dissociate and incorporate subunits from labeled WT TTR homotetramers at a rate equivalent to that exhibited by WT TTR. In contrast, the (TTR-L-TTR)(2) construct is unable to exchange any subunits, even after 180 h. All of the data presented herein and elsewhere demonstrate that the pathway of TTR tetramer dissociation occurs by scission of the tetramer along the crystallographic C(2) axis affording AB and CD dimers that rapidly dissociate into monomers. Determination of the mechanism of dissociation provides an explanation for why small molecules that bind at the AB/CD dimer-dimer interface impose kinetic stabilization upon TTR and disease-associated variants thereof.     

Models of Replicator Proliferation Involving Differential Replicator Subunit Stability

Several models for the origin of life involve molecules that are capable of self-replication, such as self-replicating polymers composed of RNA or DNA or amino acids. Here we consider a hypothetical replicator (AB) composed of two subunits, A and B. Programs written in Python and C programming languages were used to model AB replicator abundance as a function of cycles of replication (iterations), under specified hypothetical conditions. Two non-exclusive models describe how a reduced stability for B relative to A can have an advantage for replicator activity and/or evolution by generating free A subunits. In model 1, free A subunits associate with AB replicators to create AAB replicators with greater activity. In simulations, reduced stability of B was beneficial when the replication activity of AAB was greater than two times the replication activity of AB. In model 2, the free A subunit is inactive for some number of iterations before it re-creates the B subunit. A re-creates the B subunit with an equal chance of creating B or B′, where B′ is a mutant that increases AB’ replicator activity relative to AB. In simulations, at moderate number of iterations (< 15), a shorter survival time for B is beneficial when the stability of B is greater than the inactive time of A. The results are consistent with the hypothesis that reduced stability for a replicator subunit can be advantageous under appropriate conditions.     

Romanowsky dyes and Romanowsky-Giemsa effect. 5. Structural investigations of the purple DNA-AB-EY dye complexes of Romanowsky-Giemsa staining.

A reproducible Romanowsky-Giemsa staining (RGS) can be carried out with standardized staining solutions containing the two dyes azure B (AB) and eosin Y (EY). After staining, cell nuclei have a purple coloration generated by DNA-AB-EY complexes. The microspectra of cell nuclei have a sharp and intense absorption band at 18,100 cm-1 (552 nm), the so called Romanowsky band (RB), which is due to the EY chromophore of the dye complexes. Other absorption bands can be assigned to the DNA-bound AB cations. Artificial DNA-AB-EY complexes can be prepared outside the cell by subsequent staining of DNA with AB and EY. In the first step of our staining experiments we prepared thin films of blue DNA-AB complexes on microslides with 1:1 composition: each anionic phosphodiester residue of the nucleic acid was occupied by one AB cation. Microspectrophotometric investigations of the dye preparations demonstrated that, besides monomers and dimers, mainly higher AB aggregates are bound to DNA by electrostatic and hydrophobic interactions. These DNA-AB complexes are insoluble in water. Therefore it was possible to stain the DNA-AB films with aqueous EY solutions and also to prepare insoluble DNA-AB-EY films in the second step of the staining experiments. After the reaction with EY, thin sites within the dye preparations were purple. The microspectra of the purple spots show a strong Romanowsky band at 18,100 cm-1. Using a special technique it was possible to estimate the composition of the purple dye complexes. The ratio of the two dyes was approximately EY:AB approximately 1:3. The EY anions are mainly bound by hydrophobic interaction to the AB framework of the electrical neutral DNA-AB complexes. The EY absorption is red shifted by the interaction of EY with the AB framework of DNA-AB-EY. We suppose that this red shift is caused by a dielectric polarization of the bound EY dianions. The DNA chains in the DNA-AB complexes can mechanically be aligned in a preferred direction k. Highly oriented dye complexes prepared on microslides were birefringent and dichroic. The orientation is maintained during subsequent staining with aqueous EY solutions. In this way we also prepared highly orientated purple DNA-AB-EY complexes on microslides. The light absorption of both types of dye complexes was studied by means of a microspectrophotometer equipped with a polarizer and an analyser. The sites of best orientation within the dye preparations were selected under crossed nicols according to the quality of birefringence.(ABSTRACT TRUNCATED AT 400 WORDS)     

Differential evolution for protein folding optimization based on a three-dimensional AB off-lattice model

This paper presents a differential evolution algorithm that is adapted for the protein folding optimization on a three-dimensional AB off-lattice model. The proposed algorithm is based on a self-adaptive differential evolution that improves the algorithm efficiency and reduces the number of control parameters. A mutation strategy for the fast convergence is used inside the algorithm. A temporal locality is used in order to speed up the algorithm convergence additionally and to find amino-acid conformations with the lowest free energy values. Within this mechanism a new vector is calculated when the trial vector is better than the corresponding vector from the population. This new vector is likely better than the trial vector and this accelerates convergence speed. Because of the fast convergence the algorithm has some chance to be trapped into the local optima. To mitigate this problem the algorithm includes reinitialization. The proposed algorithm was tested on amino-acid sequences that are used frequently in literature. The obtained results show that the proposed algorithm is superior to the algorithms from the literature and the obtained amino-acid sequences have significantly lower free energy values.

Open image in new window

Graphical Abstract Protein folding optimization on a three-dimensional AB off-lattice model using the differential evolution algorithm.

     

A simple peak detection and label-free quantitation algorithm for chromatography-mass spectrometry

Background

Label-free quantitation of mass spectrometric data is one of the simplest and least expensive methods for differential expression profiling of proteins and metabolites. The need for high accuracy and performance computational label-free quantitation methods is still high in the biomarker and drug discovery research field. However, recent most advanced types of LC-MS generate huge amounts of analytical data with high scan speed, high accuracy and resolution, which is often impossible to interpret manually. Moreover, there are still issues to be improved for recent label-free methods, such as how to reduce false positive/negatives of the candidate peaks, how to expand scalability and how to enhance and automate data processing. AB3D (A simple label-free quantitation algorithm for Biomarker Discovery in Diagnostics and Drug discovery using LC-MS) has addressed these issues and has the capability to perform label-free quantitation using MS1 for proteomics study.

Results

We developed an algorithm called AB3D, a label free peak detection and quantitative algorithm using MS1 spectral data. To test our algorithm, practical applications of AB3D for LC-MS data sets were evaluated using 3 datasets. Comparisons were then carried out between widely used software tools such as MZmine 2, MSight, SuperHirn, OpenMS and our algorithm AB3D, using the same LC-MS datasets. All quantitative results were confirmed manually, and we found that AB3D could properly identify and quantify known peptides with fewer false positives and false negatives compared to four other existing software tools using either the standard peptide mixture or the real complex biological samples of Bartonella quintana (strain JK31). Moreover, AB3D showed the best reliability by comparing the variability between two technical replicates using a complex peptide mixture of HeLa and BSA samples. For performance, the AB3D algorithm is about 1.2 - 15 times faster than the four other existing software tools.

Conclusions

AB3D is a simple and fast algorithm for label-free quantitation using MS1 mass spectrometry data for large scale LC-MS data analysis with higher true positive and reasonable false positive rates. Furthermore, AB3D demonstrated the best reproducibility and is about 1.2- 15 times faster than those of existing 4 software tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0376-0) contains supplementary material, which is available to authorized users.     

A Fast Algorithm to Estimate the Deepest Points of Lakes for Regional Lake Registration

When conducting image registration in the U.S. state of Alaska, it is very difficult to locate satisfactory ground control points because ice, snow, and lakes cover much of the ground. However, GCPs can be located by seeking stable points from the extracted lake data. This paper defines a process to estimate the deepest points of lakes as the most stable ground control points for registration. We estimate the deepest point of a lake by computing the center point of the largest inner circle (LIC) of the polygon representing the lake. An LIC-seeking method based on Voronoi diagrams is proposed, and an algorithm based on medial axis simplification (MAS) is introduced. The proposed design also incorporates parallel data computing. A key issue of selecting a policy for partitioning vector data is carefully studied, the selected policy that equalize the algorithm complexity is proved the most optimized policy for vector parallel processing. Using several experimental applications, we conclude that the presented approach accurately estimates the deepest points in Alaskan lakes; furthermore, we gain perfect efficiency using MAS and a policy of algorithm complexity equalization.     

A Two-Stage Algorithm for Origin-Destination Matrices Estimation Considering Dynamic Dispersion Parameter for Route Choice

This paper proposes a two-stage algorithm to simultaneously estimate origin-destination (OD) matrix, link choice proportion, and dispersion parameter using partial traffic counts in a congested network. A non-linear optimization model is developed which incorporates a dynamic dispersion parameter, followed by a two-stage algorithm in which Generalized Least Squares (GLS) estimation and a Stochastic User Equilibrium (SUE) assignment model are iteratively applied until the convergence is reached. To evaluate the performance of the algorithm, the proposed approach is implemented in a hypothetical network using input data with high error, and tested under a range of variation coefficients. The root mean squared error (RMSE) of the estimated OD demand and link flows are used to evaluate the model estimation results. The results indicate that the estimated dispersion parameter theta is insensitive to the choice of variation coefficients. The proposed approach is shown to outperform two established OD estimation methods and produce parameter estimates that are close to the ground truth. In addition, the proposed approach is applied to an empirical network in Seattle, WA to validate the robustness and practicality of this methodology. In summary, this study proposes and evaluates an innovative computational approach to accurately estimate OD matrices using link-level traffic flow data, and provides useful insight for optimal parameter selection in modeling travelers’ route choice behavior.     

Phenotypic mixing between different hepadnavirus nucleocapsid proteins reveals C protein dimerization to be cis preferential.

Hepadnaviruses encode a single core (C) protein which assembles into a nucleocapsid containing the polymerase (P) protein and pregenomic RNA during viral replication in hepatocytes. We examined the ability of heterologous hepadnavirus C proteins to cross-oligomerize. Using a two-hybrid assay in HepG2 cells, we observed cross-oligomerization among the core proteins from hepatitis B virus (HBV), woodchuck hepatitis virus, and ground squirrel hepatitis virus. When expressed in Xenopus oocytes, in which hepadnavirus C proteins form capsids, the C polypeptides from woodchuck hepatitis virus and ground squirrel hepatitis virus, but not duck hepatitis B virus, can efficiently coassemble with an epitope-tagged HBV core polypeptide to form mixed capsids. However, when two different core mRNAs are coexpressed in oocytes the core monomers show a strong preference for forming homodimers rather than heterodimers. This holds true even for coexpression of two HBV C proteins differing only by an epitope tag, suggesting that core monomers are not free to diffuse and associate with other monomers. Thus, mixed capsids result from aggregation of different species of homodimers.     

Structures of the two 3D domain-swapped RNase A trimers

When concentrated in mildly acidic solutions, bovine pancreatic ribonuclease (RNase A) forms long-lived oligomers including two types of dimer, two types of trimer, and higher oligomers. In previous crystallographic work, we found that the major dimeric component forms by a swapping of the C-terminal beta-strands between the monomers, and that the minor dimeric component forms by swapping the N-terminal alpha-helices of the monomers. On the basis of these structures, we proposed that a linear RNase A trimer can form from a central molecule that simultaneously swaps its N-terminal helix with a second RNase A molecule and its C-terminal strand with a third molecule. Studies by dissociation are consistent with this model for the major trimeric component: the major trimer dissociates into both the major and the minor dimers, as well as monomers. In contrast, the minor trimer component dissociates into the monomer and the major dimer. This suggests that the minor trimer is cyclic, formed from three monomers that swap their C-terminal beta-strands into identical molecules. These conclusions are supported by cross-linking of lysyl residues, showing that the major trimer swaps its N-terminal helix, and the minor trimer does not. We verified by X-ray crystallography the proposed cyclic structure for the minor trimer, with swapping of the C-terminal beta-strands. This study thus expands the variety of domain-swapped oligomers by revealing the first example of a protein that can form both a linear and a cyclic domain-swapped oligomer. These structures permit interpretation of the enzymatic activities of the RNase A oligomers on double-stranded RNA.     

Lattice simulations of protein crystal formation

A new algorithm is presented for the lattice simulation of protein crystal growth. The algorithm allows the calculation of the size distribution of microcrystals in the volume and timescale of experiments and within the framework of the previously-published microscopic model [A.M. Kierzek, W.M. Wolf, P. Zielenkiewicz, Biophys. J. 73 (1997) 571-580]. Simulations for the tetragonal lysozyme crystal show that there are two critical sizes in the development of ordered phase. The first one corresponds to the size of the smallest stable complex which, in the case of the tetragonal lysozyme crystal, is the particular tetramer. In a volume of 5 mul the tetramer appears in the millisecond timescale. The second critical radius of approximately 100 monomers is only reached by a few of all the smallest stable complexes formed in the solution. The model predicts that out of 10(7) tetramers which appear in solution, only eight reach the size of 100 monomers within 8 h. After exceeding the second critical radius the microcrystals grow to the size of 10(4) monomers in the minute timescale and are thus assumed to quickly lead to macroscopic crystals. The predicted number of crystals formed during 8 h of nucleation is in qualitative agreement with arrested nucleation experiments.     

Frequency characteristics of airway and tissue impedances in respiratory diseases

We measured the frequency characteristics (at 10-40 Hz) of airway (Za) and tissue (Zt) impedances in cases of chronic obstructive pulmonary disease [asthmatic bronchitis (AB), chronic pulmonary emphysema (CPE)] and interstitial pneumonitis (IP) by use of an improved random noise oscillation and body box method. The results were then compared with those obtained for normal subjects. The real part of Za was markedly elevated in patients with AB but only slightly elevated in those with CPE. To interpret these data we used an electromechanical analogue including serial inhomogeneity with shunt impedance. From this model we concluded that AB causes both the central and peripheral airway resistances to increase, while CPE brings about a rise mainly in peripheral resistance. In IP patients, only the imaginary part of Zt decreased, which might reflect the decrease in both lung and chest wall compliance. In CPE patients, but not in AB patients, the real part of Zt fell. These data were consistent with the assumption that the decrease in mass per unit volume of lung tissue and hyperinflation of the chest wall in CPE patients might lower the tissue resistances.     

Membrane topography of the hydrophobic anchor sequence of poliovirus 3A and 3AB proteins and the functional effect of 3A/3AB membrane association upon RNA replication

Replication of poliovirus RNA takes place on the cytoplasmic surface of membranous vesicles that form after infection of the host cell. It is generally accepted that RNA polymerase 3D(pol) interacts with membranes in a complex with viral protein 3AB, which binds to membranes by means of a hydrophobic anchor sequence that is located near the C-terminus of the 3A domain. In this study, we used fluorescence and fluorescence quenching methods to define the topography of the anchor sequence in the context of 3A and 3AB proteins inserted in model membranes. Mutants with a single tryptophan near the center of the anchor sequence but lacking Trp elsewhere in 3A/3AB were constructed which, after the emergence of suppressor mutations, replicated well in HeLa cells. When a peptide containing the mutant anchor sequence was incorporated in model membrane vesicles, measurements of Trp depth within the lipid bilayer indicated formation of a transmembrane topography. However, rather than the 22-residue length predicted from hydrophobicity considerations, the transmembrane segment had an effective length of 16 residues, such that Gln64 likely formed the N-terminal boundary. Analogous experiments using full-length proteins bound to preformed model membrane vesicles showed that the anchor sequence formed a mixture of transmembrane and nontransmembrane topographies in the 3A protein but adopted only the nontransmembrane configuration in the context of 3AB protein. Studies of the function of 3A/3AB inserted into model membrane vesicles showed that membrane-bound 3AB is highly efficient in stimulating the activity of 3D(pol) in vitro while membrane-bound 3A totally lacks this activity. Moreover, in vitro uridylylation reactions showed that membrane-bound 3AB is not a substrate for 3D(pol), but free VPg released by cleavage of 3AB with proteinase 3CD(pro) could be uridylylated.