Scoring functions for fragment-based drug discovery.

Fragment-based drug design represents a challenge for computational drug design because almost inevitably fragments will be weak binders to the biomolecular targets of a specific disease, and the performances of the scoring functions for weak binders are usually poorer than those for the stronger binders. This protocol describes how to predict the binding modes and binding affinities of fragments towards their binding partner with our refined AutoDock scoring function incorporating a quantum chemical charge model, namely, the restrained electrostatic potential (RESP) model. This scoring function was calibrated by robust regression analysis and has been demonstrated to perform well for general classes of protein-ligand interactions and for weak binders (with root-mean square of error of about 2.1 kcal/mol).

A low-cost, multiplexed electrophysiology system for chronic µECoG recordings in rodents.

Micro-Electrocorticography (µECoG) offers a minimally invasive, high resolution interface with large areas of cortex. However, large arrays of electrodes with many contacts that are individually wired to external recording systems are cumbersome and make chronic recording in freely behaving small animals challenging. Multiplexed headstages overcome this limitation by combining the signals from many electrodes to a smaller number of connections directly on the animal's head. Commercially available multiplexed headstages provide high performance integrated amplification, multiplexing and analog to digital conversion. However, the cost of these systems can be prohibitive for small labs or for experiments that require a large number of animals to be continuously recorded at the same time. Here we have developed a multiplexed 60-channel headstage amplifier optimized to chronically record electrophysiological signals from high-density µECoG electrode arrays. A single, ultraflexible (2 mm thickness) microHDMI cable provided the data interface. Using low cost components, we have reduced the cost of the multiplexed headstage to ~$125. Paired with a custom interface printed circuit board (PCB) and a general purpose data acquisition system (M-series DAQ, National Instruments), an inexpensive and customizable electrophysiology system is assembled. Open source LabVIEW software that we have previously released controlled the system. It can also be used with other open source neural data acquisition packages. Combined, we have presented a scalable, low-cost platform for high-channel count electrophysiology.

Scoring functions for prediction of protein-ligand interactions.

The scoring functions for protein-ligand interactions plays central roles in computational drug design, virtual screening of chemical libraries for new lead identification, and prediction of possible binding targets of small chemical molecules. An ideal scoring function for protein-ligand interactions is expected to be able to recognize the native binding pose of a ligand on the protein surface among decoy poses, and to accurately predict the binding affinity (or binding free energy) so that the active molecules can be discriminated from the non-active ones. Due to the empirical nature of most, if not all, scoring functions for protein-ligand interactions, the general applicability of empirical scoring functions, especially to domains far outside training sets, is a major concern. In this review article, we will explore the foundations of different classes of scoring functions, their possible limitations, and their suitable application domains. We also provide assessments of several scoring functions on weakly-interacting protein-ligand complexes, which will be useful information in computational fragment-based drug design or virtual screening.

idTarget: a web server for identifying protein targets of small chemical molecules with robust scoring functions and a divide-and-conquer docking approach.

Identification of possible protein targets of small chemical molecules is an important step for unravelling their underlying causes of actions at the molecular level. To this end, we construct a web server, idTarget, which can predict possible binding targets of a small chemical molecule via a divide-and-conquer docking approach, in combination with our recently developed scoring functions based on robust regression analysis and quantum chemical charge models. Affinity profiles of the protein targets are used to provide the confidence levels of prediction. The divide-and-conquer docking approach uses adaptively constructed small overlapping grids to constrain the searching space, thereby achieving better docking efficiency. Unlike previous approaches that screen against a specific class of targets or a limited number of targets, idTarget screen against nearly all protein structures deposited in the Protein Data Bank (PDB). We show that idTarget is able to reproduce known off-targets of drugs or drug-like compounds, and the suggested new targets could be prioritized for further investigation. idTarget is freely available as a web-based server at http://idtarget.rcas.sinica.edu.tw.

Enhancement of initial equivalency for protein structure alignment based on encoded local structures.

Most alignment algorithms find an initial equivalent residue pair followed by an iterative optimization process to explore better near-optimal alignments in the surrounding solution space of the initial alignment. It plays a decisive role in determining the alignment quality since a poor initial alignment may make the final alignment trapped in an undesirable local optimum even with an iterative optimization. We proposed a vector-based alignment algorithm with a new initial alignment approach accounting for local structure features called MIRAGE-align. The new idea is to enhance the quality of the initial alignment based on encoded local structural alphabets to identify the protein structure pair whose sequence identity falls in or below twilight zone. The statistical analysis of alignment quality based on Match Index (MI) and computation time demonstrated that MIRAGE-align algorithm outperformed four previously published algorithms, i.e., the residue-based algorithm (CE), the vector-based algorithm (SSM), TM-align, and Fr-TM-align. MIRAGE-align yields a better estimate of initial solution to enhance the quality of initial alignment and enable the employment of a non-iterative optimization process to achieve a better alignment.

Robust scoring functions for protein-ligand interactions with quantum chemical charge models.

Ordinary least-squares (OLS) regression has been used widely for constructing the scoring functions for protein-ligand interactions. However, OLS is very sensitive to the existence of outliers, and models constructed using it are easily affected by the outliers or even the choice of the data set. On the other hand, determination of atomic charges is regarded as of central importance, because the electrostatic interaction is known to be a key contributing factor for biomolecular association. In the development of the AutoDock4 scoring function, only OLS was conducted, and the simple Gasteiger method was adopted. It is therefore of considerable interest to see whether more rigorous charge models could improve the statistical performance of the AutoDock4 scoring function. In this study, we have employed two well-established quantum chemical approaches, namely the restrained electrostatic potential (RESP) and the Austin-model 1-bond charge correction (AM1-BCC) methods, to obtain atomic partial charges, and we have compared how different charge models affect the performance of AutoDock4 scoring functions. In combination with robust regression analysis and outlier exclusion, our new protein-ligand free energy regression model with AM1-BCC charges for ligands and Amber99SB charges for proteins achieve lowest root-mean-squared error of 1.637 kcal/mol for the training set of 147 complexes and 2.176 kcal/mol for the external test set of 1427 complexes. The assessment for binding pose prediction with the 100 external decoy sets indicates very high success rate of 87% with the criteria of predicted root-mean-squared deviation of less than 2 Å. The success rates and statistical performance of our robust scoring functions are only weakly class-dependent (hydrophobic, hydrophilic, or mixed).