Team:UANL/Simulations

Modeling and Molecular Dynamics

Molecular Dynamics (MD) is a computational (numerical) method of simulation, in which the equations of motion are solved numerically to study the physical movement of atoms and molecules in a system. There are several programs that can perform Molecular Dynamics (MD) Simulations: MacroModel, LAMMPS, NAMD, etc. are some examples, but in this work GROMACS (gmx), a high performance MD program, was used to solve the equations of motion and Visual Molecular Dynamics (VMD) was used to construct the systems.

Modeling

As explained in Design, protein engineering was needed to modify one of the proteins of the Bacterial Microcompartment (BMC) to increase the permeability of substrate and enable it to pass through the proteinic shell. Specifically, an amino-acidic deletion in the n-terminal of the pduU protein (one of the structural proteins of pdu BMC) was desired.

Figure 1 - "u" subunit of the pdu BMC system (pduU). QuickSurf representation on top of NewCartoon from VMD tool. We can appreciate that the center of the protein is not a porous through which molecules can pas freely.


So, we took the FASTA format of the protein (Figure 2)...

Figure 2. FASTA: wild type pduU monomer. Underlined amino acids are the ones deleted.


... and deleted the first 17 amino-acids. We chose to delete those specifically because that modification had already been reported in an article by J. Jorda et. al. (3) for other purposes.
We then proceeded to evaluate it’s tertiary conformation with SwissModel, which uses homollogy modeling to predict the tertiary conformation of amino acidic sequences, and we obtained the predicted 3D conformation (1). It's easy to see that the porous is bigger now:). And thats VERY cool. Very nice.

Figure 3 - The "U" subunit of pdu BMC system with the amino acidic deletion (pduUmod). Surface representation on top of NewCartoon. From Visual Molecular Dynamics (VMD) (4). 3D conformation obtained from Swiss Model (1)


MD Simulatoins

MD simulations are divided in 2 parts:

  • Stability of modified protein verification
  • Permeability of new porous

Stability

We estimated that this 3D array would be stable because, as explained in Design, most of that remaining domain is a very conserved one, but we still wanted to prove it through Molecular Dynamics simulations, which we performed using GROMACS package (5) through

2 methods:

  • 1) Comparing Radius of Gyration (Rg) of the WT and modified proteins through ~90 ns simulations (Figure 4).
  • 2) Comparing Root-Mean-Square deviation (RMSD) of the WT and modified proteins through ~90 ns simulations (Figure 6).
  • 3) Comparing free energy of dimmer formation of both proteins using Umbrella Sampling (US) technique (2). See Figure 8.


Radius of Gyration

Figure 4. Change in Radius of Gyration between initial configuration of pduUWT (black) and of pduUmod (red) through +90,000ps.


Figure 5. Change in radius of Gyration between initial configuration of pduUWT (black) and of pduUmod (red) through 20,000ns.


It's easy to see that the average folding of both, wild type and modified proteins, is maintained throughout the whole simulation. The fact that the Rg is bigger for pduUmod means that the modification allowed the whole protein to increase in volume:

Root-Squared-Mean distance

Figure 6. Root-mean-squared distance between initial configuration of pduUWT (black) and of pduUmod (red) through +90,000ps.


Figure 7. Root-mean-squared distance between initial configuration of pduUWT (black) and of pduUmod (red) through 20,000ns.


We can observe that, despite the fact that pduUmod RMSD is a little bigger than pduUWT RMSD, both are still pretty stable during the simulation.

Free Energy: dimmer formation

Figure 8. Molecular Dynamics pulling simulation of one dimmer of pduU BMC subunit. From this simulations, snapshot are taken to evaluate interaction energy through the distance: Umbrella Sampling.


Permeability

We also needed to verify that the new porous allowed the substrate of the metabolic pathway to pass through the proteinic shell of the BMC. Which we also did using US technique.

Figure 9. Vanillate smoothly passing through pduUmod porus (front).
Figure 10. Vanillate smoothly passing through pduUmod porus (side).


References

¹ Waterhouse, A., Bertoni, M., Bienert, S., Studer, G., Tauriello, G., Gumienny, R., Heer, F.T., de Beer, T.A.P., Rempfer, C., Bordoli, L., Lepore, R., Schwede, T. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 46(W1), W296-W303 (2018).

² Hub JS, Groot BLD, Spoel DVD. g_wham—A Free Weighted Histogram Analysis Implementation Including Robust Error and Autocorrelation Estimates. Journal of Chemical Theory and Computation. 2010;6(12):3713–20.

³ Jorda, J., Liu, Y., Bobik, T. A., & Yeates, T. O. (2015). Exploring bacterial organelle interactomes: a model of the protein-protein interaction network in the Pdu microcompartment. PLoS computational biology, 11(2), e1004067. doi:10.1371/journal.pcbi.1004067

⁴ VMD

⁵ GROMACS