lunes, 29 de febrero de 2016

orac

================================================================================
=                                                                              =
=               oooooo       rrrrrrr        aaaaa          ccccccc             =
=             oooooooooo     rrrrrrrr     aaaaaaaaa       ccccccccc            =
=            ooo      ooo    rr    rrr    aaa   aaa      ccc                   =
=            oo        oo    rr    rr    aaa     aaa     cc                    =
=            oo        oo    rr rrr      aaaaaaaaaaa     cc                    =
=            ooo      ooo    rr  rrr     aaaaaaaaaaa     ccc                   =
=             oooooooooo     rr    rr    aa       aa      ccccccccc            =
=               oooooo       rr    rrr   aa       aa       ccccccc             =
=                                                                              =
=                              (Version 4.0.1)                                 =
=                                                                              =
=     "A Molecular Dynamics Program to Simulate Complex Molecular Systems"     =
=                                                                              =
=                     Copyright(C) 1989 - 1997                                 =
=                                                                              =
=                        All Right Reserved                                    =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=   ORAC is provided "as is" and without any warranty express or implied.      =
=   The user assumes all risks of using ORAC.                                  =
=                                                                              =
=   The user may make copies of ORAC for his/her own use, and modify those     =
=   copies. The user MAY NOT distribute any copy of the original or            =
=   modified source code to any users at any sites other than his/her own.     =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                      Centre  d'Etudes Saclay                                 =
=                       Gif sur Yvette, FRANCE                                 =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
================================================================================


=========================== INPUT TO THE RUN ===================================
=                                                                              =
=###############################################################               =
=#  Minimize Crystallographic structure form PDBank                            =
=###############################################################               =
=                                                                              =
=#                                                                             =
=# Set MD cell and read pdb coordinates                                        =
=#                                                                             =
=&SETUP                                                                        =
=   CRYSTAL  50.0 50.0 50.0 90.0 90.0 90.0                                     =
=   READ_PDB  ../pdb/bpti0.pdb                                                 =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=# read ASCII databases and build up solute                                    =
=#                                                                             =
=&PARAMETERS                                                                   =
=   READ_TPG_ASCII  ../lib/amber95.tpg                                         =
=   READ_PRM_ASCII  ../lib/amber95.prm                                         =
=   ADD_TPG SOLUTE                                                             =
=     bond 1sg 2sg residue  5 55                                               =
=     bond 1sg 2sg residue 14 38                                               =
=     bond 1sg 2sg residue 30 51                                               =
=   END                                                                        =
=   JOIN SOLUTE                                                                =
=       arg-h pro asp phe cys leu glu pro pro tyr thr gly pro                  =
=       cys lys ala arg ile ile arg tyr phe tyr asn ala lys                    =
=       ala gly leu cys gln thr phe val tyr gly gly cys arg                    =
=       ala lys arg asn asn phe lys ser ala glu asp cys met                    =
=       arg thr cys gly gly ala-o po4 hoh x 167                                =
=   END                                                                        =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  Simulation  Commands:  Minimize the structure                              =
=#                                                                             =
=#                                                                             =
=&SIMULATION                                                                   =
=   MINIMIZE                                                                   =
=      CG  0.01                                                                =
=   END                                                                        =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  Cutoff for minimize is 7.0 A.                                              =
=#                                                                             =
=&POTENTIAL                                                                    =
=   UPDATE      20.0   1.5                                                     =
=   CUTOFF  7.0                                                                =
=   STRETCHING                                                                 =
=   QQ-FUDGE  0.83333                                                          =
=   LJ-FUDGE  0.50                                                             =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  do 20 minimization step and intermediate printout every 5                  =
=#                                                                             =
=&RUN                                                                          =
=   CONTROL      0                                                             =
=   TIME         20.0                                                          =
=   PRINT         5.0                                                          =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=# write final pdb file to 1.pdb                                               =
=#                                                                             =
=&INOUT                                                                        =
=   ASCII    20.0 OPEN 1.pdb                                                   =
=&END                                                                          =
=                                                                              =
================================================================================


     Checking Input ......              ---->

  * * * W A R N I N G * * *
     PROPERTY(&RUN) is O or less than TIMESTEP(&VARIABLES); NPROP set to 1000      


 ********* W A R N I N G  W A R N I N G ************
 *  1 ERRORS WHILE EXECUTING VERIFY_INPUT                                        
 ******** Recoverable Error Program Continues ******



   ========================================================================
   =                                                                      =
   =                 SOLUTE:   Assembling Molecules                       =
   =                                                                      =
   ========================================================================

     SOLUTE:  Linking the molecule together      ---->
     SOLUTE:  Adding extra links                 ---->
     SOLUTE:  Searching for Bonds                ---->
     SOLUTE:  Searching for Bendings             ---->
     SOLUTE:  Searching for Bendings to Omit     ---->
     SOLUTE:  Assigning a Type to Each Atom      ---->
     SOLUTE:  Find All Possible 1-4 Interactions ---->
     SOLUTE:  Searching for Bond Parameters  ---->
     SOLUTE:  Searching for Bendings Parameters  ---->
     SOLUTE:  Pick rigid fragments that Match Input  ---->
     SOLUTE:  Add extra P-Torsions               ---->
     SOLUTE:  Match P-Torsions with Parameters   ---->
     SOLUTE:  Pick I-Torsions that Match Input   ---->
     SOLUTE:  Add extra I-Torsions               ---->
     SOLUTE:  Match I-Torsions with Parameters   ---->

  ************************************************************************
  *                   SOLUTE:    Initial TOPOLOGY  List                  *
  *                                                                      *
  *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
  *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
  *             199 I-Torsions     2347 1-4 Inter.                       *
  *                                                                      *
  ************************************************************************

     Found Charge   6.00000 on Solute Molecule      1
     Found Charge  -3.00000 on Solute Molecule      2
 These parameters are not compatible :
