1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132
|
/* Copyright (C) 2006 Damir Zucic */
/*=============================================================================
hyphob_function4.c
Purpose:
Draw the hydrophobicity function F4. The sequence stored to the
main sequence buffer is used to calculate the function value.
Input:
(1) Pointer to the storage where the minimal function value will
be stored.
(2) Pointer to the storage where the maximal function value will
be stored.
(3) Pointer to RuntimeS structure.
Output:
(1) Function F4 calculated and stored.
(2) Return value.
Return value:
(1) Positive on success.
(2) Negative on failure.
Notes:
(1) The function F4 may be modified and used for many purposes.
Originally, it was introduced while searching for the method
which will be suitable for prediction of the porin secondary
structure.
========includes:============================================================*/
#include <stdio.h>
#include <X11/Xlib.h>
#include <X11/Xutil.h>
#include <X11/Xos.h>
#include <X11/Xatom.h>
#include "defines.h"
#include "typedefs.h"
/*======calculate the hydrophobicity function F4:============================*/
int HyphobFunction4_ (double *min_functionP, double *max_functionP,
RuntimeS *runtimeSP)
{
int residuesN, residueI; /* Do not use size_t instead of int ! */
int max_length;
int count;
int windowI, combinedI;
char *residue_nameP;
char residue_nameA[RESNAMESIZE];
double function_value;
/*------prepare some parameters:---------------------------------------------*/
/* The number of residues in the sequence buffer: */
residuesN = (int) runtimeSP->residuesN;
if (residuesN == 0) return -1;
/* The maximal residue name length: */
max_length = RESNAMESIZE - 1;
/*------calculate the function F4:-------------------------------------------*/
/* Initialize the extreme values: */
*min_functionP = +999999.0;
*max_functionP = -999999.0;
/* 20060608.1158: */
/* Scan the window of 7 residues. Count */
/* LEU, ILE, VAL, ALA, MET, PHE and TYR. */
/* Scan the sequence: */
for (residueI = 0; residueI < residuesN; residueI++)
{
/* Reset the function value, it might be initialized before: */
*(runtimeSP->function4P + residueI) = 0.0;
/* Reset the counter: */
count = 0;
/* Scan the sliding window: */
for (windowI = -3; windowI <= 3; windowI++)
{
/* Prepare and check the combined index: */
combinedI = residueI + windowI;
if (combinedI < 0) continue;
if (combinedI >= residuesN) continue;
/* Pointer to the residue name of */
/* the residue defined by combinedI: */
residue_nameP = runtimeSP->sequenceP + combinedI * max_length;
/* Copy the residue name: */
strncpy (residue_nameA, residue_nameP, max_length);
residue_nameA[max_length] = '\0';
/* Add +1 for selected hydrophobic residues: */
if ((strcmp (residue_nameA, "LEU") == 0) ||
(strcmp (residue_nameA, "ILE") == 0) ||
(strcmp (residue_nameA, "VAL") == 0) ||
(strcmp (residue_nameA, "MET") == 0) ||
(strcmp (residue_nameA, "ALA") == 0) ||
(strcmp (residue_nameA, "PHE") == 0) ||
(strcmp (residue_nameA, "TYR") == 0))
{
count++;
}
}
/* Scale and store the function value: */
function_value = 1.0 * (double) count;
*(runtimeSP->function4P + residueI) = function_value;
/* Find the extreme values: */
if (function_value < *min_functionP) *min_functionP = function_value;
if (function_value > *max_functionP) *max_functionP = function_value;
/* End of residueI loop: */
}
/*---------------------------------------------------------------------------*/
return 1;
}
/*===========================================================================*/
|