[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [modeller_usage] Modeller9.10 on Centos 6.2



Thanks so much for the quick response. I've tried the suggestions & get similar results to before:
[newhoir@hals cpt6short]$ python cps6.py

                         MODELLER 9.10, 2011/09/28, r8346

     PROTEIN STRUCTURE MODELLING BY SATISFACTION OF SPATIAL RESTRAINTS


                     Copyright(c) 1989-2011 Andrej Sali
                            All Rights Reserved

                             Written by A. Sali
                               with help from
           B. Webb, M.S. Madhusudhan, M-Y. Shen, M.A. Marti-Renom,
                N. Eswar, F. Alber, M. Topf, B. Oliva, A. Fiser,
                    R. Sanchez, B. Yerkovich, A. Badretdinov,
                      F. Melo, J.P. Overington, E. Feyfant
                 University of California, San Francisco, USA
                    Rockefeller University, New York, USA
                      Harvard University, Cambridge, USA
                   Imperial Cancer Research Fund, London, UK
              Birkbeck College, University of London, London, UK


Kind, OS, HostName, Kernel, Processor: 4, Linux hals 2.6.32-220.4.1.el6.x86_64 x86_64
Date and time of compilation         : 2011/09/28 18:39:17
MODELLER executable type             : x86_64-intel8
Job starting time (YY/MM/DD HH:MM:SS): 2012/03/07 06:52:00

openf___224_> Open           $(LIB)/restyp.lib
openf___224_> Open           ${MODINSTALL9v10}/modlib/resgrp.lib
rdresgr_266_> Number of residue groups:        2
openf___224_> Open           ${MODINSTALL9v10}/modlib/sstruc.lib

Dynamically allocated memory at   amaxlibraries [B,KiB,MiB]:      3234076    3158.277     3.084

Dynamically allocated memory at   amaxlibraries [B,KiB,MiB]:      3234604    3158.793     3.085
openf___224_> Open           ${MODINSTALL9v10}/modlib/resdih.lib

Dynamically allocated memory at   amaxlibraries [B,KiB,MiB]:      3283204    3206.254     3.131
rdrdih__263_> Number of dihedral angle types         :        9
              Maximal number of dihedral angle optima:        3
              Dihedral angle names                   :  Alph Phi Psi Omeg chi1 chi2 chi3 chi4 chi5
openf___224_> Open           ${MODINSTALL9v10}/modlib/radii.lib

Dynamically allocated memory at   amaxlibraries [B,KiB,MiB]:      3292444    3215.277     3.140
openf___224_> Open           ${MODINSTALL9v10}/modlib/radii14.lib
openf___224_> Open           ${MODINSTALL9v10}/modlib/af_mnchdef.lib
rdwilmo_274_> Mainchain residue conformation classes:  APBLE
openf___224_> Open           ${MODINSTALL9v10}/modlib/mnch.lib
rdclass_257_> Number of classes:        5
openf___224_> Open           ${MODINSTALL9v10}/modlib/mnch1.lib
openf___224_> Open           ${MODINSTALL9v10}/modlib/mnch2.lib
openf___224_> Open           ${MODINSTALL9v10}/modlib/mnch3.lib
openf___224_> Open           ${MODINSTALL9v10}/modlib/xs4.mat
rdrrwgh_268_> Number of residue types:       21
openf______E> Cannot open file cp6s.ali: No such file or directory

Traceback (most recent call last):
  File "cps6.py", line 42, in <module>
    a.make()                            # do the actual homology modeling
  File "/usr/lib64/python2.6/site-packages/modeller/automodel/automodel.py", line 96, in make
    self.homcsr(exit_stage)
  File "/usr/lib64/python2.6/site-packages/modeller/automodel/automodel.py", line 408, in homcsr
    aln = self.read_alignment()
  File "/usr/lib64/python2.6/site-packages/modeller/automodel/automodel.py", line 398, in read_alignment
    aln.append(file=self.alnfile, align_codes=self.knowns+[self.sequence])
  File "/usr/lib64/python2.6/site-packages/modeller/alignment.py", line 78, in append
    allow_alternates)