i        1     0.00000     0.00000
j       54     2.00000     0.11000

  * * * E R R O R  * * *
     Wrong Format in 14 LJ parameter file                                          
 These parameters are not compatible :
i        1     0.00000     0.00000
j       55     2.00000     0.11000
 These parameters are not compatible :
i        1     0.00000     0.00000
j       56     2.00000     0.11000
 These parameters are not compatible :
i        1     0.00000     0.00000
j       57     1.32000     0.02200
 .... and many other error like these..

   ========================================================================
   =                                                                      =
   =                 Setting up the Simulation Box                        =
   =                                                                      =
   ========================================================================



                    ==========================================
                    =       Reading System Coordinates       =
                    ==========================================


     Reading the PDB file               ---->
     PDB file parsed:  1398 atoms; -    772 hydrogens to be assigned

                    Solute Center of Mass coordinates at
                          28.546       9.649       0.354



  ************************************************************************
  *                        New TOPOLOGY  List                            *
  *                                                                      *
  *               0 Solvent Mol       0 Atoms Each   169 Solute Mol.     *
  *                                                                      *
  *                            For the Sysyem                            *
  *                                                                      *
  *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
  *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
  *             199 I-Torsions     2347 1-4 Inter.                       *
  *                                                                      *
  ************************************************************************




********************************************************************************
*                                                                              *
*              I n p u t   O p e r a t i o n s   C o m p l e t e d             *
*                                                                              *
*              I n p u t   O p e r a t i o n s   C o m p l e t e d             *
*                                                                              *
*                                                                              *
********************************************************************************







*******************************************************************************
*                                                                             *
*            Minimization for an aggregate composed of 1398 atoms             *
*            Minimization for an aggregate composed of 1398 atoms             *
*            Minimization for an aggregate composed of 1398 atoms             *
*                                                                             *
*******************************************************************************


==============  P A R A M E T E R S    F O R   T H E    R U N  =================


           The Minimization will be run from scratch.
                     Initial atomic coordinates will be generated.
           Reject =     0 Steps     Run =    20 Steps
           Maxrun =     0 Steps

           Save Subaverages every   1000 Steps
           Save Restart file every     0 Steps
           Save PLOT file every        0 Steps
           Save PDB  file every       20 Steps
           Save DUMP file every        0 Steps

          potential truncation distance = 7.0000 ang

          time-step = 1.0000E+00 fs.

           Number of atoms   =  1398

                             CO-matrix

                         25.00000   0.00000   0.00000
                          0.00000  25.00000   0.00000
                          0.00000   0.00000  25.00000

           Volume            =     125000.000  Ang**3

----------------------------- Neighbor List ------------------------------------

           UPDATE every    20 Steps

           Dispersion+Electrostatics  rspcut  =     8.50



     ************************************************************************
     *                        System  TOPOLOGY  List                        *
     *                                                                      *
     *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
     *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
     *             199 I-Torsions     2347 1-4 Inter.   524 Atomic Groups   *
     *             226 Units           169 Molecules                        *
     *                                                                      *
     ************************************************************************



================================================================================





********************************************************************************
*                                                                              *
*                   I n t e r m e d i a t e   R e s u l t s                    *
*                                                                              *
*                   I n t e r m e d i a t e   R e s u l t s                    *
*                                                                              *
*                                                                              *
********************************************************************************



