Merge branch 'prerelease-3.2.1' into czarek
[unres.git] / doc / CLUSTER.TXT
diff --git a/doc/CLUSTER.TXT b/doc/CLUSTER.TXT
deleted file mode 100644 (file)
index f952c76..0000000
+++ /dev/null
@@ -1,671 +0,0 @@
-                             CLUSTER 
-          Cluster analysis of UNRES simulation results 
-          ---------------------------------------------
-
-TABLE OF CONTENTS
------------------
-
-1. License terms
-
-2. References
-
-3. Functions of the program
-
-4. Installation
-
-5. Running the program
-
-6. Input and output files 
-   6.1. Summary of files
-   6.2. The main input file
-        6.2.1. Title
-        6.2.2. General data
-        6.2.3. Energy-term weights and parameter files
-        6.2.4 Molecule data
-              6.2.4.1. Sequence information
-              6.2.4.2. Dihedral angle restraint information
-              6.2.4.3. Disulfide-bridge data
-        6.2.5. Reference structure
-   6.3. Main output file (out)
-   6.4. Output coordinate files
-        6.4.1. The internal coordinate (int) files
-        6.4.2. The Cartesian coordinate (x) files
-        6.4.3. The PDB files
-               6.4.3.1. CLUST-UNRES runs
-               6.4.3.2. CLUST-WHAM runs
-                        6.4.3.2.1. Conformation family files
-                        6.4.3.2.2. Average-structure file
-   6.5. The conformation-distance file
-   6.6. The clustering-tree PicTeX file
-
-7. Support
-
-1. LICENSE TERMS
-----------------
-
-* This software is provided free of charge to academic users, subject to the
-  condition that no part of it be sold or used otherwise for commercial
-  purposes, including, but not limited to its incorporation into commercial
-  software packages, without written consent from the authors. For permission
-  contact Prof. H. A. Scheraga, Cornell University.
-
-* This software package is provided on an "as is" basis. We in no way warrant
-  either this software or results it may produce.
-
-* Reports or publications using this software package must contain an
-  acknowledgment to the authors and the NIH Resource in the form commonly
-used
-  in academic research.
-
-2. REFERENCES
--------------
-
-The program incorporates the hierarchical-clustering subroutine, hc.f written
-by G. Murtagh (refs 1 and 2). The subroutine contains seven methods of 
-hierarchical clustering.
-
-[1]  F. Murtagh. Multidimensional clustering algorithms; Physica-Verlag: 
-     Vienna, Austria, 1985.
-[2]  F. Murtagh, A. Heck. MultiVariate data analysis; Kluwer Academic:
-     Dordrecht, Holland, 1987.
-[3]  A. Liwo, M. Khalili, C. Czaplewski, S. Kalinowski, S. Oldziej, K. Wachucik,
-     H.A. Scheraga.
-     Modification and optimization of the united-residue (UNRES) potential
-     energy function for canonical simulations. I. Temperature dependence of the
-     effective energy function and tests of the optimization method with single
-     training proteins. J. Phys. Chem. B, 2007, 111, 260-285.
-[4]  S. Oldziej, A. Liwo, C. Czaplewski, J. Pillardy, H.A. Scheraga.
-     Optimization of the UNRES force field by hierarchical design of the
-     potential-energy landscape. 2. Off-lattice tests of the method with single
-     proteins.  J. Phys. Chem. B., 2004, 108, 16934-16949.
-
-3. FUNCTIONS OF THE PROGRAM
----------------------------
-
-The program runs cluster analysis of UNRES simulation results. There are two
-versions of the program depending on the origin of input conformation: 
-
-1) CLUST-UNRES: performs cluster analysis of conformations that are obtained 
-   directly from UNRES runs (CSA, MCM, MD, (M)REMD, multiple-conformation 
-   energy minimization). The source code and other important files are
-   deposited in CLUST-UNRES subdirectory
-
-   The source code of this version is deposited in clust-unres/src
-
-2) CLUST-WHAM: performs cluster analysis of conformations obtained in UNRES
-   MREMD simulations and then processed with WHAM (weighted histogram analysis 
-   method). This enables the user to obtain clusters as conformational 
-   ensembles at a given temperature and to compute their probabilities 
-   (section 2.5 of ref 3). This version is deposited in the CLUST-WHAM
-   subdirectory. This version has single- and multichain variants, whose
-   source codes are deposited in the following subdirectories:
-
-   a) clust-wham/src    single-chain proteins
-
-   b) clust-wham/src-M  oligomeric proteins
-
-The version developed for oligomeric proteins treats whole system as a single
-chain with dummy residues inserted. It also works for single chains but is
-not fully checked and it is recommended to use single-chain version for 
-single-chain proteins.
-
-4. INSTALLATION
----------------
-
-Customize Makefile to your system. See section 7 of the description of UNRES
-for compiler flags that are used to created executables for a particular
-force field. There are already several Makefiles prepared for various
-systems and force fields.
-
-Run make in the appropriate source directory version. CLUST-UNRES runs
-only in single-processor mode an CLUST-WHAM runs in both serial and parallel 
-mode [only conformation-distance (rmsd) calculations are parallelized]. 
-The parallel version uses MPI.
-
-5. RUNNING THE PROGRAM
-----------------------
-
-The program requires a parallel system to run. Depending on system,
-either the wham.csh C-shell script (in WHAM/bin directory) can be started
-using mpirun or the binary in the C-shell script must be executed through
-mpirun. See the wham.csh C-shell script and section 6 for the files 
-processed by the program.
-
-6. INPUT AND OUTPUT FILES
--------------------------
-
-6.1. SUMMARY OF THE FILES
--------------------------
-
-The C-shell script wham.csh is used to run the program (see the 
-bin/WHAM directory). The data files that the script needs are mostly the same as 
-for UNRES (see section 6 of UNRES description). In addition, the environmental
-variable CONTFUN specifies the method to assess whether two side chains
-are at contact; if EONTFUN=GB, the criterion defined by eq 8 of ref 4 is
-used to assess whether two side chains are at contact. Also, the parameter
-files from the C-shell scripts are overridden if the data from Hamiltonian
-MREMD are processed; if so, the parameter files are defined in the main
-input file.
-
-The main input file must have inp extension. If it is INPUT.inp, the output
-files are as follows:
-
-Coordinate input file COORD.ext, where ext denotes file extension in one of the 
-following formats:
-
-INT (extension int; UNRES angles theta, gamma, alpha, and beta),
-X   (extension x; UNRES Cartesian coordinate format; from MD),
-PDB (extension pdb; Protein Data Bank format; fro MD),
-CX  (extension cx; xdrf format; from WHAM).
-
-INPUT_clust.out (single-processor mode) or INPUT_clust.out_xxx (parallel mode) -
-     output file(s) (INPUT.out_000 is the main output file for parallel mode).
-
-COORD_clust.int: leading (lowest-energy) members of the families
-    in internal-coordinate format.
-COORD_clust.x: leading members of the families in UNRES Cartesian coordinate
-    format.
-COORD_xxxx.pdb or COORD_xxxx_yyy.pdb (CLUST-UNRES): PDB file of member yyy
-    of family xxxx; yyy is omitted if the family contains only one member
-    within a given energy cut-off.
-COORD_TxxxK_yyyy.pdb: concatenated conformations in PDB format of the 
-    members of family yyyy clustered at T=xxxK ranked by probabilities in
-    descending order at this temperature (CLUST-WHAM).
-COORD_T_xxxK_ave.pdb: cluster-averaged coordinates and coordinates of a 
-    member of each family that is closest to the cluster average in PDB
-    format, concatenated in a single file (CLUST-WHAM).
-INPUT_clust.tex: PicTeX code of the cluster tree.
-
-INPUT.rms: rmsds between conformations.
-
-6.2. MAIN INPUT FILE
---------------------
-
-This file has the same structure as the UNRES input file; most of the data are 
-input in a keyword-based form (see section 7.1 of UNRES description). The data
-are grouped into records, referred to as lines. Each record, except for the 
-records that are input in non-keyword based form, can be continued by placing
-an ampersand (&) in column 80. Such a format is referred to as the data list
-format.
-
-In the following description, the default values are given in parentheses.
-
-6.2.1. Title (80-character string)
-----------------------------------
-
-6.2.2. General data (data list format)
---------------------------------------
-
-NRES (0) - the number of residues
-
-ONE_LETTER - if present, the sequence is input in one-letter code.
-
-SYM (1) - number of chains with same sequence (for oligomeric proteins only),
-
-WITH_DIHED_CONSTR - if present, dihedral-angle restraints were imposed in the
-    processed MREMD simulations
-
-RESCALE (1) - Choice of the type of temperature dependence of the force field.
-0  - no temperature dependence
-1  - homographic dependence (not implemented yet with any force field)
-2  - hyperbolic tangent dependence [3].
-
-DISTCHAINMAX (50.0) - for oligomeric proteins, distance between the chains
-     above which restraints will be switched on to keep the chains at a
-     reasonable distance.
-
-PDBOUT - clusters will be printed in PDB format.
-
-ECUT - energy cut-off criterion to print conformations (UNRES-CLUST runs).
-     Only those families will be output the energy of the lowest-energy
-     conformation of which is within ECUT kcal/mol above that of the 
-     lowest-energy conformation and for a family only those members will be
-     output which have energy within ECUT kcal/mol above the energy of the 
-     lowest-energy member of the family.
-
-PRINT_CART - output leading members of the families in UNRES x format.
-
-PRINT_INT - output leading members of the families in UNRES int format.
-
-REF_STR - if present, reference structure is input and rmsd will be computed 
-      with respect to it (CLUST-UNRES only; rmsd is provided in the cx file 
-      from WHAM for CLUST-WHAM runs).
-
-PDBREF - if present, reference structure will be read in from a pdb file.
-
-SIDE - side chains will be considered in superposition when calculating rmsd
-
-CA_ONLY - only the Calpha atoms will be used in rmsd calculation
-
-NSTART (0) - first residue to superpose
-
-NEND (0) - last residue to superpose
-
-NTEMP (1) - number of temperatures at which probabilities will be calculated
-         and clustering performed (CLUST-WHAM)
-
-TEMPER (NTEMP tiles) - temperatures at which clustering will be performed
-        (CLUST-WHAM)
-
-EFREE - if present, conformation entropy factor is read if the conformation
-        is input from an x or pdb file
-
-PROB (0.99) - cut-off on the summary probability of the conformations that
-     are clustered at a given temperature (CLUST-WHAM)
-
-IOPT (2) - clustering algorithm:
-
-1 - Ward's minimum variance method
-2 - single link method
-3 - complete link method
-4 - average link (or group average) method
-5 - McQuitty's method
-6 - Median (Gower's) method
-7 - centroid method
-
-Instead of IOPT=1, MINTREE and instead of IOPT=2 MINVAR can be specified
-
-NCUT (1) - number of cut-offs in clustering
-
-CUTOFF (-1.0; NCUT values) cut-offs at which clustering will be performed;
-    at the cut-off flagged by a "-" sign clustering will be performed with
-    cutoff value=abs(cutoff(i)) and conformations corresponding to clusters
-    will be output in the desired format.
-
-MAKE_TREE - if present, produce a clustering-tree graph
-
-PLOT_TREE - if present, the tree is written in PicTeX format to a file
-
-PRINT_DIST - if present, distance (rmsd) matrix is printed to main output
-    file
-PUNCH_DIST - if present, the upper-triangle of the distance matrix will be 
-    printed to a file
-
-6.2.3. Energy-term weights and parameter files
-----------------------------------------------
-
-WSC (1.0) -    side-chain-side-chain interaction energy
-
-WSCP (1.0) -   side chain-peptide group interaction energy
-
-WELEC (1.0) -  peptide-group-peptide group interaction energy
-
-WEL_LOC (1.0)- third-order backbone-local correlation energy
-
-WCORR (1.0) -  fourth-order backbone-local correlation energy
-
-WCORR5 (1.0) - fifth-order backbone-local correlation energy
-
-WCORR6 (1.0) - sixth-order backbone-local correlation energy
-
-WTURN3 (1.0) - third-order backbone-local correlation energy of pairs of 
-               peptide groups separated by a single peptide group
-
-WTURN4 (1.0) - fourth-order backbone-local correlation energy of pairs of 
-               peptide groups separated by two peptide groups
-
-WTURN6 (1.0) - sixth-order backbone-local correlation energy for pairs of 
-               peptide groups separated by four peptide groups
-
-WBOND (1.0)  - virtual-bond-stretching energy
-
-WANG (1.0) -   virtual-bond-angle-bending energy
-
-WTOR (1.0) -   virtual-bond-torsional energy
-
-WTORD (1.0) -  virtual-bond-double-torsional energy
-
-WSCCOR (1.0) - sequence-specific virtual-bond-torsional energy
-
-WDIHC (0.0)  - dihedral-angle-restraint energy
-
-WHPB (1.0)   - distance-restraint energy
-
-SCAL14 (0.4) - scaling factor of 1,4-interactions
-
-6.2.4. Molecule information
------------------------------
-
-6.2.4.1. Sequence information
------------------------------
-
-Amino-acid sequence
-
-3-letter code: Sequence is input in format 20(1X,A3) 
-
-1-letter code: Sequence is input in format 80A1
-
-6.2.4.2. Dihedral angle restraint information
----------------------------------------------
-
-This is the information about dihedral-angle restraints, if any are present.
-It is specified only when WITH_DIHED_CONSTR is present in the first record.
-
-1st line: ndih_constr - number of restraints (free format)
-
-2nd line: ftors - force constant (free format)
-
-Each of the following ndih_constr lines:
-
-idih_constr(i),phi0(i),drange(i)  (free format)
-
-idih_constr(i) - the number of the dihedral angle gamma corresponding to the
-ith restraint
-
-phi0(i) - center of dihedral-angle restraint
-
-drange(i) - range of flat well (no restraints for phi0(i) +/- drange(i))
-
-6.2.4.3. Disulfide-bridge data
-------------------------------
-
-1st line: NS, (ISS(I),I=1,NS)    (free format)
-
-NS - number of cystine residues forming disulfide bridges
-
-ISS(I) - the number of the Ith disulfide-bonding cystine in the sequence
-
-2nd line: NSS, (IHPB(I),JHPB(I),I=1,NSS) (free format)
-
-NSS - number of disulfide bridges
-
-IHPB(I),JHPB(I) - the first and the second residue of ith disulfide link
-
-Because the input is in free format, each line can be split
-
-6.2.5. Reference structure
---------------------------
-
-If PDBREF is specified, filename with reference (experimental) structure,
-otherwise UNRES internal coordinates as the theta, gamma, alpha, and beta 
-angles.
-
-6.3. Main output file (out)
-------------------------------------------------
-
-The main (with name INPUT_clust.out or INPUT_clust.out_000 for parallel runs)
-output file contains the results of clustering (numbers of families
-at different cut-off values, probabilities of clusters, composition of 
-families, and rmsd values corresponding to families (0 if rmsd was not
-computed or read from WHAM-generated cx file).
-
-The output files corresponding to non-master processors
-(INPUT_clust.out_xxx where xxx>0 contain only the information up to the
-clustering protocol. These files can be deleted right after the run.
-
-Excerpts from the a sample output file are given below:
-
-CLUST-UNRES:
-
-THERE ARE   20 FAMILIES OF CONFORMATIONS
-
-FAMILY    1 CONTAINS    2 CONFORMATION(S):
-  42 -2.9384E+03  50 -2.9134E+03
-
-
-Max. distance in the family:    14.0; average distance in the family:    14.0
-
-FAMILY    2 CONTAINS    3 CONFORMATION(S):
-  13 -2.9342E+03   7 -2.8827E+03  10 -2.8682E+03
-
-CLUST-WHAM:
-
-AT CUTOFF: 200.00000
-Maximum distance found:  137.82
-Free energies and probabilities of clusters at 325.0 K
-clust   efree    prob sumprob
-    1   -76.5 0.25035 0.25035
-    2   -76.5 0.24449 0.49484
-    3   -76.4 0.21645 0.71129
-    4   -76.4 0.20045 0.91174
-    5   -75.8 0.08826 1.00000
-
-
-THERE ARE    5 FAMILIES OF CONFORMATIONS
-
-FAMILY    1 WITH TOTAL FREE ENERGY   -7.65228E+01 CONTAINS  548 CONFORMATION(S):
-8363  -7.332E+013939  -7.332E+012583  -7.332E+017395  -7.332E+019932  -7.332E+01
-5816  -7.332E+013096  -7.332E+012663  -7.332E+014099  -7.332E+016822  -7.332E+01
-3176  -7.332E+017542  -7.332E+018933  -7.332E+017315  -7.332E+01 200  -7.332E+01.
-.
-5637  -7.062E+018060  -7.061E+013797  -7.060E+018800  -7.057E+016295  -7.057E+01
-6298  -7.057E+012332  -7.057E+012709  -7.057E+01
-
-Max. distance in the family:    16.5; average distance in the family:     8.8
-Average RMSD 8.22 A
-
-6.4. Output coordinate files
-----------------------------
-
-6.4.1. The internal coordinate (int) files
-------------------------------------------
-
-The file with name COORD_clust.int contains the angles theta, gamma, alpha,
-and beta of all residues of the leaders (lowest UNRES energy conformations
-from consecutive families for CLUST-UNRES runs and lowest free energy 
-conformations for CLUST-WHAM runs). The format is the same as that of the 
-file output by UNRES; see section 9.1.1 of UNRES description.
-
-For CLUST-WHAM runs, the first line contains more items:
-
-number of family                              (format i5)
-UNRES free energy of the conformation         (format f12.3)
-Free energy of the entire family              (format f12.3)
-number of disulfide bonds                     (format i2)
-list disulfide-bonded pairs                   (format 2i3)
-conformation class number (0 if not provided) (format i10)
-
-6.4.2. The Cartesian coordinate (x) files
------------------------------------------
-
-The file with name COORD_clust.x contains the Cartesian coordinates of the 
-alpha-carbon and side-chain-center coordinates. The coordinate format is
-as in section 9.1.2 of UNRES description and the first line contains the
-following items:
-
-Number of the family                          (format I5)
-UNRES free energy of the conformation         (format f12.3)
-Free energy of the entire family              (format f12.3)
-number of disulfide bonds                     (format i2)
-list disulfide-bonded pairs                   (format 2i3)
-conformation class number (0 if not provided) (format i10)
-
-6.4.3. The PDB files
---------------------
-
-The PDB files are in standard format (see 
-ftp://ftp.wwpdb.org/pub/pdb/doc/format_descriptions/Format_v33_Letter.pdf).
-The ATOM records contain Calpha coordinates (CA) or UNRES side-chain-center
-coordinates (CB). For oligomeric proteins chain identifiers are present
-(A, B, ..., etc.) and each chain ends with a TER record. Coordinates of a 
-single conformation or multiple conformations  The header (REMARK) records 
-and the contents depends on cluster run type. The next subsections are devoted 
-to different run types. 
-
-6.4.3.1. CLUST-UNRES runs
----------------------------
-
-The files contain the members of the families obtained from clustering such
-that the lowest-energy conformation of a family is within ECUT kcal/mol higher
-in energy than the lowest-energy conformation. Again, within a family, only
-those conformations are output whose energy is within ECUT kcal/mol above
-that of the lowest-energy member of the family. Families and the members 
-of a family within a family are ranked by increasing energy. The file names are:
-
-COORD_xxxx.pdb  where xxxx is the number of the family, if the family contains
-    only one member of if only one member is output.
-
-COORD_xxxx_yyy.pdb where xxxx is the number of the family and yyy is the number
-    of the member of this family.
-
-An example is the following:
-
-REMARK R0001                            ENERGY    -2.93843E+03
-ATOM      1  CA  GLY     1       0.000   0.000   0.000
-ATOM      2  CA  HIS     2       3.800   0.000   0.000
-ATOM      3  CB  HIS     2       5.113   1.656   0.015
-ATOM      4  CA  VAL     3       5.927  -3.149   0.000
-.
-.
-.
-ATOM    346  CB  GLU   183     -43.669 -32.853  -7.320
-TER
-CONECT    1    2
-CONECT    2    4    3
-.
-.
-.
-CONECT  341  343  342
-CONECT  343  344
-CONECT  345  346
-
-where ENERGY is the UNRES energy. The CONECT records defined the Calpha-Calpha
-and Calpha-SC connection.
-
-6.4.3.2. CLUST-WHAM runs
---------------------------
-
-The program generates a file for each family with its members and a summary
-file with ensemble-averaged conformations for all families. These are described
-in the two next sections.
-
-6.4.3.2.1. Conformation family files
-------------------------------------
-
-For each family, the file name is COORD_TxxxK_yyyy.pdb, where yyyy is the
-number of the family and xxx is the integer part of the temperature (K).
-The first REMARK line in the file contains the information about the free
-energy and average rmsd of the entire cluster and, for each conformation,
-the initial REMARK line contains these quantities for this conformation.
-Same applies to oligomeric proteins, for which the TER records separate the 
-chains and the ENDMDL record separates conformations.
-An example is given below. 
-
-REMARK CLUSTER    1 FREE ENERGY  -7.65228E+01 AVE RMSD 8.22
-REMARK 1BDD L18G full clust ENERGY    -7.33241E+01 RMS  10.40
-ATOM      1  CA  VAL     1      18.059 -33.585   4.616  1.00  5.00
-ATOM      2  CB  VAL     1      18.720 -32.797   3.592  1.00  5.00
-.
-.
-.
-ATOM    115  CA  LYS    58      29.641 -44.596  -8.159  1.00  5.00
-ATOM    116  CB  LYS    58      27.593 -45.927  -8.930  1.00  5.00
-TER
-CONECT    1    3    2
-CONECT    3    5    4
-.
-.
-CONECT  113  114
-CONECT  115  116
-TER
-REMARK 1BDD L18G full clust ENERGY    -7.33240E+01 RMS  10.04
-ATOM      1  CA  VAL     1       3.174   2.833 -34.386  1.00  5.00
-ATOM      2  CB  VAL     1       3.887   2.811 -33.168  1.00  5.00
-.
-.
-ATOM    115  CA  LYS    58      16.682   6.695 -20.438  1.00  5.00
-ATOM    116  CB  LYS    58      18.925   5.540 -20.776  1.00  5.00
-TER
-CONECT    1    3    2
-CONECT    3    5    4
-CONECT  113  114
-CONECT  115  116
-TER
-
-6.4.3.2.2. Average-structure file
----------------------------------
-
-The file name is COORD_T_xxxK_ave.pdb. The entries are in pairs; the first
-one is cluster-averaged conformation and the second is a family member which
-has the lowest rmsd from this average conformation. Computing average 
-conformations is explained in section 2.5 of ref 3. Example excerpts from
-an entry corresponding to a given family are shown below. The last 
-number in each ATOM record is the rmsd of the mean coordinate of a given
-atom averaged over the cluster.
-
-REMAR AVERAGE CONFORMATIONS AT TEMPERATURE  300.00
-REMARK CLUSTER    1
-REMARK 2HEP clustering 300K ENERGY    -8.22572E+01 RMS   3.29
-ATOM      1  CA  MET     1     -17.748  48.148 -19.284  1.00  5.96
-ATOM      2  CB  MET     1     -17.373  47.911 -19.294  1.00  6.34
-ATOM      3  CA  ILE     2     -18.770  49.138 -18.133  1.00  3.98
-.
-.
-.
-ATOM     80  CB  PHE    41     -14.353  44.680 -15.642  1.00  2.62
-ATOM     81  CA  ARG    42     -11.619  41.645 -13.117  1.00  4.06
-ATOM     82  CB  ARG    42     -11.330  40.378 -13.313  1.00  5.19
-TER
-CONECT    1    3    2
-CONECT    3    5    4
-.
-.
-.
-CONECT   76   78   77
-CONECT   78   79
-CONECT   79   80
-CONECT   81   82
-TER
-REMARK 2HEP clustering 300K ENERGY    -8.22572E+01 RMS   3.29
-ATOM      1  CA  MET     1     -37.698  40.489 -32.408  1.00  5.96
-ATOM      2  CB  MET     1     -38.477  39.426 -34.159  1.00  6.34
-.
-.
-.
-ATOM     80  CB  PHE    41     -35.345  50.342 -31.371  1.00  2.62
-ATOM     81  CA  ARG    42     -33.603  54.332 -27.130  1.00  4.06
-ATOM     82  CB  ARG    42     -33.832  53.074 -24.415  1.00  5.19
-TER
-CONECT    1    3    2
-CONECT    3    5    4
-.
-.
-.
-CONECT   76   78   77
-CONECT   78   79
-CONECT   79   80
-CONECT   81   82
-TER
-
-
-6.5. The conformation-distance file
------------------------------------
-
-The file name is INPUT_clust.rms. It contains the upper-diagonal part of 
-the matrix of rmsds between conformations and differences between their
-energies:
-
-i,j,rmsd,energy(j)-energy(i) (format 2i5,2f10.5)
-
-where i and j, j>i are the numbers of the conformations, rmsd is the rmsd
-between conformation i and conformation j and energy(i) and energy(j) are
-the UNRES energies of conformations i and j, respectively.
-
-6.6. The clustering-tree PicTeX file
-------------------------------------
-
-This file contains the PicTeX code of the clustering tree. The file name is
-INPUT_clust.tex. It should be supplemented with LaTeX preamble and final 
-commands or incorporated into a LaTeX source and compiled with LaTeX. The 
-picture is produced by running LaTeX followed by dvips, dvipdf or other command
-to convert LaTeX-generated dvi files into a human-readable files.
-
-7. SUPPORT
-----------
-
-   Dr. Adam Liwo
-   Faculty of Chemistry, University of Gdansk
-   ul. Sobieskiego 18, 80-952 Gdansk Poland.
-   phone: +48 58 523 5430
-   fax: +48 58 523 5472
-   e-mail: adam@chem.univ.gda.pl
-
-   Dr. Cezary Czaplewski
-   Faculty of Chemistry, University of Gdansk
-   ul. Sobieskiego 18, 80-952 Gdansk Poland.
-   phone: +48 58 523 5430
-   fax: +48 58 523 5472
-   e-mail: czarek@chem.univ.gda.pl
-
-Prepared by Adam Liwo, 02/19/12