Jocelyn Sunseri edited section_Query_Definition_subsection_Inputs__.tex  about 8 years ago

Commit id: 0c825a7774088b87ae16f5825f0fe0a39be8293f

deletions | additions      

       

\subsection{Inputs}  The typical starting point for a Pharmit session is a ligand-receptor complex structure. A Pharmit session can be automatically initialized using any complex in the PDB by inputing the corresponding PDB ascension code on the Pharmit main page  and selecting how active site water molecules should be treated (ignored, as part of the receptor, or as part of the ligand). The dropdown menu next to the box where a PDB code may be entered allows the user to select which ligand found in the PDB file should be used to generate a pharmacophore.  Alternatively, a user can upload their own complex, in which case the receptor and ligand structures must be in separate files. files; these are uploaded by clicking "Enter Pharmit Search" on the main page and then choosing "Load receptor" and "Load features" on the resulting page.  Any file format supported by OpenBabel \cite{O_Boyle_2011} may be used. Pharmit prepares the receptor by protonating it with OpenBabel, but no other modifications are made. Thus the user must decide whether there are missing residues that should be included, the histidine protonation state is correct, or any other structural changes to the receptor are necessary. Note that a partial charge calculation is irrelevant as it is not used by the AutoDock Vina scoring function. The user may choose to provide the ligand structure as the file utilized when loading features, in which case Pharmit will automatically generate a set of relevant features from which a pharmacophore may be produced, or external software (e.g. MOE or PharmaGIST) may be used to derive a ligand-based pharmacophore and the result may be uploaded to Pharmit.  Note that the query ligand must be pre-positioned within the binding site of the receptor - Pharmit does not perform docking. It uses the pharmacophore and shape features of a known ligand to screen for novel compounds. \subsection{Pharmacophore Queries}  A pharmacophore \cite{Koes_2015rev,Yang_2010,Leach_2010} defines the essential features of an interaction. Importantly, a pharmacophore includes the spatial arrangement of these features.  

\subsection{Shape Queries}  Similarity of molecular shape is a common method of structure-based virtual screening \cite{Nicholls_2010}. Pharmit uses the Volumetric Aligned Molecular Shapes (VAMS) \cite{vams} method of shape search, which uses inclusive and exclusive shape constraints to identify matching molecular shapes. In Pharmit, inclusive constraints are specified using the shape of the provided ligand or by manually specified inclusion spheres. Inclusive constraints specify a minimum bound on the desired molecular shape; matching compounds will overlap these constraints. Exclusive constraints are specified using the shape of the provided receptor or by manually specified exclusion spheres. Exclusive constraints are used to limit the desired molecular shape; matching compounds are prohibited from overlapping these constraints. Both constraints are represented using voxelized volumes, as shown in Figure~\ref{shapefig}, and can be adjusted by growing or shrinking the constraint volume. A shape-first search with pharmacophore restraints is similar to color ROCS, but it does not optimize the position with respect to the pharmacophores.