Neighbor Lists Dimensions     *neighbor(  18459)*


     CPUtime for update          =     0.030


                     ***************************************
                     *        Not yet at minimum !!        *
                     *        GradMax     0.69843E+06      *
                     ***************************************




           Tstep    =       5.000 Total    =   57126.724 TotPot   =   57126.724
           Coulomb  =   -7652.840 NonBond  =     403.263 Ener14   =   10117.494
           Bonded   =   56723.462 Stretch  =   45011.605 Angle    =    9592.190
           I-Tors   =      25.776 P-Tors   =    2093.891 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.26772E+01      *
                     ***************************************



           Tstep    =      10.000 Total    =   -3787.220 TotPot   =   -3787.220
           Coulomb  =   -8853.156 NonBond  =   -7784.520 Ener14   =    9909.408
           Bonded   =    3997.301 Stretch  =     440.460 Angle    =    1458.502
           I-Tors   =      61.854 P-Tors   =    2036.484 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.93677E-01      *
                     ***************************************



           Tstep    =      15.000 Total    =   -8170.596 TotPot   =   -8170.596
           Coulomb  =  -11344.386 NonBond  =  -11368.233 Ener14   =    9629.103
           Bonded   =    3197.638 Stretch  =     332.359 Angle    =     935.345
           I-Tors   =     101.311 P-Tors   =    1828.623 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.40184E-01      *
                     ***************************************



           Tstep    =      20.000 Total    =  -10481.638 TotPot   =  -10481.638
           Coulomb  =  -13032.542 NonBond  =  -13275.800 Ener14   =    9445.498
           Bonded   =    2794.163 Stretch  =     294.359 Angle    =     835.734
           I-Tors   =      76.857 P-Tors   =    1587.213 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.30852E-01      *
                     ***************************************



                     ***************************************
                     *                                     *
                     *        M I N I M I Z A T I O N      *
                     *                                     *
                     *       F I N A L   R E S U L T S     *
                     *                                     *
                     ***************************************




           Tstep    =      20.000 Total    =  -10481.638 TotPot   =  -10481.638
           Coulomb  =  -13032.542 NonBond  =  -13275.800 Ener14   =    9445.498
           Bonded   =    2794.163 Stretch  =     294.359 Angle    =     835.734
           I-Tors   =      76.857 P-Tors   =    1587.213 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.29430E-01      *
                     ***************************************



          ===========================================

                Total cpu time for the run     =     11.840
                Total elapsed time for the run =      0.000
                Averaged time per step         =      0.592
                Averaged elapsed per step      =      0.000



          ===========================================



     ************************************************************************
     *                                                                      *
     *                        PROGRAM COMPLETED                             *
     *                                                                      *
     ************************************************************************



          |--------------------------------------------------------------|
          |    ``They're unfriendly.  Which is fortunate really;         |
          |      they'd be difficult to love.''                          |
          |                              . . . Avon    Blake's Seven     |
          |                                                              |
          |--------------------------------------------------------------|

================================================================================
=                                                                              =
=               oooooo       rrrrrrr        aaaaa          ccccccc             =
=             oooooooooo     rrrrrrrr     aaaaaaaaa       ccccccccc            =
=            ooo      ooo    rr    rrr    aaa   aaa      ccc                   =
=            oo        oo    rr    rr    aaa     aaa     cc                    =
=            oo        oo    rr rrr      aaaaaaaaaaa     cc                    =
=            ooo      ooo    rr  rrr     aaaaaaaaaaa     ccc                   =
=             oooooooooo     rr    rr    aa       aa      ccccccccc            =
=               oooooo       rr    rrr   aa       aa       ccccccc             =
=                                                                              =
=                              (Version 4.0.1)                                 =
=                                                                              =
=     "A Molecular Dynamics Program to Simulate Complex Molecular Systems"     =
=                                                                              =
=                     Copyright(C) 1989 - 1997                                 =
=                                                                              =
=                        All Right Reserved                                    =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=   ORAC is provided "as is" and without any warranty express or implied.      =
=   The user assumes all risks of using ORAC.                                  =
=                                                                              =
=   The user may make copies of ORAC for his/her own use, and modify those     =
=   copies. The user MAY NOT distribute any copy of the original or            =
=   modified source code to any users at any sites other than his/her own.     =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                      Centre  d'Etudes Saclay                                 =
=                       Gif sur Yvette, FRANCE                                 =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
================================================================================


