 Software
 Open Access
 Published:
CPSPtools – Exact and complete algorithms for highthroughput 3D lattice protein studies
BMC Bioinformatics volume 9, Article number: 230 (2008)
Abstract
Background
The principles of protein folding and evolution pose problems of very high inherent complexity. Often these problems are tackled using simplified protein models, e.g. lattice proteins. The CPSPtools package provides programs to solve exactly and completely the problems typical of studies using 3D lattice protein models. Among the tasks addressed are the prediction of (all) globally optimal and/or suboptimal structures as well as sequence design and neutral network exploration.
Results
In contrast to stochastic approaches, which are not capable of answering many fundamental questions, our methods are based on fast, nonheuristic techniques. The resulting tools are designed for highthroughput studies of 3Dlattice proteins utilising the HydrophobicPolar (HP) model. The source bundle is freely available [1].
Conclusion
The CPSPtools package is the first set of exact and complete methods for extensive, highthroughput studies of nonrestricted 3Dlattice protein models. In particular, our package deals with cubic and face centered cubic (FCC) lattices.
Background
The organisation of biomolecules, in particular proteins, in the sequence and structure space has recently been attracting increased attention. Particularly questions concerning finding the native structure or investigating the kinetics and evolution of proteins have been widely studied. These problems are often tackled using simplified models such as the HydrophobicPolar (HP) model (e.g. Jacob et al. [2]). Though abstract, these models are computationally feasible and do allow for deeper insights into fundamental and general principles [2–4].
Several recurring tasks can be identified in such studies using simplified models. Namely, predicting the native structure, classifying whether a sequence is proteinlike, calculating its degeneracy and stability, or the design of sequences that optimally fold to a given structure. The problems associated with these tasks are computationally very hard (NPcomplete) [5–7]. Nevertheless, these tasks demand for exact and complete (i.e. nonheuristic) methods. It is important to note that stochastic methods cannot be used for proving optimality and in particular proving that a sequence has a unique lowest energy (proteinlike) fold [8].
Consequently, with the exception of Yue and Dill [9], all studies requiring complete and exact answers to optimal structure prediction were based on exhaustive enumeration. These studies were, hence, confined to small sequence lengths. In other approaches, structures are artificially restricted to be maximally compact (e.g. filling a 3 × 3 × 3 cube) [10]. This allows for complete enumeration but artificially biases the energy function towards overall hydrophobicity.
Furthermore, many studies are confined to extremely simplified models on the 2Dsquare or 3Ddiamondlattice [2, 11]. The coordination number, a measurement of lattice complexity, is four in both cases. The use of lattices with such a low complexity may lead to oversimplified models that are not able to reproduce real world properties. Park and Levitt [12] have shown that lattices with higher coordination number provide a much better fit to real protein structures. A further hint toward the simplicity of the 2Dlattice is the low computational complexity of inverse folding when compared to the 3Dcubic lattice [7]. The Constraintbased Protein Structure Prediction (CPSP) approach by Backofen and Will [13] provides a way to overcome the aforementioned obstacles. The method is tailored to the HP model introduced by Lau and Dill [14]. This model is widely used in the literature [15, 16]. CPSP supports complex 3D lattices (currently cubic and face centered cubic) without artificial restrictions (e.g. to be maximally compact). The approach predicts all globally optimal structures together with a proof of optimality. No naive, exhaustive enumeration of all structures is performed and it is as fast as stochastic methods that cannot prove optimality. Backofen and Will [13] showed that the CPSPapproach could fold even sequences of length 200 to optimality within seconds. In contrast, exhaustive structure enumeration as e.g. done by Blackburne and Hirst [17] is restricted to short sequence lengths. For instance, on a 3Dcubic lattice it is only viable to enumerate up to about length 20. In fact, the exact number of structures is only known up to length 23 where there are already more than 5 × 10^{15} [18]. CPSP uses constraint programming that is commonly applied to hard (NPcomplete) problems and, thus, avoids the complete expansion of the whole search space. Hence, constraintprogramming techniques are a powerful tool to handle the high complexity that typifies problems related to protein structure. Constraintprogramming techniques have successfully been applied to structure prediction with given secondary structure information [19], analysis of NMR data [20], and modeling of protein complexes [21].
Currently, we are not aware of any other complete approach that ensures optimality of the predicted structures in different lattices. There is an alternative to CPSP for the 3Dcubic lattice, the constraintbased hydrophobic core construction method by Yue and Dill [9]. This allows the prediction of optimal structures and proves their optimality. However, using the CPSPapproach, Backofen and Will showed that the method developed by Yue and Dill is not always complete in enumerating all optimal structures [13].
Complex Lattices
As mentioned before, complete structure enumeration is only applicable to simple, low coordination number lattices. In contrast, the CPSPapproach is built for the more complex 3Dcubic and 3Dfacecenteredcubic (FCC) lattices with higher coordination numbers of 6 and 12, respectively. A main feature of the CPSPtools is their applicability to the unrestricted FCC lattice. The FCC lattice lacks one of the main problems of the 3Dcubic lattice, namely that only sequence positions with different parities form contacts; the parity problem [22]. Modeling protein structures on a FCC lattice, Park and Levitt [12] demonstrated that a good approximation of real protein structures is possible. They achieved a coordinate root mean square deviation of 1.78 Å, whereas a deviation of 2.84 Å was obtained in the 3Dcubic lattice. Recently, Bagci et al. [23] have shown that the neighborhood of amino acids in proteins closely resembles a distorted FCC lattice, and that the FCC is best suited for modeling proteins. The CPSPapproach is the first exact method that allows the prediction of provable optimal structures in the FCC lattice. An example is given in Figure 1.
Implementation
CPSPtools provides a set of programs that enable typical, modern research tasks to be calculated efficiently and accurately. Here we list the programs each with a typical example application. HPSTRUCT predicts (all) optimal and suboptimal structures as required for investigating properties of low energy conformations, as e.g. studied by Jacob and Unger [16]. The statistical analysis of proteinlike sequences, see Blackburne and Hirst [11], requires a degeneracybased classification of sequences. This is possible with HPDEG. For the exploration of protein evolution, similar to Wroe and Chan [24], one needs to investigate the sequencestructure space. We provide HPDESIGN for sequence design and HPNNET for neutral network computation.
All methods can be applied to HPsequences in the cubic and the more complex face centered cubic lattice model. Before giving a detailed description of the tools, we first introduce the idea of Hcores, central to these methods.
Hcore database
In the HP lattice models, two monomers form a contact if they occupy neighboring positions in the lattice. The energy of a structure is defined by the number of contacts between Hmonomers, i.e. HHcontacts. Thus, an optimal (minimum energy) conformation maximizes the number of HHcontacts. An important observation is that optimal structures show an almost optimal (maximally compact) packing of the Hmonomers. Such dispersions of Hmonomers without any chain connectivity are called Hcores. The compactness of the Hcores is a basic feature that can be used for structure prediction and sequence design. Note that optimal Hcores are independent of a particular sequence and depend only on the number of Hmonomers. Hence, compact and nearly compact Hcores can be precalculated and stored in a database. HPSTRUCT and HPDESIGN use this database as a starting point for their calculations (details later). Thereby, redundant computation is avoided, which significantly speeds up the CPSPapproach and related applications.
The enumeration of all optimal Hcores in complex lattice models such as FCC is a computationally hard problem by itself and was solved by Backofen and Will using constraintprogramming techniques [25]. Firstly an upper bound on the number of possible contacts for a given number of monomers is calculated via dynamic programming. Subsequently, this information is used to enumerate all compact optimal and almost optimal (suboptimal) Hcores for a given number of Hmonomers using constraintprogramming. Some statistics on the number of Hcores in the 3Dcubic lattice are given in Fig. 4. It shows that the number of Hcores grows exponentially in Hcore size but still much slower than the number of structures for a corresponding sequence length.
HPstruct
Motivation
HPSTRUCT implements the CPSP approach, as introduced by Backofen and Will [13], to predict provably optimal structures of 3D lattice proteins in the HPmodel. For a given HPsequence S and a given lattice type (cubic or face centered cubic), (all) optimal structures are calculated. The CPSP approach computes the global minimal energy for S.
Methods
The CPSPapproach is based on the Hcore database as described before. For a concrete sequence S the approach systematically examines the list of Hcores compatible with S in decreasing maximal contact number. For each core, it attempts to thread the sequence through the core. Threading means to find a placement of the monomers of S in a selfavoiding walk such that all Hmonomers are elements of the given Hcore and all Pmonomers are outside of the core. Since the Hcores are considered in the order of decreasing contacts, the first successful threading results in a structure with global minimal energy. Note that at this point the algorithm has proven that there is no structure of S that forms more HHcontacts.
Technically, the threading of a sequence through a core is performed by a constraint program. For this purpose, we formulate the threading problem as a constraint satisfaction problem (CSP) [26]. It constrains the Hmonomers of the sequence to the positions in the Hcore. Further, it enforces successive monomers along the sequence to be neighbored in the lattice and prohibits the multiple use of a single position. The constraintprogramming machinery allows for the enumeration of all valid placements according to the given constraints. In this way, all (sub)optimal structures for a given sequence can be calculated. For a more detailed description of the CSP definition and the mechanisms for solving it see [13].
Advanced Features
All resulting structures of HPSTRUCT are returned in absolute move string representation. This compactly encodes the lattice position vectors between successive monomers in the structure and reduces the space consumption for huge data sets.
To handle the common case of highly degenerated sequences (with many optima), HPSTRUCT offers the possibility to limit the number of predicted structures or to generate only a representing subset. Such a subset only contains structures that are separated by at least (a user defined) distance k. The distance measure is the hamming distance on the absolute move strings.
HPdeg
Motivation
The degeneracy of an HPsequence S is the number of optimal structures S can adopt. It can be calculated using HPDEG and is the base to determine the stability of structures [27]. HPDEG specializes HPSTRUCT and completely counts all optimal structures.
An important application of HPDEG is the classification of sequences as proteinlike or not. A sequence is proteinlike if it can adopt only one optimal structure (degeneracy 1), a definition applied by Li et al. [10] and Huard et al. [4] among others.
Methods
HPDEG is directly based on the CPSPapproach to compute the degeneracy. Here, all solutions for all arbitrary Hcores/CSPs are calculated. In addition, a significant acceleration of the process can be achieved by the search decomposition methods we introduced in [28]. This is done by identifying subchains of the sequence that can be placed independently from each other. Their placements are calculated separately and the resulting numbers are multiplied to the overall structure number of the whole chain. This decomposition strategy results in a speedup of 3times and higher on average.
HPdesign
Motivation
HPDESIGN solves the inverse folding problem, i.e. the design of sequences that form a given structure X as their unique optimum. It allows deeper investigations of sequencestructure relations and a better understanding of general properties of protein folding [29].
The inverse folding problem (IFP) in 3D lattices has been shown by Berman et al. [7] to be NPcomplete, i.e. it is, as the protein folding problem, a hard computational problem. In contrast, as the same authors show, the IFP in the simple 2D lattice is solvable in polynomial time. This indicates once more the higher complexity of threedimensional lattice models. To our knowledge, HPDESIGN is the only method applicable to a 3Dmodel that calculates the desired sequence properties without exhaustive sequence space enumeration.
Methods
The approach is based on the CPSP Hcore database in order to get a set of good candidate sequences C. First, using Hcores ordered by decreasing size and optimality, a matching of the core and the structure is done. For each match a candidate sequence is derived and added to C. Afterwards, each c ∈ C is evaluated concerning degeneracy and checked if X is its optimal structure.
The candidate set C, produced by the filtering step using the Hcores, consists of sequences that can adopt X with an optimal or slightly suboptimal Hcore. Therefore, their probability to form X as their unique optimum is very high and the size of C very small compared to the whole sequence space. The latter is of high importance for the performance of the method.
Advanced Features
Often sequences with a special ratio of H/P occurrences or with only limited degeneracy are of interest. Both can be specified using HPDESIGN.
Furthermore, the number of evaluated Hcores is selectable to allow a balancing between runtime and completeness. This is done by adjusting their allowed level of optimality used in the filtering step.
HPnnet
Motivation
The organisation of sequence space in neutral networks provides insights into evolutionary principles [15, 30]. Such networks can be expanded using HPNNET. A neutral network for a given structure X is an undirected binary graph, where each node represents a sequence that forms X as its unique optimal structure. Edges connect evolutionary related sequences, i.e. sequences that differ only in one sequence position, a point mutation. HPNNET expands a neutral network starting from an initial sequence (or a set of sequences) S that folds into the structure X.
Methods
The method follows the generateandtest paradigm. Recursively, all neighboring sequences of S are tested if they adopt X as their unique optimum. If so, they are added to the network and their neighbors are checked. Therefore, HPNNET is capable of detecting and expanding connected neutral networks of different structures.
Advanced Features
Running HPNNET with S as the only start sequence results in the connected component of the network S belongs to. However, Blackburne and Hirst [17] have shown by exhaustive enumeration in restricted models that neutral networks may consist of several connected components. To find and study them in complex threedimensional lattices a combination of HPDESIGN and HPNNET can be used. The independently designed sequences resulting from HPDESIGN have a high chance to belong to different components. HPNNET supports as input such a set of sequences and expands all corresponding connected components. An example is later shown in the results section.
Utility tools
In addition to those described above, CPSPtools provides a set of utility programs helpful for lattice protein studies. For instance using HPCONVERT, it is possible to convert between absolute move strings, the 3Dposition data in XYZ, Protein Data Bank (PDB) and Chemical Markup Language (CML) format. A move string normalization, as well as a conversion into an orientation independent relative move string, is available for a symmetry independent structure comparison.
HPVIEW interactively visualizes structures in 2Dsquare, 3Dcubic, and 3DFCC lattices using the Jmol interface [31].
Installation and Usage
The package supplies standard installation procedures for Linux based on common tools (GNU automake) and can be compiled and installed easily on current 32 and 64bit Linux systems (including Cygwin for Microsoft Windows™). The programs are written in C++ for highest performance and provide a slim textbased user interface for efficient pipelining as required for highthroughput experiments. A web front end is under development.
All constraint programming based algorithms utilize the open source Gecode system [32].
The validity of the algorithms has been tested and confirmed on a large set of benchmark problems. The functionality of Hcore database access, structure prediction, and degeneracy computation are collected in the C++ CPSPlibrary. A complete API is included which allows the embedding, extension, and use of the CPSP approach in new programs.
To reduce package size, only a small fraction of the Hcore database is included in the source package. This already enables the use of CPSPtools for short sequences. The complete database is available on request.
Results and Discussion
For illustration, we provide some scenarios that exemplify the use of CPSPtools in extending known or enabling new studies. All examples are performed in the unrestricted 3Dcubic lattice with HPsequences of length 27. Note that for this length there are already more than 10^{19} possible structures, which makes an exhaustive enumeration inapplicable. Table 1 outlines the performance of programs from CPSPtools. Table 2 shows the sequences used for Table 1, their optimal energy (E), and degeneracy (deg). All tasks were performed on an Intel P4 3 GHz (using CPSP2.0.0).

