Periodic Reporting for period 1 - HYPERBIO (Hybrid Particle-Field Approach Including Electrostatics for Large-Scale Simulations of Biological Systems)
Reporting period: 2016-09-01 to 2018-08-31
1.The first objective of the proposal was to implement electrostatics into OCCAM software (WP1). The new code would be tested for its efficiency and performance in both serial and parallel implementations.
2.The second objective was the application of hPF-MD to realistic soft matter models (WP2). We aimed to validate the models for charged amphiphile systems and investigate the ability of hPF-MD in describing both the structural and dynamic properties of those charged systems. Our goal was to address the structure and assembly of complex ipopolysaccharide (LPS) membranes. In fact, we started our study with coarse grained models for simple charged amphiphile systems such as palmitoyloleoylphosphatidylglycerol (POPG) lipid bilayer and sodium dodecyl sulfate (SDS) surfactant in aqueous environment.
3.Our original plan for WP3 was to extend the charge-field formalism to incorporate dipole-field and dipole-dipole interactions for the description of proteins in a field approach compatible with CG method proposed by Michele Cascella and coworkers. On the contrary, we preferred to concentrate our efforts to introduce spatially resolved dielectric which depends on the local density of the different molecular species. This step is crucial for the better treatment of electrostatics in complex phase separated systems like biological membranes comprised of complex polyelectrolytes. In WP3, we also aimed at the introduction of protocols to perform simulations of tensionless bio membranes by considering an interaction energy due to surface tension in particle-field formalism.
Conclusion:
We introduced a series of methods in particle-field formalism to compute the electrostatic energy and forces for a mesoscopic system in the condensed phase, described with molecular resolution. The methods are inexpensive, robust and are able to reproduce the correct electrostatic features for dipolar/charged lipid membranes. A pressure tensor is derived in particle field formalism to perform simulations on tensionless bio membranes to study major conformational changes like membrane curvature, disruption and pore formation.
WP2 (01/05/2017 - 31/12/2017): By using the OCCAM code, we tested and developed models for POPG and SDS. Our work showed how our new implementation of grid based electrostatic solvers for particle field simulation can be used to investigate realistic models of polyelectrolyte soft matter systems. We demonstrated that the hPF-MD approach can treat amphiphilic systems either in extended aggregates, like lipid bilayers, or in finite micellar forms. . This work has been published in the Journal of Chemical Theory and Computation and was selected for the the cover page of the issue. As part of the continuation work on SDS self-assembly, we provided a theoretical explanation for the transition of spherical micelles to cylindrical aggregates (shown in Fig.1) in presence of ions. This manuscript is under preparation now.
WP3(01/01/2018 - 31/08/2018): A theoretical framework has been developed to compute the polarization forces through local dielectric permittivity. We do this by introducing a spatially resolved dielectric that depends on the local density of the different molecular species. We employed a finite difference scheme to solve generalized Poisson equation iteratively using the Successive over relaxation method (SOR). This has been implemented into OCCAM and tested for ion permeability in oil/water system and lipid membranes. This work is crucial for the better treatment of electrostatics in complex phase separated systems. The manuscript on this work is finished and about to be submitted. Alongside, the derivation of the pressure tensor for amphiphile mono/bi layers in hPF formalism has been performed and implemented in OCCAM.
Overview:
A series of methods are developed for rigorous and accurate treatment of electrostatics. This work opens the possibility of using hPF schemes to investigate major biological processes that are dominated by such interactions-for example ion-membrane permeation, membrane electroporation, protein/protein and protein/membrane interactions. This project has delivered two publications in reputed journals, one manuscript to be submitted and two are under preparation.