Knowledge-Based Methods To Train and Optimize Virtual Screening Ensembles.

Ensemble docking can be a successful virtual screening technique that addresses the innate conformational heterogeneity of macromolecular drug targets. Yet, lacking a method to identify a subset of conformational states that effectively segregates active and inactive small molecules, ensemble docking may result in the recommendation of a large number of false positives. Here, three knowledge-based methods that construct structural ensembles for virtual screening are presented. Each method selects ensembles by optimizing an objective function calculated using the receiver operating characteristic (ROC) curve: either the area under the ROC curve (AUC) or a ROC enrichment factor (EF). As the number of receptor conformations, N, becomes large, the methods differ in their asymptotic scaling. Given a set of small molecules with known activities and a collection of target conformations, the most resource intense method is guaranteed to find the optimal ensemble but scales as O(2(N)). A recursive approximation to the optimal solution scales as O(N(2)), and a more severe approximation leads to a faster method that scales linearly, O(N). The techniques are generally applicable to any system, and we demonstrate their effectiveness on the androgen nuclear hormone receptor (AR), cyclin-dependent kinase 2 (CDK2), and the peroxisome proliferator-activated receptor δ (PPAR-δ) drug targets. Conformations that consisted of a crystal structure and molecular dynamics simulation cluster centroids were used to form AR and CDK2 ensembles. Multiple available crystal structures were used to form PPAR-δ ensembles. For each target, we show that the three methods perform similarly to one another on both the training and test sets.

Conserved gating elements in TRPC4 and TRPC5 channels.

TRPC4 and TRPC5 proteins share 65% amino acid sequence identity and form Ca(2+)-permeable nonselective cation channels. They are activated by stimulation of receptors coupled to the phosphoinositide signaling cascade. Replacing a conserved glycine residue within the cytosolic S4-S5 linker of both proteins by a serine residue forces the channels into an open conformation. Expression of the TRPC4G503S and TRPC5G504S mutants causes cell death, which could be prevented by buffering the Ca(2+) of the culture medium. Current-voltage relationships of the TRPC4G503S and TRPC5G504S mutant ion channels resemble that of fully activated TRPC4 and TRPC5 wild-type channels, respectively. Modeling the structure of the transmembrane domains and the pore region (S4-S6) of TRPC4 predicts a conserved serine residue within the C-terminal sequence of the predicted S6 helix as a potential interaction site. Introduction of a second mutation (S623A) into TRPC4G503S suppressed the constitutive activation and partially rescued its function. These results indicate that the S4-S5 linker is a critical constituent of TRPC4/C5 channel gating and that disturbance of its sequence allows channel opening independent of any sensor domain.