Assessment of mutation probabilities of KRAS G12 missense mutants and their long-timescale dynamics by atomistic molecular simulations and Markov state modeling.

A mutated KRAS protein is frequently observed in human cancers. Traditionally, the oncogenic properties of KRAS missense mutants at position 12 (G12X) have been considered as equal. Here, by assessing the probabilities of occurrence of all KRAS G12X mutations and KRAS dynamics we show that this assumption does not hold true. Instead, our findings revealed an outstanding mutational bias. We conducted a thorough mutational analysis of KRAS G12X mutations and assessed to what extent the observed mutation frequencies follow a random distribution. Unique tissue-specific frequencies are displayed with specific mutations, especially with G12R, which cannot be explained by random probabilities. To clarify the underlying causes for the nonrandom probabilities, we conducted extensive atomistic molecular dynamics simulations (170 µs) to study the differences of G12X mutations on a molecular level. The simulations revealed an allosteric hydrophobic signaling network in KRAS, and that protein dynamics is altered among the G12X mutants and as such differs from the wild-type and is mutation-specific. The shift in long-timescale conformational dynamics was confirmed with Markov state modeling. A G12X mutation was found to modify KRAS dynamics in an allosteric way, which is especially manifested in the switch regions that are responsible for the effector protein binding. The findings provide a basis to understand better the oncogenic properties of KRAS G12X mutants and the consequences of the observed nonrandom frequencies of specific G12X mutations.

Calcium Directly Regulates Phosphatidylinositol 4,5-Bisphosphate Headgroup Conformation and Recognition.

The orchestrated recognition of phosphoinositides and concomitant intracellular release of Ca2+ is pivotal to almost every aspect of cellular processes, including membrane homeostasis, cell division and growth, vesicle trafficking, as well as secretion. Although Ca2+ is known to directly impact phosphoinositide clustering, little is known about the molecular basis for this or its significance in cellular signaling. Here, we study the direct interaction of Ca2+ with phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2), the main lipid marker of the plasma membrane. Electrokinetic potential measurements of PI(4,5)P2 containing liposomes reveal that Ca2+ as well as Mg2+ reduce the zeta potential of liposomes to nearly background levels of pure phosphatidylcholine membranes. Strikingly, lipid recognition by the default PI(4,5)P2 lipid sensor, phospholipase C delta 1 pleckstrin homology domain (PLC δ1-PH), is completely inhibited in the presence of Ca2+, while Mg2+ has no effect with 100 nm liposomes and modest effect with giant unilamellar vesicles. Consistent with biochemical data, vibrational sum frequency spectroscopy and atomistic molecular dynamics simulations reveal how Ca2+ binding to the PI(4,5)P2 headgroup and carbonyl regions leads to confined lipid headgroup tilting and conformational rearrangements. We rationalize these findings by the ability of calcium to block a highly specific interaction between PLC δ1-PH and PI(4,5)P2, encoded within the conformational properties of the lipid itself. Our studies demonstrate the possibility that switchable phosphoinositide conformational states can serve as lipid recognition and controlled cell signaling mechanisms.

doGlycans-Tools for Preparing Carbohydrate Structures for Atomistic Simulations of Glycoproteins, Glycolipids, and Carbohydrate Polymers for GROMACS.

Carbohydrates constitute a structurally and functionally diverse group of biological molecules and macromolecules. In cells they are involved in, e.g., energy storage, signaling, and cell-cell recognition. All of these phenomena take place in atomistic scales, thus atomistic simulation would be the method of choice to explore how carbohydrates function. However, the progress in the field is limited by the lack of appropriate tools for preparing carbohydrate structures and related topology files for the simulation models. Here we present tools that fill this gap. Applications where the tools discussed in this paper are particularly useful include, among others, the preparation of structures for glycolipids, nanocellulose, and glycans linked to glycoproteins. The molecular structures and simulation files generated by the tools are compatible with GROMACS.

Apolipoprotein A-I mimetic peptide 4F blocks sphingomyelinase-induced LDL aggregation.

Lipolytic modification of LDL particles by SMase generates LDL aggregates with a strong affinity for human arterial proteoglycans and may so enhance LDL retention in the arterial wall. Here, we evaluated the effects of apoA-I mimetic peptide 4F on structural and functional properties of the SMase-modified LDL particles. LDL particles with and without 4F were incubated with SMase, after which their aggregation, structure, and proteoglycan binding were analyzed. At a molar ratio of L-4F to apoB-100 of 2.5 to 20:1, 4F dose-dependently inhibited SMase-induced LDL aggregation. At a molar ratio of 20:1, SMase-induced aggregation was fully blocked. Binding of 4F to LDL particles inhibited SMase-induced hydrolysis of LDL by 10% and prevented SMase-induced LDL aggregation. In addition, the binding of the SMase-modified LDL particles to human aortic proteoglycans was dose-dependently inhibited by pretreating LDL with 4F. The 4F stabilized apoB-100 conformation and inhibited SMase-induced conformational changes of apoB-100. Molecular dynamic simulations showed that upon binding to protein-free LDL surface, 4F locally alters membrane order and fluidity and induces structural changes to the lipid layer. Collectively, 4F stabilizes LDL particles by preventing the SMase-induced conformational changes in apoB-100 and so blocks SMase-induced LDL aggregation and the resulting increase in LDL retention.

How anacetrapib inhibits the activity of the cholesteryl ester transfer protein? Perspective through atomistic simulations.

Cholesteryl ester transfer protein (CETP) mediates the reciprocal transfer of neutral lipids (cholesteryl esters, triglycerides) and phospholipids between different lipoprotein fractions in human blood plasma. A novel molecular agent known as anacetrapib has been shown to inhibit CETP activity and thereby raise high density lipoprotein (HDL)-cholesterol and decrease low density lipoprotein (LDL)-cholesterol, thus rendering CETP inhibition an attractive target to prevent and treat the development of various cardiovascular diseases. Our objective in this work is to use atomistic molecular dynamics simulations to shed light on the inhibitory mechanism of anacetrapib and unlock the interactions between the drug and CETP. The results show an evident affinity of anacetrapib towards the concave surface of CETP, and especially towards the region of the N-terminal tunnel opening. The primary binding site of anacetrapib turns out to reside in the tunnel inside CETP, near the residues surrounding the N-terminal opening. Free energy calculations show that when anacetrapib resides in this area, it hinders the ability of cholesteryl ester to diffuse out from CETP. The simulations further bring out the ability of anacetrapib to regulate the structure-function relationships of phospholipids and helix X, the latter representing the structural region of CETP important to the process of neutral lipid exchange with lipoproteins. Altogether, the simulations propose CETP inhibition to be realized when anacetrapib is transferred into the lipid binding pocket. The novel insight gained in this study has potential use in the development of new molecular agents capable of preventing the progression of cardiovascular diseases.

Membrane-Dependent Binding and Entry Mechanism of Dopamine into Its Receptor.

Synaptic neurotransmission has recently been proposed to function via either a membrane-independent or a membrane-dependent mechanism, depending on the neurotransmitter type. In the membrane-dependent mechanism, amphipathic neurotransmitters first partition to the lipid headgroup region and then diffuse along the membrane plane to their membrane-buried receptors. However, to date, this mechanism has not been demonstrated for any neurotransmitter-receptor complex. Here, we combined isothermal calorimetry measurements with a diverse set of molecular dynamics simulation methods to investigate the partitioning of an amphipathic neurotransmitter (dopamine) and the mechanism of its entry into the ligand-binding site. Our results show that the binding of dopamine to its receptor is consistent with the membrane-dependent binding and entry mechanism. Both experimental and simulation results showed that dopamine favors binding to lipid membranes especially in the headgroup region. Moreover, our simulations revealed a ligand-entry pathway from the membrane to the binding site. This pathway passes through a lateral gate between transmembrane alpha-helices 5 and 6 on the membrane-facing side of the protein. All in all, our results demonstrate that dopamine binds to its receptor by a membrane-dependent mechanism, and this is complemented by the more traditional binding mechanism directly through the aqueous phase. The results suggest that the membrane-dependent mechanism is common in other synaptic receptors, too.

Negatively Charged Gangliosides Promote Membrane Association of Amphipathic Neurotransmitters.

