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.
View Article and Find Full Text PDFIn 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.
View Article and Find Full Text PDFBiochim Biophys Acta Gen Subj
July 2017
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.
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.
View Article and Find Full Text PDFExperiments 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.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
March 2012
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.
View Article and Find Full Text PDFProc Natl Acad Sci U S A
May 2012
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.
View Article and Find Full Text PDFGxxxG 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.
View Article and Find Full Text PDFMost 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.
View Article and Find Full Text PDFThe 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.
View Article and Find Full Text PDFAssociation 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.
View Article and Find Full Text PDFThe 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.
View Article and Find Full Text PDFJ Chem Theory Comput
October 2010
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.
View Article and Find Full Text PDFProtein 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.
View Article and Find Full Text PDFPotential 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.
View Article and Find Full Text PDFWe 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.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
November 2008
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.
View Article and Find Full Text PDFA 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.
View Article and Find Full Text PDFIn 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.
View Article and Find Full Text PDF