=========================== INPUT TO THE RUN ===================================
=                                                                              =
=###############################################################               =
=#  Minimize Crystallographic structure form PDBank                            =
=###############################################################               =
=                                                                              =
=#                                                                             =
=# Set MD cell and read pdb coordinates                                        =
=#                                                                             =
=&SETUP                                                                        =
=   CRYSTAL  50.0 50.0 50.0 90.0 90.0 90.0                                     =
=   READ_PDB  ../pdb/bpti0.pdb                                                 =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=# read ASCII databases and build up solute                                    =
=#                                                                             =
=&PARAMETERS                                                                   =
=   READ_TPG_ASCII  ../lib/amber95.tpg                                         =
=   READ_PRM_ASCII  ../lib/amber95.prm                                         =
=   ADD_TPG SOLUTE                                                             =
=     bond 1sg 2sg residue  5 55                                               =
=     bond 1sg 2sg residue 14 38                                               =
=     bond 1sg 2sg residue 30 51                                               =
=   END                                                                        =
=   JOIN SOLUTE                                                                =
=       arg-h pro asp phe cys leu glu pro pro tyr thr gly pro                  =
=       cys lys ala arg ile ile arg tyr phe tyr asn ala lys                    =
=       ala gly leu cys gln thr phe val tyr gly gly cys arg                    =
=       ala lys arg asn asn phe lys ser ala glu asp cys met                    =
=       arg thr cys gly gly ala-o po4 hoh x 167                                =
=   END                                                                        =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  Simulation  Commands:  Minimize the structure                              =
=#                                                                             =
=#                                                                             =
=&SIMULATION                                                                   =
=   MINIMIZE                                                                   =
=      CG  0.01                                                                =
=   END                                                                        =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  Cutoff for minimize is 7.0 A.                                              =
=#                                                                             =
=&POTENTIAL                                                                    =
=   UPDATE      20.0   1.5                                                     =
=   CUTOFF  7.0                                                                =
=   STRETCHING                                                                 =
=   QQ-FUDGE  0.83333                                                          =
=   LJ-FUDGE  0.50                                                             =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  do 20 minimization step and intermediate printout every 5                  =
=#                                                                             =
=&RUN                                                                          =
=   CONTROL      0                                                             =
=   TIME         20.0                                                          =
=   PRINT         5.0                                                          =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=# write final pdb file to 1.pdb                                               =
=#                                                                             =
=&INOUT                                                                        =
=   ASCII    20.0 OPEN 1.pdb                                                   =
=&END                                                                          =
=                                                                              =
================================================================================


     Checking Input ......              ---->

  * * * W A R N I N G * * *
     PROPERTY(&RUN) is O or less than TIMESTEP(&VARIABLES); NPROP set to 1000      


 ********* W A R N I N G  W A R N I N G ************
 *  1 ERRORS WHILE EXECUTING VERIFY_INPUT                                        
 ******** Recoverable Error Program Continues ******



   ========================================================================
   =                                                                      =
   =                 SOLUTE:   Assembling Molecules                       =
   =                                                                      =
   ========================================================================

     SOLUTE:  Linking the molecule together      ---->
     SOLUTE:  Adding extra links                 ---->
     SOLUTE:  Searching for Bonds                ---->
     SOLUTE:  Searching for Bendings             ---->
     SOLUTE:  Searching for Bendings to Omit     ---->
     SOLUTE:  Assigning a Type to Each Atom      ---->
     SOLUTE:  Find All Possible 1-4 Interactions ---->
     SOLUTE:  Searching for Bond Parameters  ---->
     SOLUTE:  Searching for Bendings Parameters  ---->
     SOLUTE:  Pick rigid fragments that Match Input  ---->
     SOLUTE:  Add extra P-Torsions               ---->
     SOLUTE:  Match P-Torsions with Parameters   ---->
     SOLUTE:  Pick I-Torsions that Match Input   ---->
     SOLUTE:  Add extra I-Torsions               ---->
     SOLUTE:  Match I-Torsions with Parameters   ---->

  ************************************************************************
  *                   SOLUTE:    Initial TOPOLOGY  List                  *
  *                                                                      *
  *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
  *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
  *             199 I-Torsions     2347 1-4 Inter.                       *
  *                                                                      *
  ************************************************************************

     Found Charge   6.00000 on Solute Molecule      1
     Found Charge  -3.00000 on Solute Molecule      2
 These parameters are not compatible :
i        1     0.00000     0.00000
j       54     2.00000     0.11000

  * * * E R R O R  * * *
     Wrong Format in 14 LJ parameter file                                          
 These parameters are not compatible :
i        1     0.00000     0.00000
j       55     2.00000     0.11000
 These parameters are not compatible :
i        1     0.00000     0.00000
j       56     2.00000     0.11000
 These parameters are not compatible :