Lipophilic neurotransmitters (NTs) such as dopamine are chemical messengers enabling neurotransmission by adhering onto the extracellular surface of the post-synaptic membrane in a synapse, followed by binding to their receptors. Previous studies have shown that the strength of the NT-membrane association is dependent on the lipid composition of the membrane. Negatively charged lipids such as phosphatidylserine, phosphatidylglycerol, and phosphatidic acid have been indicated to promote NT-membrane binding, however these anionic lipids reside almost exclusively in the intracellular leaflet of the post-synaptic membrane instead of the extracellular leaflet facing the synaptic cleft. Meanwhile, the extracellular leaflet is relatively rich in biologically relevant anionic gangliosides such as monosialotetrahexosylganglioside (GM1), yet the role of gangliosides in NT-membrane association is not clear. Here, we explored the role of GM1 in modulating the binding of dopamine and histamine (as amphipathic/cationic NTs) as well as acetylcholine (as a hydrophilic/cationic NT) with the post-synaptic membrane surface. Atomistic molecular dynamics simulations and free energy calculations indicated that GM1 fosters membrane association of histamine and dopamine. For acetylcholine, this effect was not observed. The in silico results suggest that gangliosides form a charge-based vestibule in front of the post-synaptic membrane, attracting amphipathic NTs to the vicinity of the membrane. The results also stress the importance to understand the significance of the structural details of NTs, as exemplified by the GM1-acetylcholine interaction. In a larger context, the NT-membrane adherence, coupled to lateral diffusion in the membrane plane, is proposed to improve neurotransmission efficiency by advancing NT entry into the membrane-embedded ligand-binding sites.

Calcium Assists Dopamine Release by Preventing Aggregation on the Inner Leaflet of Presynaptic Vesicles.

In this study, the dopamine-lipid bilayer interactions were probed with three physiologically relevant ion compositions using atomistic molecular dynamics simulations and free energy calculations. The in silico results indicate that calcium is able to decrease significantly the binding of dopamine to a neutral (zwitterionic) phosphatidylcholine lipid bilayer model mimicking the inner leaflet of a presynaptic vesicle. We argue that the observed calcium-induced effect is likely in crucial role in the neurotransmitter release from the presynaptic vesicles docked in the active zone of nerve terminals. The inner leaflets of presynaptic vesicles, which are responsible for releasing neurotransmitters into the synaptic cleft, are mainly composed of neutral lipids such as phosphatidylcholine and phosphatidylethanolamine. The neutrality of the lipid head group region, enhanced by a low pH level, should limit membrane aggregation of transmitters. In addition, the simulations suggest that the high calcium levels inside presynaptic vesicles prevent even the most lipophilic transmitters such as dopamine from adhering to the inner leaflet surface, thus rendering unhindered neurotransmitter release feasible.

Phase Partitioning of GM1 and Its Bodipy-Labeled Analog Determine Their Different Binding to Cholera Toxin.

Driven by interactions between lipids and proteins, biological membranes display lateral heterogeneity that manifests itself in a mosaic of liquid-ordered (Lo) or raft, and liquid-disordered (Ld) or non-raft domains with a wide range of different properties and compositions. In giant plasma membrane vesicles and giant unilamellar vesicles, specific binding of Cholera Toxin (CTxB) to GM1 glycolipids is a commonly used strategy to label raft domains or Lo membrane environments. However, these studies often use acyl-chain labeled bodipy-GM1 (bdGM1), whose headgroup accessibility and membrane order or phase partitioning may differ from those of GM1, rendering the interpretation of CTxB binding data quite problematic. To unravel the molecular basis of CTxB binding to GM1 and bdGM1, we explored the partitioning and the headgroup presentation of these gangliosides in the Lo and Ld phases using atomistic molecular dynamics simulations complemented by CTxB binding experiments. The conformation of both GM1 and bdGM1 was shown to be largely similar in the Lo and Ld phases. However, bdGM1 showed reduction in receptor availability when reconstituted into synthetic bilayer mixtures, highlighting that membrane phase partitioning of the gangliosides plays a considerable role in CTxB binding. Our results suggest that the CTxB binding is predominately modulated by the partitioning of the receptor to an appropriate membrane phase. Further, given that the Lo and Ld partitioning of bdGM1 differs from those of GM1, usage of bdGM1 for studying GM1 behavior in cells can lead to invalid interpretation of experimental data.

Effect of piroxicam on lipid membranes: Drug encapsulation and gastric toxicity aspects.

