curvesplus package

Submodules

curvesplus.biobb_curves module

Module containing the Curves class and the command line interface.

class curvesplus.biobb_curves.Curves(input_struc_path, output_lis_path, output_cda_path, output_zip_path=None, input_top_path=None, properties=None, **kwargs)[source]

Bases: biobb_common.generic.biobb_object.BiobbObject

biobb_dna Curves
Wrapper for the Cur+ executable that is part of the Curves+ software suite.
Parameters
  • input_struc_path (str) – Trajectory or PDB input file. File type: input. Sample file. Accepted formats: trj (edam:format_3910), pdb (edam:format_1476).

  • input_top_path (str) (Optional) –

    Topology file, needed along with .trj file (optional). File type: input. Sample file. Accepted formats: top (edam:format_3881).

  • output_cda_path (str) –

    Filename for Curves+ output .cda file. File type: output. Sample file. Accepted formats: cda (edam:format_2330).

  • output_lis_path (str) –

    Filename for Curves+ output .lis file. File type: output. Sample file. Accepted formats: lis (edam:format_2330).

  • output_zip_path (str) (Optional) – Filename for .zip files containing Curves+ output that is not .cda or .lis files. File type: output. Accepted formats: zip (edam:format_3987).

  • properties (dict) –

    • s1range (str) - (None) Range of first strand. Must be specified in the form “start:end”.

    • s2range (str) - (None) Range of second strand. Must be specified in the form “start:end”.

    • stdlib_path (str) - (‘standard’) Path to Curves’ standard library files for nucleotides. If not specified will look for ‘standard’ files in current directory.

    • itst (int) - (0) Iteration start index.

    • itnd (int) - (0) Iteration end index.

    • itdel (int) - (1) Iteration delimiter.

    • ions (bool) - (False) If True, helicoidal analysis of ions (or solvent molecules) around solute is carried out.

    • test (bool) - (False) If True, provide addition output in .lis file on fitting and axis generation.

    • line (bool) - (False) if True, find the best linear helical axis.

    • fit (bool) - (False) if True, fit a standard bases to the input coordinates (important for MD snapshots to avoid base distortions leading to noisy helical parameters).

    • axfrm (bool) - (False) if True, generates closely spaced helical axis frames as input for Canal and Canion.

    • curves_exec (str) - (Cur+) Path to Curves+ executable, otherwise the program wil look for Cur+ executable in the binaries folder.

Examples

This is a use example of how to use the building block from Python:

from biobb_dna.curvesplus.biobb_curves import biobb_curves
prop = {
    's1range': '1:12',
    's2range': '24:13',
}
biobb_curves(
    input_struc_path='/path/to/structure/file.trj',
    input_top_path='/path/to/topology/file.top',
    output_cda_path='/path/to/output/file.cda',
    output_lis_path='/path/to/output/file.lis',
    properties=prop)
Info:
launch() int[source]

Execute the Curves object.

curvesplus.biobb_curves.biobb_curves(input_struc_path: str, output_lis_path: str, output_cda_path: str, input_top_path: Optional[str] = None, output_zip_path: Optional[str] = None, properties: Optional[dict] = None, **kwargs) int[source]

Create Curves class and execute the launch() method.

curvesplus.biobb_curves.main()[source]

Command line execution of this building block. Please check the command line documentation.

curvesplus.biobb_canal module

Module containing the Canal class and the command line interface.

class curvesplus.biobb_canal.Canal(input_cda_file, input_lis_file=None, output_zip_path=None, properties=None, **kwargs)[source]

Bases: biobb_common.generic.biobb_object.BiobbObject

biobb_dna Canal
Wrapper for the Canal executable that is part of the Curves+ software suite.
Parameters
  • input_cda_file (str) –

    Input cda file, from Cur+ output. File type: input. Sample file. Accepted formats: cda (edam:format_2330).

  • input_lis_file (str) (Optional) – Input lis file, from Cur+ output. File type: input. Accepted formats: lis (edam:format_2330).

  • output_zip_path (str) –

    zip filename for output files. File type: output. Sample file. Accepted formats: zip (edam:format_3987).

  • properties (dic) –

    • bases (str) - (None) sequence of bases to be searched for in the I/P data (default is blank, meaning no specified sequence).

    • itst (int) - (0) Iteration start index.

    • itnd (int) - (0) Iteration end index.

    • itdel (int) - (1) Iteration delimiter.

    • lev1 (int) - (0) Lower base level limit (i.e. base pairs) used for analysis.

    • lev2 (int) - (0) Upper base level limit used for analysis. If lev1 > 0 and lev2 = 0, lev2 is set to lev1 (i.e. analyze lev1 only). If lev1=lev2=0, lev1 is set to 1 and lev2 is set to the length of the oligmer (i.e. analyze all levels).

    • nastr (str) - (‘NA’) character string used to indicate missing data in .ser files.

    • cormin (float) - (0.6) minimal absolute value for printing linear correlation coefficients between pairs of analyzed variables.

    • series (str) - (False) if True then output spatial or time series data. Only possible for the analysis of single structures or single trajectories.

    • histo (str) - (False) if True then output histogram data.

    • corr (str) - (False) if True than output linear correlation coefficients between all variables.

    • sequence (str) - (Optional) sequence of the first strand of the corresponding DNA fragment, for each .cda file. If not given it will be parsed from .lis file.

    • canal_exec (str) - (‘Canal’) Path to Canal executable, otherwise the program wil look for Canal executable in the binaries folder.

