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).

No country for old antibiotics! Antimicrobial peptides (AMPs) as next-generation treatment for skin and soft tissue infection.

In recent years, the unprecedented rise of bacterial antibiotic resistance together with the lack of adequate therapies have made the treatment of skin infections and chronic wounds challenging, urging the scientific community to focus on the development of new and more efficient treatment strategies. In this context, there is a growing interest in the use of natural molecules with antimicrobial features, capable of supporting wound healing i.e., antimicrobial peptides (AMPs), for the treatment of skin and soft tissue infections. In this review, we give a short overview of the bacterial skin infections as well as some of the classic treatments used for topical application. We then summarize the AMPs classes, stressing the importance of the appropriate selection of the peptides based on their characteristics and physicochemical properties in order to maximize the antibacterial efficacy of the therapeutic systems against multi-drug resistant pathogens. Additionally, the present paper provides a comprehensive and rigorous assessment of the latest clinical trials investigating the efficacy of AMPs in the treatment of skin and soft tissue infections, highlighting the relevant outcomes. Seeking to obtain novel and improved compounds with synergistic activity, while also decreasing some of the known side effects of AMPs, we present two employed strategies using AMPs: (i) AMPs-conjugated nanosystems for systemic and topical drug delivery systems and (ii) antibiotics-peptide conjugates as a strategy to overcome antibiotics resistance. Finally, an important property of some of the AMPs used in wound treatment is highlighted: their ability to help in wound healing by generally promoting cell proliferation and migration, and in some cases re-epithelialization and angiogenesis among others. Thus, as the pursuit of improvement is an ongoing effort, this work presents the advances made in the treatment of skin and soft tissue infections along with their advantages and limitations, while the still remaining challenges are addressed by providing future prospects and strategies to overcome them.

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.

The gating mechanism of the human aquaporin 5 revealed by molecular dynamics simulations.

Aquaporins are protein channels located across the cell membrane with the role of conducting water or other small sugar alcohol molecules (aquaglyceroporins). The high-resolution X-ray structure of the human aquaporin 5 (HsAQP5) shows that HsAQP5, as all the other known aquaporins, exhibits tetrameric structure. By means of molecular dynamics simulations we analyzed the role of spontaneous fluctuations on the structural behavior of the human AQP5. We found that different conformations within the tetramer lead to a distribution of monomeric channel structures, which can be characterized as open or closed. The switch between the two states of a channel is a tap-like mechanism at the cytoplasmic end which regulates the water passage through the pore. The channel is closed by a translation of the His67 residue inside the pore. Moreover, water permeation rate calculations revealed that the selectivity filter, located at the other end of the channel, regulates the flow rate of water molecules when the channel is open, by locally modifying the orientation of His173. Furthermore, the calculated permeation rates of a fully open channel are in good agreement with the reported experimental value.

Formation and domain partitioning of H-ras peptide nanoclusters: effects of peptide concentration and lipid composition.

Experiments have shown that homologous Ras proteins containing different lipid modification, which is required for membrane binding, form nonoverlapping nanoclusters on the plasma membrane. However, the physical basis for clustering and lateral organization remains poorly understood. We have begun to tackle this issue using coarse-grained molecular dynamics simulations of the H-ras lipid anchor (tH), a triply lipid-modified heptapeptide embedded in a domain-forming mixed lipid bilayer [Janosi L. et al. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 8097]. Here we use the same simulation approach to investigate the effect of peptide concentration and bilayer composition on the clustering and lateral distribution of tH. We found no major difference in the clustering behavior of tH above a certain concentration. However, the simulations predict the existence of a critical concentration below which tH does not form nanoclusters. Moreover, our data demonstrate that cholesterol enhances the stability of tH nanoclusters but is not required for their formation. Finally, analyses of peptide distributions and partition free energies allowed us to quantitatively describe how clustering facilitates the accumulation of tH at the interface between ordered and disordered domains of the simulated bilayer systems. These thermodynamic insights represent some of the key elements for a comprehensive understanding of the molecular basis for the formation and stability of Ras signaling platforms.

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.

Organization, dynamics, and segregation of Ras nanoclusters in membrane domains.

Recent experiments have shown that membrane-bound Ras proteins form transient, nanoscale signaling platforms that play a crucial role in high-fidelity signal transmission. However, a detailed characterization of these dynamic proteolipid substructures by high-resolution experimental techniques remains elusive. Here we use extensive semiatomic simulations to reveal the molecular basis for the formation and domain-specific distribution of Ras nanoclusters. As model systems, we chose the triply lipidated membrane targeting motif of H-ras (tH) and a large bilayer made up of di160-PC (DPPC), di182-PC (DLiPC), and cholesterol. We found that 4-10 tH molecules assemble into clusters that undergo molecular exchange in the sub-µs to µs time scale, depending on the simulation temperature and hence the stability of lipid domains. Driven by the opposite preference of tH palmitoyls and farnesyl for ordered and disordered membrane domains, clustered tH molecules segregate to the boundary of lipid domains. Additionally, a systematic analysis of depalmitoylated and defarnesylated tH variants allowed us to decipher the role of individual lipid modifications in domain-specific nanocluster localization and thereby explain why homologous Ras isoforms form nonoverlapping nanoclusters. Moreover, the localization of tH nanoclusters at domain boundaries resulted in a significantly lower line tension and increased membrane curvature. Taken together, these results provide a unique mechanistic insight into how protein assembly promoted by lipid-modification modulates bilayer shape to generate functional signaling platforms.

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.

In silico predictions of LH2 ring sizes from the crystal structure of a single subunit using molecular dynamics simulations.

Most of the currently known light-harvesting complexes 2 (LH2) rings are formed by 8 or 9 subunits. As of now, questions like "what factors govern the LH2 ring size?" and "are there other ring sizes possible?" remain largely unanswered. Here, we investigate by means of molecular dynamics (MD) simulations and stochastic modeling the possibility of predicting the size of an LH2 ring from the sole knowledge of the high resolution crystal structure of a single subunit. Starting with single subunits of two LH2 rings with known size, that is, an 8-ring from Rs. moliscianum (MOLI) and a 9-ring from Rps. acidophila (ACI), and one with unknown size (referred to as X), we build atomic models of subunit dimers corresponding to assumed 8-, 9-, and 10-ring geometries. After inserting each of the dimers into a lipid-water environment, we determine the preferred angle between the corresponding subunits by three methods: (1) energy minimization, (2) free MD simulations, and (3) potential of mean force calculations. We find that the results from all three methods are consistent with each other, and when taken together, it allows one to predict with reasonable level of confidence the sizes of the corresponding ring structures. One finds that X and ACI very likely form a 9-ring, while MOLI is more likely to form an 8-ring than a 9-ring. Finally, we discuss both the merits and limitations of all three prediction methods.

Importance of the sphingosine base double-bond geometry for the structural and thermodynamic properties of sphingomyelin bilayers.

The precise role of the sphingosine base trans double bond for the unique properties of sphingomyelins (SMs), one of the main lipid components in raftlike structures of biological membranes, has not been fully explored. Several reports comparing the hydration, lipid packing, and hydrogen-bonding behaviors of SM and glycerophospholipid bilayers found remarkable differences overall. However, the atomic interactions linking the double-bond geometry with these thermodynamic and structural changes remained elusive. A recent report on ceramides, which differ from SMs only by their hydroxyl headgroup, has shown that replacing the trans double bond of the sphingosine base by cis weakens the hydrogen-bonding potential of these lipids and thereby alters their biological activity. Based on data from extensive (a total 0.75 µs) atomistic molecular dynamics simulations of bilayers composed of all-trans, all-cis, and a trans/cis (4:1 ratio) racemic mixture of sphingomyelin lipids, here we show that the trans configuration allows for the formation of significantly more hydrogen bonds than the cis. The extra hydrogen bonds enabled tighter packing of lipids in the all-trans and trans/cis bilayers, thus reducing the average area per lipid while increasing the chain order and the bilayer thickness. Moreover, fewer water molecules access the lipid-water interface of the all-trans bilayer than of the all-cis bilayer. These results provide the atomic basis for the importance of the natural sphingomyelin trans double-bond conformation for the formation of ordered membrane domains.

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.

Lipid-modulated sequence-specific association of glycophorin A in membranes.

