Continuum Electromechanical Modeling of Protein-Membrane Interaction

A continuum electromechanical model is proposed to describe the membrane curvature induced by electrostatic interactions in a solvated protein-membrane system. The model couples the macroscopic strain energy of membrane and the electrostatic solvation energy of the system, and equilibrium membrane deformation is obtained by minimizing the electro-elastic energy functional with respect to the dielectric interface. The model is illustrated with the systems with increasing geometry complexity and captures the sensitivity of membrane curvature to the permanent and mobile charge distributions.Comment: 5 pages, 12 figure

Water-Membrane Partition Thermodynamics of an Amphiphilic Lipopeptide: An Enthalpy-Driven Hydrophobic Effect

To shed light on the driving force for the hydrophobic effect that partitions amphiphilic lipoproteins between water and membrane, we carried out an atomically detailed thermodynamic analysis of a triply lipid modified H-ras heptapeptide anchor (ANCH) in water and in a DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) bilayer. Combining molecular mechanical and continuum solvent approaches with an improved technique for solute entropy calculation, we obtained an overall transfer free energy of ∼−13 kcal mol−1. This value is in qualitative agreement with free energy changes derived from a potential of mean force calculation and indirect experimental observations. Changes in free energies of solvation and ANCH conformational reorganization are unfavorable, whereas ANCH-DMPC interactions—especially van der Waals—favor insertion. These results are consistent with an enthalpy-driven hydrophobic effect, in accord with earlier calorimetric data on the membrane partition of other amphiphiles. Furthermore, structural and entropic analysis of molecular dynamics-generated ensembles suggests that conformational selection may play a hitherto unappreciated role in membrane insertion of lipid-modified peptides and proteins

Non-bisphosphonate inhibitors of isoprenoid biosynthesis identified via computer-aided drug design.

The relaxed complex scheme, a virtual-screening methodology that accounts for protein receptor flexibility, was used to identify a low-micromolar, non-bisphosphonate inhibitor of farnesyl diphosphate synthase. Serendipitously, we also found that several predicted farnesyl diphosphate synthase inhibitors were low-micromolar inhibitors of undecaprenyl diphosphate synthase. These results are of interest because farnesyl diphosphate synthase inhibitors are being pursued as both anti-infective and anticancer agents, and undecaprenyl diphosphate synthase inhibitors are antibacterial drug leads

KRAS Ubiquitination at Lysine 104 Retains Exchange Factor Regulation by Dynamically Modulating the Conformation of the Interface

RAS proteins function as highly regulated molecular switches that control cellular growth. In addition to regulatory proteins, RAS undergoes a number of posttranslational modifications (PTMs) that regulate its activity. Lysine 104, a hot spot for multiple PTMs, is a highly conserved residue that forms key interactions that stabilize the RAS helix-2(H2)/helix-3(H3) interface. Mutation at 104 attenuates interaction with guanine nucleotide exchange factors (GEFs), whereas ubiquitination at lysine 104 retains GEF regulation. To elucidate how ubiquitination modulates RAS function, we generated monoubiquitinated KRAS at 104 using chemical biology approaches and conducted biochemical, NMR, and computational analyses. We find that ubiquitination promotes a new dynamic interaction network and alters RAS conformational dynamics to retain GEF function. These findings reveal a mechanism by which ubiquitination can regulate protein function

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

Gated Diffusion-controlled Reactions

The binding and active sites of proteins are often dynamically occluded by motion of the nearby polypeptide. A variety of theoretical and computational methods have been developed to predict rates of ligand binding and reactivity in such cases. Two general approaches exist, "protein centric" approaches that explicitly treat only the protein target, and more detailed dynamical simulation approaches in which target and ligand are both treated explicitly. This mini-review describes recent work in this area and some of the biological implications

Dynamics of the Acetylcholinesterase Tetramer

Acetylcholinesterase rapidly hydrolyzes the neurotransmitter acetylcholine in cholinergic synapses, including the neuromuscular junction. The tetramer is the most important functional form of the enzyme. Two low-resolution crystal structures have been solved. One is compact with two of its four peripheral anionic sites (PAS) sterically blocked by complementary subunits. The other is a loose tetramer with all four subunits accessible to solvent. These structures lacked the C-terminal amphipathic t-peptide (WAT domain) that interacts with the proline-rich attachment domain (PRAD). A complete tetramer model (AChEt) was built based on the structure of the PRAD/WAT complex and the compact tetramer. Normal mode analysis suggested that AChEt could exist in several conformations with subunits fluctuating relative to one another. Here, a multiscale simulation involving all-atom molecular dynamics and Cα-based coarse-grained Brownian dynamics simulations was carried out to investigate the large-scale intersubunit dynamics in AChEt. We sampled the ns-μs timescale motions and found that the tetramer indeed constitutes a dynamic assembly of monomers. The intersubunit fluctuation is correlated with the occlusion of the PAS. Such motions of the subunits “gate” ligand-protein association. The gates are open more than 80% of the time on average, which suggests a small reduction in ligand-protein binding. Despite the limitations in the starting model and approximations inherent in coarse graining, these results are consistent with experiments which suggest that binding of a substrate to the PAS is only somewhat hindered by the association of the subunits

Palmitoylated ras proteins traffic through recycling endosomes to the plasma membrane during exocytosis

Ras proteins regulate cell growth, death, and differentiation, and it is well established that this functional versatility is accomplished through their different subcellular localizations. Palmitoylated H- and N-Ras are believed to localize at the perinuclear Golgi and plasma membrane (PM). Notably, however, recycling endosomes (REs) also localize to a perinuclear region, which is often indistinguishable from the Golgi. In this study, we show that active palmitoylated Ras proteins mainly localize intracellularly at REs and that REs act as a way station along the post-Golgi exocytic pathway to the PM. H-Ras requires two palmitoyl groups for RE targeting. The lack of either or both palmitoyl groups leads to the mislocalization of the mutant proteins to the endoplasmic reticulum, Golgi apparatus, or the PM. Therefore, we demonstrate that palmitoylation directs Ras proteins to the correct intracellular organelles for trafficking and activity. © 2010 Misaki et al

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

The role of tyrosine 71 in modulating the flap conformations of BACE1

β‐site amyloid precursor protein cleaving enzyme 1 (BACE1) is a potential target for treating Alzheimer's disease. BACE1's binding site is partially covered by a flexible loop on its N‐terminal domain, known as the “flap,” which has been found in several conformations in crystal structures of BACE1 and other aspartyl proteases. The side chain of the invariant residue Tyr71 on the flap adopts several rotameric orientations, leading to our hypothesis that the orientation of this residue dictates the movement and conformations available to the flap. We investigated this hypothesis by performing 220 ns of molecular dynamics simulations of bound and unbound wild‐type BACE1 as well as the unbound Y71A mutant. Our findings indicate that the flap exhibits various degrees of mobility and adopts different conformations depending on the Tyr71 orientation. Surprisingly, the “self‐inhibited” form is stable in our simulations, making it a reasonable target for drug design. The alanine mutant, lacking a large side chain at position 71, displays significant differences in flap dynamics from wild type, freely sampling very open and closed conformations. Our simulations show that Tyr71, in addition to its previously determined functions in catalysis and substrate binding, has the important role of modulating flap conformations in BACE1. Proteins 2011; © 2011 Wiley‐Liss, Inc.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/87074/1/PROT_23050_sm_suppinfoCoverArt.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/87074/2/PROT_23050_sm_suppinfo.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/87074/3/23050_ftp.pd