(1)
Studies of sequence or structure features of proteins as done by Huard et al. [4] require a classification of sequences as proteinlike. One way is to classify them by the number of optimal structures, i.e. their degeneracy. The fast calculation of this sequence property by HPDEG allows production of sufficiently large benchmark sets for detailed studies. To illustrate this, we run HPDEG for a random HPsequence S_{0} revealing an enormous degeneracy, which is a frequent finding in the HPmodel. As a starting point for the following scenarios, we evaluate the degeneracy of S_{1}, a sequence with a single optimal structure. The very short runtimes for both checks are given in Table 1.

(2)
Calculating the globally optimal structure for a given sequence is the main task in many studies, e.g. see Jacob and Unger [16]. Furthermore, in stochastic folding simulation approaches knowing the minimal possible energy is favorable. Both can be calculated extremely rapidly using HPSTRUCT. Again, We demonstrate this with sequences S_{0} and S_{1}. This results in an energy of 13 and 22 and the optimal structures X_{0} and X_{1}, respectively. Both structures are visualized in Figure 2.

(3)
To study protein evolution on the sequence level, neutral networks are widely utilized [17]. Using HPNNET we can span the connected component of the neutral network for a given sequence with a unique optimal structure. Applied to S_{1} with X_{1} we find four sequences S_{2} .. S_{4} sharing X_{1} as their unique optimal structure. Note, this can be done without exhaustive sequence enumeration for a given structure.

