Barry J. Grant

Bio3d: an R package for the comparative analysis of protein structures

UNLABELLED An automated procedure for the analysis of homologous protein structures has been developed. The method facilitates the characterization of internal conformational differences and inter-conformer relationships and provides a framework for the analysis of protein structural evolution. The method is implemented in bio3d, an R package for the exploratory analysis of structure and sequence data. AVAILABILITY 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/

Drug hypersensitivity caused by alteration of the MHC-presented self-peptide repertoire

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 HLA-linked 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.

Ras conformational switching: simulating nucleotide-dependent conformational transitions with accelerated molecular dynamics.

Ras mediates signaling pathways controlling cell proliferation and development by cycling between GTP- and GDP-bound active and inactive conformational states. Understanding the complete reaction path of this conformational change and its intermediary structures is critical to understanding Ras signaling. We characterize nucleotide-dependent conformational transition using multiple-barrier-crossing accelerated molecular dynamics (aMD) simulations. These transitions, achieved for the first time for wild-type Ras, are impossible to observe with classical molecular dynamics (cMD) simulations due to the large energetic barrier between end states. Mapping the reaction path onto a conformer plot describing the distribution of the crystallographic structures enabled identification of highly populated intermediate structures. These structures have unique switch orientations (residues 25–40 and 57–75) intermediate between GTP and GDP states, or distinct loop3 (46–49), loop7 (105–110), and α5 C-terminus (159–166) conformations distal from the nucleotide-binding site. In addition, these barrier-crossing trajectories predict novel nucleotide-dependent correlated motions, including correlations of α2 (residues 66–74) with α3-loop7 (93–110), loop2 (26–37) with loop10 (145–151), and loop3 (46–49) with α5 (152–167). The interconversion between newly identified Ras conformations revealed by this study advances our mechanistic understanding of Ras function. In addition, the pattern of correlated motions provides new evidence for a dynamic linkage between the nucleotide-binding site and the membrane interacting C-terminus critical for the signaling function of Ras. Furthermore, normal mode analysis indicates that the dominant collective motion that occurs during nucleotide-dependent conformational exchange, and captured in aMD (but absent in cMD) simulations, is a low-frequency motion intrinsic to the structure.

Mapping the Nucleotide and Isoform-Dependent Structural and Dynamical Features of Ras Proteins

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.

Targeted molecular dynamics study of C-loop closure and channel gating in nicotinic receptors

The initial coupling between ligand binding and channel gating in the human α7 nicotinic acetylcholine receptor (nAChR) has been investigated with targeted molecular dynamics (TMD) simulation. During the simulation, eight residues at the tip of the C-loop in two alternating subunits were forced to move toward a ligand-bound conformation as captured in the crystallographic structure of acetylcholine binding protein (AChBP) in complex with carbamoylcholine. Comparison of apo- and ligand-bound AChBP structures shows only minor rearrangements distal from the ligand-binding site. In contrast, comparison of apo and TMD simulation structures of the nAChR reveals significant changes toward the bottom of the ligand-binding domain. These structural rearrangements are subsequently translated to the pore domain, leading to a partly open channel within 4 ns of TMD simulation. Furthermore, we confirmed that two highly conserved residue pairs, one located near the ligand-binding pocket (Lys145 and Tyr188), and the other located toward the bottom of the ligand-binding domain (Arg206 and Glu45), are likely to play important roles in coupling agonist binding to channel gating. Overall, our simulations suggest that gating movements of the α7 receptor may involve relatively small structural changes within the ligand-binding domain, implying that the gating transition is energy-efficient and can be easily modulated by agonist binding/unbinding.

Novel allosteric sites on ras for lead generation

Aberrant Ras activity is a hallmark of diverse cancers and developmental diseases. Unfortunately, conventional efforts to develop effective small molecule Ras inhibitors have met with limited success. We have developed a novel multi-level computational approach to discover potential inhibitors of previously uncharacterized allosteric sites. Our approach couples bioinformatics analysis, advanced molecular simulations, ensemble docking and initial experimental testing of potential inhibitors. Molecular dynamics simulation highlighted conserved allosteric coupling of the nucleotide-binding switch region with distal regions, including loop 7 and helix 5. Bioinformatics methods identified novel transient small molecule binding pockets close to these regions and in the vicinity of the conformationally responsive switch region. Candidate binders for these pockets were selected through ensemble docking of ZINC and NCI compound libraries. Finally, cell-based assays confirmed our hypothesis that the chosen binders can inhibit the downstream signaling activity of Ras. We thus propose that the predicted allosteric sites are viable targets for the development and optimization of new drugs.

Accessing a hidden conformation of the maltose binding protein using accelerated molecular dynamics.

