Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

The results ligand is different from SMILES I have entered #189

Open
MinsikBae opened this issue Feb 20, 2025 · 1 comment
Open

The results ligand is different from SMILES I have entered #189

MinsikBae opened this issue Feb 20, 2025 · 1 comment

Comments

@MinsikBae
Copy link

I am experiencing a discrepancy between the ligand structure I expect based on my input SMILES and the 3D structure generated by Boltz. My SMILES string specifies an alkane chain (CCCCCCCN) as a linker, but the output structure has an alkene (C=C) in that linker(1 to 3 alkene generated), and the molecule forms an unexpected macrocycle.

I have already attempted the following troubleshooting steps:

Used the MMFF94 force field instead of UFF.
Modified the compute_3d_conformer function in schema.py to check for and convert unexpected double bonds to single bonds.
Verified that useRandomCoords is effectively disabled.
Checked input SMILES string
Adjust recycling_steps , sampling_steps, step_scale

@JesusSeco
Copy link

I had a very similar problem. Using canonical SMILES worked, but the most effective solution, is to turn off in schema.py the Sanitize option to False, to avoid rdkit to change the SMILES accodingly to sanitization rules

Additionally, you can try to use the implicit smiles syntax as below. I tested your smiles with random sequence and the output worked for me

  • ligand:
    id: B
    smiles: '[CH3][CH2][CH2][CH2][CH2][CH2][CH2]N'

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants