GTP-Bound N-Ras Conformational States and Substates Are Modulated by Membrane and Point Mutation.

Oncogenic Ras proteins are known to present multiple conformational states, as reported by the great variety of crystallographic structures. The GTP-bound states are grouped into two main states: the "inactive" state 1 and the "active" state 2. Recent reports on H-Ras have shown that state 2 exhibits two substates, directly related to the orientation of Tyr32: toward the GTP-bound pocket and outwards. In this paper, we show that N-Ras exhibits another substate of state 2, related to a third orientation of Tyr32, toward Ala18 and parallel to the GTP-bound pocket. We also show that this substate is highly sampled in the G12V mutation of N-Ras and barely present in its wild-type form, and that the G12V mutation prohibits the sampling of the GTPase-activating protein (GAP) binding substate, rendering this mutation oncogenic. Furthermore, using molecular dynamics simulations, we explore the importance of the membrane on N-Ras' conformational state dynamics and its strong influence on Ras protein stability. Moreover, the membrane has a significant influence on the conformational (sub)states sampling of Ras. This, in turn, is of crucial importance in the activation/deactivation cycle of Ras, due to the binding of guanine nucleotide exchange factor proteins (GEFs)/GTPase-activating proteins (GAPs).

Modulating short tryptophan- and arginine-rich peptides activity by substitution with histidine.

BACKGROUND: High antimicrobial efficacy of short tryptophan-and arginine-rich peptides makes them good candidates in the fight against pathogens. Substitution of tryptophan and arginine by histidine could be used to modulate the peptides efficacy by optimizing their structures. METHODS: The peptide (RRWWRWWRR), reported to showed good antimicrobial efficacy, was used as template, seven new analogs being designed substituting tryptophan or arginine with histidine. The peptides' efficacy was tested against E. coli, B. subtilis and S. aureus. The cytotoxicity and hemolytic effect were evaluated and the therapeutic index was inferred for each peptide. Atomic force microscopy and molecular simulation were used to analyze the effects of peptides on bacterial membrane. RESULTS: The substitution of tryptophan by histidine proved to strongly modulate the antimicrobial activity, mainly by changing the peptide-to-membrane binding energy. The substitution of arginine has low effect on the antimicrobial efficacy. The presence of histidine residue reduced the cytotoxic and hemolytic activity of the peptides in some cases maintaining the same efficacy against bacteria. The peptides' antimicrobial activity was correlated to the 3D-hydrophobic moment and to a simple structure-based packing parameter. CONCLUSION: The results show that some of these peptides have the potential to become good candidates to fight against bacteria. The substitution by histidine proved to fine tune the therapeutic index allowing the optimization of the peptide structure mainly by changing its binding energy and 3D-hydrophobic moment. GENERAL SIGNIFICANCE: The short tryptophan reach peptides therapeutic index can be maximized using the histidine substitution to optimize their structure.

Kinetic Monte Carlo and cellular particle dynamics simulations of multicellular systems.

Computer modeling of multicellular systems has been a valuable tool for interpreting and guiding in vitro experiments relevant to embryonic morphogenesis, tumor growth, angiogenesis and, lately, structure formation following the printing of cell aggregates as bioink particles. Here we formulate two computer simulation methods: (1) a kinetic Monte Carlo (KMC) and (2) a cellular particle dynamics (CPD) method, which are capable of describing and predicting the shape evolution in time of three-dimensional multicellular systems during their biomechanical relaxation. Our work is motivated by the need of developing quantitative methods for optimizing postprinting structure formation in bioprinting-assisted tissue engineering. The KMC and CPD model parameters are determined and calibrated by using an original computational-theoretical-experimental framework applied to the fusion of two spherical cell aggregates. The two methods are used to predict the (1) formation of a toroidal structure through fusion of spherical aggregates and (2) cell sorting within an aggregate formed by two types of cells with different adhesivities.

GxxxG motifs, phenylalanine, and cholesterol guide the self-association of transmembrane domains of ErbB2 receptors.