Uptake of piroxicam, a non-steroidal anti-inflammatory drug, from the intestines after oral intake is limited due to its low solubility and its wide use is associated with several side effects related to the gastrointestinal tract. In this study, all-atom molecular dynamics (MD) simulations and fluorescent spectroscopy were employed to investigate the interaction of piroxicam in neutral, zwitterionic, and cationic forms with lipid bilayers composed of phosphatidylcholine, cholesterol, and PEGylated lipids. Our study was aimed to assess the potential for encapsulation of piroxicam in liposomal carriers and to shed more light on the process of gastrointestinal tract injury by the drug. Through both the MD simulations and laser scanning confocal microscopy, we have demonstrated that all forms of piroxicam can associate with the lipid bilayers and locate close to the water-membrane interface. Conventional liposomes used in drug delivery are usually stabilized by the addition of cholesterol and have their bloodstream lifetime extended through the inclusion of PEGylated lipids in the formulation to create a protective polymer corona. For this reason, we tested the effect of these two modifications on the behavior of piroxicam in the membrane. When the bilayer was PEGylated, piroxicam localize to the PEG layer and within the lipid headgroup region. This suggests that PEGylated liposomes are capable of carrying a larger quantity of piroxicam than the conventional ones.

PIP2 and Talin Join Forces to Activate Integrin.

Integrins are major players in cell adhesion and migration, and malfunctions in controlling their activity are associated with various diseases. Nevertheless, the details of integrin activation are not completely understood, and the role of lipids in the process is largely unknown. Herein, we show using atomistic molecular dynamics simulations that the interplay of phosphatidylinositol 4,5-bisphosphate (PIP2) and talin may directly alter the conformation of integrin αIIbß3. Our results provide a new perspective on the role of PIP2 in integrin activation and indicate that the charged PIP2 lipid headgroup can perturb a clasp at the cytoplasmic face of the integrin heterodimer.

Mutually Exclusive Roles of SHARPIN in Integrin Inactivation and NF-κB Signaling.

SHANK-associated RH domain interactor (SHARPIN) inhibits integrins through interaction with the integrin α-subunit. In addition, SHARPIN enhances nuclear factor-kappaB (NF-κB) activity as a component of the linear ubiquitin chain assembly complex (LUBAC). However, it is currently unclear how regulation of these seemingly different roles is coordinated. Here, we show that SHARPIN binds integrin and LUBAC in a mutually exclusive manner. We map the integrin binding site on SHARPIN to the ubiquitin-like (UBL) domain, the same domain implicated in SHARPIN interaction with LUBAC component RNF31 (ring finger protein 31), and identify two SHARPIN residues (V267, L276) required for both integrin and RNF31 regulation. Accordingly, the integrin α-tail is capable of competing with RNF31 for SHARPIN binding in vitro. Importantly, the full SHARPIN RNF31-binding site contains residues (F263A/I272A) that are dispensable for SHARPIN-integrin interaction. Importantly, disrupting SHARPIN interaction with integrin or RNF31 abolishes SHARPIN-mediated regulation of integrin or NF-κB activity, respectively. Altogether these data suggest that the roles of SHARPIN in inhibiting integrin activity and supporting linear ubiquitination are (molecularly) distinct.

PEGylated Liposomes as Carriers of Hydrophobic Porphyrins.

Sterically stabilized liposomes (SSLs) (PEGylated liposomes) are applied as effective drug delivery vehicles. Understanding the interactions between hydrophobic compounds and PEGylated membranes is therefore important to determine the effectiveness of PEGylated liposomes for delivery of drugs or other bioactive substances. In this study, we have combined fluorescence quenching analysis (FQA) experiments and all-atom molecular dynamics (MD) simulations to study the effect of membrane PEGylation on the location and orientation of 5,10,15,20-tetrakis(4-hydroxyphenyl)porphyrin (p-THPP) that has been used in our study as a model hydrophobic compound. First, we consider the properties of p-THPP in the presence of different fluid phosphatidylcholine bilayers that we use as model systems for protein-free cell membranes. Next, we studied the interaction between PEGylated membranes and p-THPP. Our MD simulation results indicated that the arrangement of p-THPP within zwitterionic membranes is dependent on their free volume, and p-THPP solubilized in PEGylated liposomes is localized in two preferred positions: deep within the membrane (close to the center of the bilayer) and in the outer PEG corona (p-THPP molecules being wrapped with the polymer chains). Fluorescence quenching methods confirmed the results of atomistic MD simulations and showed two populations of p-THPP molecules as in MD simulations. Our results provide both an explanation for the experimental observation that PEGylation improves the drug-loading efficiency of membranes and also a more detailed molecular-level description of the interactions between porphyrins and lipid membranes.

