- modeller_usage - salilab.org

(no subject)
by Fred Pontes 14 Apr '25

14 Apr '25

Hi! I'm using MODELLER to build a homology model of my protein. In my first try, I create a model without any restraint, and then I compared the secondary structure of this model (Using STRIDE and VMD) with the secondary structure predict by PSIPRED webserver. For the first 26 residues, I have these predictions: PSIPRED: CHHHHHHHHHHHHHHHHHHHHHHHHC Homology model: CCCCCCCCCCCCHHHHHHHHHHHHHC My idea is restraint the region where the predictions are coincidences (13:A to 26:A) and applying … [View More]

2 1

Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES( 1)
by Fred Pontes 01 Apr '25

01 Apr '25

Hi! Sorry for the very basic question, but again I'm facing issues with MODELLER. When I run MODELLER I got the following error: "*_modeller.ModellerError: read_al_373E> Protein specified in ALIGN_CODES(i) was not found in the alignment file; ALIGN_CODES( 1) = 6KG7.pdb*" But my structure file is present in my directory: *(base) [medfpo@uol-pc-301925 hPIEZO2_structure]$ ls -l 6KG7.pdb -rw-r--r--. 1 medfpo Domain Users 1177728 Mar 31 13:25 6KG7.pdb* I don't have any idea what's … [View More]

3 3

(no subject)
by Fred Pontes 27 Mar '25

27 Mar '25

Hi everyone, I'm a beginner in using MODELLER and I need help to fix my files in my modelling. I'm trying to build a homology model using two templates. I'd like to select different regions of each template to build my model. I attached the script below, In my alignment.pir file I include the sequence of my template in the following way: >P1;template1 structureX:protA.pdb:8:A:2822:A:::-1.00:-1.00 --------------------------------------------------------------------------- ... >P1;… [View More]

2 1

what is non-bonded restraints in modeller molpdf
by Li Xue 17 Feb '25

17 Feb '25

Hello all, I am trying to understand the components of modeller molpdf. Does anyone know what the non-bonded restraints are, since modeller by default does not calculate L-J nor electrostatics? Best, Li

2 1

how to calculate energy using modeller
by Li Xue 21 Jan '25

21 Jan '25

Dear Modeller Team, I am currently looking for a quick method for calculating the energy of protein structures. I understand that Modeller utilizes a simplified version of the CHARMM force field, which suits my needs, as I’m looking for a fast and approximate energy calculation. I would like to print out each energy term, but with the current code I’m using, the Lennard-Jones and Coulomb terms are returning zero. Could you please advise on how to resolve this issue? The code that I used: … [View More]

2 1

Reconstruction of CA position only files with all-atom template
by jingkaizeng 13 Nov '24

13 Nov '24