GxxxG motifs are common in transmembrane domains of membrane proteins and are often introduced to artificial peptides to inhibit or promote association to stable structures. The transmembrane domain of ErbB2 presents two separate such motifs that are proposed to be connected to stability and activity of the dimer. Using molecular simulations, we show that these sequences play a critical role during the recognition stage, forming transient complexes that lead to stable dimers. In pure phospholipid bilayers association occurs by contacts formed at the C-terminus promoted by the presence of phenylalanine residues. Helices subsequently rotate to eventually pack at short separations favored by lipid entropic contributions. In contrast, at intermediate cholesterol concentrations, a different pathway is followed that involves dimers with a weaker interface toward the N-terminus. However, at high cholesterol content, a switch toward the C-terminus is observed with an overall nonmonotonic change of the dimerization affinity. This conformational switch modulated by cholesterol has important implications on the thermodynamic, structural, and kinetic characteristics of helix-helix association in lipid membranes.

Self-association of models of transmembrane domains of ErbB receptors in a lipid bilayer.

Association of transmembrane (TM) helices is facilitated by the close packing of small residues present along the amino-acid sequence. Extensive studies have established the role of such small residue motifs (GxxxG) in the dimerization of Glycophorin A (GpA) and helped to elucidate the association of TM domains in the epidermal growth factor family of receptors (ErbBs). Although membrane-mediated interactions are known to contribute under certain conditions to the dimerization of proteins, their effect is often considered nonspecific, and any potential dependence on protein sequence has not been thoroughly investigated. We recently reported that the association of GpA is significantly assisted by membrane-induced contributions as quantified in different lipid bilayers. Herein we extend our studies to explore the origin of these effects and quantify their magnitude using different amino-acid sequences in the same lipid environment. Using a coarse-grained model that accounts for amino-acid specificity, we perform extensive parallel Monte Carlo simulations of ErbB homodimerization in dipalmitoyl-phosphatidylcholine lipid bilayers. A detailed characterization of dimer formation and estimates of the free energy of association reveal that the TM domains show a significant affinity to self-associate in lipid bilayers, in qualitative agreement with experimental findings. The presence of GxxxG motifs enhances favorable protein-protein interactions at short separations. However, the lipid-induced attraction presents a more complex character than anticipated. Depending on the interfacial residues, lipid-entropic contributions support a decrease of separation or a parallel orientation to the membrane normal, with important implications for protein function.

Segregation of negatively charged phospholipids by the polycationic and farnesylated membrane anchor of Kras.

The Kras protein, a member of the Ras family of bio-switches that are frequently mutated in cancer and developmental disorders, becomes functional when anchored to the inner surface of the plasma membrane. It is well known that membrane attachment involves the farnesylated and poylcationic C-terminus of the protein. However, little is known about the structure of the complex and the specific protein-lipid interactions that are responsible for the binding. On the basis of data from extensive (>0.55 µs) molecular dynamics simulations of multiple Kras anchors in bilayers of POPC/POPG lipids (4:1 ratio), we show that, as expected, Kras is tethered to the bilayer surface by specific lysine-POPG salt bridges and by nonspecific farnesyl-phospholipid van der Waals interactions. Unexpectedly, however, only the C-terminal five of the eight Kras Lys side chains were found to directly interact with the bilayer, with the N-terminal ones staying in water. Furthermore, the positively charged Kras anchors pull the negatively charged POPG lipids together, leading to the clustering of the POPG lipids around the proteins. This selective Kras-POPG interaction is directly related to the specific geometry of the backbone, which exists in two major conformational states: 1), a stable native-like ensemble of structures characterized by an extended geometry with a pseudohelical turn; and 2), less stable nonnative ensembles of conformers characterized by severely bent geometries. Finally, although the interface-bound anchor has little effect on the overall structure of the bilayer, it induces local thinning within a persistence length of ∼12 Å. Our results thus go beyond documenting how Kras attaches to a mixed bilayer of charged and neutral lipids; they highlight a fascinating process of protein-induced lipid sorting coupled with the (re)shaping of a surface-bound protein by the host lipids.