Recent studies found that membrane-bound K-Ras dimers are important for biological function. However, the structure and thermodynamic stability of these complexes remained unknown because they are hard to probe by conventional approaches. Combining data from a wide range of computational and experimental approaches, here we describe the structure, dynamics, energetics and mechanism of assembly of multiple K-Ras dimers. Utilizing a range of techniques for the detection of reactive surfaces, protein-protein docking and molecular simulations, we found that two largely polar and partially overlapping surfaces underlie the formation of multiple K-Ras dimers. For validation we used mutagenesis, electron microscopy and biochemical assays under non-denaturing conditions. We show that partial disruption of a predicted interface through charge reversal mutation of apposed residues reduces oligomerization while introduction of cysteines at these positions enhanced dimerization likely through the formation of an intermolecular disulfide bond. Free energy calculations indicated that K-Ras dimerization involves direct but weak protein-protein interactions in solution, consistent with the notion that dimerization is facilitated by membrane binding. Taken together, our atomically detailed analyses provide unique mechanistic insights into K-Ras dimer formation and membrane organization as well as the conformational fluctuations and equilibrium thermodynamics underlying these processes.

Mounting evidence suggests that the neuronal cell membrane is the main site of oligomer-mediated neuronal toxicity of amyloid-β peptides in Alzheimer’s disease. To gain a detailed understanding of the mutual interference of amyloid-β oligomers and the neuronal membrane, we carried out microseconds of all-atom molecular dynamics (MD) simulations on the dimerization of amyloid-β (Aβ)42 in the aqueous phase and in the presence of a lipid bilayer mimicking the in vivo composition of neuronal membranes. The dimerization in solution is characterized by a random coil to β-sheet transition that seems on pathway to amyloid aggregation, while the interactions with the neuronal membrane decrease the order of the Aβ42 dimer by attenuating its propensity to form a β-sheet structure. The main lipid interaction partners of Aβ42 are the surface-exposed sugar groups of the gangliosides GM1. As the neurotoxic activity of amyloid oligomers increases with oligomer order, these results suggest that GM1 is neuroprotective against Aβ-mediated toxicity.

We conducted over 150 ns of simulation of a protegrin-1 octamer pore in a lipid bilayer composed of palmitoyloleoyl-phosphatidylethanolamine (POPE) and palmitoyloleoyl-phosphatidylglycerol (POPG) lipids mimicking the inner membrane of a bacterial cell. The simulations improve on a model of a pore proposed from recent NMR experiments and provide a coherent understanding of the molecular mechanism of antimicrobial activity. Although lipids tilt somewhat toward the peptides, the simulated protegrin-1 pore more closely follows the barrel-stave model than the toroidal-pore model. The movement of ions is investigated through the pore. The pore selectively allows negatively charged chloride ions to pass through at an average rate of one ion every two nanoseconds. Only two events are observed of sodium ions crossing through the pore. The potential of mean force is calculated for the water and both ion types. It is determined that the chloride ions move through the pore with ease, similarly to the water molecules with the exception of a zone of restricted movement midway through the pore. In bacteria, ions moving through the pore will compromise the integrity of the transmembrane potential. Without the transmembrane potential as a countermeasure, water will readily flow inside the higher osmolality cytoplasm. We determine that the diffusivity of water through a single PG-1 pore is sufficient to cause fast cell death by osmotic lysis.

Self-assembly of plasma membrane-associated Ras GTPases has major implications to the regulation of cell signaling. However, the structural basis of homo-oligomerization and the fractional distribution of oligomeric states remained undetermined. We have addressed these issues by deciphering the distribution of dimers and higher-order oligomers of K-Ras4B, the most frequently mutated Ras isoform in human cancers. We focused on the constitutively active G12V K-Ras and two of its variants, K101E and K101C/E107C, which respectively destabilize and stabilize oligomers. Using raster image correlation spectroscopy and number and brightness analysis combined with fluorescence recovery after photobleaching, fluorescence correlation spectroscopy and electron microscopy in live cells, we show that G12V K-Ras exists as a mixture of monomers, dimers and larger oligomers, while the K101E mutant is predominantly monomeric and K101C/E107C is dominated by oligomers. This observation demonstrates the ability of K-Ras to exist in multiple oligomeric states whose population can be altered by interfacial mutations. Using molecular modeling and simulations we further show that K-Ras uses two partially overlapping interfaces to form compositionally and topologically diverse oligomers. Our results thus provide the first detailed insight into the multiplicity, structure, and membrane organization of K-Ras homomers.