Hi all, Here is my task: I have some CA position only files of a protein of interested that underlying some important biological pathways. We also have the protein full Cryo-EM structures (pdb:8io4, the exact same protein as well as sequence as the CA only files). Then I would like to fill the keep the CA position only files' CA position (don't move at all) then reconstruct the detail like side chain/secondary structure information from template. So this task is not homology modeling but … [View More]reconstruction. After getting the reconstructed files, I would like to use these into MD simulations. I searched from previous mail list, and I found one did this task but not based on template. https://www.salilab.org/archives/modeller_usage/2008/msg00285.html I tried with " env = Environ() # Load the CA-only seed model seed_model = Model(env, file=seed_file) # Load the template structure (with full atomistic detail) template_model = Model(env, file=template_file) # Create an alignment between the seed and itself (1:1 alignment) aln = Alignment(env) aln.append_model(seed_model, align_codes='seed') aln.append_model(seed_model, align_codes='seed') # Write the alignment for debugging purposes aln.write(file='alignment.ali', alignment_format='PIR') # Define a custom model class to fix the Cα positions class MyModel(automodel): def select_atoms(self): # Select all atoms except the Cα atoms to allow refinement of everything but Cα s = selection(self) return s - s.only_atom_types('CA') # Set up environment for refinement env.edat.nonbonded_sel_atoms = 2 # Disable interactions between selected and unselected atoms a = MyModel(env, alnfile = aln, # Alignment file created earlier knowns = 'seed', # Seed model (CA-only) sequence = 'seed') # Target is also the seed (1:1 alignment) a.starting_model = 1 a.ending_model = 1 a.make() " But this code didn't use the side chain atom from the Cryo-EM structure, and just use CA position to reconstruct the detail. It seems the reconstruction by this code is not reliable enough that cause infinite energy error in MD simulation, it may indicate some overlap or incorrect side-chain orientations, or improper hydrogen placements. Therefore, I would need to try to reconstruct with the side chain information from template Cryo-EM structure's pdb. I tried some methods that use the template Cryo-EM structure, but it seems we just do the homology modeling but NOT reconstruction. The production files we get is almost same as the template without make use of the CA position from CA only files. The code I used like : " # Load the CA-only seed model seed_model = Model(env, file=seed_file) # Load the template structure (with full atomistic detail) template_model = Model(env, file=template_file) # Create an alignment between the seed and the template aln = Alignment(env) aln.append_model(template_model, align_codes='template') aln.append_model(seed_model, align_codes='seed') # Perform the alignment aln.align() # Now defining a custom AutoModel class to apply Cα restraints class MyModel(AutoModel): def special_restraints(self, aln): # Restrain all Cα atoms based on the seed model rsr = self.restraints atmsel = selection(self).only_atom_types('CA') # Apply restraints to the selected Cα atoms for atom in atmsel: rsr.add(forms.Gaussian(group=physical.xy_distance, feature=features.Distance(atom, atom), mean=0.0, # No change in CA positions stdev=0.01)) # Instantiate the custom model class a = MyModel(env, alnfile=aln, knowns='template', sequence='seed') a.starting_model = 1 a.ending_model = 1 a.make() " I think this task can be done since I found a paper (https://doi.org/10.7554/eLife.68369) that do the similar task, they claim to use Modeller to perform this: "First, side chain atoms from the template X-ray structure (PDB ID 4HFI) were added to each model, followed by a cycle of refinement with all Ca atoms restrained. Restraints on Ca atoms were then substituted with restraints on backbone hydrogen bonds, taken from helix and sheet annotations in the template PDB file, for another cycle of refinement to ensure proper secondary structure." Can anyone give some suggestions about this reconstruction task? Thanks. Best, Jingkai [View Less]

2 1

Missing residues modeling in Modeller
by bbmpresi＠gmail.com 04 Nov '24

04 Nov '24

