added:
Catagories of scripts
<a href="#MarkScriptsToPlotData">Scripts to plot data</a>
changed:
-Scripts to Plot Data:
<a name="MarkScriptsToPlotData">Scripts to Plot Data:</a>
Return to Main CALUA page
Catagories of scripts
Scripts to plot data
Scripts to modify repository globally:
- ResidueTypes:
- CorrectAtomMagnitudes.lua : Adjusts all methyl group atom magnitudes to 3 (this fixes an error in the standard template Start1.2.cara)
- LoadBmrbStats.lua : Updates ResidueType Ave,Dev values from a BMRB statistics file. See header of script for details.
- UPDATE Aug.29 16: LoadBmrbStats.lua BMRB changed their datestamp format from MO-DAY-YEAR to MO/DAY/YEAR, updated script to read this. Also added check for IUPAC nomenclature and script now stores the value of SDmult as Attribute of each modified ResidueType?. See header for details.
Scripts to write Analysis Reports to Script Window or to File:
- AssignmentReport.lua : Writes out a report of chemical shift assignments to either terminal window or to a file. Identifies unusual shifts. It is possible to filter for only missing assignments (H atoms only or also heteroatoms).
- UPDATE 20.May 05: Fixed a bug which occured if NOT all residues were assigned. Error message was: "line 83: attempt to index field
Sys
(a nil value)".
- UPDATE 4.6.06: Includes a nice user-interface which remembers settings. Gives average offsets from random coil values. Allows user to set minimum SD to report outlier shifts. This is the version included in the latest template Start1.2.cara.
- UPDATE 19.6.07: Fixed a bug which occured when there is only one project Error message was "No projects found in repository". Added a summary of pseudoatom conflicts (assignment of both the real atom HB2 and pseudoatom HB is reported in a warning).
- CalcMassOfLabeledMolecule.lua : Calculates the molecular weight of a sequence in a project assuming the isotope-enrichment and percent d2o entered by the user. 2H, 13C, and 15N are available by default. Other nuclei can be accessed by commenting in an option at the top of the script.
- UPDATED It now includes a checkbox to calculate the expected mass for mass-spectrometry and can also calculate the mass in solution for a given pH.
- ChemShiftDeviationsFile.lua : Writes a file containing Chemical Shift Deviations from the average values stored in the repository library. See also ChemShiftDeviationsPlot.lua
- NEW ChemShiftMapFromTwoPeakLists.lua : Writes combined chemical shift changes to script terminal window using two peaklists as input. The two peaklists should have the same peak Ids for same peaks and ideally peaks should be numbered by residue (You can use PeakListNumberedByResidue.lua to generate them). Combined chemical shift differences are calculated using the difference in chemical shift of the peaks with same Id in the two peaklists. Script also outputs atom selections in MOLMOL format for those residues having shifts larger that StdevThreshold x standard deviation of the dataset as well as a MOLMOL selection corresponding to all peaks with color = PeakColor in PeakList1.
- CompareProjects.lua : Compares the assigned resonances and shift values in two projects - reports assignments unique to one project and shift values that differ signficantly. See also CopyAssignmentsToProject.lua
- ExportAllSpins.lua : Writes a file containing an atom list in XEASY format. All spins in the selected project are included (also CA-1 etc)
- GetLineWidthsFromPeakLists.lua : Writes out a file with one line for each peak in a peaklist reporting peaknumber and linewidth at half height along a one dimension of a 2D spectrum (see header for details)
- HetNoe.lua : Uses a batchlist consisting of two spectra (noe & ref) to calculate the heteronuclear nOe (ratio of peaks) and stores the result in the Peak attribute HetNoe. It also writes the results to a file in format suitable for excel. See the beginning of the script for instructions.
- ReportSpinsWithColorNinPeakList.lua : Writes out a text file containing spins which occur more than N times in a peak with color C in peaklist P. (User enters N,C & P).
- ShiftDiffBetweenTwoProjects.lua : Writes out a report of the combined chemical shift differences between two projects from each residue. Optionally the alias shifts of two spectra can be used. This can be useful for chemical shift mapping.
- SideChainDeviations.lua : Writes out a report of average deviation of chemical shifts from each residue normalized to the standard deviation of the database versus residue number (H atoms only or also heteroatoms).
- WriteShiftsInColumns.lua : Writes out shifts for selected atoms into columns, one line per residue. The user can edit a table at the beginning of the script to define which atoms shifts are written out. Currently: "H", "HA", "CA", "CB".
Scripts to Export Data:
- Analysis Reports to Script Window or to File:
- AssignmentReport.lua : Writes out a report of chemical shift assignments to either terminal window or to a file. Identifies unusual shifts. It is possible to filter for only missing assignments (H atoms only or also heteroatoms).
- UPDATE 20.May 05: Fixed a bug which occured if NOT all residues were assigned. Error message was: "line 83: attempt to index field
Sys
(a nil value)".
- UPDATE 4.6.06: Includes a nice user-interface which remembers settings. Gives average offsets from random coil values. Allows user to set minimum SD to report outlier shifts. This is the version included in the latest template Start1.2.cara.
- CalcMassOfLabeledMolecule.lua : Calculates the molecular weight of a sequence in a project assuming the isotope-enrichment and percent d2o entered by the user. 2H, 13C, and 15N are available by default. Other nuclei can be accessed by commenting in an option at the top of the script.
- UPDATED 8.11.07: It now includes a checkbox to calculate the expected mass for mass-spectrometry and can also calculate the mass in solution for a given pH.
- ChemShiftDeviationsFile.lua : Writes a file containing Chemical Shift Deviations from the average values stored in the repository library. See also ChemShiftDeviationsPlot.lua
- ExportAllSpins.lua : Writes a file containing an atom list in XEASY format. All spins in the selected project are included (also CA-1 etc)
- ReportSpinsWithColorNinPeakList.lua : Writes out a text file containing spins which occur more than N times in a peak with color C in peaklist P. (User enters N,C & P).
- SideChainDeviations.lua : Writes out a report of average deviation of chemical shifts from each residue normalized to the standard deviation of the database versus residue number (H atoms only or also heteroatoms).
- WriteAssignments.lua : Writes out assignments for BMRB, DYANA, CYANA or XEASY (see "Export File for a Program" for details)
- WriteShiftsInColumns.lua : Writes out shifts for selected atoms into columns, one line per residue. The user can edit a table at the beginning of the script to define which atoms shifts are written out. Currently: "H", "HA", "CA", "CB".
- Export File for a Program:
- ATNOS/CANDID & DYANA or CYANA2.0
- ExportSpinLinksToUpls.lua : Writes a UPL file in CYANA format containing one UPL for each spinlink connecting the atoms corresponding to the linked spins. The upper distance limits are all set to 5.0 Angstroms. This script can be used together with UplsToSpinLinks.lua in order to transfer spinlinks from one project to another. Thanks to D. Frueh for this contribution.
- WriteAssignments.lua : Writes out assignments in a nomenclature consistent for structure calculations using ATNOS/CANDID together with DYANA or CYANA2.0 (seq and prot files). Can also write out assignments for submission to BMRB. This is the successor to WriteAssignmentsInFormatX.lua. The main new feature is that WriteAssignments.lua can write out alias shifts for a selected spectrum.
- UPDATED 30.10.07: BMRB format now writes out one line for each atom with a degenerate shift.
- WriteAssignmentsInFormatX.lua : Writes out assignments in a nomenclature consistent for structure calculations using ATNOS/CANDID together with DYANA or CYANA2.0 (seq and prot files). Can also write out assignments for submission to BMRB.
- EXCEL
- WriteBatchTable.lua : Write out a table from a batch integration into a text file (for EXCEL).
- HetNoe.lua: Uses a batchlist consisting of two spectra (noe & ref) to calculate the heteronuclear nOe (ratio of peaks) and stores the result in the Peak attribute HetNoe. It also writes the results to a file in format suitable for excel. See the beginning of the script for instructions. See the beginning of the script for instructions.
- UPDATE: Improved User interface, output format and introduced the option to write out the spectrum-dependent parameter (e.g. MixTime, ExchTime)
- WriteTitrationCurveDataToFile.lua : Writes chemical shift values of spins with a user-selected label for a series of spectra into a spreadsheet file. (shifts are taken from the alias values in each spectrum). See header for details.
- UPDATE 2.Jan.05: Fixed an error in line "t.Line = LineTable..." and added dialog for Project selection and entry of ProtonListName
- MARS
- MOLMOL
- ShowNoesInMolmol.lua : Writes out a molmol macro which highlights a selected target atom in cyan and all candidate atoms within a tolerance range of a user-determined chemical shift. Useful for NOE assignment during structure refinement. See header for details
- PACES
- ExportToPaces.lua : Writes out a comma-separated list of HN,N,CA,CB chemical shifts for import into Excel and use with PACES.
- UPDATE 20. May 05: Adjusted formatting to more closely adhere to PACES. Included CO, CO-1, HA, HA-1.
- SPARKY
- WriteSparkyPeakList.lua : Writes out a peaklist in sparky format from a Peaklist in a CARA project.
- UPDATED (Ver. 2) Script now writes out sparky format for 2D, 3D and 4D peaklists.
- TALOS+
- XEASY
- WriteProtonList.lua
: Writes out an xeasy-formatted proton list for a project (no unassigned "999.000" atoms included)
- PeakListNumberedByResidue.lua : Write out a peaklist where the peaknumber corresponds to the residue number for a specific atom in each dimension. E.g. you want to write out a peaklist for an HSQC15N where each residue gives rise to a peak (HN/N) whose peaknumber corresponds to the residue number. This is useful for plotting the results of NH dynamics or exchange vs. sequence.
- X-PLOR
- ExportToXplor.lua : Converts spinlinks from a project into upper distance limits with the maximum distance (5.5A) Thanks to S. Ilin for this contribution.
- PeakListToXplorConstraints.lua : Generates a XPLOR-NIH constraints file from a PeakList? in a project. The intensities are divided into four catagories and converted into distance constraints which are written to an external file. Thanks to C. Hilty for this contribution.
Scripts to convert external files:
- OneLetterFileToSeqFile.lua : Converts a file containing an amino acid sequence in one letter code to an XEASY seq file for easy read in to CARA. It ignores case and all characters which do not code for one of the 20 standard amino acids. Copy-paste from ExPASy or NCBI website directly into a .aa file and convert to .seq file!
- UPDATE 8.Apr.05: Modified so that last line of seq file is not empty.
Scripts to modify Projects (always back-up Repository before running these scripts - there is no undo!)