The bio3d package is distributed with full source code as a platform-independent R package under a GPL2 license from: http://mccammon.ucsd.edu/~bgrant/bio3d/

Idiosyncratic adverse drug reactions are unpredictable, dose-independent and potentially life threatening; this makes them a major factor contributing to the cost and uncertainty of drug development. Clinical data suggest that many such reactions involve immune mechanisms, and genetic association studies have identified strong linkages between drug hypersensitivity reactions to several drugs and specific HLA alleles. One of the strongest such genetic associations found has been for the antiviral drug abacavir, which causes severe adverse reactions exclusively in patients expressing the HLA molecular variant B*57:01. Abacavir adverse reactions were recently shown to be driven by drug-specific activation of cytokine-producing, cytotoxic CD8 + T cells that required HLA-B*57:01 molecules for their function; however, the mechanism by which abacavir induces this pathologic T-cell response remains unclear. Here we show that abacavir can bind within the F pocket of the peptide-binding groove of HLA-B*57:01, thereby altering its specificity. This provides an explanation for HLA-linked idiosyncratic adverse drug reactions, namely that drugs can alter the repertoire of self-peptides presented to T cells, thus causing the equivalent of an alloreactive T-cell response. Indeed, we identified specific self-peptides that are presented only in the presence of abacavir and that were recognized by T cells of hypersensitive patients. The assays that we have established can be applied to test additional compounds with suspected HLAlinked hypersensitivities in vitro. Where successful, these assays could speed up the discovery and mechanistic understanding of HLA-linked hypersensitivities, and guide the development of safer drugs.3D structure | small molecule | binding site A bacavir is a nucleoside analog that suppresses HIV replication. In approximately 8% of recipients, abacavir is associated with significant immune-mediated drug hypersensitivity, which is strongly associated with the presence of the HLA-B*57:01 allele (1, 2). Three complementary models for the mechanism of immune-mediated severe adverse drug reactions have traditionally been discussed (3, 4). The hapten (or prohapten) model states that drugs and their metabolites are too small to be immunogenic on their own, but rather act like haptens and modify certain self-proteins in the host that lead to immune recognition of the resulting hapten-self-peptide complexes as de novo antigens (5-7). The pharmacologic interaction with immune receptors (p-i) model states that drugs can induce the formation of HLA-drug complexes that can activate T-cell immune responses directly without requiring a specific peptide ligand (8). The danger model, which is in principle compatible with other models, states that danger signals other than the drug itself (e.g., chemical, physical, or viral stress) are required to overcome immune tolerance barriers that otherwise suppress drug hypersensitivity reactions (7).None of these existing models provides a convincing mechanism explaining how abacav...

BackgroundPopular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution.ResultsHere, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case.ConclusionsThe integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-014-0399-6) contains supplementary material, which is available to authorized users.

Ras GTPases are conformational switches controlling cell proliferation, differentiation, and development. Despite their prominent role in many forms of cancer, the mechanism of conformational transition between inactive GDP-bound and active GTP-bound states remains unclear. Here we describe a detailed analysis of available experimental structures and molecular dynamics simulations to quantitatively assess the structural and dynamical features of active and inactive states and their interconversion. We demonstrate that GTP-bound and nucleotide-free G12V H-ras sample a wide region of conformational space, and show that the inactive-to-active transition is a multiphase process defined by the relative rearrangement of the two switches and the orientation of Tyr32. We also modeled and simulated N- and K-ras proteins and found that K-ras is more flexible than N- and H-ras. We identified a number of isoform-specific, long-range side chain interactions that define unique pathways of communication between the nucleotide binding site and the C terminus.