Summary: We will use structure-based ultra-large virtual screenings using VirtualFlow. Step 1: Protein preparation Protein structures will be prepared with Maestro from Schrödinger (protonation state assignment, assignment of missing atoms/side chains, hydrogen atoms, ...). MD simulations of the target protein will be carried out using Amber 18. Conformations will be clustered, and representative structures of the clusters will be used for the virtual screens. Step 2: Hit identification The hit identification step will consist of two virtual screening stages. - Stage 1: We will use structure-based ultra-large virtual screenings using physics-based docking methods AutoDock Vina, QuickVina, Smina, and PLANTS. We will screen multi-billion ligand libraries with VirtuaFlow, and open-source platform for ultra-large virtual screens. The libraries we are using are the ZINC and the Enamine REAL libraries. Due to the large-scale computations required for this approach, we will use the supercomputer of our university, as well as the cloud if additional computation time is required. We have extensive experience using the cloud and how to use over 100,000 CPUs in parallel. The protein will be held rigid in stage 1. The ligand libraries which we will be using are the libraries part of the VirtualFlow project. These libraries contain ligands which have been protonated, tautomerized, the 3D conformation has been computed, and the ligands are the ready-to-dock PDBQT format. - Stage 2: We will rescreen the top 1 million compounds of stage 1 in stage 2, and will allow the protein side chains at the binding site to be flexible. Multiple protein backbone conformations might be used to carry out ensemble dockings in addition, based on the results of the MD simulations in the protein preparation step (see section above). Step 3: Hit optimization In the second round of the challenge (hit optimization), we will search the chemical space (our available libraries) for the most similar analogs, and screen them with again VirtualFlow as described above in step 2. Step 4: Postprocessing of the results The screened compounds will be ranked by their docking score. Of the top 1000 compounds, biophysical and pharmacokinetic properties will be computed. Compounds with unfavorable properties (e.g. too high logP) will be filtered out.