i        1     0.00000     0.00000
j       57     1.32000     0.02200
 .... and many other error like these..

   ========================================================================
   =                                                                      =
   =                 Setting up the Simulation Box                        =
   =                                                                      =
   ========================================================================



                    ==========================================
                    =       Reading System Coordinates       =
                    ==========================================


     Reading the PDB file               ---->
     PDB file parsed:  1398 atoms; -    772 hydrogens to be assigned

                    Solute Center of Mass coordinates at
                          28.546       9.649       0.354



  ************************************************************************
  *                        New TOPOLOGY  List                            *
  *                                                                      *
  *               0 Solvent Mol       0 Atoms Each   169 Solute Mol.     *
  *                                                                      *
  *                            For the Sysyem                            *
  *                                                                      *
  *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
  *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
  *             199 I-Torsions     2347 1-4 Inter.                       *
  *                                                                      *
  ************************************************************************




********************************************************************************
*                                                                              *
*              I n p u t   O p e r a t i o n s   C o m p l e t e d             *
*                                                                              *
*              I n p u t   O p e r a t i o n s   C o m p l e t e d             *
*                                                                              *
*                                                                              *
********************************************************************************







*******************************************************************************
*                                                                             *
*            Minimization for an aggregate composed of 1398 atoms             *
*            Minimization for an aggregate composed of 1398 atoms             *
*            Minimization for an aggregate composed of 1398 atoms             *
*                                                                             *
*******************************************************************************


==============  P A R A M E T E R S    F O R   T H E    R U N  =================


           The Minimization will be run from scratch.
                     Initial atomic coordinates will be generated.
           Reject =     0 Steps     Run =    20 Steps
           Maxrun =     0 Steps

           Save Subaverages every   1000 Steps
           Save Restart file every     0 Steps
           Save PLOT file every        0 Steps
           Save PDB  file every       20 Steps
           Save DUMP file every        0 Steps

          potential truncation distance = 7.0000 ang

          time-step = 1.0000E+00 fs.

           Number of atoms   =  1398

                             CO-matrix

                         25.00000   0.00000   0.00000
                          0.00000  25.00000   0.00000
                          0.00000   0.00000  25.00000

           Volume            =     125000.000  Ang**3

----------------------------- Neighbor List ------------------------------------

           UPDATE every    20 Steps

           Dispersion+Electrostatics  rspcut  =     8.50



     ************************************************************************
     *                        System  TOPOLOGY  List                        *
     *                                                                      *
     *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
     *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
     *             199 I-Torsions     2347 1-4 Inter.   524 Atomic Groups   *
     *             226 Units           169 Molecules                        *
     *                                                                      *
     ************************************************************************



================================================================================





********************************************************************************
*                                                                              *
*                   I n t e r m e d i a t e   R e s u l t s                    *
*                                                                              *
*                   I n t e r m e d i a t e   R e s u l t s                    *
*                                                                              *
*                                                                              *
********************************************************************************



Neighbor Lists Dimensions     *neighbor(  18459)*


     CPUtime for update          =     0.030


                     ***************************************
                     *        Not yet at minimum !!        *
                     *        GradMax     0.69843E+06      *
                     ***************************************




           Tstep    =       5.000 Total    =   57126.724 TotPot   =   57126.724
           Coulomb  =   -7652.840 NonBond  =     403.263 Ener14   =   10117.494
           Bonded   =   56723.462 Stretch  =   45011.605 Angle    =    9592.190
           I-Tors   =      25.776 P-Tors   =    2093.891 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.26772E+01      *
                     ***************************************



           Tstep    =      10.000 Total    =   -3787.220 TotPot   =   -3787.220
           Coulomb  =   -8853.156 NonBond  =   -7784.520 Ener14   =    9909.408
           Bonded   =    3997.301 Stretch  =     440.460 Angle    =    1458.502
           I-Tors   =      61.854 P-Tors   =    2036.484 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.93677E-01      *
                     ***************************************



           Tstep    =      15.000 Total    =   -8170.596 TotPot   =   -8170.596
           Coulomb  =  -11344.386 NonBond  =  -11368.233 Ener14   =    9629.103
           Bonded   =    3197.638 Stretch  =     332.359 Angle    =     935.345
           I-Tors   =     101.311 P-Tors   =    1828.623 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.40184E-01      *
                     ***************************************



           Tstep    =      20.000 Total    =  -10481.638 TotPot   =  -10481.638
           Coulomb  =  -13032.542 NonBond  =  -13275.800 Ener14   =    9445.498
           Bonded   =    2794.163 Stretch  =     294.359 Angle    =     835.734
           I-Tors   =      76.857 P-Tors   =    1587.213 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.30852E-01      *
                     ***************************************



                     ***************************************
                     *                                     *
                     *        M I N I M I Z A T I O N      *
                     *                                     *
                     *       F I N A L   R E S U L T S     *
                     *                                     *
                     ***************************************




           Tstep    =      20.000 Total    =  -10481.638 TotPot   =  -10481.638
           Coulomb  =  -13032.542 NonBond  =  -13275.800 Ener14   =    9445.498
           Bonded   =    2794.163 Stretch  =     294.359 Angle    =     835.734
           I-Tors   =      76.857 P-Tors   =    1587.213 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.29430E-01      *
                     ***************************************



          ===========================================

                Total cpu time for the run     =     11.880
                Total elapsed time for the run =      0.000
                Averaged time per step         =      0.594
                Averaged elapsed per step      =      0.000



          ===========================================



     ************************************************************************
     *                                                                      *
     *                        PROGRAM COMPLETED                             *
     *                                                                      *
     ************************************************************************



          |--------------------------------------------------------------|
          |    ``They're unfriendly.  Which is fortunate really;         |
          |      they'd be difficult to love.''                          |
          |                              . . . Avon    Blake's Seven     |
          |                                                              |
          |--------------------------------------------------------------|