Dear all, I am trying to model the missing residues in a beta barrel shaped protein that contains 12 missing residues in the original pdb. As instructed for the 1qg8 tutorial, I also first obtained the full amino acid sequence of the protein using: from modeller import * # Get the sequence of the 1qg8 PDB file, and write to an alignment file code = 'pdbid' e = Environ() m = Model(e, file=code) aln = Alignment(e) aln.append_model(m, align_codes=code) aln.write(file=code+'.seq') and I am … [View More]obtaining an output file named pdbid.seq The actual protein sequence contains 324 amino acid residues and 12 residues are missing in the pdbid.pdb file This is the alignment.ali file I am preparing for missing residue addition: >P1;pdbid structureX:pdbid:1:A:+312:A:MOL_ID 1; ASDQRGYKP------------GGHVGTSVEYEDKVTRGFNNTDKKEKTITNEVFNFFYNNPQWNFMGFYSFKIENREQKEPGYYENEDGIKQLFSLNKGHDLGNGWATGLIYELEYTRSKVYSPDVSGLRKNLAEHSIRPYLTYWNNDYNMGFYSNLEYLLSKEDRNAWGKRQEQGYSALFKPYKRFGNWEVGVEFYYQIKTNDEKQPDGTINEKSDFNERYIEPIVQYSFDDAGTLYTRVRVGKNETKNTDRSGGGNAGINYFKDIRKATVGYEQSIGESWVAKAEYEYANEVEKKSRLSGWEARNKSELTQHTFYAQALYRF* >P1;pdb_fill sequence::::::::: ASDQRGYKPEDVAFDESFFSFGGHVGTSVEYEDKVTRGFNNTDKKEKTITNEVFNFFYNNPQWNFMGFYSFKIENREQKEPGYYENEDGIKQLFSLNKGHDLGNGWATGLIYELEYTRSKVYSPDVSGLRKNLAEHSIRPYLTYWNNDYNMGFYSNLEYLLSKEDRNAWGKRQEQGYSALFKPYKRFGNWEVGVEFYYQIKTNDEKQPDGTINEKSDFNERYIEPIVQYSFDDAGTLYTRVRVGKNETKNTDRSGGGNAGINYFKDIRKATVGYEQSIGESWVAKAEYEYANEVEKKSRLSGWEARNKSELTQHTFYAQALYRF* But in the template I get residue number written as 312 instead of 324, I directly copied the entry from the pdbid.seq file that was provided by Modeller at the first step that I explained at the beginning. Due to 12 missing residues Modeller is reporting 312 residues in the second line of the alignment.ali file: structureX:pdbid:1:A:+312:A:MOL_ID 1; However, I checked other examples on Modeller tutorial and it seems in some examples the total residue numbers (including the missing residue numbers) are written in the template section for the alignment.ali file, so I could have used structureX:pdbid:1:A:+324:A:MOL_ID 1; writing 324 as total number of residues for the template; so my question is that does it matter if write structureX:pdbid:1:A:+312:A:MOL_ID 1; or structureX:pdbid:1:A:+324:A:MOL_ID 1; I generated a new pdb with all missing residues modelled by Modeller using: from modeller import * from modeller.automodel import * # Load the AutoModel class log.verbose() env = Environ() # directories for input atom files env.io.atom_files_directory = ['.', '../atom_files'] class MyModel(AutoModel): def select_atoms(self): return Selection(self.residue_range('10:A', '21:A')) #a = MyModel(env, alnfile = 'alignment.ali', # knowns = 'pdbid', sequence = 'pdbid_fill') a = AutoModel(env, alnfile = 'alignment.ali', knowns = 'pdbid', sequence = 'pdbid_fill') a.starting_model= 1 a.ending_model = 5 a.make() I checked that using structureX:pdbid:1:A:+312:A:MOL_ID 1; or structureX:pdbid:1:A:+324:A:MOL_ID 1; in the second line of alignment.ali for the template, I am obtaining a very similar or more or less identical modeled structure. Can anybody tell me if the number of residues written in the second line of alignment.ali file for template pdb matters or not? Or only the template sequence and the gaps really matter? Any help would be much appreciated, thank you [View Less]

2 1

Fwd: Average the models
by Marco Sette 24 Oct '24

24 Oct '24

Thanks Ben, is there a way to avoid that an insertion from becoming a curious loop like the one shown at the top in the figure? I’m trying to relax the system using MD but maybe I can avoid it during modeling. Best, Marco Il 22/10/2024 01:26, Modeller Caretaker ha scritto: > On 10/21/24 2:28 AM, Marco Sette via modeller_usage wrote: >> I have obtained five models of a protein with Modeller and I am >> analyzing them in relation to the wt. There are of course small >&… [View More]

2 1

D-Ala
by Xiaosa Wu 22 Oct '24

22 Oct '24

Hi all, I tried to use modeller10.4 to build a homology model for a peptide, which contains a D-Ala. Here the following steps I did: 1: in modlib/top.lib, Copy topology entry of L-Ala, then change From IC N C *CA CB 1.4592 114.4400 123.2300 111.0900 1.5461 IC N C *CA HA 1.4592 114.4400 -120.4500 106.3900 1.0840 To IC N C *CA CB 1.4592 114.4400 -123.2300 111.0900 1.5461 IC N C *CA HA 1.4592 114.4400 … [View More]

2 1

Modeller 10.6 release
by Modeller Caretaker 22 Oct '24

22 Oct '24

A new version of Modeller, 10.6, is now available for download! Please see the download page at https://salilab.org/modeller/ for more information. If you have a license key for Modeller 9 or 10, there is no need to reregister for Modeller 10.6 - the same license key will work. (It won't do any harm to reregister if you want to, though!) 10.6 is a minor update that adds support for Python 3.13, and modifies the behavior of Alignment.append_sequence() to make it easier to construct … [View More]

1 0