IOError: openf______E> Cannot open file cp6s.ali: No such file or directory

Yes, the ali file is there:
[newhoir@hals cpt6short]$ ls *ali

My .py file is:

cps6.alifrom modeller import *             
from modeller.automodel import *  

# build 20 copies of a dimer from a single template

log.verbose()    # request verbose output
# Override the 'special_restraints' and 'user_after_single_model' methods:
class MyModel(automodel):
    def special_restraints(self, aln):
        # Constrain the A and B chains to be identical (but only restrain
        # the C-alpha atoms, to reduce the number of interatomic distances
        # that need to be calculated):
        s1 = selection(self.chains['A']).only_atom_types('CA')
        s2 = selection(self.chains['B']).only_atom_types('CA')
        self.restraints.symmetry.append(symmetry(s1, s2, 1.0))
    def user_after_single_model(self):
        # Report on symmetry violations greater than 1A after building
        # each model:
        self.restraints.symmetry.report(1.0)

env = environ()  # create a new MODELLER environment to build this model in

# directories for input atom files
env.io.atom_files_directory = ['.', '/home/newhoir/rub-atfiles/']

#read HETATMS in template
env.io.hetatm = True

# Be sure to use 'MyModel' rather than 'automodel' here!
a = MyModel(env,
            alnfile  = 'cp6s.ali' ,     # alignment filename
            knowns   = ('2x06'),   # codes of the templates
            sequence = 'cps6',              # code of the target
            assess_methods=assess.GA341)     # request model assessment
a.starting_model= 1                 # index of the first model
a.ending_model  = 20               # index of the last model
                                    # (determines how many models to calculate)
a.library_schedule = autosched.slow # thorough VTFM opt
a.max_var_iterations = 300
a.md_level = refine.very_slow
a.make()                            # do the actual homology modeling

It makes no difference if I insert a first line  #!/usr/bin/python or not.

Thanks again!
Irene

----- Original Message -----
From: Modeller Caretaker <>
Date: Tuesday, March 6, 2012 12:03 pm
Subject: Re: [modeller_usage] Modeller9.10 on Centos 6.2
To: Irene Newhouse <>
Cc:

> On 3/6/12 1:31 PM, Irene Newhouse wrote:
> > I've installed Modeller9.10 on my Centos 62. linux station. I
> had an
> > earlier version working under Centos 5.x. Now I find that if I
> set an
> > alias I call mod9v10 to point to modeller9.10_x86_64-intel8 in
> > analogy to the alias I successfully set to its analog previously
>
> There's no need to do that - there's already a script called
> mod9.10
> which does basically that, which ships as part of Modeller. But
> we
> generally recommend to use Python instead if you can.
>
> > Furthermore, I get an error indicating that my
> > .ali file can't be found&  I've triple-checked to make
> sure I've
> > entered the name correctly.
>
> Are you sure it's complaining about the .ali file and not
> something
> else? Or maybe you're running in the wrong directory.
>
> > There is no modpy.sh in the directory to which the web pages point.
>
> If you installed the RPM, there is no modpy.sh because it isn't
> necessary. You can just run your regular system Python, e.g.
> "python
> foo.py".
>
> > it cannot find the equivalent of foo.py unless I give the full
> path to foo.py
>
> If you're trying to run "foo.py" in the current directory by
> just typing
> "foo.py", that typically won't work on a Unix system since the
> current
> directory isn't usually in your path. Try "./foo.py" instead.
>
> > When I specify the full path, I get error
> > messages such as "from:  command not found"&  
> a crash, indicative of
> > a massive failure to access python.
>
> This sounds like you have tried to run the Python script as a
> shell
> script. Either prefix it with "python", as above, or add
> "#!/usr/bin/python" to the first line of your Python script to
> tell the
> system to use Python to run it.
>
> Ben Webb, Modeller Caretaker
> --
> modeller-
>              http://www.salilab.org/modeller/
> Modeller mail list: http://salilab.org/mailman/listinfo/modeller_usage