May be you try Marvin Sketch or Chemsketch. But why ligand -ligand docking?
@nurnadiah44905 күн бұрын
Good day Dr. I want to ask if multiple ligand can interact with receptor in autodock 4?
@BioinfoCopilot5 күн бұрын
No, only possible with Vina.
@nurnadiah44905 күн бұрын
@@BioinfoCopilot thank you Dr. may I ask for docking between protein-protein interaction, which autodock is suitable?
@BioinfoCopilot5 күн бұрын
Cluspro is suitable for protein-protein docking
@nurnadiah44905 күн бұрын
@@BioinfoCopilot Ok thank you for the respons
@pabitra.835 күн бұрын
please make a video to implement hse06. Also please let me know whether HONPAS is compatible to siesta 4.1 b4 or not.
@BioinfoCopilot5 күн бұрын
I will try to. Thanks
@SirwanMansoori5 күн бұрын
Thanks. ☘
@BioinfoCopilot5 күн бұрын
My pleasure
@pabitra.836 күн бұрын
How to implement hse06 here?
@BioinfoCopilot5 күн бұрын
Check the latest siesta documentation
@Scienceiirwn7 күн бұрын
Would it be possible to show a tutorial on how to make heatmap or UMAP/TSNE using the data? That would be very helpful
@BioinfoCopilot7 күн бұрын
Yes that’s a great suggestion. I will try to make one
@ES-yd1ze8 күн бұрын
So diff dock can achieve by Python shell script or by gui of nvidia ? and same result ?
@BioinfoCopilot8 күн бұрын
Both are same!
@ES-yd1ze8 күн бұрын
What is the difference between rmsd of ligand with and without h heavy atom ?
@BioinfoCopilot5 күн бұрын
The calculation is done for ligand without having H atoms called non-hydrogen atoms. Heavy atoms is nothing but the ligand without H atoms. Difference is that when you include H atoms your RMSD values go crazy because of H atoms.
@skyshines901610 күн бұрын
Hello Sir I want to ask that in papers the values are written as for example: -7.8 ± 2.3 , so I would like to know that how can we write the values from the results like this? Do we use SD or SEM? and there are two columns of each so which values we should use. Kindly help me
@BioinfoCopilot10 күн бұрын
Repeat for many intervals and then you get the values. Do some statistical analyses to get the SD values.
@sethuraman887211 күн бұрын
I have a huge interest in structural bioinformatics and i am planning to pursue a PhD in this field along with generative drug design. this channel is very helpful and i came across all of your videos, very useful resources. please post regularly so that a beginner can learn from a pro in this field. thank you!
@BioinfoCopilot10 күн бұрын
Much appreciated and thank you so much 😊
@olgalanzetta12 күн бұрын
Hello, how can I save the flowset with the cells that are only includes in the gates?
@BioinfoCopilot5 күн бұрын
Use cyto_save rdrr.io/github/DillonHammill/CytoExploreR/man/cyto_save.html
@pabitra.8314 күн бұрын
What if i use two types of atoms and use hybrid functional for both, how to define the species then?
@BioinfoCopilot14 күн бұрын
You can define different atoms in hybrid functionals in the fdf file. Check my previous tutorials how to define a fdf file
@pabitra.8314 күн бұрын
@@BioinfoCopilot ok
@hoangmaihuy897916 күн бұрын
I don't think there is a parameter called --log in autodockvina, when i ran this sh script it showed "Command line parse error: unrecognised option '--log' ". I use autodock vina version 1.2.5
@BioinfoCopilot16 күн бұрын
Check the manual “Summary Usage” section vina.scripps.edu/manual/
@mohamedmarzouk253718 күн бұрын
thank you which mac processor you are using ? Is it M2 or other M series or intel versions? I am trying to create the environment but it is not working !
@BioinfoCopilot18 күн бұрын
I use M2 and intel both. Works fine for me.
@mdshakilhossen275519 күн бұрын
Very clear and well understood explanation. How to run simulation with 10 chunks for 100ns simulation. Please make a video or guide me. I tried it for several days but coudn't get over the line.😢
@syarifuddinhusain-x3m21 күн бұрын
Thank you for the video. I have a question if i run a 100 ns simulation but i want to calculate only the last 10 ns, what is the startframe and endframe should i use?
@BioinfoCopilot20 күн бұрын
Thank you 😊. 100ns if you get 1000 frames, then 900 to 1000 will be your frames for last 10ns.
@simmim31021 күн бұрын
How to make size x y z .. is it standard value
@BioinfoCopilot21 күн бұрын
No you have to set it using MGLtools
@simmim31019 күн бұрын
@@BioinfoCopilot I have many different ligands and one protein.. each ligand will interact at some different place.. there are many active sites present on protein,,,,,so how can I define coordinates .. xyz
@roberthill537322 күн бұрын
Nice presentation! These days, only programs that harness the powerful GPU could be popular since over 98% of world compute power comes from GPU in the current days. Could you introduce another one with use of GPU?
@BioinfoCopilot22 күн бұрын
Thank you 😊. Yes, I will try my best to showcase the power of GPU using maybe any online resources. Thanks
@rimabd609623 күн бұрын
Thanks 🙏 you are the best 🎉 I have a problem.. after docking the ligand was fragmented.. what I have to do
@BioinfoCopilot23 күн бұрын
Thank you 😊. You have carefully optimized the ligand using Avogadro software and design it using Marvin Sketch. Then do docking
@MuhannadDehni24 күн бұрын
thank you for great work, I notice from your result Gibbes energy landscape start from 0 Kj/mol. and I try to do that I get the same ,why we didn't get negative value even MMPBSA is = -9 kcal/mol
@BioinfoCopilot22 күн бұрын
The tutorial shown was for demonstration purposes. It might vary from specifics and your simulation environment.
@zechariaholuwadamilola308925 күн бұрын
Thank you very much for this session sir. However, i would like to ask, is it possible to use files (xtc, tpr , ndx) generated from version 2024 of gromacs in this MMPBSA analysis?
@BioinfoCopilot25 күн бұрын
Thank you 😊. Yes you can use those files for MMPBSA analysis.
@evalopez817626 күн бұрын
Good afternoon, I enter the mgltools page but I do not get the download links, is there any alternative to download it? I don't know if they have been removed or it's my problem. Thanks
@BioinfoCopilot26 күн бұрын
ccsb.scripps.edu/mgltools/downloads/
@evalopez817624 күн бұрын
I can install it now, i think it was a problem with my computer. Thanks
@misaelbermudez283927 күн бұрын
Thank you so much for the video. One question, blind docking can be used also to see covalent interactions between a ligand and a protein ?.
@BioinfoCopilot27 күн бұрын
😊 Thank you. The main driving forces between protein and ligand is non-covalent interactions. For covalent interactions you might need to consider Gromacs to run your simulation (if there are any)
@muhriddinmahkamov839228 күн бұрын
Sir, CGenFF server is no longer freely available . I recommend ATB, Acpype or such kind of other servers.
@yunheekim886229 күн бұрын
Thank you so much for your great teaching. It is very helpful for me.
@BioinfoCopilot28 күн бұрын
@@yunheekim8862 My pleasure 😇 and thanks for your support
@nhimcaycu61729 күн бұрын
Hi Sir, Your videos are very useful and engaging. However, using the mdp file from OPLS to simulate a protein in water with the CHARMM36 force field seems inappropriate, as the OPLS mdp file lacks parameters for the CHARMM36 force field. Below are the missing parameters: constraints = h-bonds cutoff-scheme = Verlet vdwtype = cutoff vdw-modifier = force-switch rlist = 1.2 rvdw = 1.2 rvdw-switch = 1.0 coulombtype = PME rcoulomb = 1.2 DispCorr = no
@BioinfoCopilot29 күн бұрын
Thanks for your insights. My tutorials are meant for demonstration purposes. Users can define their own parameters and tune in to adjust to their specific requirements. However thanks for pointing it out!
@emanueleraggi272Ай бұрын
Thank you very much for sharing your knowledge on siesta! In this tutorial you mentioned that other than the molecule (NO2) you could also put other molecules or even DNA. How do you put DNA? What would that be as a representation and how would you do that? Thank you very much!!
@BioinfoCopilotАй бұрын
Thank you 😊. Yes you can put individual DNA nucleobases or pairs based on your setup. You need to optimize the structures first using Siesta and then place on top of your desired material using ASE. If you want to find the optimum configuration for your DNA bases, you can also perform NEB using Siesta and or VASP.
@emanueleraggi272Ай бұрын
@@BioinfoCopilot Thank you very much for the quick reply and for the suggestions. I guess that my biggest problem now is how do I represent a DNA nucleotide base (e.g. ATTAGCCG) in Siesta instead of using a molecule of NO2 like you did on the video? Have you ever done that before? And could you point me to any specific direction to do that? On the Siesta manual there is no clue and on the Siesta website they don't have specific procedure for that...
@arijitsamanta3660Ай бұрын
Sometimes the automation is stopped due to maxwarn warning. How can i resume the automation?
@BioinfoCopilotАй бұрын
You can set maxwarn= 10 in parfile.par
@BioinfoCopilotАй бұрын
And then you can use chaperong options to resume your automation. Explore the options in chaperong.
@AYESHASOHAIL-l9bАй бұрын
I want to know how you install windows on mac book and use both side by side
@BioinfoCopilotАй бұрын
You can use Parallels to install Windows and Ubuntu as well and use it alongside Mac
@amisupriyo1Ай бұрын
Hello Pritam in my case sh script.sh command not working. IT show bad for loop variable. Please help
@BioinfoCopilotАй бұрын
First do chmod +x script.sh Then try to execute it from the parent folder.
@erenkasimfirtina7870Ай бұрын
I've been using the old Avogadro version for a while, according to the one you're using, and I've been encountering too many mistakes. So I need a new version of Avogadro; could you help me find the download file?
@BioinfoCopilotАй бұрын
Check this one: two.avogadro.cc/install/index.html Or sourceforge.net/projects/avogadro/files/latest/download
@KamalSingh-dn7gvАй бұрын
Hello Pritam, Introducing Nvidia's DiffDock, which uses generative AI for the docking of small molecules is indeed a promising idea. However, there are some inaccuracies in your video that should be addressed. It's important to ensure the information shared is scientifically sound, as misleading content can misinform viewers who may trust your statements. Specifically, there are two points of concern regarding your docking of remdesivir into the structure of SARS-CoV-2 Mpro: (i) Remdesivir is not a substrate of Mpro. (ii) The form of remdesivir you docked is not its active form. It is a prodrug form that requires metabolic conversion to a form that is accepted as a substrate by the target enzyme. I recommend reviewing these aspects and conducting further research to ensure accuracy in future content. It's essential to have a solid understanding of the science behind such AI techniques to maximize the potential of tools like Nvidia DiffDock and ensure their effective use. Otherwise, these AI tools can also be considered useless! Thank you - Kamal
@BioinfoCopilotАй бұрын
Hello Kamal, Thank you for your insights. The tutorials I made so far related to docking or dynamics are for demonstration purposes. The tutorials must not be taken into consideration for replicating or to publish etc. I think people who watch my videos know it very well. I don’t have to put a disclaimer or something in the beginning of the videos. Thank you 😊
@thebestusernameevrАй бұрын
How necessary do you think a Ph.D. is in the US? Or do you think a master's is sufficient?
@BioinfoCopilotАй бұрын
PhD from Europe and Job in US is the perfect combination.
@er.gautamkumarjaiswal8664Ай бұрын
Dear Why you are not describe how to make .mpr file?
@BioinfoCopilotАй бұрын
@@er.gautamkumarjaiswal8664 I have described .mpr file at 6:13
@ES-yd1zeАй бұрын
I noticed that you have linux system but you are using windows version of swiss PDB viewer Why not not using Linux version? Excuse me sir but I have same issue
@BioinfoCopilotАй бұрын
I have Mac not Linux. That’s why I use wine to simulate windows inside Mac.
@penscriber9320Ай бұрын
Its really cool. Thanks
@BioinfoCopilotАй бұрын
My pleasure
@erenkasimfirtina7870Ай бұрын
When I try to open the pdbqt file of ligand conformations using Discovery Studio, the conformations of the ligand are looking destroyed. So, the main chain of the ligand is not looking truly. As if molecule is containing Siyano, etc. groups. I just prepared the ligand using Autodock 4. The same problem is popping up again. Whenever I want to save my ligand as pdbqt, Discovery Studio is showing a destroyed structure, not the real molecule structure. I've done what you've done. If I write as complex, I can see the ligand structure as it should be.
@BioinfoCopilotАй бұрын
I see the problem. Please open the ligand in Marvin Sketch or Chemsketch. Do a quick check for the structure. Optimize it. And then use Autodock to create pdbqt file. After docking if you open in DS Visualizer it will be shown correctly.
@erenkasimfirtina7870Ай бұрын
Normally, dock.dpf file containing AD4.1_bound.dat isn't running because the program can't find AD4.1_bound.dat file. When I encountered this problem, I looked at your dpf file and saw outlev 1 instead of AD4.1_bound.dat. Then I changed as outlev 1. Now it's running without any mistake, but according to what I've read, outlev just ignore the mistake. I'm looking forward to hearing from you. 🤔
@BioinfoCopilotАй бұрын
It should be named as exactly what I have shown. Sometimes it throws error but as you said outlev ignores the mistakes. Generally if the ligand is made up of generic atoms it will consider it for docking. If metal ions are present then it will throw an error.
@arijitsamanta3660Ай бұрын
Sir if I want to change the simulation run time length, do I change the steps only in the mp.mdp file or other mdp files too?
@BioinfoCopilotАй бұрын
You have to change in the mdp file and then you can start it from the point where the simulation is done. Remember that you have to use grompp to generate new files and then run the simulation.
@sinhasuman1944Ай бұрын
Is this for just single molecule docking or it can do virtual screening?
@BioinfoCopilotАй бұрын
I have shown an example for single docking but you can perform virtual screening as well.
@imanebensahbane6318Ай бұрын
Thanks, your channel is very useful !
@BioinfoCopilotАй бұрын
Thank you 😊
@MuhannadDehniАй бұрын
IF I want hydrogen bond occupancy percentage for one residue during the simulation what cod should I use, Thank you.
@BioinfoCopilotАй бұрын
gmx_hbond analysis
@은진-e8nАй бұрын
Hi, I have some questions. I finished docking and opened the results. but, my results have 23. The only top 10 results are not showing
@BioinfoCopilotАй бұрын
If you have selected binding poses to 23 then it will show 23 binding poses. The top binding poses will be at the end of the file. Check the histogram that you obtained in your .dlg file.
@은진-e8nАй бұрын
Where can I specify the number of binding poses?
@은진-e8nАй бұрын
Is it okay if I send the file via email?
@AshrafwazirkhanАй бұрын
Another masterpiece of work from you. We are learning a lot from you. Sir one request if possible can you upload video on 3D graph of PCA analysis. Further if possible can we correlate the Gibbs free energy with the RMSD. Thanks for wonderful work❤
@BioinfoCopilotАй бұрын
Thank you 😊. I will try to do so soon!
@effyc9843Ай бұрын
Thank you sooo much for the video. This is really helpful for me! I followed through the tutorial and obtained my analysis results, but I'm having a bit difficulty understanding the plots. Could you make another one on explaining how to read them?
@BioinfoCopilotАй бұрын
Thank you very much. Regarding the plots, only one value that you should consider is the total free energy value. Rest of the plots are for publication (mostly supplementary to support your analysis). Thanks
@penscriber9320Ай бұрын
What is liked in this video is "Correcting the Errors". Like you really understood the errors and then solved them. Now, I will be doing MD simulation for my protein-ligand complex. Subscribed
@BioinfoCopilotАй бұрын
Thank you very much 😊
@LucasPaul-d9qАй бұрын
Do have any of your videos with docking validation and statistical analysis? If possible I would recommend to do that
@BioinfoCopilotАй бұрын
Sure will do that. Thanks for the recommendation.
@varuntk1183Ай бұрын
Hi, Can we extract the most stable comformation from FEL?
@BioinfoCopilotАй бұрын
There is no guarantee of this, and MD simulations are not driven by free energy, which is what is truly necessary to say one has the most “stable” conformation. The total potential energy in the MD simulation is not something that force fields are generally parametrized to ascribe any meaningful value, so there is no way to know if you are in a global energy minimum at any point during the trajectory.
@varuntk1183Ай бұрын
@@BioinfoCopilot certain publication have done like extract the ensemble having lowest energy minima from fel plot, I was also not able to understand how it can be done
@BioinfoCopilotАй бұрын
Use trjconv to dump pdb files at several intervals. Then check the binding orientations. For each pdb file at certain interval check for free energy like I did for 0-100ns.