Adverse drug reaction prediction using scores produced by large-scale drug-protein target docking on high-performance computing machines.

Late-stage or post-market identification of adverse drug reactions (ADRs) is a significant public health issue and a source of major economic liability for drug development. Thus, reliable in silico screening of drug candidates for possible ADRs would be advantageous. In this work, we introduce a computational approach that predicts ADRs by combining the results of molecular docking and leverages known ADR information from DrugBank and SIDER. We employed a recently parallelized version of AutoDock Vina (VinaLC) to dock 906 small molecule drugs to a virtual panel of 409 DrugBank protein targets. L1-regularized logistic regression models were trained on the resulting docking scores of a 560 compound subset from the initial 906 compounds to predict 85 side effects, grouped into 10 ADR phenotype groups. Only 21% (87 out of 409) of the drug-protein binding features involve known targets of the drug subset, providing a significant probe of off-target effects. As a control, associations of this drug subset with the 555 annotated targets of these compounds, as reported in DrugBank, were used as features to train a separate group of models. The Vina off-target models and the DrugBank on-target models yielded comparable median area-under-the-receiver-operating-characteristic-curves (AUCs) during 10-fold cross-validation (0.60-0.69 and 0.61-0.74, respectively). Evidence was found in the PubMed literature to support several putative ADR-protein associations identified by our analysis. Among them, several associations between neoplasm-related ADRs and known tumor suppressor and tumor invasiveness marker proteins were found. A dual role for interstitial collagenase in both neoplasms and aneurysm formation was also identified. These associations all involve off-target proteins and could not have been found using available drug/on-target interaction data. This study illustrates a path forward to comprehensive ADR virtual screening that can potentially scale with increasing number of CPUs to tens of thousands of protein targets and millions of potential drug candidates.

Mutually exclusive T-cell receptor induction and differential susceptibility to human immunodeficiency virus type 1 mutational escape associated with a two-amino-acid difference between HLA class I subtypes.

The relative contributions of HLA alleles and T-cell receptors (TCRs) to the prevention of mutational viral escape are unclear. Here, we examined human immunodeficiency virus type 1 (HIV-1)-specific CD8(+) T-cell responses restricted by two closely related HLA class I alleles, B*5701 and B*5703, that differ by two amino acids but are both associated with a dominant response to the same HIV-1 Gag epitope KF11 (KAFSPEVIPMF). When this epitope is presented by HLA-B*5701, it induces a TCR repertoire that is highly conserved among individuals, cross-recognizes viral epitope variants, and is rarely associated with mutational escape. In contrast, KF11 presented by HLA-B*5703 induces an entirely different, more heterogeneous TCR beta-chain repertoire that fails to recognize specific KF11 escape variants which frequently arise in clade C-infected HLA-B*5703(+) individuals. These data show the influence of HLA allele subtypes on TCR selection and indicate that extensive TCR diversity is not a prerequisite to prevention of allowable viral mutations.

Conformational dependence of a protein kinase phosphate transfer reaction.

Atomic motions and energetics for a phosphate transfer reaction catalyzed by the cAMP-dependent protein kinase are calculated by plane-wave density functional theory, starting from structures of proteins crystallized in both the reactant conformation (RC) and the transition-state conformation (TC). In TC, we calculate that the reactants and products are nearly isoenergetic with a 20-kJ/mol barrier, whereas phosphate transfer is unfavorable by 120 kJ/mol in the RC, with an even higher barrier. With the protein in TC, the motions involved in reaction are small, with only P(gamma) and the catalytic proton moving >0.5 A. Examination of the structures reveals that in the RC the active site cleft is not completely closed and there is insufficient space for the phosphorylated serine residue in the product state. Together, these observations imply that the phosphate transfer reaction occurs rapidly and reversibly in a particular conformation of the protein, and that the reaction can be gated by changes of a few tenths of an angstrom in the catalytic site.