(4)
The detailed study of neutral networks by Blackburne and Hirst [17] has shown that neutral networks may decompose into connected components. Their results are based on full enumeration of sequences and structures in the diamond lattice. This approach does not extend to complex lattice models due to the enormous size of the structure space as discussed above.
HPDESIGN can overcome that problem by directly designing sequences of the neutral network. Recall that the neutral network contains only sequences with the same unique optimal structure. The described design approach allows one to generate sequences of independent components in the neutral network without exhaustive enumeration. Afterwards, the full components can be expanded via HPNNET.
We apply this approach to the neutral network of the structure X_{1}. HPDESIGN calculates 12 members of the network (S_{1} .. S_{12}), including the four sequences S_{1} .. S_{4} known from scenario (3). Expanding the network N from these sequences via HPNNET reveals two further sequences S_{13}, S_{14} and two independent connected components as shown in Figure 3.
Preliminary studies performed with CPSPtools indicate that neutral networks as large as N with several large independent components are rare in the unrestricted 3Dcubic model.
Conclusion
For complex 3D models, mainly heuristic and/or stochastic approaches to search for optimal structures of a given sequence are available [8, 33]. However, these methods are (a) incomplete and (b) cannot ensure the global optimality of the predicted structures. In consequence, the investigation of problems requiring this information was only possible using exhaustive enumeration, which is not possible for longer sequence lengths.
The CPSP approach is as fast as common stochastic methods but ensures that all predicted structures are globally optimal, and that none are missing. This is done without exhaustive structure space exploration applying constraintprogramming techniques. Therefore, it is well suited to many studies in complex 3D models; especially for finding proteinlike sequences, the investigation of neutral networks or sequence design. Further applications range from the generation of candidate sets to the validation of results of folding simulations and stochastic optimization methods.
The CPSPtools package combines several applications in the field of bioinformatics concerning 3D lattice proteins. It allows advanced investigation of problems related to protein structure prediction, sequence evolution, inverse folding, and energy landscapes.
Availability and requirements
Project name: CPSPtools
Project home page: http://www.bioinf.unifreiburg.de/sw/cpsp/
Operating system(s): all Linux based systems (including Cygwin for MS Windows™)
Programming language: C++
Other requirements: Gecode and BIU library (a source bundle is provided)
License: BSDstyle license
Any restrictions to use by nonacademics: none
References
 1.