Periplasmic binding proteins (PBPs) are a large family of molecular transporters that play a key role in nutrient uptake and chemotaxis in Gram-negative bacteria. All PBPs have characteristic two-domain architecture with a central interdomain ligand-binding cleft. Upon binding to their respective ligands, PBPs undergo a large conformational change that effectively closes the binding cleft. This conformational change is traditionally viewed as a ligand induced-fit process; however, the intrinsic dynamics of the protein may also be crucial for ligand recognition. Recent NMR paramagnetic relaxation enhancement (PRE) experiments have shown that the maltose binding protein (MBP) - a prototypical member of the PBP superfamily - exists in a rapidly exchanging (ns to µs regime) mixture comprising an open state (approx 95%), and a minor partially closed state (approx 5%). Here we describe accelerated MD simulations that provide a detailed picture of the transition between the open and partially closed states, and confirm the existence of a dynamical equilibrium between these two states in apo MBP. We find that a flexible part of the protein called the balancing interface motif (residues 175–184) is displaced during the transformation. Continuum electrostatic calculations indicate that the repacking of non-polar residues near the hinge region plays an important role in driving the conformational change. Oscillations between open and partially closed states create variations in the shape and size of the binding site. The study provides a detailed description of the conformational space available to ligand-free MBP, and has implications for understanding ligand recognition and allostery in related proteins.

The Distinct Conformational Dynamics of K-Ras and H-Ras A59G

Ras proteins regulate signaling cascades crucial for cell proliferation and differentiation by switching between GTP- and GDP-bound conformations. Distinct Ras isoforms have unique physiological functions with individual isoforms associated with different cancers and developmental diseases. Given the small structural differences among isoforms and mutants, it is currently unclear how these functional differences and aberrant properties arise. Here we investigate whether the subtle differences among isoforms and mutants are associated with detectable dynamical differences. Extensive molecular dynamics simulations reveal that wild-type K-Ras and mutant H-Ras A59G are intrinsically more dynamic than wild-type H-Ras. The crucial switch 1 and switch 2 regions along with loop 3, helix 3, and loop 7 contribute to this enhanced flexibility. Removing the gamma-phosphate of the bound GTP from the structure of A59G led to a spontaneous GTP-to-GDP conformational transition in a 20-ns unbiased simulation. The switch 1 and 2 regions exhibit enhanced flexibility and correlated motion when compared to non-transitioning wild-type H-Ras over a similar timeframe. Correlated motions between loop 3 and helix 5 of wild-type H-Ras are absent in the mutant A59G reflecting the enhanced dynamics of the loop 3 region. Taken together with earlier findings, these results suggest the existence of a lower energetic barrier between GTP and GDP states of the mutant. Molecular dynamics simulations combined with principal component analysis of available Ras crystallographic structures can be used to discriminate ligand- and sequence-based dynamic perturbations with potential functional implications. Furthermore, the identification of specific conformations associated with distinct Ras isoforms and mutants provides useful information for efforts that attempt to selectively interfere with the aberrant functions of these species.

Induced Fit or Conformational Selection? The Role of the Semi-closed State in the Maltose Binding Protein

A full characterization of the thermodynamic forces underlying ligand-associated conformational changes in proteins is essential for understanding and manipulating diverse biological processes, including transport, signaling, and enzymatic activity. Recent experiments on the maltose binding protein (MBP) have provided valuable data about the different conformational states implicated in the ligand recognition process; however, a complete picture of the accessible pathways and the associated changes in free energy remains elusive. Here we describe results from advanced accelerated molecular dynamics (aMD) simulations, coupled with adaptively biased force (ABF) and thermodynamic integration (TI) free energy methods. The combination of approaches allows us to track the ligand recognition process on the microsecond time scale and provides a detailed characterization of the protein’s dynamic and the relative energy of stable states. We find that an induced-fit (IF) mechanism is most likely and that a mechanism involving both a conformational selection (CS) step and an IF step is also possible. The complete recognition process is best viewed as a “Pac Man” type action where the ligand is initially localized to one domain and naturally occurring hinge-bending vibrations in the protein are able to assist the recognition process by increasing the chances of a favorable encounter with side chains on the other domain, leading to a population shift. This interpretation is consistent with experiments and provides new insight into the complex recognition mechanism. The methods employed here are able to describe IF and CS effects and provide formally rigorous means of computing free energy changes. As such, they are superior to conventional MD and flexible docking alone and hold great promise for future development and applications to drug discovery.

Conformational Selection in G-Proteins: Lessons from Ras and Rho

The induced fit model has traditionally been invoked to describe the activating conformational change of the monomeric G-proteins, such as Ras and Rho. With this scheme, the presence or absence of the γ-phosphate of GTP leads to an instantaneous switch in conformation. Here we describe atomistic molecular simulations that demonstrate that both Ras and Rho superfamily members harbor an intrinsic susceptibility to sample multiple conformational states in the absence of nucleotide ligand. By comparing the distribution of conformers in the presence and absence of nucleotide, we show that conformational selection is the dominant mechanism by which Ras and Rho undergo nucleotide-dependent conformational changes. Furthermore, the pattern of correlated motions revealed by these simulations predicts a preserved allosteric coupling of the nucleotide-binding site with the membrane interacting C-terminus in both Rho and Ras.