![]() |
Rosetta
3.8
|
Application-level code for simple_cycpep_predict application. More...
#include <SimpleCycpepPredictApplication.hh>
Public Member Functions | |
SimpleCycpepPredictApplication (bool const allow_file_read=true) | |
Constructor. More... | |
~SimpleCycpepPredictApplication () override | |
Explicit virtual destructor. More... | |
SimpleCycpepPredictApplication (SimpleCycpepPredictApplication const &src) | |
Explicit copy constructor. More... | |
void | initialize_from_options () |
Initialize the application. More... | |
void | set_scorefxn (core::scoring::ScoreFunctionCOP sfxn_in) |
Sets the default scorefunction to use. More... | |
void | set_native (core::pose::PoseCOP native) |
Allows external code to provide a native, so that the SimpleCycpepPredictApplication doesn't have to read directly from disk. More... | |
void | set_sequence (std::string const &seq) |
Allows external code to provide a sequence, so that the SimpleCycpepPredictApplication doesn't have to read directly from disk. More... | |
void | set_allowed_residues_by_position (std::map< core::Size, utility::vector1< std::string > > const &allowed_canonicals, std::map< core::Size, utility::vector1< std::string > > const &allowed_noncanonicals) |
Allows external code to set the allowed residues by position, so that this needn't be read directly from disk. More... | |
void | set_silentstructure_outputlist (utility::vector1< core::io::silent::SilentStructOP > *silentlist, utility::vector1< SimpleCycpepPredictApplication_MPI_JobResultsSummaryOP > *summarylist) |
Allows external code to specify that output should be appended to a list of SilentStructureOPs, so that the SimpleCycpepPredictApplication doesn't have to write directly to disk. More... | |
void | set_suppress_checkpoints (bool const suppress_checkpoints) |
Allows external code to suppress checkpointing, to prevent direct file I/O from disk. More... | |
void | set_my_rank (int const rank_in) |
If called by MPI code, the rank of the current process can be stored here. More... | |
void | set_already_completed_job_count (core::Size const count_in) |
Set the number of jobs that this process has already completed. More... | |
void | set_nstruct (core::Size const nstruct_in) |
Allows external code to override the number of structures that this should generate (otherwise set by options system. More... | |
void | set_L_alpha_compfile_contents (std::string const &contents_in) |
Allows external code to set the file contents for the L-alpha aa_composition file. More... | |
void | set_D_alpha_compfile_contents (std::string const &contents_in) |
Allows external code to set the file contents for the D-alpha aa_composition file. More... | |
void | set_L_beta_compfile_contents (std::string const &contents_in) |
Allows external code to set the file contents for the L-beta aa_composition file. More... | |
void | set_D_beta_compfile_contents (std::string const &contents_in) |
Allows external code to set the file contents for the D-beta aa_composition file. More... | |
void | set_abba_bins_binfile_contents (std::string const &contents_in) |
Allows external code to set the ABBA bin parameters without having to read the bin params file directly from disk. More... | |
void | set_sample_cis_pro_frequency (core::Real const &freq_in) |
Set the frequency with which we sample cis proline. More... | |
void | disable_cis_pro_sampling () |
Set cis proline sampling OFF. More... | |
void | set_total_energy_cutoff (core::Real const &value_in) |
Set the total energy cutoff. More... | |
void | disable_total_energy_cutoff () |
Sets use_total_energy_cutoff_ to 'false'. More... | |
void | set_angle_relax_rounds (core::Size const rounds_in) |
Set the number of rounds of relaxation with flexible bond angles. More... | |
void | set_angle_length_relax_rounds (core::Size const rounds_in) |
Set the number of rounds of relaxation with flexible bond angles and bond lengths. More... | |
void | set_cartesian_relax_rounds (core::Size const rounds_in) |
Set the number of rounds of Cartesian relaxation. More... | |
void | set_use_rama_prepro_for_sampling (bool const setting) |
Set whether we're using RamaPrePro tables for sampling. More... | |
void | run () const |
Actually run the application. More... | |
Static Public Member Functions | |
static void | register_options () |
Register the set of options that this application uses (for the help menu). More... | |
Private Member Functions | |
core::Size | count_cis_peptide_bonds (core::pose::PoseCOP pose) const |
-------—— Methods ----------------------—— More... | |
void | do_final_fastrelax (core::pose::PoseOP pose, core::scoring::ScoreFunctionOP sfxn, core::Size const relax_rounds, bool const angle_min, bool const length_min, bool const cartesian_min) const |
Carry out the final FastRelax. More... | |
void | build_polymer (core::pose::PoseOP pose, utility::vector1< std::string > const &restypes) const |
Actually build the geometry that we'll be working with. More... | |
void | add_n_methylation (core::pose::PoseOP pose, core::Size const cyclic_offset) const |
Add N-methylation. More... | |
void | set_default_rama_table_type (std::string const &type_name) |
Given the name of a Rama_Table_Type, set the default Rama_Table_Type. More... | |
void | set_rama_table_type_by_res (utility::vector1< std::string > const &type_name_vector) |
Given a string vector that we need to parse, populate the rama_table_type_by_res_ map. More... | |
core::scoring::Rama_Table_Type | get_rama_table_type_from_name (std::string const &type_name) const |
Given a Rama_Table_Type name, return the Rama_Table_Type, or an informative error message on failure. More... | |
void | read_sequence (std::string const &seqfile, utility::vector1< std::string > &resnames) const |
Read a sequence (as a series of full names, separated by whitespace) and store it in a string vector. More... | |
void | set_up_termini_mover (protocols::cyclic_peptide::DeclareBondOP termini, core::pose::PoseCOP pose, bool const native=false, core::Size const last_res=0) const |
Set up the DeclareBond mover used to connect the termini. More... | |
core::Size | do_cyclic_permutation (utility::vector1< std::string > const &resnames, utility::vector1< std::string > &resnames_copy) const |
Takes a vector of residue names, chooses a random number for cyclic offset, and does a cyclic permutation. More... | |
void | import_and_set_up_native (std::string const &native_file, core::pose::PoseOP native_pose, core::Size const expected_residue_count) const |
Imports the native pose and sets up a terminial peptide bond. More... | |
void | set_up_native (core::pose::PoseOP native_pose, core::Size const expected_residue_count) const |
Sets up a terminial peptide bond and does some checks. More... | |
void | add_cyclic_constraints (core::pose::PoseOP pose) const |
Function to add cyclic constraints to a pose. More... | |
void | set_mainchain_torsions (core::pose::PoseOP pose, core::Size const cyclic_offset) const |
Sets all omega values to 180, and randomizes mainchain torsions. More... | |
void | set_up_hbond_filter (protocols::filters::CombinedFilterOP total_hbond, core::Size const nres, core::scoring::ScoreFunctionOP sfxn, core::Real const &min_hbonds) const |
Set up the filters for the mainchain hydrogen bonds that will be used to discard solutions with too little mainchain hydrogen bonding. More... | |
bool | genkic_close (core::pose::PoseOP pose, core::scoring::ScoreFunctionOP sfxn_highhbond, core::scoring::ScoreFunctionOP sfxn_highhbond_cart, core::scoring::ScoreFunctionCOP sfxn_default, protocols::filters::CombinedFilterOP total_hbond, core::Size const cyclic_offset) const |
Use GeneralizedKIC to close the pose. More... | |
void | set_up_design_taskoperations (protocols::denovo_design::movers::FastDesignOP fdes, core::Size const cyclic_offset, core::Size const nres) const |
Set up the TaskOperations that conrol the design process, given user inputs. More... | |
std::string | get_oneletter_codes (utility::vector1< std::string > const &fullnames) const |
Given a vector of full residue names of canonical residues, give me a concatenated list of one-letter codes. More... | |
std::string | get_nc_threeletter_codes (utility::vector1< std::string > const &fullnames) const |
Given a vector of full residue names, give me a string of the form "NC <3-letter code> NC <3-letter code> NC <3-letter code> ...". More... | |
void | store_disulfides (core::pose::PoseCOP pose, utility::vector1< std::pair< core::Size, core::Size > > &old_disulfides) const |
Given a pose, store a list of the disulfides in the pose. More... | |
void | break_disulfides (core::pose::PoseOP pose, utility::vector1< std::pair< core::Size, core::Size > > const &disulfides) const |
Given a pose and a list of the disulfides in the pose, break the disulfides. More... | |
void | rebuild_disulfides (core::pose::PoseOP pose, utility::vector1< std::pair< core::Size, core::Size > > const &disulfides) const |
Given a pose and a list of the disulfides that should be in the pose, form the disulfides. More... | |
void | depermute_disulfide_list (utility::vector1< std::pair< core::Size, core::Size > > const &old_disulfides, utility::vector1< std::pair< core::Size, core::Size > > &new_disulfides, core::Size const offset, core::Size const nres) const |
Given a list of old disulfide positions, generate a list of new disulfide positions based on the offset. More... | |
void | depermute (core::pose::PoseOP pose, core::Size const offset) const |
Given a pose that has undergone an N-residue cyclic permutation, restore the original pose, without the permutation. More... | |
core::Real | align_and_calculate_rmsd (core::pose::PoseOP pose, core::pose::PoseCOP native_pose) const |
Align pose to native_pose, and return the RMSD between the two poses. More... | |
void | new_checkpoint_file () const |
Create a new checkpoint file. More... | |
void | initialize_checkpointing (core::Size &lastjob, core::Size &successes) const |
Initialize checkpointing for this run. More... | |
void | checkpoint (core::Size const curjob, core::Size const successes) const |
Add a checkpoint to the checkpoint file. More... | |
void | end_checkpointing () const |
End checkpointing and delete the checkpoint file. More... | |
void | get_random_seed_info () const |
Restore the state of the random generator from a previous run. More... | |
void | store_random_seed_info () const |
Store the state of the random generator from a previous run. More... | |
void | erase_random_seed_info () const |
Erase the stored state of the random generator from a previous run. More... | |
core::Size | current_position (core::Size const absolute_position, core::Size const permutation_offset, core::Size const nresidue) const |
Given an absolute position in the sequence and the current cyclic permuation offset, return the position in the current pose corresponding to that absolute postion. More... | |
core::Size | original_position (core::Size const curr_position, core::Size const permutation_offset, core::Size const nresidue) const |
Given an position in the current (perturbed) pose, return the position in the original (unperturbed) pose. More... | |
bool | custom_rama_table_defined (core::Size const absolute_position) const |
Does a position have a custom Rama table defined? More... | |
core::scoring::Rama_Table_Type | rama_table_type_by_res (core::Size const absolute_position) const |
Custom Rama table for a position. More... | |
core::scoring::Rama_Table_Type | default_rama_table_type () const |
Get the default custom Rama table type. More... | |
bool | use_rama_filter () const |
Are we using a rama filter? Defaults to true. More... | |
bool | sample_cis_pro () const |
Are we sampling cis-prolines? More... | |
core::Real const & | sample_cis_pro_frequency () const |
Frequency for sampling cis prolines. More... | |
core::Size | angle_relax_rounds () const |
Bond angle relax rounds. More... | |
core::Size | angle_length_relax_rounds () const |
Bond angle / bond length relax rounds. More... | |
core::Size | cartesian_relax_rounds () const |
Cartesian relax rounds. More... | |
bool | use_rama_prepro_for_sampling () const |
Are we using RamaPrePro tables for sampling? More... | |
core::Size | sequence_length () const |
The length of the sequence, excluding crosslinkers. More... | |
void | re_append_tbmb_residues (core::pose::PoseCOP pose, core::pose::PoseOP newpose, core::Size const offset) const |
Given a pose with TBMB in it and another pose without TBMB, copy the TBMB residues from the first to the second, and add back covalent bonds. More... | |
Private Attributes | |
int | my_rank_ |
If this is called by MPI code, this can store the rank of the current process. Zero otherwise. More... | |
core::Size | already_completed_job_count_ |
The number of jobs that this slave process has already completed. Only used in MPI mode; zero otherwise. More... | |
core::scoring::ScoreFunctionOP | scorefxn_ |
The default ScoreFunction to use. The high h-bond version is constructed from this. More... | |
bool | suppress_checkpoints_ |
Allows external code to suppress checkpointing, so that the SimpleCycpepPredictApplication doesn't write directly to disk. More... | |
bool | silent_out_ |
Should this application produce silent file output? More... | |
bool | silentlist_out_ |
Should this application produce silent structure OP output, appending to a list? More... | |
utility::vector1 < core::io::silent::SilentStructOP > * | silentlist_ |
If silentlist_out_ is used, this is the list to append SilentStructurOPs to. More... | |
utility::vector1 < SimpleCycpepPredictApplication_MPI_JobResultsSummaryOP > * | summarylist_ |
If silentlist_out_ is used, this is the list of job summaries to which summaries should be appended. More... | |
core::pose::PoseCOP | native_pose_ |
A native pose, provided by external code. More... | |
std::string | out_filename_ |
The prefix for the output filename. More... | |
std::string | out_scorefilename_ |
The output score file name. More... | |
std::string | sequence_file_ |
Filename for the text file containing the sequence of the peptide. More... | |
std::string | sequence_string_ |
The string that would be read from a sequence file. More... | |
core::Size | sequence_length_ |
The length of the sequence, excluding any cross-linkers. More... | |
core::Size | genkic_closure_attempts_ |
The number of attempts that will be made by the generalized kinematic closure machinery. More... | |
core::Size | genkic_min_solution_count_ |
The minimum number of solutions that must be found by the generalized kinematic closure machinery before a single solution is chosen. More... | |
bool | cyclic_permutations_ |
Should we consider cyclic permutations? More... | |
bool | use_rama_filter_ |
Should solutions be filtered by Ramachandran energy (in the case of alpha-amino acids)? More... | |
core::Real | rama_cutoff_ |
The energy cutoff for the rama filter. More... | |
core::Real | high_hbond_weight_multiplier_ |
For parts of the protocol that upweight the hydrogen bond terms, this is the factor by which these terms are upweighted. More... | |
core::Real | min_genkic_hbonds_ |
The minimum number of mainchain hydrogen bonds that a generalized kinematic closure solution must have. More... | |
core::Real | min_final_hbonds_ |
The minimum number of mainchain hydrogen bonds that a final solution must have. More... | |
core::Real | total_energy_cutoff_ |
The total energy cutoff above which solutions are discarded. More... | |
bool | use_total_energy_cutoff_ |
Determines whether the total energy cutoff should be used. More... | |
core::Real | hbond_energy_cutoff_ |
The hbond energy cutoff above which a hbond is no longer counted. More... | |
core::Size | fast_relax_rounds_ |
The number of FastReleax rounds to be applied whenever FastRelax is used. More... | |
bool | count_sc_hbonds_ |
Should sidechain-mainchain hydrogen bonds be counted when counting hydrogen bonds? More... | |
bool | native_exists_ |
Has the user specified a native structure? More... | |
std::string | native_filename_ |
Filename for the native, if one is provided. More... | |
core::Size | nstruct_ |
Max number of structures to generate. More... | |
std::string | checkpoint_job_identifier_ |
A unique job name for checkpointing. More... | |
std::string | checkpoint_filename_ |
The name of the checkpoint file. More... | |
core::scoring::Rama_Table_Type | default_rama_table_type_ |
A default custom Ramachandran table to be used for sampling. More... | |
std::map< core::Size, core::scoring::Rama_Table_Type > | rama_table_type_by_res_ |
Custom Ramachandran tables to be used for sampling, listed by residue. More... | |
std::string | rand_checkpoint_file_ |
The name of the checkpoint file for the random number generator. More... | |
bool | try_all_disulfides_ |
Should we try all disulfide combinations during structure prediction? More... | |
core::Real | disulf_energy_cutoff_prerelax_ |
The cutoff dslf_fa13 energy, pre-relaxation, above which closure solutions are rejected. More... | |
core::Real | disulf_energy_cutoff_postrelax_ |
The cutoff dslf_fa13 energy, post-relaxation, above which closure solutions are rejected. More... | |
std::map< core::Size, utility::vector1< core::Real > > | user_set_alpha_dihedrals_ |
Map of (seqpos -> phi/psi/omega triplet). More... | |
core::Real | user_set_dihedral_perturbation_ |
A small random value added to all user-set dihedral values. Defaults to 0. More... | |
bool | filter_oversaturated_hbond_acceptors_ |
Should we filter out solutions that have more than the allowed number of hydrogen bonds to an acceptor? More... | |
core::Real | oversaturated_hbond_cutoff_energy_ |
The energy above which we no longer count a hydrogen bond to an acceptor for the filter. More... | |
bool | sample_cis_pro_ |
At residues that are more likely to have cis omega values, sample cis some fraction of the time. More... | |
core::Real | sample_cis_pro_frequency_ |
At residues that are more likely to have cis omega values, sample cis this fraction of the time. More... | |
bool | design_peptide_ |
Should we design the peptide after sampling the conformation, or just relax it? More... | |
std::string | design_filename_ |
If we're designing, this is the filename for the file that tells us what residues are allowed at what postion. More... | |
bool | prevent_design_file_read_ |
Used to prevent this process from directly reading the input file if it has been read by another process and the information has been transmitted to this process. More... | |
std::map< core::Size, utility::vector1< std::string > > | allowed_canonicals_by_position_ |
Allowed canonical residues at each position. More... | |
std::map< core::Size, utility::vector1< std::string > > | allowed_noncanonicals_by_position_ |
Allowed noncanonical residues at each position. More... | |
bool | prohibit_D_at_negative_phi_ |
Should D-residues be prohibited at positions with negative phi values? More... | |
bool | prohibit_L_at_positive_phi_ |
Should L-residues be prohibited at positions with positive phi values? More... | |
bool | use_aa_comp_ |
Should we use the aa_composition score term in design? More... | |
bool | L_alpha_comp_file_exists_ |
Has an aa_composition setup file ben provided for residues in the L-alpha helix region of Ramachadran space? More... | |
bool | D_alpha_comp_file_exists_ |
Has an aa_composition setup file ben provided for residues in the D-alpha helix region of Ramachadran space? More... | |
bool | L_beta_comp_file_exists_ |
Has an aa_composition setup file ben provided for residues in the L-beta strand region of Ramachadran space? More... | |
bool | D_beta_comp_file_exists_ |
Has an aa_composition setup file ben provided for residues in the D-beta strand region of Ramachadran space? More... | |
std::string | comp_file_contents_L_alpha_ |
Storage for the composition constraint setup for the L-alpha helix region of Ramachandran space. More... | |
std::string | comp_file_contents_D_alpha_ |
Storage for the composition constraint setup for the D-alpha helix region of Ramachandran space. More... | |
std::string | comp_file_contents_L_beta_ |
Storage for the composition constraint setup for the L-beta strand region of Ramachandran space. More... | |
std::string | comp_file_contents_D_beta_ |
Storage for the composition constraint setup for the D-beta strand region of Ramachandran space. More... | |
std::string | abba_bins_ |
Storage for a bin definition file. More... | |
bool | do_not_count_adjacent_res_hbonds_ |
Prevent counting of hydrogen bonds to adjacent residues? More... | |
core::Size | angle_relax_rounds_ |
Number of rounds of relaxation with flexible bond angles. More... | |
core::Size | angle_length_relax_rounds_ |
Number of rounds of relaxation with flexible bond angles and bond lengths. More... | |
core::Size | cartesian_relax_rounds_ |
Number of rounds of Cartesian-space relaxation. More... | |
bool | use_rama_prepro_for_sampling_ |
Should we use rama_prepro tables for sampling? More... | |
utility::vector1< core::Size > | n_methyl_positions_ |
List of positions (in original sequence indexing – not permuted) that are N-methylated. More... | |
utility::vector1 < utility::vector1< core::Size > > | tbmb_positions_ |
List of positions linked by 1,3,5-tris(bromomethyl)benzene. More... | |
bool | link_all_cys_with_tbmb_ |
Should all cysteine residues be linked with 1,3,5-tris(bromomethyl)benzene? More... | |
bool | use_tbmb_filters_ |
If true, filters are applied based on distance between TBMB cysteines and on constraints to discard GenKIC solutions that can't be crosslinked easily. More... | |
core::Real | tbmb_sidechain_distance_filter_multiplier_ |
Multiplier to make the TBMB distance filter more permissive. More... | |
core::Real | tbmb_constraints_energy_filter_multiplier_ |
Multiplier to make the TBMB constraints energy filter more permissive. More... | |
core::Size | required_symmetry_repeats_ |
If this option is used, then only backbones that are cN (or cN/m, if mirror symmetry is required) symmetric will be accepted. For example, if set to 2, then only c2-symmetric backbones will be accepted. More... | |
bool | required_symmetry_mirroring_ |
If this option is used, then only backbones with mirror symmetry are accepted. Must be used with the -cyclic_peptide:require_symmetry_repeats flag. More... | |
core::Real | required_symmetry_angle_threshold_ |
The cutoff, in degrees, to use when comparing mainchain torsion values to determine whether symmetry repeats are truely symmetric. Defaults to 10 degrees. More... | |
core::Real | required_symmetry_perturbation_ |
The random perturbation, in degrees, to be applied when copying mainchain torsion values to produce symmetric conformations. More... | |
Application-level code for simple_cycpep_predict application.
Also called by the BOINC minirosetta app.
protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::SimpleCycpepPredictApplication | ( | bool const | allow_file_read = true | ) |
Constructor.
If allow_file_read is true, initialization triggers reads from the filesystem.
References initialize_from_options().
|
overridedefault |
Explicit virtual destructor.
protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::SimpleCycpepPredictApplication | ( | SimpleCycpepPredictApplication const & | src | ) |
Explicit copy constructor.
References scorefxn_.
|
private |
Function to add cyclic constraints to a pose.
References core::id::AtomID::atomno(), core::id::AtomID::rsd(), sequence_length(), SimpleCycpepPredictApplication_PEPBOND_C_ANGLE, SimpleCycpepPredictApplication_PEPBOND_LENGTH, SimpleCycpepPredictApplication_PEPBOND_N_ANGLE, and protocols::TR().
Referenced by run().
|
private |
Add N-methylation.
Must be called after pose is cyclized.
References protocols::simple_moves::ModifyVariantTypeMover::apply(), n_methyl_positions_, sequence_length(), protocols::simple_moves::ModifyVariantTypeMover::set_additional_type_to_add(), protocols::simple_moves::ModifyVariantTypeMover::set_residue_selector(), and protocols::TR().
Referenced by run().
|
private |
Align pose to native_pose, and return the RMSD between the two poses.
Assumes that the pose has already been de-permuted (i.e. the native and the pose line up). Only uses alpha-amino acids for the alignment, currently.
References core::id::BOGUS_ATOM_ID, core::pose::initialize_atomid_map(), sequence_length(), and core::scoring::superimpose_pose().
Referenced by run().
|
inlineprivate |
Bond angle / bond length relax rounds.
References angle_length_relax_rounds_.
Referenced by genkic_close(), and run().
|
inlineprivate |
|
private |
Given a pose and a list of the disulfides in the pose, break the disulfides.
References core::conformation::break_disulfide().
Referenced by depermute().
|
private |
Actually build the geometry that we'll be working with.
References protocols::cyclic_peptide::PeptideStubMover::set_reset_mode(), and protocols::TR().
Referenced by run().
|
inlineprivate |
Cartesian relax rounds.
References cartesian_relax_rounds_.
Referenced by genkic_close(), and run().
|
private |
Add a checkpoint to the checkpoint file.
The checkpoint file must already exist. Does nothing if checkpointing is disabled.
[in] | curjob | The index of the current job just run, for writing to the checkpoint file. |
[in] | successes | The number of successes so far, for writing to the checkpoint file. |
References checkpoint_filename_, checkpoint_job_identifier_, store_random_seed_info(), suppress_checkpoints_, and protocols::boinc::Boinc::update_pct_complete().
Referenced by run().
|
private |
-------—— Methods ----------------------——
Count the number of cis-peptide bonds in the pose.
Counts as cis if in the range (-90,90].
References sequence_length(), and protocols::TR().
Referenced by run().
|
inlineprivate |
Given an absolute position in the sequence and the current cyclic permuation offset, return the position in the current pose corresponding to that absolute postion.
Referenced by genkic_close().
|
inlineprivate |
Does a position have a custom Rama table defined?
Does not include a default custom Rama table – only position-specific Rama tables are checked.
References rama_table_type_by_res_.
Referenced by genkic_close(), rama_table_type_by_res(), and set_mainchain_torsions().
|
inlineprivate |
Get the default custom Rama table type.
Defaults to unknown_ramatable_type if not set.
References default_rama_table_type_.
Referenced by genkic_close(), and set_mainchain_torsions().
|
private |
Given a pose that has undergone an N-residue cyclic permutation, restore the original pose, without the permutation.
References break_disulfides(), depermute_disulfide_list(), re_append_tbmb_residues(), rebuild_disulfides(), sequence_length(), set_up_termini_mover(), store_disulfides(), and tbmb_positions_.
Referenced by run().
|
private |
Given a list of old disulfide positions, generate a list of new disulfide positions based on the offset.
Replaces the new_disulfides list.
Referenced by depermute().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::disable_cis_pro_sampling | ( | ) |
Set cis proline sampling OFF.
References sample_cis_pro_, and sample_cis_pro_frequency_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::disable_total_energy_cutoff | ( | ) |
Sets use_total_energy_cutoff_ to 'false'.
References use_total_energy_cutoff_.
|
private |
Takes a vector of residue names, chooses a random number for cyclic offset, and does a cyclic permutation.
Returns the offset and stores the new string vector in resnames_copy.
References core::scoring::rg, and protocols::TR().
Referenced by run().
|
private |
Carry out the final FastRelax.
References set_up_termini_mover(), and protocols::TR().
Referenced by run().
|
private |
End checkpointing and delete the checkpoint file.
Does nothing if checkpointing is disabled.
References checkpoint_filename_, checkpoint_job_identifier_, erase_random_seed_info(), and suppress_checkpoints_.
Referenced by run().
|
private |
Erase the stored state of the random generator from a previous run.
References rand_checkpoint_file_, and suppress_checkpoints_.
Referenced by end_checkpointing(), and new_checkpoint_file().
|
private |
Use GeneralizedKIC to close the pose.
References core::chemical::aa_dpr, core::chemical::aa_pro, abba_bins_, protocols::generalized_kinematic_closure::filter::alpha_aa_rama_check, angle_length_relax_rounds(), angle_relax_rounds(), cartesian_relax_rounds(), comp_file_contents_D_alpha_, comp_file_contents_D_beta_, comp_file_contents_L_alpha_, comp_file_contents_L_beta_, protocols::generalized_kinematic_closure::perturber::copy_backbone_dihedrals, current_position(), custom_rama_table_defined(), D_alpha_comp_file_exists_, D_beta_comp_file_exists_, default_rama_table_type(), design_peptide_, disulf_energy_cutoff_postrelax_, disulf_energy_cutoff_prerelax_, core::scoring::dslf_fa13, fast_relax_rounds_, filter_oversaturated_hbond_acceptors_, genkic_closure_attempts_, genkic_min_solution_count_, L_alpha_comp_file_exists_, L_beta_comp_file_exists_, protocols::generalized_kinematic_closure::filter::loop_bump_check, protocols::generalized_kinematic_closure::selector::lowest_energy_selector, min_genkic_hbonds_, protocols::generalized_kinematic_closure::perturber::mirror_backbone_dihedrals, original_position(), oversaturated_hbond_cutoff_energy_, protocols::generalized_kinematic_closure::perturber::perturb_dihedral, rama_cutoff_, protocols::generalized_kinematic_closure::filter::rama_prepro_check, rama_table_type_by_res(), protocols::generalized_kinematic_closure::perturber::randomize_alpha_backbone_by_rama, protocols::generalized_kinematic_closure::perturber::randomize_backbone_by_rama_prepro, required_symmetry_angle_threshold_, required_symmetry_mirroring_, required_symmetry_perturbation_, required_symmetry_repeats_, core::scoring::rg, protocols::generalized_kinematic_closure::perturber::sample_cis_peptide_bond, sample_cis_pro(), sample_cis_pro_frequency(), sequence_length(), protocols::generalized_kinematic_closure::perturber::set_dihedral, set_up_design_taskoperations(), set_up_termini_mover(), SimpleCycpepPredictApplication_PEPBOND_C_ANGLE, SimpleCycpepPredictApplication_PEPBOND_LENGTH, SimpleCycpepPredictApplication_PEPBOND_N_ANGLE, tbmb_constraints_energy_filter_multiplier_, tbmb_positions_, tbmb_sidechain_distance_filter_multiplier_, protocols::TR(), try_all_disulfides_, core::scoring::unknown_ramatable_type, use_rama_filter(), use_rama_prepro_for_sampling(), use_tbmb_filters_, user_set_alpha_dihedrals_, and user_set_dihedral_perturbation_.
Referenced by run().
|
private |
Given a vector of full residue names, give me a string of the form "NC <3-letter code> NC <3-letter code> NC <3-letter code> ...".
Does no checking for duplicates. Will fail gracelessly with invalid names.
References core::chemical::FA_STANDARD, core::chemical::ResidueTypeFinder::get_representative_type(), and core::chemical::ResidueTypeFinder::residue_base_name().
Referenced by set_up_design_taskoperations().
|
private |
Given a vector of full residue names of canonical residues, give me a concatenated list of one-letter codes.
Does no checking for duplicates.
References core::chemical::aa_from_name(), and core::chemical::oneletter_code_from_aa().
Referenced by set_up_design_taskoperations().
|
private |
Given a Rama_Table_Type name, return the Rama_Table_Type, or an informative error message on failure.
References core::scoring::Ramachandran::get_ramatable_name_by_type(), core::scoring::Ramachandran::get_ramatable_type_by_name(), core::scoring::rama, and core::scoring::unknown_ramatable_type.
Referenced by set_default_rama_table_type(), and set_rama_table_type_by_res().
|
private |
Restore the state of the random generator from a previous run.
References rand_checkpoint_file_, and core::scoring::rg.
Referenced by initialize_checkpointing().
|
private |
Imports the native pose and sets up a terminial peptide bond.
References core::import_pose::PDB_file, core::import_pose::pose_from_file(), set_up_native(), and protocols::TR().
Referenced by run().
|
private |
Initialize checkpointing for this run.
This function does several things. First, it checks for an existing checkpoint file. If one exists, it checks whether the unique job name in the file matches the current job. If it does, then this job has already been attempted, and we're somewhere in the middle of it. The function reads the last attempt number and success count from the checkpoint file, and returns these values. Otherwise, it creates a new checkpoint file with the current job name and returns (0,0). If checkpointing is disabled, this function does nothing, and returns (0,0).
[out] | lastjob | The index of the last job run. Set to zero if checkpointing is disabled or if we're creating a new checkpoint file (first job run). |
[out] | successes | The number of successes so far. Set to zero if checkpointing is disabled or if we're creating a new checkpoint file (first job run). |
References checkpoint_filename_, checkpoint_job_identifier_, get_random_seed_info(), new_checkpoint_file(), suppress_checkpoints_, and protocols::TR().
Referenced by run().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::initialize_from_options | ( | ) |
Initialize the application.
Initializes using the option system.
References abba_bins_, allowed_canonicals_by_position_, allowed_noncanonicals_by_position_, checkpoint_filename_, checkpoint_job_identifier_, comp_file_contents_D_alpha_, comp_file_contents_D_beta_, comp_file_contents_L_alpha_, comp_file_contents_L_beta_, count_sc_hbonds_, cyclic_permutations_, D_alpha_comp_file_exists_, D_beta_comp_file_exists_, design_filename_, design_peptide_, disulf_energy_cutoff_postrelax_, disulf_energy_cutoff_prerelax_, do_not_count_adjacent_res_hbonds_, fast_relax_rounds_, filter_oversaturated_hbond_acceptors_, genkic_closure_attempts_, genkic_min_solution_count_, core::scoring::get_score_function(), hbond_energy_cutoff_, high_hbond_weight_multiplier_, L_alpha_comp_file_exists_, L_beta_comp_file_exists_, link_all_cys_with_tbmb_, min_final_hbonds_, min_genkic_hbonds_, n_methyl_positions_, native_exists_, native_filename_, nstruct_, out_filename_, out_scorefilename_, oversaturated_hbond_cutoff_energy_, prevent_design_file_read_, prohibit_D_at_negative_phi_, prohibit_L_at_positive_phi_, rama_cutoff_, rand_checkpoint_file_, protocols::cyclic_peptide_predict::read_file_into_string(), protocols::cyclic_peptide_predict::read_peptide_design_file(), required_symmetry_angle_threshold_, required_symmetry_mirroring_, required_symmetry_perturbation_, required_symmetry_repeats_, scorefxn_, sequence_file_, set_angle_length_relax_rounds(), set_angle_relax_rounds(), set_cartesian_relax_rounds(), set_default_rama_table_type(), set_rama_table_type_by_res(), set_sample_cis_pro_frequency(), set_total_energy_cutoff(), set_use_rama_prepro_for_sampling(), silent_out_, tbmb_constraints_energy_filter_multiplier_, tbmb_positions_, tbmb_sidechain_distance_filter_multiplier_, protocols::TR(), try_all_disulfides_, use_aa_comp_, use_rama_filter_, use_tbmb_filters_, user_set_alpha_dihedrals_, and user_set_dihedral_perturbation_.
Referenced by SimpleCycpepPredictApplication().
|
private |
Create a new checkpoint file.
References checkpoint_filename_, checkpoint_job_identifier_, erase_random_seed_info(), store_random_seed_info(), and suppress_checkpoints_.
Referenced by initialize_checkpointing().
|
inlineprivate |
Given an position in the current (perturbed) pose, return the position in the original (unperturbed) pose.
Referenced by genkic_close(), and set_mainchain_torsions().
|
inlineprivate |
Custom Rama table for a position.
Returns unknown_ramatable_type if none defined for the position.
References custom_rama_table_defined(), rama_table_type_by_res_, and core::scoring::unknown_ramatable_type.
Referenced by genkic_close(), and set_mainchain_torsions().
|
private |
Given a pose with TBMB in it and another pose without TBMB, copy the TBMB residues from the first to the second, and add back covalent bonds.
This function is called at the end of the protocol, and therefore doesn't bother to add back constraints.
References protocols::cyclic_peptide::threefold_linker::TBMB_Helper::add_linker_bonds_asymmetric(), and tbmb_positions_.
Referenced by depermute().
|
private |
Read a sequence (as a series of full names, separated by whitespace) and store it in a string vector.
References sequence_string_, and protocols::TR().
Referenced by run().
|
private |
Given a pose and a list of the disulfides that should be in the pose, form the disulfides.
References core::conformation::form_disulfide().
Referenced by depermute().
|
static |
Register the set of options that this application uses (for the help menu).
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::run | ( | ) | const |
Actually run the application.
The initialize_from_options() function must be called before calling this. (Called by default constructor.)
References core::scoring::aa_composition, add_cyclic_constraints(), add_n_methylation(), align_and_calculate_rmsd(), already_completed_job_count_, core::scoring::angle_constraint, angle_length_relax_rounds(), angle_relax_rounds(), core::scoring::atom_pair_constraint, build_polymer(), core::scoring::cart_bonded, cartesian_relax_rounds(), checkpoint(), protocols::fldsgn::topology::compare(), count_cis_peptide_bonds(), cyclic_permutations_, depermute(), core::scoring::dihedral_constraint, do_cyclic_permutation(), do_final_fastrelax(), end_checkpointing(), fast_relax_rounds_, genkic_close(), core::scoring::hbond_lr_bb, core::scoring::hbond_sr_bb, high_hbond_weight_multiplier_, import_and_set_up_native(), initialize_checkpointing(), link_all_cys_with_tbmb_, min_final_hbonds_, min_genkic_hbonds_, protocols::boinc::BoincSharedMemory::model_count, my_rank_, native_exists_, native_filename_, native_pose_, nstruct_, out_filename_, core::scoring::pro_close, read_sequence(), required_symmetry_angle_threshold_, required_symmetry_mirroring_, required_symmetry_repeats_, protocols::noesy_assign::round(), scorefxn_, sequence_file_, sequence_length(), sequence_length_, set_mainchain_torsions(), set_up_hbond_filter(), set_up_termini_mover(), silent_out_, silentlist_, silentlist_out_, summarylist_, tbmb_positions_, total_energy_cutoff_, protocols::TR(), use_total_energy_cutoff_, and protocols::boinc::Boinc::worker_is_finished().
|
inlineprivate |
Are we sampling cis-prolines?
Defaults to false.
References sample_cis_pro_.
Referenced by genkic_close().
|
inlineprivate |
Frequency for sampling cis prolines.
References sample_cis_pro_frequency_.
Referenced by genkic_close().
|
inlineprivate |
The length of the sequence, excluding crosslinkers.
This is also the index of the last sequence residue. Crosslinker residues follow this position.
References sequence_length_.
Referenced by add_cyclic_constraints(), add_n_methylation(), align_and_calculate_rmsd(), count_cis_peptide_bonds(), depermute(), genkic_close(), run(), set_mainchain_torsions(), and set_up_termini_mover().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_abba_bins_binfile_contents | ( | std::string const & | contents_in | ) |
Allows external code to set the ABBA bin parameters without having to read the bin params file directly from disk.
References abba_bins_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_allowed_residues_by_position | ( | std::map< core::Size, utility::vector1< std::string > > const & | allowed_canonicals, |
std::map< core::Size, utility::vector1< std::string > > const & | allowed_noncanonicals | ||
) |
Allows external code to set the allowed residues by position, so that this needn't be read directly from disk.
References allowed_canonicals_by_position_, and allowed_noncanonicals_by_position_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_already_completed_job_count | ( | core::Size const | count_in | ) |
Set the number of jobs that this process has already completed.
References already_completed_job_count_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_angle_length_relax_rounds | ( | core::Size const | rounds_in | ) |
Set the number of rounds of relaxation with flexible bond angles and bond lengths.
References angle_length_relax_rounds_.
Referenced by initialize_from_options().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_angle_relax_rounds | ( | core::Size const | rounds_in | ) |
Set the number of rounds of relaxation with flexible bond angles.
References angle_relax_rounds_.
Referenced by initialize_from_options().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_cartesian_relax_rounds | ( | core::Size const | rounds_in | ) |
Set the number of rounds of Cartesian relaxation.
References cartesian_relax_rounds_.
Referenced by initialize_from_options().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_D_alpha_compfile_contents | ( | std::string const & | contents_in | ) |
Allows external code to set the file contents for the D-alpha aa_composition file.
References comp_file_contents_D_alpha_, and D_alpha_comp_file_exists_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_D_beta_compfile_contents | ( | std::string const & | contents_in | ) |
Allows external code to set the file contents for the D-beta aa_composition file.
References comp_file_contents_D_beta_, and D_beta_comp_file_exists_.
|
private |
Given the name of a Rama_Table_Type, set the default Rama_Table_Type.
Error if unknown type.
References default_rama_table_type_, get_rama_table_type_from_name(), and protocols::TR().
Referenced by initialize_from_options().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_L_alpha_compfile_contents | ( | std::string const & | contents_in | ) |
Allows external code to set the file contents for the L-alpha aa_composition file.
References comp_file_contents_L_alpha_, and L_alpha_comp_file_exists_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_L_beta_compfile_contents | ( | std::string const & | contents_in | ) |
Allows external code to set the file contents for the L-beta aa_composition file.
References comp_file_contents_L_beta_, and L_beta_comp_file_exists_.
|
private |
Sets all omega values to 180, and randomizes mainchain torsions.
For alpha-amino acids, mainchain torsions are randomized by the Ramachandran plot. For other residue types, just randomizes mainchain torsions other than peptide bonds.
References core::chemical::aa_unk, core::id::BB, custom_rama_table_defined(), default_rama_table_type(), core::scoring::Ramachandran::draw_random_phi_psi_from_extra_cdf(), original_position(), core::scoring::rama, rama_table_type_by_res(), core::scoring::RamaPrePro::random_mainchain_torsions(), core::scoring::Ramachandran::random_phipsi_from_rama(), core::scoring::rg, sequence_length(), protocols::TR(), core::scoring::unknown_ramatable_type, and use_rama_prepro_for_sampling().
Referenced by run().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_my_rank | ( | int const | rank_in | ) |
If called by MPI code, the rank of the current process can be stored here.
Used for output of job summaries.
References my_rank_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_native | ( | core::pose::PoseCOP | native | ) |
Allows external code to provide a native, so that the SimpleCycpepPredictApplication doesn't have to read directly from disk.
References native_exists_, native_pose_, and set_up_native().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_nstruct | ( | core::Size const | nstruct_in | ) |
Allows external code to override the number of structures that this should generate (otherwise set by options system.
References nstruct_.
|
private |
Given a string vector that we need to parse, populate the rama_table_type_by_res_ map.
The string vector must be of the format: [integer] [rama type name] [integer] [rama type name] etc. Throws error if could not parse.
References get_rama_table_type_from_name(), rama_table_type_by_res_, and protocols::TR().
Referenced by initialize_from_options().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_sample_cis_pro_frequency | ( | core::Real const & | freq_in | ) |
Set the frequency with which we sample cis proline.
Implicitly sets sample_cis_pro_ to "true" if freq_in is not 0.0, "false" if it is.
References sample_cis_pro_, sample_cis_pro_frequency_, and protocols::TR().
Referenced by initialize_from_options().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_scorefxn | ( | core::scoring::ScoreFunctionCOP | sfxn_in | ) |
Sets the default scorefunction to use.
The scorefunction is cloned. The high-hbond version is constructed from this one. If necessary, the aa_composition score term will be turned on in that one; it needn't be turned on in this one.
References scorefxn_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_sequence | ( | std::string const & | seq | ) |
Allows external code to provide a sequence, so that the SimpleCycpepPredictApplication doesn't have to read directly from disk.
References sequence_string_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_silentstructure_outputlist | ( | utility::vector1< core::io::silent::SilentStructOP > * | silentlist, |
utility::vector1< SimpleCycpepPredictApplication_MPI_JobResultsSummaryOP > * | summarylist | ||
) |
Allows external code to specify that output should be appended to a list of SilentStructureOPs, so that the SimpleCycpepPredictApplication doesn't have to write directly to disk.
References silent_out_, silentlist_, silentlist_out_, and summarylist_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_suppress_checkpoints | ( | bool const | suppress_checkpoints | ) |
Allows external code to suppress checkpointing, to prevent direct file I/O from disk.
Useful on Blue Gene.
References suppress_checkpoints_.
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_total_energy_cutoff | ( | core::Real const & | value_in | ) |
Set the total energy cutoff.
Also sets use_total_energy_cutoff_ to 'true'.
References total_energy_cutoff_, and use_total_energy_cutoff_.
Referenced by initialize_from_options().
|
private |
Set up the TaskOperations that conrol the design process, given user inputs.
Default behaviour is designing all positions with L-canonicals and their D-equivalents EXCEPT cys and met (and gly), unless the user overrides this.
References allowed_canonicals_by_position_, allowed_noncanonicals_by_position_, get_nc_threeletter_codes(), get_oneletter_codes(), prohibit_D_at_negative_phi_, prohibit_L_at_positive_phi_, and protocols::TR().
Referenced by genkic_close().
|
private |
Set up the filters for the mainchain hydrogen bonds that will be used to discard solutions with too little mainchain hydrogen bonding.
References count_sc_hbonds_, do_not_count_adjacent_res_hbonds_, and hbond_energy_cutoff_.
Referenced by run().
|
private |
Sets up a terminial peptide bond and does some checks.
References core::pose::remove_lower_terminus_type_from_pose_residue(), core::pose::remove_upper_terminus_type_from_pose_residue(), set_up_termini_mover(), and protocols::TR().
Referenced by import_and_set_up_native(), and set_native().
|
private |
Set up the DeclareBond mover used to connect the termini.
References sequence_length(), and protocols::TR().
Referenced by depermute(), do_final_fastrelax(), genkic_close(), run(), and set_up_native().
void protocols::cyclic_peptide_predict::SimpleCycpepPredictApplication::set_use_rama_prepro_for_sampling | ( | bool const | setting | ) |
Set whether we're using RamaPrePro tables for sampling.
Setting this to "false" lets us use classic rama tables. True by default.
References use_rama_prepro_for_sampling_.
Referenced by initialize_from_options().
|
private |
Given a pose, store a list of the disulfides in the pose.
Clears the old_disulfides list and repopulates it.
References core::conformation::disulfide_bonds().
Referenced by depermute().
|
private |
Store the state of the random generator from a previous run.
References rand_checkpoint_file_, core::scoring::rg, and suppress_checkpoints_.
Referenced by checkpoint(), and new_checkpoint_file().
|
inlineprivate |
Are we using a rama filter? Defaults to true.
References use_rama_filter_.
Referenced by genkic_close().
|
inlineprivate |
Are we using RamaPrePro tables for sampling?
True if we are, false if we're using classic rama tables.
References use_rama_prepro_for_sampling_.
Referenced by genkic_close(), and set_mainchain_torsions().
|
private |
Storage for a bin definition file.
Cached to prevent repeated read from disk.
Referenced by genkic_close(), initialize_from_options(), and set_abba_bins_binfile_contents().
|
private |
Allowed canonical residues at each position.
Only used for design. Map key "0" is used for default settings applied to any position not explicitly specified.
Referenced by initialize_from_options(), set_allowed_residues_by_position(), and set_up_design_taskoperations().
|
private |
Allowed noncanonical residues at each position.
Only used for design. Map key "0" is used for default settings applied to any position not explicitly specified.
Referenced by initialize_from_options(), set_allowed_residues_by_position(), and set_up_design_taskoperations().
|
private |
The number of jobs that this slave process has already completed. Only used in MPI mode; zero otherwise.
Referenced by run(), and set_already_completed_job_count().
|
private |
Number of rounds of relaxation with flexible bond angles and bond lengths.
Default 0.
Referenced by angle_length_relax_rounds(), and set_angle_length_relax_rounds().
|
private |
Number of rounds of relaxation with flexible bond angles.
Default 0.
Referenced by angle_relax_rounds(), and set_angle_relax_rounds().
|
private |
Number of rounds of Cartesian-space relaxation.
Default 0.
Referenced by cartesian_relax_rounds(), and set_cartesian_relax_rounds().
|
private |
The name of the checkpoint file.
Defaults to "checkpoint.txt". Read from options.
Referenced by checkpoint(), end_checkpointing(), initialize_checkpointing(), initialize_from_options(), and new_checkpoint_file().
|
private |
A unique job name for checkpointing.
If none is provided, job is not checkpointed.
Referenced by checkpoint(), end_checkpointing(), initialize_checkpointing(), initialize_from_options(), and new_checkpoint_file().
|
private |
Storage for the composition constraint setup for the D-alpha helix region of Ramachandran space.
Cached to prevent repeated read from disk.
Referenced by genkic_close(), initialize_from_options(), and set_D_alpha_compfile_contents().
|
private |
Storage for the composition constraint setup for the D-beta strand region of Ramachandran space.
Cached to prevent repeated read from disk.
Referenced by genkic_close(), initialize_from_options(), and set_D_beta_compfile_contents().
|
private |
Storage for the composition constraint setup for the L-alpha helix region of Ramachandran space.
Cached to prevent repeated read from disk.
Referenced by genkic_close(), initialize_from_options(), and set_L_alpha_compfile_contents().
|
private |
Storage for the composition constraint setup for the L-beta strand region of Ramachandran space.
Cached to prevent repeated read from disk.
Referenced by genkic_close(), initialize_from_options(), and set_L_beta_compfile_contents().
|
private |
Should sidechain-mainchain hydrogen bonds be counted when counting hydrogen bonds?
Defaults to false.
Referenced by initialize_from_options(), and set_up_hbond_filter().
|
private |
Should we consider cyclic permutations?
Defaults to true.
Referenced by initialize_from_options(), and run().
|
private |
Has an aa_composition setup file ben provided for residues in the D-alpha helix region of Ramachadran space?
Referenced by genkic_close(), initialize_from_options(), and set_D_alpha_compfile_contents().
|
private |
Has an aa_composition setup file ben provided for residues in the D-beta strand region of Ramachadran space?
Referenced by genkic_close(), initialize_from_options(), and set_D_beta_compfile_contents().
|
private |
A default custom Ramachandran table to be used for sampling.
Defaults to unknown_ramatable_type (which means that no custom Rama table is applied by default).
Referenced by default_rama_table_type(), and set_default_rama_table_type().
|
private |
If we're designing, this is the filename for the file that tells us what residues are allowed at what postion.
If left blank, no file is read.
Referenced by initialize_from_options().
|
private |
Should we design the peptide after sampling the conformation, or just relax it?
Default "false" (relax only; no design).
Referenced by genkic_close(), and initialize_from_options().
|
private |
The cutoff dslf_fa13 energy, post-relaxation, above which closure solutions are rejected.
Read from options.
Referenced by genkic_close(), and initialize_from_options().
|
private |
The cutoff dslf_fa13 energy, pre-relaxation, above which closure solutions are rejected.
Read from options.
Referenced by genkic_close(), and initialize_from_options().
|
private |
Prevent counting of hydrogen bonds to adjacent residues?
Default true.
Referenced by initialize_from_options(), and set_up_hbond_filter().
|
private |
The number of FastReleax rounds to be applied whenever FastRelax is used.
Defaults to 3.
Referenced by genkic_close(), initialize_from_options(), and run().
|
private |
Should we filter out solutions that have more than the allowed number of hydrogen bonds to an acceptor?
Default true.
Referenced by genkic_close(), and initialize_from_options().
|
private |
The number of attempts that will be made by the generalized kinematic closure machinery.
Defaults to 1.
Referenced by genkic_close(), and initialize_from_options().
|
private |
The minimum number of solutions that must be found by the generalized kinematic closure machinery before a single solution is chosen.
Defaults to 1.
Referenced by genkic_close(), and initialize_from_options().
|
private |
The hbond energy cutoff above which a hbond is no longer counted.
Defaults to -0.25.
Referenced by initialize_from_options(), and set_up_hbond_filter().
|
private |
For parts of the protocol that upweight the hydrogen bond terms, this is the factor by which these terms are upweighted.
Defaults to 10.0.
Referenced by initialize_from_options(), and run().
|
private |
Has an aa_composition setup file ben provided for residues in the L-alpha helix region of Ramachadran space?
Referenced by genkic_close(), initialize_from_options(), and set_L_alpha_compfile_contents().
|
private |
Has an aa_composition setup file ben provided for residues in the L-beta strand region of Ramachadran space?
Referenced by genkic_close(), initialize_from_options(), and set_L_beta_compfile_contents().
|
private |
Should all cysteine residues be linked with 1,3,5-tris(bromomethyl)benzene?
False by default.
Referenced by initialize_from_options(), and run().
|
private |
The minimum number of mainchain hydrogen bonds that a final solution must have.
Defaults to 0 (report only).
Referenced by initialize_from_options(), and run().
|
private |
The minimum number of mainchain hydrogen bonds that a generalized kinematic closure solution must have.
Defaults to 3.
Referenced by genkic_close(), initialize_from_options(), and run().
|
private |
If this is called by MPI code, this can store the rank of the current process. Zero otherwise.
-------—— Data -------------------------—— --—— When you add new data to this class, -— --—— you must update the copy constructor -—
Referenced by run(), and set_my_rank().
|
private |
List of positions (in original sequence indexing – not permuted) that are N-methylated.
Defaults to empty list.
Referenced by add_n_methylation(), and initialize_from_options().
|
private |
Has the user specified a native structure?
Referenced by initialize_from_options(), run(), and set_native().
|
private |
Filename for the native, if one is provided.
Referenced by initialize_from_options(), and run().
|
private |
A native pose, provided by external code.
If provided, this prevents the app from reading from the filesystem in the run() function.
Referenced by run(), and set_native().
|
private |
Max number of structures to generate.
Referenced by initialize_from_options(), run(), and set_nstruct().
|
private |
The prefix for the output filename.
Defaults to "S_".
Referenced by initialize_from_options(), and run().
|
private |
|
private |
The energy above which we no longer count a hydrogen bond to an acceptor for the filter.
Default -0.1.
Referenced by genkic_close(), and initialize_from_options().
|
private |
Used to prevent this process from directly reading the input file if it has been read by another process and the information has been transmitted to this process.
Default false; set to true by set_allowed_residues_by_position().
Referenced by initialize_from_options().
|
private |
Should D-residues be prohibited at positions with negative phi values?
Default true.
Referenced by initialize_from_options(), and set_up_design_taskoperations().
|
private |
Should L-residues be prohibited at positions with positive phi values?
Default true.
Referenced by initialize_from_options(), and set_up_design_taskoperations().
|
private |
The energy cutoff for the rama filter.
Defaults to 0.3.
Referenced by genkic_close(), and initialize_from_options().
|
private |
Custom Ramachandran tables to be used for sampling, listed by residue.
Defaults to an empty map. Only mapped indicies have custom rama tables applied. Overrides default_rama_table_type_ (if used) at the residue indices in the map.
Referenced by custom_rama_table_defined(), rama_table_type_by_res(), and set_rama_table_type_by_res().
|
private |
The name of the checkpoint file for the random number generator.
Defaults to "rng.state.gz". Read from options.
Referenced by erase_random_seed_info(), get_random_seed_info(), initialize_from_options(), and store_random_seed_info().
|
private |
The cutoff, in degrees, to use when comparing mainchain torsion values to determine whether symmetry repeats are truely symmetric. Defaults to 10 degrees.
Referenced by genkic_close(), initialize_from_options(), and run().
|
private |
If this option is used, then only backbones with mirror symmetry are accepted. Must be used with the -cyclic_peptide:require_symmetry_repeats flag.
Referenced by genkic_close(), initialize_from_options(), and run().
|
private |
The random perturbation, in degrees, to be applied when copying mainchain torsion values to produce symmetric conformations.
Defaults to 0.0 (no perturbation).
Referenced by genkic_close(), and initialize_from_options().
|
private |
If this option is used, then only backbones that are cN (or cN/m, if mirror symmetry is required) symmetric will be accepted. For example, if set to 2, then only c2-symmetric backbones will be accepted.
Referenced by genkic_close(), initialize_from_options(), and run().
|
private |
At residues that are more likely to have cis omega values, sample cis some fraction of the time.
Currently applies only to L-proline and D-proline. Could be extended to beta-3-proline and to N-methylated amino acids, once those are added. Read from database.
Referenced by disable_cis_pro_sampling(), sample_cis_pro(), and set_sample_cis_pro_frequency().
|
private |
At residues that are more likely to have cis omega values, sample cis this fraction of the time.
Read from database.
Referenced by disable_cis_pro_sampling(), sample_cis_pro_frequency(), and set_sample_cis_pro_frequency().
|
private |
The default ScoreFunction to use. The high h-bond version is constructed from this.
Referenced by initialize_from_options(), run(), set_scorefxn(), and SimpleCycpepPredictApplication().
|
private |
Filename for the text file containing the sequence of the peptide.
Must be provided with the -cyclic_peptide:sequence_file flag.
Referenced by initialize_from_options(), and run().
|
mutableprivate |
The length of the sequence, excluding any cross-linkers.
Computed internally and stored for reference. This is also the index of the last sequence residue (with cross-linkers following in linear sequence).
Referenced by run(), and sequence_length().
|
private |
The string that would be read from a sequence file.
If provided by external code, prevents filesystem read in run() function.
Referenced by read_sequence(), and set_sequence().
|
private |
Should this application produce silent file output?
Referenced by initialize_from_options(), run(), and set_silentstructure_outputlist().
|
private |
If silentlist_out_ is used, this is the list to append SilentStructurOPs to.
Note that this is a non-const pointer to a vector.
Referenced by run(), and set_silentstructure_outputlist().
|
private |
Should this application produce silent structure OP output, appending to a list?
Referenced by run(), and set_silentstructure_outputlist().
|
private |
If silentlist_out_ is used, this is the list of job summaries to which summaries should be appended.
Note that this is a non-const pointer to a vector.
Referenced by run(), and set_silentstructure_outputlist().
|
private |
Allows external code to suppress checkpointing, so that the SimpleCycpepPredictApplication doesn't write directly to disk.
Referenced by checkpoint(), end_checkpointing(), erase_random_seed_info(), initialize_checkpointing(), new_checkpoint_file(), set_suppress_checkpoints(), and store_random_seed_info().
|
private |
Multiplier to make the TBMB constraints energy filter more permissive.
Default 1.0.
Referenced by genkic_close(), and initialize_from_options().
|
mutableprivate |
List of positions linked by 1,3,5-tris(bromomethyl)benzene.
This is a vector of lists of three residues.
Referenced by depermute(), genkic_close(), initialize_from_options(), re_append_tbmb_residues(), and run().
|
private |
Multiplier to make the TBMB distance filter more permissive.
Default 1.0.
Referenced by genkic_close(), and initialize_from_options().
|
private |
The total energy cutoff above which solutions are discarded.
Defaults to 0.0.
Referenced by run(), and set_total_energy_cutoff().
|
private |
Should we try all disulfide combinations during structure prediction?
Default false. Read from options.
Referenced by genkic_close(), and initialize_from_options().
|
private |
Should we use the aa_composition score term in design?
Default false.
Referenced by initialize_from_options().
|
private |
Should solutions be filtered by Ramachandran energy (in the case of alpha-amino acids)?
Defaults to true.
Referenced by initialize_from_options(), and use_rama_filter().
|
private |
Should we use rama_prepro tables for sampling?
Default true.
Referenced by set_use_rama_prepro_for_sampling(), and use_rama_prepro_for_sampling().
|
private |
If true, filters are applied based on distance between TBMB cysteines and on constraints to discard GenKIC solutions that can't be crosslinked easily.
True by default.
Referenced by genkic_close(), and initialize_from_options().
|
private |
Determines whether the total energy cutoff should be used.
Defaults to 'false'.
Referenced by disable_total_energy_cutoff(), run(), and set_total_energy_cutoff().
|
private |
Map of (seqpos -> phi/psi/omega triplet).
The user can set certain alpha-amino acid mainchain dihedral values, if he or she so wishes.
Referenced by genkic_close(), and initialize_from_options().
|
private |
A small random value added to all user-set dihedral values. Defaults to 0.
Referenced by genkic_close(), and initialize_from_options().