Protegrin peptides are potent antimicrobial agents believed to act against a variety of pathogens by forming nonselective transmembrane pores in the bacterial cell membrane. We have employed 3D Poisson-Nernst-Planck (PNP) calculations to determine the steady-state ion conduction characteristics of such pores at applied voltages in the range of −100 to +100 mV in 0.1 M KCl bath solutions. We have tested a variety of pore structures extracted from molecular dynamics (MD) simulations based on an experimentally proposed octomeric pore structure. The computed single-channel conductance values were in the range of 290–680 pS. Better agreement with the experimental range of 40–360 pS was obtained using structures from the last 40 ns of the MD simulation, where conductance values range from 280 to 430 pS. We observed no significant variation of the conductance with applied voltage in any of the structures that we tested, suggesting that the voltage dependence observed experimentally is a result of voltage-dependent channel formation rather than an inherent feature of the open pore structure. We have found the pore to be highly selective for anions, with anionic to cationic current ratios (ICl−/IK+) on the order of 103. This is consistent with the highly cationic nature of the pore but surprisingly in disagreement with the experimental finding of only slight anionic selectivity. We have additionally tested the sensitivity of our PNP model to several parameters and found the ion diffusion coefficients to have a significant influence on conductance characteristics. The best agreement with experimental data was obtained using a diffusion coefficient for each ion set to 10% of the bulk literature value everywhere inside the channel, a scaling used by several other studies employing PNP calculations. Overall, this work presents a useful link between previous work focused on the structure of protegrin pores and experimental efforts aimed at investigating their conductance characteristics.

Probe-based or mixed solvent molecular dynamics simulation is a useful approach for the identification and characterization of druggable sites in drug targets. However, thus far the method has been applied only to soluble proteins. A major reason for this is the potential effect of the probe molecules on membrane structure. We have developed a technique to overcome this limitation that entails modification of force field parameters to reduce a few pairwise non-bonded interactions between selected atoms of the probe molecules and bilayer lipids. We used the resulting technique, termed pMD-membrane, to identify allosteric ligand binding sites on the G12D and G13D oncogenic mutants of the K-Ras protein bound to a negatively charged lipid bilayer. In addition, we show that differences in probe occupancy can be used to quantify changes in the accessibility of druggable sites due to conformational changes induced by membrane binding or mutation.

To investigate the stability and functional role of long-residence water molecules in the Q61H variant of the signaling protein K-ras, we analyzed all available Ras crystal structures and conformers derived from a series of independent explicit solvent molecular dynamics (MD) simulations totaling 1.76 µs. We show that the protein samples a different region of phase space in the presence and absence of several crystallographically conserved and buried water molecules. The dynamics of these waters is coupled with the local as well as the global motions of the protein, in contrast to less buried waters whose exchange with bulk is only loosely coupled with the motion of loops in their vicinity. Aided by two novel reaction coordinates involving the distance (d) between the Cα atoms of G60 at switch 2 and G10 at the P-loop and the N-Cα-C-O dihedral (ξ) of G60, we further show that three water molecules located in lobe1, at the interface between the lobes and at lobe2, are involved in the relative motion of residues at the two lobes of Q61H K-ras. Moreover, a d/ξ plot classifies the available Ras x-ray structures and MD-derived K-ras conformers into active GTP-, intermediate GTP-, inactive GDP-bound, and nucleotide-free conformational states. The population of these states and the transition between them is modulated by water-mediated correlated motions involving the functionally critical switch 2, P-loop and helix 3. These results suggest that water molecules act as allosteric ligands to induce a population shift among distinct switch 2 conformations that differ in effector recognition.

Bile acids are powerful detergents that emulsify and solubilize lipids, vitamins, cholesterol and other molecules in the biliary tract and intestines. It has long been known that bile acids form soluble mixed micelles with lipids. However, the detailed thermodynamic and structural properties of these micelles are not fully understood. This study sheds light on this issue based on results from multiple molecular dynamics simulations of cholic acid (CA) and dodecylphosphocholine (DPC) mixed micelles. We found that CA molecules form aggregates of up to 12 monomers with a mean size of 5–6. In agreement with several previous simulations and earlier predictions, the overall shape of these CA clusters is oblate disk-like such that the methyl groups point toward the core of the aggregate and the hydroxyl groups point away from it. The self-aggregation behavior of the CA clusters in the DPC-CA mixture is similar to the pure CA. Furthermore, the sizes and aggregation numbers of the DPC-CA mixed micelles are linearly dependent on CA molarity. In agreement with the radial shell model of Nichols and Ozarowski [Nichols, J. W.; Ozarowski, J. Biochemistry 1990, 29, 4600]1, our results demonstrate that CA molecules form a wedge between the DPC molecules with their hydroxyl and carboxyl groups facing the aqueous phase while their methyl groups are buried in and face the hydrocarbon core of the DPC micelle. The DPC-CA micelles simulated here tend to be spherical to prolate in shape, with the deviation from spherical geometry significantly increasing with increasing CA:DPC ratio.