================================================================================
=                                                                              =
=               oooooo       rrrrrrr        aaaaa          ccccccc             =
=             oooooooooo     rrrrrrrr     aaaaaaaaa       ccccccccc            =
=            ooo      ooo    rr    rrr    aaa   aaa      ccc                   =
=            oo        oo    rr    rr    aaa     aaa     cc                    =
=            oo        oo    rr rrr      aaaaaaaaaaa     cc                    =
=            ooo      ooo    rr  rrr     aaaaaaaaaaa     ccc                   =
=             oooooooooo     rr    rr    aa       aa      ccccccccc            =
=               oooooo       rr    rrr   aa       aa       ccccccc             =
=                                                                              =
=                              (Version 4.0.1)                                 =
=                                                                              =
=     "A Molecular Dynamics Program to Simulate Complex Molecular Systems"     =
=                                                                              =
=                     Copyright(C) 1989 - 1997                                 =
=                                                                              =
=                        All Right Reserved                                    =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=   ORAC is provided "as is" and without any warranty express or implied.      =
=   The user assumes all risks of using ORAC.                                  =
=                                                                              =
=   The user may make copies of ORAC for his/her own use, and modify those     =
=   copies. The user MAY NOT distribute any copy of the original or            =
=   modified source code to any users at any sites other than his/her own.     =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                      Centre  d'Etudes Saclay                                 =
=                       Gif sur Yvette, FRANCE                                 =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
=                                                                              =
================================================================================