Protein association in lipid membranes is a complex process with thermodynamics directed by a multitude of different factors. Amino-acid sequence is a molecular parameter that affects dimerization as shown by limited directed mutations along the transmembrane domains. Membrane-mediated interactions are also important although details of such contributions remain largely unclear. In this study, we probe directly the free energy of association of Glycophorin A by means of extensive parallel Monte Carlo simulations with recently developed methods and a model that accounts for sequence-specificity while representing lipid membranes faithfully. We find that lipid-induced interactions are significant both at short and intermediate separations. The ability of molecules to tilt in a specific hydrophobic environment extends their accessible interfaces, leading to intermittent contacts during protein recognition. The dimer with the lowest free energy is largely determined by the favorable lipid-induced attractive interactions at the closest distance. Finally, the coarse-grained model employed herein, together with the extensive sampling performed, provides estimates of the free energy of association that are in excellent agreement with existing data.

Simulating POPC and POPC/POPG Bilayers: Conserved Packing and Altered Surface Reactivity.

Molecular dynamics (MD) simulation is a popular technique to study bilayer structural properties, but it has not been widely used in mixed bilayers of neutral and charged lipids. Here, we present results from constant temperature and pressure MD simulations of a 2-oleoyl-1-pamlitoyl-sn-glyecro-3-phosphocholine (POPC) bilayer containing 23% 2-oleoyl-1-pamlitoyl-sn-glyecro-3-glycerol (POPG). The simulations were performed using the recently updated CHARMM force field and involved two bilayers of 104 and 416 lipids. A control simulation of a pure POPC bilayer of 128 lipids yielded equilibrium structural properties that compare very well with experimental data. The average equilibrium properties of the mixed bilayer systems were very similar to those of the pure POPC. However, nearly one-half of all the POPG lipids were found to be involved in hydrogen bonding with POPC lipids. Furthermore, the hydration of the mixed bilayer is different from that of the pure POPC, with the former inducing ordering of water molecules at longer distances. Thus, a phospholipid bilayer with ∼23% negative charge content in the liquid crystalline phase differs from its neutral counterpart only at the headgroup.

Accelerating flat-histogram methods for potential of mean force calculations.

Potential of mean force calculations along a reaction coordinate (RC) demand exhaustive sampling, which often leads to prohibitively long computational times. The expanded ensemble density of states (EXEDOS) [E. B. Kim, R. Faller, Q. Yan et al., J. Chem. Phys. 117, 7781 (2002)] is a simple flat-histogram Monte Carlo method based on the density of states algorithm proposed by Wang and Landau [Phys. Rev. Lett. 86, 2050 (2001)]. EXEDOS offers the advantage of continuous uniform sampling of the RC with no a priori knowledge of the free energy profile. However, the method is not certain to converge within accessible simulation time. Furthermore, the strongly asymmetric distribution of tunneling times inherent in flat-histogram sampling imposes additional limitations. We propose several improvements that accelerate the EXEDOS method and can be generally applicable in free energy calculations. First, we propose an asynchronous parallel implementation of the density of states algorithm in a multiple-walkers multiple-windows scheme and extend the algorithm in an expanded ensemble [(MW)(2)-XDOS] for PMF calculations as the original EXEDOS. Despite the nonideal scaling over a number of processors this technique overcomes limitations by extreme values of tunneling times and allows consistent evaluations of performance. The second set of improvements addresses the dependence of convergence times on system size, density, and sampling rate of the RC. At low densities, the coupling of (MW)(2)-XDOS with the rejection-free geometric cluster move provides impressive performance that overshadows any other technique. However, the limited applicability of cluster moves at high densities requires an alternative approach. We propose the coupling of (MW)(2)-XDOS with preferential sampling methods. In the systems studied, single displacements in the proximity of particles defining the RC accelerate calculations significantly and render the simulation nearly size-independent. A further modification of preferential sampling involves collective displacements of particles performed in a "smart Monte Carlo" scheme. This "local Brownian dynamics" algorithm can be generally applicable to many free energy simulation methods and would be particularly beneficial at high densities and molecular systems with strong intramolecular potentials.

Using stochastic models calibrated from nanosecond nonequilibrium simulations to approximate mesoscale information.

