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

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

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

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

In silico and saturation transfer difference NMR approaches to unravel the binding mode of an andrographolide derivative to K-Ras oncoprotein

Background: Andrographolide and its benzylidene derivatives, SRJ09 and SRJ23, potentially bind oncogenic K-Ras to exert anticancer activity. Their molecular interactions with K-Ras oncoproteins that lead to effective biological activity are of major interest. Methods & results: In silico docking and molecular dynamics simulation were performed using Glide and Desmond, respectively; while saturation transfer difference NMR was performed using GDP-bound K-RasG12V. SRJ23 was found to bind strongly and selectively to K-RasG12V, by anchoring to a binding pocket (namely p2) principally via hydrogen bond and hydrophobic interactions. The saturation transfer difference NMR analysis revealed the proximity of protons of functional moieties in SRJ23 to K-RasG12V, suggesting positive binding. Conclusion: SRJ23 binds strongly and interacts stably with K-RasG12V to exhibit its inhibitory activity

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

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

AMPK and endothelial nitric oxide synthase signaling regulates K-Ras plasma membrane interactions via cyclic GMP-dependent protein kinase 2

K-Ras must localize to the plasma membrane and be arrayed in nanoclusters for biological activity. We show here that K-Ras is a substrate for cyclic GMP-dependent protein kinases (PKGs). In intact cells, activated PKG2 selectively colocalizes with K-Ras on the plasma membrane and phosphorylates K-Ras at Ser181 in the C-terminal polybasic domain. K-Ras phosphorylation by PKG2 is triggered by activation of AMP-activated protein kinase (AMPK) and requires endothelial nitric oxide synthase and soluble guanylyl cyclase. Phosphorylated K-Ras reorganizes into distinct nanoclusters that retune the signal output. Phosphorylation acutely enhances K-Ras plasma membrane affinity, but phosphorylated K-Ras is progressively lost from the plasma membrane via endocytic recycling. Concordantly, chronic pharmacological activation of AMPK → PKG2 signaling with mitochondrial inhibitors, nitric oxide, or sildenafil inhibits proliferation of K-Ras-positive non-small cell lung cancer cells. The study shows that K-Ras is a target of a metabolic stress-signaling pathway that can be leveraged to inhibit oncogenic K-Ras function

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

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