CPSPtools[http://www.bioinf.unifreiburg.de/sw/cpsp/]
 2.
Jacob E, Horovitz A, Unger R: Different mechanistic requirements for prokaryotic and eukaryotic chaperonins: a lattice study. Bioinformatics 2007, 23: 240–248. 10.1093/bioinformatics/btm180
 3.
Wolfinger MT, Will S, Hofacker IL, Backofen R, Stadler PF: Exploring the lower part of discrete polymer model energy landscapes. Europhysics Lett 2006, 74: 725–732. 10.1209/epl/i2005105770
 4.
Huard FP, Deane CM, Woo GR: Modelling sequential protein folding under kinetic control. Bioinformatics 2006, 22: 202–210. 10.1093/bioinformatics/btl248
 5.
Unger R, Moult J: Finding the lowest free energy conformation of a protein is an NPhard problem: proof and implications. Bull Math Biol 1993, 55: 1183–1198.
 6.
Berger B, Leighton T: Protein folding in the hydrophobichydrophilic (HP) model is NPcomplete. J Comp Biol 1998, 5: 27–40.
 7.
Berman P, DasGupta B, Mubayi D, Sloan R, Turán G, Zhang Y: The protein sequence design problem in canonical model on 2D and 3D lattices. In Combinatorial Pattern Matching. Volume 3109. Springer; 2004:244–253.
 8.
Hsu HP, Mehra V, Nadler W, Grassberger P: Growthbased optimization algorithm for lattice heteropolymers. Phys Rev E 2003, 68: 021113. 10.1103/PhysRevE.68.021113
 9.
Yue K, Dill KA: Forces of tertiary structural organization in globular proteins. Proc Natl Acad Sci 1995, 92: 146–150. 10.1073/pnas.92.1.146
 10.
Li H, Tang C, Wingreen NS: Designability of protein structures: a latticemodel study using the MiyazawaJernigan matrix. Proteins 2002, 49: 403–412. 10.1002/prot.10239
 11.
Blackburne BP, Hirst JD: Population dynamics simulations of functional model proteins. J Chem Phys 2005, 123: 154907–9. 10.1063/1.2056545
 12.
Park BH, Levitt M: The complexity and accuracy of discrete state models of protein structure. J Mol Biol 1995, 249: 493–507. 10.1006/jmbi.1995.0311
 13.
Backofen R, Will S: A constraintbased approach to fast and exact structure prediction in threedimensional protein models. Constraints 2006, 11: 5–30. 10.1007/s1060100668488
 14.
Lau KF, Dill KA: A lattice statistical mechanics model of the conformational and sequence spaces of proteins. Macromolecules 1989, 22: 3986–3997. 10.1021/ma00200a030
 15.
Cui Y, Wong WH, BornbergBauer E, Chan HS: Recombinatoric exploration of novel folded structures: a heteropolymerbased model of protein evolutionary landscapes. Proc Natl Acad Sci 2002, 99: 809–814. 10.1073/pnas.022240299
 16.
Jacob E, Unger R: A tale of two tails: why are terminal residues of proteins exposed? Bioinformatics 2007, 23: 225–230. 10.1093/bioinformatics/btl318
 17.
Blackburne BP, Hirst JD: Threedimensional functional model proteins: structure function and evolution. J Chem Phys 2003, 119: 3453–3460. 10.1063/1.1590310
 18.
Sloane NJA: Number of nstep selfavoiding walks on cubic lattice. OnLine Encyclopedia of Integer Sequences 2007. [http://www.research.att.com/~njas/sequences/A001412]
 19.
Dal Palu A, Dovier A, Fogolari F: Constraint Logic Programming approach to protein structure prediction. BMC Bioinformatics 2004, 5: 186. 10.1186/147121055186
 20.
Krippahl L, Barahona P: PSICO: Solving protein structures with constraint programming and optimization. Constraints 2002, 7: 317–331. 10.1023/A:1020577603762
 21.
Krippahl L, Moura JJ, Palma PN: Modeling protein complexes with BiGGER. Proteins 2003, 52: 19–23. 10.1002/prot.10387
 22.
Decatur SE: Protein folding in the generalized hydrophobicpolar model on the triangular lattice. 1996.
 23.
Bagci Z, Jernigan RL, Bahar I: Residue coordination in proteins conforms to the closest packing of spheres. Polymer 2002, 43: 451–459. 10.1016/S00323861(01)00427X
 24.
Wroe R, Chan HS, BornbergBauer E: A structural model of latent evolutionary potentials underlying neutral networks in proteins. HFSP J 2007, 1: 79–87. 10.2976/1.2739116
 25.
Backofen R, Will S: Optimally Compact Finite Sphere Packings – Hydrophobic Cores in the FCC. In Proc of the 12th Annual Symposium on Combinatorial Pattern Matching. Volume 2089. Springer; 2001:257–272.
 26.
Marriott K, Stuckey PJ: Programming with Constraints: an Introduction. The MIT Press; 1998.
 27.
Shortle D, Chan HS, Dill KA: Modeling the effects of mutations on the denatured states of proteins. Prot Sci 1992, 1: 201–215.
 28.
Will S, Mann M: Counting protein structures by DFS with dynamic decomposition. Proc of Workshop on Constraint Based Methods for Bioinformatics 2006, 83–90.
 29.
Gupta A, Manuch J, Stacho L: Structureapproximating inverse protein folding problem in the 2D HP model. J Comp Biol 2005, 12: 1328–1345. 10.1089/cmb.2005.12.1328
 30.
Schuster P, Stadler PF: Networks in molecular evolution. Complexity 2002, 8: 34–42. 10.1002/cplx.10052
 31.
Jmol: an opensource Java viewer for chemical structures in 3D[http://jmol.sourceforge.net/]
 32.
Gecode – generic constraint development environment[http://www.gecode.org]
 33.
Shmygelska A, Hoos HH: An ant colony optimisation algorithm for the 2D and 3D hydrophobic polar protein folding problem. BMC Bioinformatics 2005, 6: 30. 10.1186/14712105630
Acknowledgements
Martin Mann is supported by the EU project EMBIO (EC contract number 012835). Sebastian Will is partially supported by the EU Network of Excellence REWERSE (project reference number 506779).
Further, thanks to the reviewers of an earlier version of the manuscript for their helpful comments and Rhodri Saunders for proofreading.
Author information
Affiliations
Corresponding author
Additional information
Authors' contributions
Implementation and software design was done by MM and SW. The CPSP method was developed by RB and SW and extended by SW and MM. The CPSP derived algorithms are designed by all authors. All authors have approved and contributed to the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Mann, M., Will, S. & Backofen, R. CPSPtools – Exact and complete algorithms for highthroughput 3D lattice protein studies. BMC Bioinformatics 9, 230 (2008). https://doi.org/10.1186/147121059230
Received:
Accepted:
Published:
Keywords
 Optimal Structure
 Constraint Satisfaction Problem
 Neutral Network
 Face Centered Cubic
 Exhaustive Enumeration