=========================== INPUT TO THE RUN ===================================
=                                                                              =
=###############################################################               =
=#  Minimize Crystallographic structure form PDBank                            =
=###############################################################               =
=                                                                              =
=#                                                                             =
=# Set MD cell and read pdb coordinates                                        =
=#                                                                             =
=&SETUP                                                                        =
=   CRYSTAL  50.0 50.0 50.0 90.0 90.0 90.0                                     =
=   READ_PDB  ../pdb/bpti0.pdb                                                 =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=# read ASCII databases and build up solute                                    =
=#                                                                             =
=&PARAMETERS                                                                   =
=   READ_TPG_ASCII  ../lib/amber95.tpg                                         =
=   READ_PRM_ASCII  ../lib/amber95.prm                                         =
=   ADD_TPG SOLUTE                                                             =
=     bond 1sg 2sg residue  5 55                                               =
=     bond 1sg 2sg residue 14 38                                               =
=     bond 1sg 2sg residue 30 51                                               =
=   END                                                                        =
=   JOIN SOLUTE                                                                =
=       arg-h pro asp phe cys leu glu pro pro tyr thr gly pro                  =
=       cys lys ala arg ile ile arg tyr phe tyr asn ala lys                    =
=       ala gly leu cys gln thr phe val tyr gly gly cys arg                    =
=       ala lys arg asn asn phe lys ser ala glu asp cys met                    =
=       arg thr cys gly gly ala-o po4 hoh x 167                                =
=   END                                                                        =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  Simulation  Commands:  Minimize the structure                              =
=#                                                                             =
=#                                                                             =
=&SIMULATION                                                                   =
=   MINIMIZE                                                                   =
=      CG  0.01                                                                =
=   END                                                                        =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  Cutoff for minimize is 7.0 A.                                              =
=#                                                                             =
=&POTENTIAL                                                                    =
=   UPDATE      20.0   1.5                                                     =
=   CUTOFF  7.0                                                                =
=   STRETCHING                                                                 =
=   QQ-FUDGE  0.83333                                                          =
=   LJ-FUDGE  0.50                                                             =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=#  do 20 minimization step and intermediate printout every 5                  =
=#                                                                             =
=&RUN                                                                          =
=   CONTROL      0                                                             =
=   TIME         20.0                                                          =
=   PRINT         5.0                                                          =
=&END                                                                          =
=                                                                              =
=#                                                                             =
=# write final pdb file to 1.pdb                                               =
=#                                                                             =
=&INOUT                                                                        =
=   ASCII    20.0 OPEN 1.pdb                                                   =
=&END                                                                          =
=                                                                              =
================================================================================


     Checking Input ......              ---->

  * * * W A R N I N G * * *
     PROPERTY(&RUN) is O or less than TIMESTEP(&VARIABLES); NPROP set to 1000      


 ********* W A R N I N G  W A R N I N G ************
 *  1 ERRORS WHILE EXECUTING VERIFY_INPUT                                        
 ******** Recoverable Error Program Continues ******



   ========================================================================
   =                                                                      =
   =                 SOLUTE:   Assembling Molecules                       =
   =                                                                      =
   ========================================================================

     SOLUTE:  Linking the molecule together      ---->
     SOLUTE:  Adding extra links                 ---->
     SOLUTE:  Searching for Bonds                ---->
     SOLUTE:  Searching for Bendings             ---->
     SOLUTE:  Searching for Bendings to Omit     ---->
     SOLUTE:  Assigning a Type to Each Atom      ---->
     SOLUTE:  Find All Possible 1-4 Interactions ---->
     SOLUTE:  Searching for Bond Parameters  ---->
     SOLUTE:  Searching for Bendings Parameters  ---->
     SOLUTE:  Pick rigid fragments that Match Input  ---->
     SOLUTE:  Add extra P-Torsions               ---->
     SOLUTE:  Match P-Torsions with Parameters   ---->
     SOLUTE:  Pick I-Torsions that Match Input   ---->
     SOLUTE:  Add extra I-Torsions               ---->
     SOLUTE:  Match I-Torsions with Parameters   ---->

  ************************************************************************
  *                   SOLUTE:    Initial TOPOLOGY  List                  *
  *                                                                      *
  *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
  *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
  *             199 I-Torsions     2347 1-4 Inter.                       *
  *                                                                      *
  ************************************************************************

     Found Charge   6.00000 on Solute Molecule      1
     Found Charge  -3.00000 on Solute Molecule      2
 These parameters are not compatible :
i        1     0.00000     0.00000
j       54     2.00000     0.11000

  * * * E R R O R  * * *
     Wrong Format in 14 LJ parameter file                                          
 These parameters are not compatible :
i        1     0.00000     0.00000
j       55     2.00000     0.11000
 These parameters are not compatible :
i        1     0.00000     0.00000
j       56     2.00000     0.11000
 These parameters are not compatible :
i        1     0.00000     0.00000
j       57     1.32000     0.02200
 .... and many other error like these..

   ========================================================================
   =                                                                      =
   =                 Setting up the Simulation Box                        =
   =                                                                      =
   ========================================================================



                    ==========================================
                    =       Reading System Coordinates       =
                    ==========================================


     Reading the PDB file               ---->
     PDB file parsed:  1398 atoms; -    772 hydrogens to be assigned

                    Solute Center of Mass coordinates at
                          28.546       9.649       0.354



  ************************************************************************
  *                        New TOPOLOGY  List                            *
  *                                                                      *
  *               0 Solvent Mol       0 Atoms Each   169 Solute Mol.     *
  *                                                                      *
  *                            For the Sysyem                            *
  *                                                                      *
  *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
  *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
  *             199 I-Torsions     2347 1-4 Inter.                       *
  *                                                                      *
  ************************************************************************




********************************************************************************
*                                                                              *
*              I n p u t   O p e r a t i o n s   C o m p l e t e d             *
*                                                                              *
*              I n p u t   O p e r a t i o n s   C o m p l e t e d             *
*                                                                              *
*                                                                              *
********************************************************************************







*******************************************************************************
*                                                                             *
*            Minimization for an aggregate composed of 1398 atoms             *
*            Minimization for an aggregate composed of 1398 atoms             *
*            Minimization for an aggregate composed of 1398 atoms             *
*                                                                             *
*******************************************************************************


==============  P A R A M E T E R S    F O R   T H E    R U N  =================


           The Minimization will be run from scratch.
                     Initial atomic coordinates will be generated.
           Reject =     0 Steps     Run =    20 Steps
           Maxrun =     0 Steps

           Save Subaverages every   1000 Steps
           Save Restart file every     0 Steps
           Save PLOT file every        0 Steps
           Save PDB  file every       20 Steps
           Save DUMP file every        0 Steps

          potential truncation distance = 7.0000 ang

          time-step = 1.0000E+00 fs.

           Number of atoms   =  1398

                             CO-matrix

                         25.00000   0.00000   0.00000
                          0.00000  25.00000   0.00000
                          0.00000   0.00000  25.00000

           Volume            =     125000.000  Ang**3