Effect of PEGylation on drug entry into lipid bilayer.

Poly(ethylene glycol) (PEG) is a polymer commonly used for functionalization of drug molecules to increase their bloodstream lifetime, hence efficacy. However, the interactions between the PEGylated drugs and biomembranes are not clearly understood. In this study, we employed atomic-scale molecular dynamics (MD) simulations to consider the behavior of two drug molecules functionalized with PEG (tetraphenylporphyrin used in cancer phototherapy and biochanin A belonging to the isoflavone family) in the presence of a lipid bilayer. The commonly held view is that functionalization of a drug molecule with a polymer acts as an entropic barrier, inhibiting the penetration of the drug molecule through a cell membrane. Our results indicate that in the bloodstream there is an additional source of electrostatic repulsive interactions between the PEGylated drugs and the lipid bilayer. Both the PEG chain and lipids can bind Na(+) ions, thus effectively becoming positively charged molecules. This leads to an extra repulsive effect resulting from the presence of salt in the bloodstream. Thus, our study sheds further light on the role of PEG in drug delivery.

Interaction of hematoporphyrin with lipid membranes.

Natural or synthetic porphyrins are being used as photosensitizers in photodiagnosis (PD) and photodynamic therapy (PDT) of malignancies and some other diseases. Understanding the interactions between porphyrins and cell membranes is therefore important to rationalize the uptake of photosensitizers and their passive transport through cell membranes. In this study, we consider the properties of hematoporphyrin (Hp), a well-known photosensitizer for PD and PDT, in the presence of a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer that we use as a model system for protein-free cell membranes. For this purpose, we employed 200 ns atomic-scale molecular dynamics (MD) simulations for five systems containing the neutral (Hp(0)) or the dianionic form (Hp(2-)) of Hp and the POPC bilayer. MD simulations allowed one to estimate the position, orientation, and dynamics of Hp molecules inside the membrane. The dye molecules were found to reside in the phospholipid headgroup area close to the carbonyl groups of the POPC acyl chains. Their orientations were dependent on the protonation state of two propionic groups. Hp(2-) was found to have a lower affinity to enter the membrane than the neutral form. The dianions, being in the aqueous phase, formed stable dimers with a strictly determined geometry. Our results fully supported the experimental data and provide a more detailed molecular-level description of the interactions of photosensitizers with lipid membranes.

Study of interaction between PEG carrier and three relevant drug molecules: piroxicam, paclitaxel, and hematoporphyrin.

Molecular dynamics simulation has been used to study the specific interactions between poly(ethylene glycol) (PEG) and three drug molecules for which PEG is used to aid delivery: paclitaxel and piroxicam, where PEG is a carrier agent, and hematoporphyrin, where PEG is covalently attached to form a "stealth shield". Simulating at physiological salt concentration, we found no evidence of any specific interaction between paclitaxel or piroxicam with PEG, but found a strong interaction for the case of hematoporphyrin. This interaction is lipophilic in nature, between the nonpolar (CH(2))(2) groups of the PEG and the porphin ring of the hematoporphyrin. This interaction was found to be strong enough that the PEG aggregated to the hematoporphyrin, independent of whether or not it was covalently bound. Interestingly, when the simulation was repeated in absence of salt we found evidence of this interaction being weakened. This led us to hypothesize a previously unforeseen mechanism: interaction with salt cations cause the PEG to coil around the salt ions, each ion binding to many PEG oxygens, increasing the exposure of the nonpolar ethylene groups, thus increasing the effective hydrophobicity of PEG. The Hydrophobic ethylene groups of the PEG chains adhere strongly to the hydrophobic porphin ring. Experiments involving absorption spectra measurements were conducted, and these results also indicated that presence of salt at physiological level increases the effective attractive interaction between PEG and hematoporphyrin. Taken together, our results demonstrate that while PEG, due to its solubility in both polar and nonpolar solvents, may act as a dissolution aid for paclitaxel and piroxicam, of the three drug molecules studied it will only have a protective role for the case of the hematoporphyrin.