tutorial update, d-amino correction
authorCezary Czaplewski <czarek@chem.univ.gda.pl>
Tue, 30 Jan 2018 15:43:33 +0000 (16:43 +0100)
committerCezary Czaplewski <czarek@chem.univ.gda.pl>
Tue, 30 Jan 2018 15:43:33 +0000 (16:43 +0100)
django_simple/todo/jobfiles.py
django_simple/todo/models.py
django_simple/todo/templates/changelog.html
django_simple/todo/templates/details.html
django_simple/todo/templates/edit.html
django_simple/todo/templates/tutorial.html
django_simple/todo/views.py
files/pbs8.csh
files/pbs8_new.csh

index 9284dd3..44cf4a4 100644 (file)
@@ -368,6 +368,13 @@ CUTOFF=7.00000 WCORR4=0.00000"""
            f2.write(line.replace('$temperatures','"'+" ".join(tmp1)+'" '+str(nreplicas)))
       f1.close()
       f2.close()
+
+      f1 = open('file.seq', 'w')
+      for i in range(0,len(seq),80):
+             f1.write(seq[i:i+80]+'\n')
+      f1.write('\n')       
+      f1.close()
+      
       os.chdir('..')  
 
 # write wham & cluster input
index 765e3c6..3c93f63 100644 (file)
@@ -96,6 +96,7 @@ class Task(models.Model):
     created_date = models.DateTimeField(default=timezone.now)
     jobdirname = models.CharField(max_length=100,default='')
     example = models.CharField(max_length=10,default='')
+    damino = models.BooleanField(default=False)
     
     def __str__(self):
         if self.done:
index 9507fc8..dd54ec5 100644 (file)
  Extension and update of the Tutorial.
  </dd>
 
+<dt>30.01.2018</dt>
+ <dd>
+ Update of the Tutorial<br>
+ (M)REMD simulations of proteins with D-aminoacids show only UNRES models,
+ conversion to all-atom models was not correct.
+ </dd>
+
  
 </dl>
               
index a64237e..06ae3db 100644 (file)
@@ -781,6 +781,15 @@ function display_model(numstr){
       }),
       stage.loadFile(
       "/myfiles/download-file/{{task.jobdirname}}/MODEL"+numstr+".pdb",{defaultRepresentation: false}) .then( function( o ){
+ {% if task.damino %}
+        o.addRepresentation( "licorice",{ name: "unres", color:
+               "residueindex",colorScale:["blue","cyan", "green",
+               "yellow","orange", "red"], sele:".CA"
+                       })
+        o.addRepresentation( "licorice",{ name: "unressc", visible: false, color:
+       "residueindex",colorScale:["blue","cyan", "green", "yellow","orange", "red"]
+        , sele:"*" })    
+ {% else %}
         o.addRepresentation( "cartoon",{ name: "unres", color:
                "residueindex",colorScale:["blue","cyan", "green",
                "yellow","orange", "red"]
@@ -788,6 +797,7 @@ function display_model(numstr){
         o.addRepresentation( "line",{ name: "unressc", visible: false, color:
        "residueindex",colorScale:["blue","cyan", "green", "yellow","orange", "red"]
         , sele:"not hydrogen and sidechainAttached" })    
+ {% endif %}
         return o
       })
   ]).then(function (ol1) {
@@ -829,6 +839,15 @@ function display_model(numstr){
   var stage = new NGL.Stage("viewport_M"+numstr);
   stage.setParameters( { backgroundColor: "white",fogFar: 60, fogNear: 45 } );
   stage.loadFile( "/myfiles/download-file/{{task.jobdirname}}/MODEL"+numstr+".pdb",{defaultRepresentation: false}) .then( function( o ){
+ {% if task.damino %}
+        o.addRepresentation( "licorice",{ name: "unres", color:
+               "residueindex",colorScale:["blue","cyan", "green",
+               "yellow","orange", "red"],sele:".CA"
+                       });
+        o.addRepresentation( "licorice",{ name: "unressc", visible: false, color:
+       "residueindex",colorScale:["blue","cyan", "green", "yellow","orange", "red"]
+        , sele:"*" });
+ {% else %}
         o.addRepresentation( "cartoon",{ name: "unres", color:
                "residueindex",colorScale:["blue","cyan", "green",
                "yellow","orange", "red"]
@@ -836,6 +855,7 @@ function display_model(numstr){
         o.addRepresentation( "line",{ name: "unressc", visible: false, color:
        "residueindex",colorScale:["blue","cyan", "green", "yellow","orange", "red"]
         , sele:"not hydrogen and sidechainAttached" });
+ {% endif %}
         o.autoView();
   })             
      var toggleSidechains =document.getElementById("toggleSidechainsM"+numstr );
index f7648b7..3423543 100644 (file)
@@ -79,9 +79,9 @@ Calculation type {{ p_type }}
                <button type="submit" value="example" name="_example" class="save btn btn-default">Load example data</button>           
                {% if p_type == "replica exchange molecular dynamics - advanced options" %}
                <button type="submit" value="example_saxs"
-               name="_example_saxs" class="save btn btn-default">Load example SAXS 1</button>          
+               name="_example_saxs" class="save btn btn-default">Load SAXS 1 example</button>          
                <button type="submit" value="example_saxs1"
-               name="_example_saxs1" class="save btn btn-default">Load example SAXS 2</button>         
+               name="_example_saxs1" class="save btn btn-default">Load SAXS 2 example</button>         
 
                {% endif %}
        </form>
index e271941..6fd03ad 100644 (file)
@@ -8,30 +8,35 @@ In order to submit a job you need to provide its name:</h4>
 <p>
 <img src="static/name.png" style="border:2px solid blueviolet">
 <p>
-UNRES server allows three types of simulations: 
+The following three types of simulations can be run with the UNRES server: 
 local minimization, molecular dynamics and replica exchange molecular dynamics:
 <p>
 <img src="static/types.png" style="border:2px solid blueviolet">
 <p>
-Use "Save & submit" button to start calculations.
+Use the "Save & submit" button to start calculations.
 <img src="static/submit.png" style="border:2px solid blueviolet">
 <p>
-Use "Refresh" button to check the status of simulation:
+Use the "Refresh" button to check the status of a simulation:
 <p>
 <img src="static/refresh.png" style="border:2px solid blueviolet">
 <p>
 until
 <img src="static/done.png" style="border:2px solid blueviolet">
 <p>
-The job status information will start from the "waiting in the queue to
-start", then "running" and finally "postprocessing" before "done".
-Postprocessing do not count to 100%. Autorefresh is active every 30 sec.
-<p>
-Any single job can by accessed later using the adress of the web page
-displayed after job submission: 
+After the server job has been submitted the "waiting in the queue to
+start" message will be displayed in the status field, which will subsequently change
+to "running",  "postprocessing", and "done". While the job has the "running" 
+status, the percentage of job accomplishment is displayed.
+It should be noted that the "100%" accomplishment refers to the production
+phase of a job, after which the postprocessing phase is triggered which also takes
+some time (up to several minutes when making a movie from an MD trajectory is 
+requested) and whose progress is not monitored. The page is autorefreshed every 30 sec.
+<p>
+Any single job can by accessed later using the address of the web page
+displayed after job submission: e.g. 
 http://unres-server.chem.ug.edu.pl/details1/570cf15fc638493893ece1f011ea0182/984/
 <br>
-Registered users can save all their jobs and access them  after login.
+Registered users can save all their jobs and access them later after logging in.
 <p>
 <h4>You can use "Load example data" button before submitting calculations 
 to try examples listed below:</h4>
@@ -45,24 +50,25 @@ to try examples listed below:</h4>
     
     <div class="hiders" style="display:none" >
 <p>
-In this example only the PDB code of selected protein (1BDD) is provided
+In this example only the PDB code of the selected protein (1BDD) is input.
 <p>
 <img src="static/min.png" style="border:2px solid blueviolet">
 
 
-The results of local minimizations show the UNRES model of 1BDD protein, 
-its total UNRES energy and overlap on starting structure:
+After energy minimization is accomplished, the UNRES model of the protein, 
+its total UNRES energy and superposition on starting structure are displayed.
+The interactive <a href="https://github.com/arose/ngl">NGL Viewer</a> is employed, which enables the users to manipulate the structure.
 <p>
 <a href="static/min1.png"><img src="static/min1.png" width="300px"></a>
 <a href="static/min2.png"><img src="static/min2.png" width="300px"></a>
 <p>
 All files generated during calculations 
-can be viewed/download by clicking on "Directory" link.
+can be viewed/download by clicking on the "Directory" link.
 </div>
 </div>
 
 <h4><li>Molecular dynamics of IGG-binding domain from streptococcal protein G
-(PDB code:1IGD) starting from the native structure</h4>
+(PDB code:1IGD) starting from the experimental structure</h4>
 <div>
     <fieldset class="majorpoints">
         <legend class="majorpointslegend"
@@ -73,10 +79,10 @@ can be viewed/download by clicking on "Directory" link.
 <p>
 <img src="static/md.png" style="border:2px solid blueviolet">
 <p>
-The results of md simulations show temperature histogram, 
-potential energy changes, movie generated from trajectory,
-evolution of radius of gyration,
-RMSD, fraction of native contacts, and CA fluctuations:
+After the MD simulation is accomplished temperature histogram, 
+plots of UNRES energy, RMSD from the experimental structure, fraction of native
+contacts, and radius of gyration vs. time, a movie of the trajectory, and 
+CA-atom fluctuations are displayed.
 <p>
 <a href="static/md1.png"><img src="static/md1.png" width="300px"></a>
 <a href="static/md2.png"><img src="static/md2.png" width="300px"></a>
@@ -89,11 +95,11 @@ RMSD, fraction of native contacts, and CA fluctuations:
 <p>
 Evolution of RMSD, fraction of native contacts and 
 comparison of CA fluctuations to Bfactor is presented
-only when reference structure is provided as in this example.
+only when the reference structure is provided, as in this example.
 </div>
 </div>
 <h4><li>
-Replica exchange molecular dynamics of Trp-Cage miniprotein (PDB
+Replica exchange molecular dynamics of the Trp-Cage miniprotein (PDB
 code:1L2Y) starting from the extended chain </h4>  
 <div>
     <fieldset class="majorpoints">
@@ -103,22 +109,32 @@ code:1L2Y) starting from the extended chain </h4>
     <div class="hiders" style="display:none" >
 <img src="static/remd.png" style="border:2px solid blueviolet">
 <p>
-
-Plots of histograms of UNRES energy for each temperature, and
-energy vs temperature are presented.
-The weighted histogram analysis (WHAM) is applied to compute the
+<p>
+After a (M)REMD run is accomplished, a table showing replica-exchange statistics
+(a), plots of histograms of UNRES energy for each temperature (b), a plot of
+energy vs. temperature (c), a temperature-colored scatter plot of energy vs. RMSD
+(d), and the temperature-colored cumulative plot of rmsd vs. step number for all replicas
+(e), and a cumulative plot of walk of the replicas in temperature space (f) are
+displayed. The colors corresponding to replica temperatures are shown in panel (c),
+while different replicas are colored from black to yellow in panel (f).
+<p>
+<a href="static/remd1.png"><img src="static/remd1.png" width="300px"></a>(a)
+<a href="static/remd2.png"><img src="static/remd2.png" width="300px"></a>(b)
+<a href="static/remd3.png"><img src="static/remd3.png" width="300px"></a>(c)
+<a href="static/remd6.png"><img src="static/remd6.png" width="300px"></a>(d)
+<a href="static/remd7.png"><img src="static/remd7.png" width="300px"></a>(e)
+<a href="static/remd8.png"><img src="static/remd8.png" width="300px"></a>(f)
+<p>
+The weighted histogram analysis method (WHAM) is applied to compute the
 probabilities of the obtained conformations to occur at particular
-temperatures, plots of heat capacity and average RMSD as a functions
-of temperature are shown. Replica exchanges are analyzed.
-<p>
-<a href="static/remd1.png"><img src="static/remd1.png" width="300px"></a>
-<a href="static/remd2.png"><img src="static/remd2.png" width="300px"></a>
-<a href="static/remd3.png"><img src="static/remd3.png" width="300px"></a>
-<a href="static/remd4.png"><img src="static/remd4.png" width="300px"></a>
-<a href="static/remd5.png"><img src="static/remd5.png" width="300px"></a>
-<a href="static/remd6.png"><img src="static/remd6.png" width="300px"></a>
-<a href="static/remd7.png"><img src="static/remd7.png" width="300px"></a>
-<a href="static/remd8.png"><img src="static/remd8.png" width="300px"></a>
+temperatures and, consequently, the plots of heat capacity (g) and ensemble-averaged RMSD (h) as a functions
+of temperature; these are shown in the respective graphs. 
+<p>
+<a href="static/remd4.png"><img src="static/remd4.png" width="300px"></a>(g)
+<a href="static/remd5.png"><img src="static/remd5.png" width="300px"></a>(h)
+<p>
+Note! Panels are labeled in this tutorial only and labels are not displayed 
+in the outputs from real server jobs.
 <p>
 Finally cluster analysis is performed to select 5 families of conformations,
 and representative model from each family is converted to all-atom
@@ -134,9 +150,9 @@ PDB files can be downloaded by clicking on the <i>Download</i> button.
 </div>
 </ol>
 <hr>
-Advanced mode allows for changes of more parameters of each type of
-simulation. Separate examples are provided (use <i>Load example data</i> button as
-in Basic mode): 
+Advanced mode enables the user to change more parameters for each type of
+simulation. Separate examples are provided (use the <i>Load example data</i> button as
+in the Basic mode): 
 <ol>
 <h4><li>
 Minimization of the P8MTCP1 disulfide-bonded helical hairpin miniprotein
@@ -148,14 +164,14 @@ Minimization of the P8MTCP1 disulfide-bonded helical hairpin miniprotein
     
     <div class="hiders" style="display:none" >
 <p>
-In this example the PDB code of selected protein (1EI0) and opt-wtfsa-2
+In this example the PDB code of the selected protein (1EI0) and the OPT-WTFSA-2
 force field are selected:  
 <p>
 <img src="static/min_adv.png" style="border:2px solid blueviolet">
 
-The results of local minimizations show the UNRES model of 1EI0 protein
-with two dissulfide bonds as read from pdb, 
-its total UNRES energy and overlap on starting structure:
+After the minimization is accomplished, the UNRES model of 1EI0 protein
+with two disulfide bonds as read from the respective PDB file, 
+its total UNRES energy and superposition on the starting structure are displayed:
 <p>
 <a href="static/min_adv1.png"><img src="static/min_adv1.png" width="300px"></a>
 <a href="static/min_adv2.png"><img src="static/min_adv2.png" width="300px"></a>
@@ -165,18 +181,18 @@ its total UNRES energy and overlap on starting structure:
 </div>
 
 <h4><li>
-Canonical MD simulations of Trp-Cage miniprotein (PDB code:1L2Y) 
-starting from extended chain.</h4>
+Canonical MD simulations of the Trp-Cage miniprotein (PDB code:1L2Y) 
+starting from the extended chain.</h4>
 <div>
     <fieldset class="majorpoints">
         <legend class="majorpointslegend"
             style="background-color:#d3d3d3;cursor: pointer;">Show</legend>
     
     <div class="hiders" style="display:none" >
-In this example the PDB code of selected protein (1L2Y) and 
-option to write trajectory as pdb are selected, saving trajectory 
-as pdb allows its 3D visualisation in NGL Viewer in addtion to
-movie which is always generated from the MD trajectory 
+In this example the PDB code of the selected protein (1L2Y) and 
+option to write trajectory as a PDB file are selected, saving trajectory 
+as a PDB file allows its 3D visualisation by NGL Viewer in addition to
+the movie which is always generated from the MD trajectory 
 <p>
 <img src="static/md_adv.png" style="border:2px solid blueviolet">
 <br>
@@ -184,11 +200,9 @@ movie which is always generated from the MD trajectory
 <br>
 <img src="static/md_adv0.png" style="border:2px solid blueviolet">
 <p>
-The results of md simulations show temperature histogram, 
-potential energy changes, movie generated from trajectory,
-trajectory visualization in NGL Viewer,
-evolution of radius of gyration,
-RMSD, fraction of native contacts:
+After the MD simulation is accomplished temperature histogram, 
+plots of UNRES energy, RMSD from the experimental structure, fraction of native
+contacts, and radius of gyration vs. time, and a movie of the trajectory are displayed.
 <p>
 <a href="static/md_adv1.png"><img src="static/md_adv1.png" width="300px"></a>
 <a href="static/md_adv2.png"><img src="static/md_adv2.png" width="300px"></a>
@@ -198,8 +212,6 @@ RMSD, fraction of native contacts:
 <a href="static/md_adv6.png"><img src="static/md_adv6.png" width="300px"></a>
 <a href="static/md_adv7.png"><img src="static/md_adv7.png" width="300px"></a>
 <p>
-Evolution of RMSD and fraction of native contacts 
-is shown only when reference structure is provided as in this example.
 The CA fluctuations are not analyzed for start from the extended structure
 because of large scale conformational changes during folding.
 
@@ -216,14 +228,15 @@ MREMD simulations of 5G3Q:B (CASP12 target T0882)</h4>
             style="background-color:#d3d3d3;cursor: pointer;">Show</legend>
     
     <div class="hiders" style="display:none" >
-This simulation starts from extended chain, 2000000 steps for each of 16
+This simulation starts from the extended chain, 2000000 steps for each of 16
 replicas (8 temperatures with multiplexing 2)
-are run and secondary structure restraints predicted by PSIPRED and 
+are run and secondary-structure restraints predicted by PSIPRED and 
 obtained at the CASP12 time when the experimental 5G3Q structure
-was not included in the PDB database are used. Berendsen thermostat
-and clustering temperature 290K are selected. 
+was not included in the PDB database are used. The Berendsen thermostat
+was employed and the resulting conformational ensemble was clustered 
+at the temperature of 290K. 
 This is a full-blown structure-prediction run, but shorter
-in comparison to UNRES runs during CASP.
+compared to the UNRES runs during the CASP exercises.
 
 <p>
 <img src="static/remd_adv.png" style="border:2px solid blueviolet">
@@ -242,29 +255,22 @@ in comparison to UNRES runs during CASP.
 <p>
 
 
-Plots of histograms of UNRES energy for each temperature, and
-energy vs temperature are presented.
-The weighted histogram analysis (WHAM) is applied to compute the
-probabilities of the obtained conformations to occur at particular
-temperatures, plots of heat capacity and average RMSD as a functions
-of temperature are shown. Replica exchanges are analyzed.
+The output items are the same as those discussed with the Basic mode Example 3.
 <p>
 <a href="static/remd_adv1.png"><img src="static/remd_adv1.png" width="300px"></a>
 <a href="static/remd_adv2.png"><img src="static/remd_adv2.png" width="300px"></a>
 <a href="static/remd_adv3.png"><img src="static/remd_adv3.png" width="300px"></a>
-<a href="static/remd_adv4.png"><img src="static/remd_adv4.png" width="300px"></a>
-<a href="static/remd_adv5.png"><img src="static/remd_adv5.png" width="300px"></a>
 <a href="static/remd_adv6.png"><img src="static/remd_adv6.png" width="300px"></a>
 <a href="static/remd_adv7.png"><img src="static/remd_adv7.png" width="300px"></a>
 <a href="static/remd_adv8.png"><img src="static/remd_adv8.png" width="300px"></a>
+<a href="static/remd_adv4.png"><img src="static/remd_adv4.png" width="300px"></a>
+<a href="static/remd_adv5.png"><img src="static/remd_adv5.png" width="300px"></a>
 <p>
-Finally cluster analysis is performed to select 5 families of conformations,
-and representative model from each family is converted to all-atom
-and refined. Models are shown using NGL Viewer.
+The final models are shown using NGL Viewer.
 PDB files can be downloaded by clicking on the <i>Download</i> button.
-RMSD of model 1 
+The RMSD of model 1 
 from the experimental structure is 3.6 A and 
-the GDT TS is 58.8 %
+the GDT_TS is 58.8 %
 
 <p>
 <a href="static/remd_adv9.png"><img src="static/remd_adv9.png" width="300px"></a>
@@ -282,7 +288,7 @@ Distance distribution restrained (simulated SAXS data)
 replica exchange molecular dynamics of 
 Bacteriocin CbnXY (PDB code:5UJQ) starting from the extended chain.
 <br>
-(Use <i>Load example SAXS 1</i> button in advanced mode)
+(Use <i>Load SAXS 1 example</i> button in advanced mode)
 </h4>  
 <div>
     <fieldset class="majorpoints">
@@ -297,31 +303,24 @@ Bacteriocin CbnXY (PDB code:5UJQ) starting from the extended chain.
 <img src="static/saxs2.png" style="border:2px solid blueviolet">
 <p>
 
-Plots of histograms of UNRES energy for each temperature, and
-energy vs temperature are presented.
-The weighted histogram analysis (WHAM) is applied to compute the
-probabilities of the obtained conformations to occur at particular
-temperatures, plots of heat capacity and average RMSD as a functions
-of temperature are shown. Replica exchanges are analyzed.
+The output items are the same as those displayed with the Basic mode example 3.
+Additionaly the input reference distance distribution and 
+distance distributions for 5 final models are plotted.
 <p>
 <a href="static/saxs3.png"><img src="static/saxs3.png" width="300px"></a>
 <a href="static/sasx4.png"><img src="static/saxs4.png" width="300px"></a>
 <a href="static/saxs5.png"><img src="static/saxs5.png" width="300px"></a>
-<a href="static/saxs6.png"><img src="static/saxs6.png" width="300px"></a>
-<a href="static/saxs7.png"><img src="static/saxs7.png" width="300px"></a>
 <a href="static/saxs8.png"><img src="static/saxs8.png" width="300px"></a>
 <a href="static/saxs9.png"><img src="static/saxs9.png" width="300px"></a>
 <a href="static/saxs10.png"><img src="static/saxs10.png" width="300px"></a>
+<a href="static/saxs6.png"><img src="static/saxs6.png" width="300px"></a>
+<a href="static/saxs7.png"><img src="static/saxs7.png" width="300px"></a>
 <p>
-Finally cluster analysis is performed to select 5 families of conformations,
-and representative model from each family is converted to all-atom
-and refined. 
+Finally models.
 <p>
 <a href="static/saxs12.png"><img src="static/saxs12.png" width="300px"></a>
-
 <p>
-Additionaly the input reference distance distribution and 
-distance distributions for 5 final models are plotted.
+Comparison of the input distance distributions with those calculated from the 5 models.
 <p>
 <a href="static/saxs11.png"><img src="static/saxs11.png" width="300px"></a>
 
@@ -333,7 +332,7 @@ Distance distribution restrained (real SAXS data)
 replica exchange molecular dynamics of 
 the central portion of Factor H (PDB code:2KMS).
 <br>
-(Use <i>Load example SAXS 2</i> button in advanced mode)
+(Use <i>Load SAXS 2 example</i> button in advanced mode)
 </h4>  
 <div>
     <fieldset class="majorpoints">
@@ -342,9 +341,9 @@ the central portion of Factor H (PDB code:2KMS).
     
     <div class="hiders" style="display:none" >
 REMD simulation starts from the experimental NMR structure. Secondary
-structure restraints are used. Only short (50 energy evaluations)
+structure restraints are used. Only a short (50 energy evaluations)
 initial minimization is requested.
-The SASDA25 entry of the SASBDB is used for SAXS distance distribution.
+The distance distribution has been downloaded from the SASBDB database (the SASDA25 entry).
 
 <img src="static/saxs2_1.png" style="border:2px solid blueviolet">
 <br>
@@ -358,25 +357,18 @@ The SASDA25 entry of the SASBDB is used for SAXS distance distribution.
 
 <p>
 
-Plots of histograms of UNRES energy for each temperature, and
-energy vs temperature are presented.
-The weighted histogram analysis (WHAM) is applied to compute the
-probabilities of the obtained conformations to occur at particular
-temperatures, plots of heat capacity and average RMSD as a functions
-of temperature are shown. Replica exchanges are analyzed.
+The input items are the same as those in the Advanced mode Example 2.
 <p>
 <a href="static/sasx2_4.png"><img src="static/saxs2_4.png" width="300px"></a>
 <a href="static/saxs2_5.png"><img src="static/saxs2_5.png" width="300px"></a>
 <a href="static/saxs2_6.png"><img src="static/saxs2_6.png" width="300px"></a>
-<a href="static/saxs2_7.png"><img src="static/saxs2_7.png" width="300px"></a>
-<a href="static/saxs2_8.png"><img src="static/saxs2_8.png" width="300px"></a>
 <a href="static/saxs2_9.png"><img src="static/saxs2_9.png" width="300px"></a>
 <a href="static/saxs2_10.png"><img src="static/saxs2_10.png" width="300px"></a>
 <a href="static/saxs2_11.png"><img src="static/saxs2_11.png" width="300px"></a>
+<a href="static/saxs2_7.png"><img src="static/saxs2_7.png" width="300px"></a>
+<a href="static/saxs2_8.png"><img src="static/saxs2_8.png" width="300px"></a>
 <p>
-Finally cluster analysis is performed to select 5 families of conformations,
-and representative model from each family is converted to all-atom
-and refined. 
+Final models.
 In the UNRES calculated
 structures, the domains are at a larger angle than
 those in the experimental NMR structure.
@@ -385,12 +377,7 @@ those in the experimental NMR structure.
 <a href="static/saxs2_13.png"><img src="static/saxs2_13.png" width="300px"></a>
 
 <p>
-Additionaly the input reference distance distribution and 
-distance distributions for 5 final models are plotted.
-The distance distribution calculated from the
-experimental NMR structure does not fully overlap with that from
-SAXS; in particular, the plot calculated from the 2KMS
-structure decays to 0 quicker than that from SAXS.
+Comparison of the input distance distributions with those calculated from the 5 models.
 Much better agreement, in
 particular in the long-distance part is obtained
 for structres from the UNRES simulation.
@@ -417,4 +404,4 @@ $('.majorpointslegend').click(function(){
 </script>
 
 
-{% endblock %}
\ No newline at end of file
+{% endblock %}
index 118de52..45741f2 100644 (file)
@@ -316,6 +316,7 @@ def add_md(request,task_id):
 
              if any(c.islower() for c in seq):
               task.md_respa=False
+              task.damino=True
 
              task.ready=True
              
@@ -387,6 +388,7 @@ def add_md_a(request,task_id):
              task.md_respa=form.cleaned_data["md_respa"]             
              if any(c.islower() for c in seq):
               task.md_respa=False
+              task.damino=True
              task.boxx=form.cleaned_data["boxx"]
              task.boxy=form.cleaned_data["boxy"]             
              task.boxz=form.cleaned_data["boxz"]             
@@ -450,6 +452,7 @@ def add_remd(request,task_id):
 
              if any(c.islower() for c in seq):
               task.md_respa=False
+              task.damino=True
 
 
              task.ready=True
@@ -675,6 +678,7 @@ def add_remd_a(request,task_id):
              task.md_respa=form.cleaned_data["md_respa"]
              if any(c.islower() for c in seq):
               task.md_respa=False
+              task.damino=True
              
 #             task.remd_traj1file=form.cleaned_data["remd_traj1file"]
 #             task.remd_rest1file=form.cleaned_data["remd_rest1file"]
index 7f8f0be..650f7f9 100755 (executable)
@@ -67,6 +67,21 @@ setenv CLUSTER_WHAM_BIN /users2/czarek/UNRES/git3_tmp/build/bin/cluster_wham-mul
 
 $MPIRUN -machinefile $PBS_NODEFILE -np 2 $CLUSTER_WHAM_BIN
 
+if `awk '{cap=toupper($0); if (cap!=$0) {print 1;exit}}' file.seq` then
+#no allatom conversion/refinementkx for proteins with D-aminoacids
+ ln -s file_wham_T*K_ave.pdb ave
+ sed -n '/ENERGY/,/ENDMDL/p' ave >tmp.pdb
+ setenv numstruc `grep ENERGY tmp.pdb|wc -l`
+ setenv allline `cat tmp.pdb|wc -l`
+ setenv onestruc `echo $allline "/" $numstruc| bc -l|sed 's/\.0*//'`
+ foreach i (2 4 6 8 10)
+   setenv headval `echo $i "*" $onestruc|bc -l`
+   setenv halfi `echo $i/2|bc`
+   head -n $headval tmp.pdb | tail -n $onestruc >MODEL${halfi}.pdb
+ end
+ rm ave tmp.pdb
+else
+# proteins with L-aminoacids only
 ln -s file_wham_T*K_ave.pdb ave
 ../files/cluster2allatom.sh ave
 rm ave
@@ -75,6 +90,7 @@ rm ave
 foreach i (`seq 1 5`)
 ../files/refine.sh MODEL$i.pdb $ssbond
 end
+endif
 
 /users2/local/pymol_1.6/pymol -c ../files/model.pml
 
index 6d9d029..ec056c1 100755 (executable)
@@ -67,6 +67,21 @@ setenv CLUSTER_WHAM_BIN /users2/czarek/UNRES/git3_tmp/build/bin/cluster_wham-mul
 
 $MPIRUN -machinefile $PBS_NODEFILE -np 2 $CLUSTER_WHAM_BIN
 
+if `awk '{cap=toupper($0); if (cap!=$0) {print 1;exit}}' file.seq` then
+#no allatom conversion/refinementkx for proteins with D-aminoacids
+ ln -s file_wham_T*K_ave.pdb ave
+ sed -n '/ENERGY/,/ENDMDL/p' ave >tmp.pdb
+ setenv numstruc `grep ENERGY tmp.pdb|wc -l`
+ setenv allline `cat tmp.pdb|wc -l`
+ setenv onestruc `echo $allline "/" $numstruc| bc -l|sed 's/\.0*//'`
+ foreach i (2 4 6 8 10)
+   setenv headval `echo $i "*" $onestruc|bc -l`
+   setenv halfi `echo $i/2|bc`
+   head -n $headval tmp.pdb | tail -n $onestruc >MODEL${halfi}.pdb
+ end
+ rm ave tmp.pdb
+else
+# proteins with L-aminoacids only
 ln -s file_wham_T*K_ave.pdb ave
 ../files/cluster2allatom.sh ave
 rm ave
@@ -75,6 +90,7 @@ rm ave
 foreach i (`seq 1 5`)
 ../files/refine.sh MODEL$i.pdb $ssbond
 end
+endif
 
 /users2/local/pymol_1.6/pymol -c ../files/model.pml