Electrostatic-Consistent Coarse-Grained Potentials for Molecular Simulations of Proteins

We present a new generation of coarse-grained (CG) potentials that account for a simplified electrostatic description of soluble proteins. The treatment of permanent electrostatic dipoles of the backbone and polar side-chains allows to simulate proteins, preserving an excellent structural and dynamic agreement with respective reference structures and all-atom molecular dynamics simulations. Moreover, multiprotein complexes can be well described maintaining their molecular interfaces thanks to the ability of this scheme to better describe the actual electrostatics at a CG level of resolution. An efficient and robust heuristic algorithm based on particle swarm optimization is used for the derivation of CG parameters via a force-matching procedure. The ability of this protocol to deal with high dimensional search spaces suggests that the extension of this optimization procedure to larger data sets may lead to the generation of a fully transferable CG force field. At the present stage, these electrostatic-consistent CG potentials are easily and efficiently parametrized, show a good degree of transferability, and can be used to simulate soluble proteins or, more interestingly, large macromolecular assemblies for which long all-atom simulations may not be easily affordable.

Multiscale simulations of protein dynamics

Enrico Spiga

Characteristic timescales associated with the function of biomolecules, like proteins, range from femtoseconds up to minutes, whereas their corresponding spatial extent ranges from few ̊A to μm when associating in large macromolecular complexes. Moreover, biomolecules are functional in a large variety of different physico-chemical conditions strongly dependent on pH, ionic strength, crowding agents, etc. This huge complexity is hard to be studied with an arbitrary level of resolution embracing all these spatial and temporal scales. Molecular simulation is a well established approach to gain mechanistic insights into the function of biomolecular systems, producing atomistically detailed models of in vitro and/or in vivo conditions. I present in this thesis two projects that aim at improving on the current limitations of multiscale molecular simulations, namely (i) the sampling of large systems, and (ii) the detailed representation and description of realistic physiological conditions. Addressing the first issue, I propose a new coarse-grained model for proteins to be used in molecular dynamics simulations. This coarse-grained model is based on a more accurate description of protein electrostatics, which accounts for dipolar contributions. The parameterization of this force field is based on force-matching methods and on the use of a particle swarm optimization heuristic algorithm. The obtained results are encouraging being structural and electrostatic properties accurately reproduced with the coarse-grained model for a variety of protein folds. Moreover, the parameterization procedure can be straightforwardly applied to any protein, and can be extended to a larger dataset to generate a fully transferable coarse-grained force field, to be applied to any protein and any large macromolecular assemblies for which long all-atom simulations are still a challenge. While the development and use of coarse-grained models are important to tackle the limited sampling of large systems, it is still important to use all-atom molecular dynamics simulation to investigate with high accuracy in physiological conditions protein dynamics. For this reason, I present here the results of state-of-art molecular dynamics simulations applied to study the influence of crowding agents on the internal dynamics of the protein ubiquitin. Their analysis allows to describe how ubiquitin dynamics is slaved by crowding agents. This work demonstrates that the description of protein dynamics should take into account its intrinsic multiscale nature. The development and applications of coarse-grained models permit to simulate proteins at low computational cost, the use of atomistic simulations allows to accurately describe proteins in absence and presence of crowding agents, and both of them permit to highlight the essence of protein dynamics.

EPFL2013

Molecular Modeling of Bacterial Nanomachineries

Matteo Thomas Degiacomi

Proteins have the ability to assemble in multimeric states to perform their specific biological function. Unfortunately, characterizing experimentally these structures at atomistic resolution is usually difficult. For this reason, in silico methodologies aiming at predicting how multiple protein copies arrange to forma multimeric complex would be desirable. We present Parallel OptimizationWorkbench (POW), a swarm intelligence based optimization framework able to deal, in principle, with any optimization problem. We show that POW can be applied to biologically relevant problems such as prediction of protein assemblies and the parameterization of a Coarse-Grained force field for proteins. By combining POW optimizations, Molecular Dynamics simulations, Poisson-Boltzmann calculations and a variety of experiments, we subsequently study two bacterial nanomachieries: Aeromonas hydrophila's pore-forming toxin aerolysin, and Yersinia enterocolitica injectisome. These structures are challenging both for their size, and for the timescales involved in their functioning. Aerolysin is a pore-forming toxin secreted as an hydrophilic monomer. By means of large conformational changes, the protein heptamerizes on the target cell's surface, and finally inserts β-barrel into its lipid bilayer, causing cell death. The main hurdle in the study of this structure is the complexity of the mode of action, which spans timescales currently unreachable by classical molecular dynamics. We show that aerolysin C-terminal region has the dual role of preventing premature oligomerization and helping the folding of tertiary structure, qualifying therefore as an intramolecular chaperone. We study the transmembrane β-barrel properties and compare them with those of the homologous protein α-hemolysin. We show that aerolysin's barrel is more rigid than α-hemolysin's, and should be anion selective. We present models for aerolysin heptamer both in prepore and, for the first time, in membrane-inserted conformation. Our results are validated experimentally, and are consistent with known biochemical and structural data. The injectisome is an example of a type III secretion system. Its most striking feature is probably its size: hundreds of proteins assemble in a unique structure spanning the Gram-negative bacterial double membrane, and protruding outside the cell as a needle for tenth of nanometers. Obtaining an atomistic representation of this massive structure, and therefore some insights about its mode of action, is one of the greatest challenges. We show that the final length of injectisome's needle is determined by the secondary structure content of a ruler protein located inside its cavity during assembly. Using POW, we also produce the first model for Yersinia injectisome's basal body, highlighting the flexibility of this region in adapting between the inner and outer membranes. As a whole, this work demonstrates that a synergy of dry and wet experiments can provide precious insights into macromolecular structure and function.

EPFL2012

Multiscale simulations of protein dynamics

Enrico Spiga

EPFL2013

Molecular Modeling of Bacterial Nanomachineries

Matteo Thomas Degiacomi

EPFL2012