You are here

params file for beta-amino acids

5 posts / 0 new
Last post
params file for beta-amino acids
#1

Hi all,

I've been trying to insert a beta-amino acid into peptide, let's say, beta-homo-leucine. It seems that a params file for beta-hLeu is needed. Unforutnately, there is no params file for beta-hLeu in database in Rosetta, and molfile_to_params.py only works for alpha-amino acid. I've also visited https://rosie.rosettacommons.org/documentation/beta_peptide_design, but did't find helpful solution. So my question is:

Is there any script for generating params file for beta-amino acid?

Any helpful information is appreciated. Thank you!

 

Post Situation: 
Mon, 2021-03-29 18:30
phonez

When you say beta-homo-leucine, do you mean "leucine, but there's an extra backbone carbon atom" (which happens to be chemically equivalent to "homoleucine, but the nitrogen is bonded to the beta carbon")? In that common case B3L is the parameter set you want.

Separately, molfile_to_params_polymer.py definitely works for betas.

Mon, 2021-03-29 20:17
everyday847

Thank you for your reply! 

The "beta-homo-leucine" I'm talking about is indeed the same as B3L in Rosetta database, and I was wrong to assume they're different things just by the way they're named.

And to learn to generate params file for beta-aa, I took B3L as an example. Compared with B3L.params in Rosetta database, my result:
1. got wrong CHIs(containing dihedral angles of backbone), so I deleted them manually. I wander if it's ok to do that.
2. had a different atom name for alpha-C, which in my result was CB1 but CM in Rosetta database. I don't know how the atom tree is spanned in the latter and if it will matter.
 
Anoter question is: I used FlexPepDock Refinement Protocol with a peptide containing B3L as input. But it got an error:
'''
core.chemical.GlobalResidueTypeSet: Loading (but possibly not actually using) 'B3L' from the PDB components dictionary for residue type 'pdb_B3L'
core.conformation.Residue: [ ERROR ] New Residue cannot be aligned with template Residue, because the numbers of main-chain atoms differ.
protocols.jd2.JobDistributor: [ ERROR ] 
'''
I tried to substitute B3L with a natural amino acid and it worked. So I guessed this FlexPepDock Refinement Protocol is not compatible with beta-aa. How should I modify this protocol for beta-aa? Or change with other Dock Protocols?

Thank you!

Wed, 2021-03-31 22:04
phonez

That might be a research question for flexpepdock! I assume they're trying to make fragment moves and using alpha amino acid fragments. That's going to be really hard.

 

I am tempted to say "write your own Monte Carlo protocol" but that is a tall order. What's your scientific use case?

Thu, 2021-04-01 11:58
everyday847

Oh and yes -- atom naming is expected to be inconsistent with the raw molfile to params output. You do have to make some manual changes to echo the rosetta params style. At least it calculates all the internal coordinates right!

Thu, 2021-04-01 11:59
everyday847