We demonstrate how the surrogate process approximation (SPA) method can be used to compute both the potential of mean force along a reaction coordinate and the associated diffusion coefficient using a relatively small number (10-20) of bidirectional nonequilibrium trajectories coming from a complex system. Our method provides confidence bands which take the variability of the initial configuration of the high-dimensional system, continuous nature of the work paths, and thermal fluctuations into account. Maximum-likelihood-type methods are used to estimate a stochastic differential equation (SDE) approximating the dynamics. For each observed time series, we estimate a new SDE resulting in a collection of SPA models. The physical significance of the collection of SPA models is discussed and methods for exploiting information in the population of estimated SPA models are demonstrated and suggested. Molecular dynamics simulations of potassium ion dynamics inside a gramicidin A channel are used to demonstrate the methodology, although SPA-type modeling has also proven useful in analyzing single-molecule experimental time series [J. Phys. Chem. B 113, 118 (2009)].

Calculating free-energy profiles in biomolecular systems from fast nonequilibrium processes.

Often gaining insight into the functioning of biomolecular systems requires to follow their dynamics along a microscopic reaction coordinate (RC) on a macroscopic time scale, which is beyond the reach of current all atom molecular dynamics (MD) simulations. A practical approach to this inherently multiscale problem is to model the system as a fictitious overdamped Brownian particle that diffuses along the RC in the presence of an effective potential of mean force (PMF) due to the rest of the system. By employing the recently proposed FR method [I. Kosztin, J. Chem. Phys. 124, 064106 (2006)], which requires only a small number of fast nonequilibrium MD simulations of the system in both forward and time reversed directions along the RC, we reconstruct the PMF: (1) of deca-alanine as a function of its end-to-end distance, and (2) that guides the motion of potassium ions through the gramicidin A channel. In both cases the computed PMFs are found to be in good agreement with previous results obtained by different methods. Our approach appears to be about one order of magnitude faster than the other PMF calculation methods and, in addition, it also provides the position-dependent diffusion coefficient along the RC. Thus, the obtained PMF and diffusion coefficient can be used in an overdamped Brownian model to estimate important characteristics of the studied systems, e.g., the mean folding time of the stretched deca-alanine and the mean diffusion time of the potassium ion through gramicidin A.

Theoretical prediction of spectral and optical properties of bacteriochlorophylls in thermally disordered LH2 antenna complexes.

A general approach for calculating spectral and optical properties of pigment-protein complexes of known atomic structure is presented. The method, that combines molecular dynamics simulations, quantum chemistry calculations, and statistical mechanical modeling, is demonstrated by calculating the absorption and circular dichroism spectra of the B800-B850 bacteriochlorophylls of the LH2 antenna complex from Rs. molischianum at room temperature. The calculated spectra are found to be in good agreement with the available experimental results. The calculations reveal that the broadening of the B800 band is mainly caused by the interactions with the polar protein environment, while the broadening of the B850 band is due to the excitonic interactions. Since it contains no fitting parameters, in principle, the proposed method can be used to predict optical spectra of arbitrary pigment-protein complexes of known structure.

Calculating potentials of mean force and diffusion coefficients from nonequilibrium processes without Jarzynski's equality.

In general, the direct application of the Jarzynski equality (JE) to reconstruct potentials of mean force (PMFs) from a small number of nonequilibrium unidirectional steered molecular-dynamics (SMD) paths is hindered by the lack of sampling of extremely rare paths with negative dissipative work. Such trajectories that transiently violate the second law of thermodynamics are crucial for the validity of JE. As a solution to this daunting problem, we propose a simple and efficient method, referred to as the FR method, for calculating simultaneously both the PMF U(z) and the corresponding diffusion coefficient D(z) along a reaction coordinate z for a classical many-particle system by employing a small number of fast SMD pullings in both forward (F) and time reverse (R) directions, without invoking JE. By employing Crooks [Phys. Rev. E 61, 2361 (2000)] transient fluctuation theorem (that is more general than JE) and the stiff-spring approximation, we show that (i) the mean dissipative work W(d) in the F and R pullings is the same, (ii) both U(z) and W(d) can be expressed in terms of the easily calculable mean work of the F and R processes, and (iii) D(z) can be expressed in terms of the slope of W(d). To test its viability, the FR method is applied to determine U(z) and D(z) of single-file water molecules in single-walled carbon nanotubes (SWNTs). The obtained U(z) is found to be in very good agreement with the results from other PMF calculation methods, e.g., umbrella sampling. Finally, U(z) and D(z) are used as input in a stochastic model, based on the Fokker-Planck equation, for describing water transport through SWNTs on a mesoscopic time scale that in general is inaccessible to MD simulations.