Open Cucco23 opened 6 months ago
Certainly! You can dock multiple ligands (proteins B, C, D, E, etc.) to a fixed receptor (protein A) using molecular docking software that supports docking multiple ligands simultaneously in a single run. Prepare your receptor and ligands, create input files specifying them, and run the docking simulation. Analyze the results afterward. This approach saves time and resources compared to docking each ligand individually.
import subprocess
def run_autodock_vina(receptor_file, ligand_files, output_dir): for ligand_file in ligand_files: output_file = f"{outputdir}/output{ligand_file.split('.')[0]}.pdbqt" command = f"vina --receptor {receptor_file} --ligand {ligand_file} --out {output_file} --log {output_dir}/log.txt" subprocess.run(command, shell=True)
if name == "main": receptor_file = "receptor.pdbqt" ligand_files = ["ligand1.pdbqt", "ligand2.pdbqt", "ligand3.pdbqt"] output_dir = "dock_output"
run_autodock_vina(receptor_file, ligand_files, output_dir)
Thank you for the answer.
I meant by using alphafold multimer.
Moreover, are you sure you can dock proteins with vina? Where is this publushed? To my knowledge it can be used only with small molecules since it test conformations.
On another note lets say af multimer provides me a model. How do you test the quality of the docking?
Thanks
From: Tejas Amol Hande @.> Sent: Friday, March 8, 2024 1:47 PM To: google-deepmind/alphafold @.> Cc: Cocucci, Ema @.>; Author @.> Subject: Re: [google-deepmind/alphafold] Run multimer with a "single receptor" and "multiple partners" (Issue #902)
Certainly! You can dock multiple ligands (proteins B, C, D, E, etc. ) to a fixed receptor (protein A) using molecular docking software that supports docking multiple ligands simultaneously in a single run. Prepare your receptor and ligands,
Certainly! You can dock multiple ligands (proteins B, C, D, E, etc.) to a fixed receptor (protein A) using molecular docking software that supports docking multiple ligands simultaneously in a single run. Prepare your receptor and ligands, create input files specifying them, and run the docking simulation. Analyze the results afterward. This approach saves time and resources compared to docking each ligand individually.
import subprocess
def run_autodock_vina(receptor_file, ligand_files, output_dir): for ligand_file in ligand_files: output_file = f"{outputdir}/output{ligand_file.split('.')[0]}.pdbqt" command = f"vina --receptor {receptor_file} --ligand {ligand_file} --out {output_file} --log {output_dir}/log.txt" subprocess.run(command, shell=True)
if name == "main": receptor_file = "receptor.pdbqt" ligand_files = ["ligand1.pdbqt", "ligand2.pdbqt", "ligand3.pdbqt"] output_dir = "dock_output"
run_autodock_vina(receptor_file, ligand_files, output_dir)
— Reply to this email directly, view it on GitHubhttps://urldefense.com/v3/__https://github.com/google-deepmind/alphafold/issues/902*issuecomment-1986228464__;Iw!!KGKeukY!2dhx2hKwyp9ZINvsfetg2Gnj7QIfGEMDfgPIYbTaK1jpxo3jIMeMmkPD7dj6QL6LVVMeZhagh8GnrTN_amjy-ooJqw$, or unsubscribehttps://urldefense.com/v3/__https://github.com/notifications/unsubscribe-auth/BEVFMI6EQMDJRIGPYOMAYLDYXIBVVAVCNFSM6AAAAABD2ZJCLCVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSOBWGIZDQNBWGQ__;!!KGKeukY!2dhx2hKwyp9ZINvsfetg2Gnj7QIfGEMDfgPIYbTaK1jpxo3jIMeMmkPD7dj6QL6LVVMeZhagh8GnrTN_ami0rxQjGw$. You are receiving this because you authored the thread.Message ID: @.***>
The Alphafold source code does not support this feature, but my solution is manually generates MSA and imports it into the corresponding output folder. The simplest solution is to first calculate A-B, copy their msas folder to the output folder of A-C, but please note that it is important to ensure that the protein A is in the first position in the fasta file and delete the B folder inside the msas folder. This will automatically read the MSA of protein A and calculate the MSA of protein C when you start with --use_precomputed_msas flag. @Cucco23
Thank you very much for the suggestion! I will try this.
From: Xiao Guo 2208401002 @.> Sent: Monday, March 11, 2024 10:30 PM To: google-deepmind/alphafold @.> Cc: Cocucci, Ema @.>; Mention @.> Subject: Re: [google-deepmind/alphafold] Run multimer with a "single receptor" and "multiple partners" (Issue #902)
The Alphafold source code does not support this feature, but my solution is manually generates MSA and imports it into the corresponding output folder. The simplest solution is to first calculate A-B, copy their msas folder to the output folder
The Alphafold source code does not support this feature, but my solution is manually generates MSA and imports it into the corresponding output folder. The simplest solution is to first calculate A-B, copy their msas folder to the output folder of A-C, but please note that it is important to ensure that the protein A is in the first position in the fasta file and delete the B folder inside the msas folder. This will automatically read the MSA of protein A and calculate the MSA of protein C when you start with --use_precomputed_msas flag. @Cucco23https://urldefense.com/v3/__https:/github.com/Cucco23__;!!KGKeukY!xXcOSCdA3x5UiTa2VjDKKVEM6SobW7oqY8RMv4qk2x_SPXSvLjWXY7edbDdPmzqzjy-1r856ZX5qCeVFdN16QBETTQ$
— Reply to this email directly, view it on GitHubhttps://urldefense.com/v3/__https:/github.com/google-deepmind/alphafold/issues/902*issuecomment-1989860542__;Iw!!KGKeukY!xXcOSCdA3x5UiTa2VjDKKVEM6SobW7oqY8RMv4qk2x_SPXSvLjWXY7edbDdPmzqzjy-1r856ZX5qCeVFdN0zsdGJ2g$, or unsubscribehttps://urldefense.com/v3/__https:/github.com/notifications/unsubscribe-auth/BEVFMI25A5H72B4NRBPRS6TYXZSCFAVCNFSM6AAAAABD2ZJCLCVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSOBZHA3DANJUGI__;!!KGKeukY!xXcOSCdA3x5UiTa2VjDKKVEM6SobW7oqY8RMv4qk2x_SPXSvLjWXY7edbDdPmzqzjy-1r856ZX5qCeVFdN0uT5eyLQ$. You are receiving this because you were mentioned.Message ID: @.**@.>>
I believe the AlphaPulldown repository (by the Kosinski lab) seems to match what you want to use AF to do
Hallo all,
Let's say that I have a protein to which I would like to dock multiple different proteins but always in a 1:1 ratio. there is a way to provide a list in the fasta file so that the first protein (receptor) does not require to be re-calculated? e.g. protein A with proteins B, C, D, E....
I hope my question was clear and that someone has a good suggestion
Thank you!