Predictive Power of Molecular Dynamics Receptor Structures in Virtual Screening

Molecular dynamics (MD) simulation is a well-established method for understanding protein dynamics. Conformations from unrestrained MD simulations have yet to be assessed for blind virtual screening (VS) by docking. This study presents a critical analysis of the predictive power of MD snapshots to this regard, evaluating two well-characterized systems of varying flexibility in ligand-bound and unbound configurations. Results from such VS predictions are discussed with respect to experimentally determined structures. In all cases, MD simulations provide snapshots that improve VS predictive power over known crystal structures, possibly due to sampling more relevant receptor conformations. Additionally, MD can move conformations previously not amenable to docking into the predictive range

Asymmetric Switching in a Homodimeric ABC Transporter: A Simulation Study

ABC transporters are a large family of membrane proteins involved in a variety of cellular processes, including multidrug and tumor resistance and ion channel regulation. Advances in the structural and functional understanding of ABC transporters have revealed that hydrolysis at the two canonical nucleotide-binding sites (NBSs) is co-operative and non-simultaneous. A conserved core architecture of bacterial and eukaryotic ABC exporters has been established, as exemplified by the crystal structure of the homodimeric multidrug exporter Sav1866. Currently, it is unclear how sequential ATP hydrolysis arises in a symmetric homodimeric transporter, since it implies at least transient asymmetry at the NBSs. We show by molecular dynamics simulation that the initially symmetric structure of Sav1866 readily undergoes asymmetric transitions at its NBSs in a pre-hydrolytic nucleotide configuration. MgATP-binding residues and a network of charged residues at the dimer interface are shown to form a sequence of putative molecular switches that allow ATP hydrolysis only at one NBS. We extend our findings to eukaryotic ABC exporters which often consist of two non-identical half-transporters, frequently with degeneracy substitutions at one of their two NBSs. Interestingly, many residues involved in asymmetric conformational switching in Sav1866 are substituted in degenerate eukaryotic NBS. This finding strengthens recent suggestions that the interplay of a consensus and a degenerate NBS in eukaroytic ABC proteins pre-determines the sequence of hydrolysis at the two NBSs

Activation of Akt by the Bacterial Inositol Phosphatase, SopB, is Wortmannin Insensitive

Salmonella enterica uses effector proteins translocated by a Type III Secretion System to invade epithelial cells. One of the invasion-associated effectors, SopB, is an inositol phosphatase that mediates sustained activation of the pro-survival kinase Akt in infected cells. Canonical activation of Akt involves membrane translocation and phosphorylation and is dependent on phosphatidyl inositide 3 kinase (PI3K). Here we have investigated these two distinct processes in Salmonella infected HeLa cells. Firstly, we found that SopB-dependent membrane translocation and phosphorylation of Akt are insensitive to the PI3K inhibitor wortmannin. Similarly, depletion of the PI3K regulatory subunits p85α and p85ß by RNAi had no inhibitory effect on SopB-dependent Akt phosphorylation. Nevertheless, SopB-dependent phosphorylation does depend on the Akt kinases, PDK1 and rictor-mTOR. Membrane translocation assays revealed a dependence on SopB for Akt recruitment to Salmonella ruffles and suggest that this is mediated by phosphoinositide (3,4) P2 rather than phosphoinositide (3,4,5) P3. Altogether these data demonstrate that Salmonella activates Akt via a wortmannin insensitive mechanism that is likely a class I PI3K-independent process that incorporates some essential elements of the canonical pathway

The intrinsic dynamics of enzymes plays a dominant role in determining the structural changes induced upon inhibitor binding

The conformational flexibility of target proteins continues to be a major challenge in accurate modeling of protein–inhibitor interactions. A fundamental issue, yet to be clarified, is whether the observed conformational changes are controlled by the protein or induced by the inhibitor. Although the concept of induced fit has been widely adopted for describing the structural changes that accompany ligand binding, there is growing evidence in support of the dominance of proteins' intrinsic dynamics which has been evolutionarily optimized to accommodate its functional interactions. The wealth of structural data for target proteins in the presence of different ligands now permits us to make a critical assessment of the balance between these two effects in selecting the bound forms. We focused on three widely studied drug targets, HIV-1 reverse transcriptase, p38 MAP kinase, and cyclin-dependent kinase 2. A total of 292 structures determined for these enzymes in the presence of different inhibitors and unbound form permitted us to perform an extensive comparative analysis of the conformational space accessed upon ligand binding, and its relation to the intrinsic dynamics before ligand binding as predicted by elastic network model analysis. Our results show that the ligand selects the conformer that best matches its structural and dynamic properties among the conformers intrinsically accessible to the protein in the unliganded form. The results suggest that simple but robust rules encoded in the protein structure play a dominant role in predefining the mechanisms of ligand binding, which may be advantageously exploited in designing inhibitors

The Mechanical Properties of PCNA: Implications for the Loading and Function of a DNA Sliding Clamp

Sliding clamps are toroidal proteins that encircle DNA and act as mobile platforms for DNA replication and repair machinery. To be loaded onto DNA, the eukaryotic sliding clamp Proliferating Cell Nuclear Antigen (PCNA) must be splayed open at one of the subunit-subunit interfaces by the ATP-dependent clamp loader, Replication Factor C, whose clamp-interacting sites form a right-handed spiral. Earlier molecular dynamics (MD) studies suggested that when PCNA opens, it preferentially adopts a right-handed spiral to match the spiral of the clamp loader. Here, analysis of considerably longer MD simulations shows that although the opened form of PCNA can achieve conformations matching the helical pitch of Replication Factor C, it is not biased toward a right-handed spiral structure. A coarse-grained elastic model was also built; its strong correspondence to the all-atom MD simulations of PCNA suggests that the behavior of the open clamp is primarily due to elastic deformation governed by the topology of the clamp domains. The elastic model was further used to construct the energy landscape of the opened PCNA clamp, including conformations that would allow binding to the clamp loader and loading onto double-stranded DNA. A picture of PCNA emerges of a rather flexible protein that, once opened, is mechanically compliant in the clamp opening process

PI3K-C2γ is a Rab5 effector selectively controlling endosomal Akt2 activation downstream of insulin signalling

In the liver, insulin-mediated activation of the phosphatidylinositol 3-kinase (PI3K)/Akt pathway is at the core of metabolic control. Multiple PI3K and Akt isoenzymes are found in hepatocytes and whether isoform-selective interplays exist is currently unclear. Here we report that insulin signalling triggers the association of the liver-specific class II PI3K isoform γ (PI3K-C2γ) with Rab5-GTP, and its recruitment to Rab5-positive early endosomes. In these vesicles, PI3K-C2γ produces a phosphatidylinositol-3,4-bisphosphate pool specifically required for delayed and sustained endosomal Akt2 stimulation. Accordingly, loss of PI3K-C2γ does not affect insulin-dependent Akt1 activation as well as S6K and FoxO1-3 phosphorylation, but selectively reduces Akt2 activation, which specifically inhibits glycogen synthase activity. As a consequence, PI3K-C2γ-deficient mice display severely reduced liver accumulation of glycogen and develop hyperlipidemia, adiposity as well as insulin resistance with age or after consumption of a high-fat diet. Our data indicate PI3K-C2γ supports an isoenzyme-specific forking of insulin-mediated signal transduction to an endosomal pool of Akt2, required for glucose homeostasis