Examples

This is a use example of how to use the building block from Python:

from biobb_dna.curvesplus.biobb_canal import biobb_canal
prop = {
    'series': '.t.',
    'histo': '.t.',
    'sequence': 'CGCGAATTCGCG'
}
biobb_canal(
    input_cda_file='/path/to/curves/output.cda',
    output_zip_path='/path/to/output.zip',
    properties=prop)
Info:
launch() int[source]

Execute the Canal object.

curvesplus.biobb_canal.biobb_canal(input_cda_file: str, output_zip_path: str, input_lis_file: Optional[str] = None, properties: Optional[dict] = None, **kwargs) int[source]

Create Canal class and execute the launch() method.

curvesplus.biobb_canal.main()[source]

Command line execution of this building block. Please check the command line documentation.

curvesplus.biobb_canion module

Module containing the Canion class and the command line interface.

class curvesplus.biobb_canion.Canion(input_cdi_path, input_afr_path, input_avg_struc_path, output_zip_path=None, properties=None, **kwargs)[source]

Bases: biobb_common.generic.biobb_object.BiobbObject

biobb_dna Canion
Wrapper for the Canion executable that is part of the Curves+ software suite.
Parameters
  • input_cdi_path (str) –

    Trajectory input file. File type: input. Sample file. Accepted formats: cdi (edam:format_2330).

  • input_afr_path (str) –

    Helical axis frames corresponding to the input conformation to be analyzed. File type: input. Sample file. Accepted formats: afr (edam:format_2330).

  • input_avg_struc_path (str) –

    Average DNA conformation. File type: input. Sample file. Accepted formats: pdb (edam:format_1476).

  • output_zip_path (str) (Optional) –

    Filename for .zip files containing Canion output files. File type: output. Sample file. Accepted formats: zip (edam:format_3987).

  • properties (dict) –

    • bases (str) - (None) Sequence of bases to be analyzed (default is blank, meaning no specified sequence).

    • type (str) - (‘*’) Ions (or atoms) to be analyzed. Options are ‘Na+’, ‘K’, ‘K+’, ‘Cl’, ‘Cl-’, ‘CL’, ‘P’, ‘C1*’, ‘NH1’, ‘NH2’, ‘NZ’, ‘1’ for all cations, ‘-1’ for all anions, ‘0’ for neutral species or ‘*’ for all available data.

    • dlow (float) - (0) Select starting segment of the oglimer to analyze. If both dhig and dlow are 0, entire oglimer is analyzed.

    • dhig (float) - (0) Select ending segment of the oglimer to analyze, being the maximum value the total number of base pairs in the oligomer. If both dhig and dlow are 0, entire oglimer is analyzed.

    • rlow (float) - (0) Minimal distances from the helical axis taken into account in the analysis.

    • rhig (float) - (0) Maximal distances from the helical axis taken into account in the analysis.

    • alow (float) - (0) Minimal angle range to analyze.

    • ahig (float) - (360) Maximal angle range to analyze.

    • itst (int) - (None) Number of first snapshot to be analyzed.

    • itnd (int) - (None) Number of last snapshot to be analyzed.

    • itdel (int) - (None) Spacing between analyzed snapshots.

    • rmsf (bool) - (False) If set to True uses the combination of the helical ion parameters and an average helical axis to map the ions into Cartesian space and then calculates their average position (pdb output) and their root mean square fluctuation values (rmsf output). A single pass rmsf algorithm to make this calculation possible with a single read of the trajectory file. This option is generally used for solute atoms and not for solvent molecules or ions.

    • circ (bool) - (False) If set to True, minicircles are analyzed.

    • canion_exec (str) - (Canion) Path to Canion executable, otherwise the program wil look for Canion executable in the binaries folder.

Examples

This is a use example of how to use the building block from Python:

from biobb_dna.curvesplus.biobb_canion import biobb_canion
prop = {
    'type': 'K+',
    'bases': 'G'
}
biobb_canion(
    input_cdi_path='/path/to/input.cdi',
    input_afr_path='/path/to/input.afr',
    input_avg_struc_path='/path/to/input.pdb',
    output_zip_path='/path/to/output.zip',
    properties=prop)
Info:
launch() int[source]

Execute the Canion object.

curvesplus.biobb_canion.biobb_canion(input_cdi_path: str, input_afr_path: str, input_avg_struc_path: str, output_zip_path: Optional[str] = None, properties: Optional[dict] = None, **kwargs) int[source]

Create Canion class and execute the launch() method.

curvesplus.biobb_canion.main()[source]

Command line execution of this building block. Please check the command line documentation.