----------------------------- Neighbor List ------------------------------------

           UPDATE every    20 Steps

           Dispersion+Electrostatics  rspcut  =     8.50



     ************************************************************************
     *                        System  TOPOLOGY  List                        *
     *                                                                      *
     *            1398 Atoms          1244 Bonds       1244 FLexible Bonds  *
     *               0 Rigid Bonds    1799 Angles      2732 P-Torsions      *
     *             199 I-Torsions     2347 1-4 Inter.   524 Atomic Groups   *
     *             226 Units           169 Molecules                        *
     *                                                                      *
     ************************************************************************



================================================================================





********************************************************************************
*                                                                              *
*                   I n t e r m e d i a t e   R e s u l t s                    *
*                                                                              *
*                   I n t e r m e d i a t e   R e s u l t s                    *
*                                                                              *
*                                                                              *
********************************************************************************



Neighbor Lists Dimensions     *neighbor(  18459)*


     CPUtime for update          =     0.030


                     ***************************************
                     *        Not yet at minimum !!        *
                     *        GradMax     0.69843E+06      *
                     ***************************************




           Tstep    =       5.000 Total    =   57126.724 TotPot   =   57126.724
           Coulomb  =   -7652.840 NonBond  =     403.263 Ener14   =   10117.494
           Bonded   =   56723.462 Stretch  =   45011.605 Angle    =    9592.190
           I-Tors   =      25.776 P-Tors   =    2093.891 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.26772E+01      *
                     ***************************************



           Tstep    =      10.000 Total    =   -3787.220 TotPot   =   -3787.220
           Coulomb  =   -8853.156 NonBond  =   -7784.520 Ener14   =    9909.408
           Bonded   =    3997.301 Stretch  =     440.460 Angle    =    1458.502
           I-Tors   =      61.854 P-Tors   =    2036.484 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.93677E-01      *
                     ***************************************



           Tstep    =      15.000 Total    =   -8170.596 TotPot   =   -8170.596
           Coulomb  =  -11344.386 NonBond  =  -11368.233 Ener14   =    9629.103
           Bonded   =    3197.638 Stretch  =     332.359 Angle    =     935.345
           I-Tors   =     101.311 P-Tors   =    1828.623 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.40184E-01      *
                     ***************************************



           Tstep    =      20.000 Total    =  -10481.638 TotPot   =  -10481.638
           Coulomb  =  -13032.542 NonBond  =  -13275.800 Ener14   =    9445.498
           Bonded   =    2794.163 Stretch  =     294.359 Angle    =     835.734
           I-Tors   =      76.857 P-Tors   =    1587.213 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.30852E-01      *
                     ***************************************



                     ***************************************
                     *                                     *
                     *        M I N I M I Z A T I O N      *
                     *                                     *
                     *       F I N A L   R E S U L T S     *
                     *                                     *
                     ***************************************




           Tstep    =      20.000 Total    =  -10481.638 TotPot   =  -10481.638
           Coulomb  =  -13032.542 NonBond  =  -13275.800 Ener14   =    9445.498
           Bonded   =    2794.163 Stretch  =     294.359 Angle    =     835.734
           I-Tors   =      76.857 P-Tors   =    1587.213 TotTemp  =       0.000
           ResTemp  =       0.000 TraTemp  =       0.000

                     ***************************************
                     *        GradMax     0.29430E-01      *
                     ***************************************



          ===========================================

                Total cpu time for the run     =     11.880
                Total elapsed time for the run =      0.000
                Averaged time per step         =      0.594
                Averaged elapsed per step      =      0.000



          ===========================================



     ************************************************************************
     *                                                                      *
     *                        PROGRAM COMPLETED                             *
     *                                                                      *
     ************************************************************************



          |--------------------------------------------------------------|
          |    ``They're unfriendly.  Which is fortunate really;         |
          |      they'd be difficult to love.''                          |
          |                              . . . Avon    Blake's Seven     |
          |                                                              |
          |--------------------------------------------------------------|


No hay comentarios:

Publicar un comentario

zen consultora

Blogger Widgets

Entrada destacada

Platzy y el payaso Freddy Vega, PLATZI APESTA, PLATZI NO SIRVE, PLATZI ES UNA ESTAFA

  Platzy y los payasos fredy vega y cvander parte 1,  PLATZI ES UNA ESTAFA Hola amigos, este post va a ir creciendo conforme vaya escribiend...