Difference between revisions of "Team:Paris Bettencourt/Modeling"

 
(4 intermediate revisions by 2 users not shown)
Line 5: Line 5:
 
     <img src="https://static.igem.org/mediawiki/2018/d/d6/T--Paris_Bettencourt--Model_Ban.png">
 
     <img src="https://static.igem.org/mediawiki/2018/d/d6/T--Paris_Bettencourt--Model_Ban.png">
 
</div>
 
</div>
 +
 
<div class='textbody h1'>
 
<div class='textbody h1'>
 
<h1>Modeling</h1>
 
<h1>Modeling</h1>
 
</div>
 
</div>
  
−
<div class='textbody h2'
+
<div class='textbody h2'>
 
<h2>INTRODUCTION</h2>
 
<h2>INTRODUCTION</h2>
 
</div>
 
</div>
  
−
<div class='textbody'>
+
<div class='textbody h3'>
−
<p>Our project is based on the consequences of conformational change of antimicrobial peptides. Results generated by the testing group showed that MIC is not a reliable criteria to understand the activity of our StarCores while it has been  
+
<h3>Our project is based on the consequences of conformational change of antimicrobial peptides. Results generated by the testing group showed that MIC is not a reliable criteria to understand the activity of our StarCores while it has been  
−
previously used for species. It was crucial to have some models to:</p>
+
previously used for species. It was crucial to have some models to:</h3>
 +
</div>
 +
 
 
</div>
 
</div>
 
<table>
 
<table>
Line 28: Line 31:
 
</tr>
 
</tr>
 
</tbody>
 
</tbody>
−
</table><hr><h2 id="Results">Results</h2><h3 id="Our-Modelling-workflow-could-be-summed-up-as-">Our Modelling workflow could be summed up as :</h3><ul>
+
 
−
<li>Two hackathons to screen the most promising cores and AMPs based on rationnal criteria</li>
+
 
−
<li>Homology modeling of the 210 constructs with Modeller on Pymol</li>
+
<div class='textbody h2'>
−
<li>Molecular Dynamic simulations to test the stability of the predicted structures using Yasara and Inserm Cluster</li>
+
<h2>RESULTS</h2>
−
</ul><h3 id="I-Homology-Modelling">I. Homology Modelling</h3><p>Followed by cores and AMPs selection, we were interested in determing the model of the fusion. This led us to the modeling expert, Antoine Tally, who guided us to standardise a protocol for our analysis.</p><ul>
+
</div>
−
<li>Aim: To construct an atomic-resolution model of Star-core monomer via comparative homology modelling</li>
+
 
−
</ul><h4 id="Making-fusion-protein-of-2-known-PDB-files-using-CHIMERA-amp-MODELLER">Making fusion protein of 2 known PDB files using CHIMERA &amp; MODELLER</h4><ul>
+
<div class="textbody h4">
−
<li>Monomer construct on Chimera using Modeller graphic interface<br>
+
<p>Our Modelling workflow could be summed up as :</p>
−
<img src="https://static.igem.org/mediawiki/2018/8/87/T--Paris_Bettencourt--modelling-construct.png" alt="image alt"></li>
+
</div>
−
<li>construct optimization on Chimera<br>
+
 
−
<img src="https://static.igem.org/mediawiki/2018/7/7f/T--Paris_Bettencourt--modelling-optimization.png" alt="image alt"></li>
+
 
−
<li>Superpozition on Pymol thanks to handmade script called <a href="http://Superpoz.py" target="_blank">Superpoz.py</a><br>
+
<div class="textbody">
−
<img src="https://static.igem.org/mediawiki/2018/b/b1/T--Paris_Bettencourt--modelling-superpoz.png" alt="image alt"></li>
+
<p>Two hackathons to screen the most promising cores and AMPs based on rationnal criteria</p>
−
</ul><h3 id="II-Molecular-dynamic-simulation">II. Molecular dynamic simulation</h3><h4 id="Objective-">Objective :</h4><h5 id="To-study-behaviour-of-fusion-molecule-in-vicinity-of-cell">To study behaviour of fusion molecule in vicinity of cell.</h5><p>
+
</div>
−
<em>Aim:</em></p>
+
 
−
<p>To assess the stability of Star cores in terms of protein folding, structural characteristics and energy.</p><p>Objectives and steps:</p><ul>
+
<div class="textbody">
−
<li>Define constraints and parameters for simulations</li>
+
<p>Homology modeling of the 210 constructs with Modeller on Pymol</p>
−
<li>Using Yasara to do MDS</li>
+
</div>
−
<li>Using Yasara to analyze MDS data</li>
+
 
−
</ul><p>
+
<div class="textbody">
−
<img src="https://static.igem.org/mediawiki/2018/4/4f/T--Paris_Bettencourt--modelling-cube.png" alt="image alt"></p>
+
<p>Molecular Dynamic simulations to test the stability of the predicted structures using Yasara and Inserm Cluster</p>
−
<hr><h2 id="Conclusion">Conclusion</h2><p>
+
</div>
−
In this project, we combinatorially fused a set of known AMPs to structurally diverse, self-assembling protein cores to produce star-shaped complexes.</p><p>We selected 14 cores that already exist and 15 AMPs, hence over 200 fusions were designed and expressed in a cell-free system, then screened for activity, biocompatibility, and membrane selectivity.</p><p>To study the behaviour of in-vitro synthesised molecules we designed the fusion molecules for all constructs using Chimera, Modeller and Pymol. We visualized the assembly and fusion monomer (core + AMP) and studied their behavioural attributes, changes in folding of alpha helix and beta sheath using Yasara as it mimics atom behaviour in real life. Using visualization and MDS studies we confirm that proteins folds well and on an average maintain constant RMSF (Root Mean Square Fluctuations) with all amino acid residues, also expressed well in cell free system.</p>
+
 
−
<hr><hr><hr><h2 id="Methods">Methods</h2><p><strong>These protocols were defined under certain parameters:</strong></p>
+
 
 +
 
 +
<div class='textbody h2'>
 +
<h2>I. Homology Modelling</h2>
 +
</div>
 +
 
 +
<div class="textbody">
 +
<p>Followed by cores and AMPs selection, we were interested in determing the model of the fusion. This led us to the modeling expert, Antoine Tally, who guided us to standardise a protocol for our analysis. Aim: To construct an atomic-resolution model of Star-core monomer via comparative homology modelling</p>
 +
</div>
 +
 
 +
<div class="textbody">
 +
<p>Making fusion protein of 2 known PDB files using CHIMERA</p>
 +
</div>
 +
 
 +
 
 +
<div class='textbody h4'>
 +
<h4>Monomer construct on Chimera using Modeller graphic interface</h4>
 +
</div>
 +
 
 +
 
 +
<div class="text3 img">
 +
<img src="https://static.igem.org/mediawiki/2018/8/87/T--Paris_Bettencourt--modelling-construct.png">
 +
</div>
 +
 
 +
<div class='textbody h4'>
 +
<h4>construct optimization on Chimera</h4>
 +
</div>
 +
 
 +
 
 +
<div class="text3 img">
 +
<img src="https://static.igem.org/mediawiki/2018/7/7f/T--Paris_Bettencourt--modelling-optimization.png">
 +
</div>
 +
 
 +
<div class='textbody h4'>
 +
<h4>Superpozition on Pymol thanks to handmade script called <a href="http://Superpoz.py" target="_blank">Superpoz.py</a><br></h4>
 +
</div>
 +
 
 +
 
 +
<div class="text3 img">
 +
<img src="https://static.igem.org/mediawiki/2018/b/b1/T--Paris_Bettencourt--modelling-superpoz.png">
 +
</div>
 +
 
 +
 
 +
<br/><br/>
 +
 
 +
<div class='textbody h2'>
 +
<h2>II. Molecular dynamic simulation</h2>
 +
</div>
 +
 
 +
 
 +
<div class="textbody">
 +
<p>OBJECTIVE: </p>
 +
</div>
 +
<div class="textbody">
 +
<p>-Study behaviour of fusion molecule in vicinity of cell.</p>
 +
</div>
 +
<div class="textbody">
 +
<p>-Assess the stability of Star cores in terms of protein folding, structural characteristics and energy.</p>
 +
</div>
 +
<div class="textbody">
 +
<p>-Define constraints and parameters for simulations.</p>
 +
</div>
 +
<div class="textbody">
 +
<p>-Use Yasara to analyze MDS data.</p>
 +
</div>
 +
 
 +
 
 +
 
 +
 
 +
 
 +
<div class="text3 img">
 +
<img src="https://static.igem.org/mediawiki/2018/4/4f/T--Paris_Bettencourt--modelling-cube.png">
 +
</div>
 +
 
 +
 
 +
<div class='textbody h2'>
 +
<h2>CONCLUSION</h2>
 +
</div>
 +
 
 +
<div class="textbody">
 +
<p>In this project, we combinatorially fused a set of known AMPs to structurally diverse, self-assembling protein cores to produce star-shaped complexes.</p><p>We selected 14 cores that already exist and 15 AMPs, hence over 200 fusions were designed and expressed in a cell-free system, then screened for activity, biocompatibility, and membrane selectivity.</p><p>To study the behaviour of in-vitro synthesised molecules we designed the fusion molecules for all constructs using Chimera, Modeller and Pymol. We visualized the assembly and fusion monomer (core + AMP) and studied their behavioural attributes, changes in folding of alpha helix and beta sheath using Yasara as it mimics atom behaviour in real life. Using visualization and MDS studies we confirm that proteins folds well and on an average maintain constant RMSF (Root Mean Square Fluctuations) with all amino acid residues, also expressed well in cell free system.</p>
 +
</div>
 +
 
 +
 
 +
<div class='textbody h2'>
 +
<h2>METHODS</h2>
 +
</div>
 +
 
 +
<div class='textbody h4'>
 +
<h4>These protocols were defined under certain parameters:</h4>
 +
</div>
 +
 
 +
 
 +
 
 
<ul>
 
<ul>
 
<li>
 
<li>
Line 62: Line 158:
 
<p>However we verified this hypothesis with our stability MD assay developped with Marc Baaden.</p>
 
<p>However we verified this hypothesis with our stability MD assay developped with Marc Baaden.</p>
 
</li>
 
</li>
−
</ul><h3 id="Homology-Modelling">Homology Modelling</h3><h4 id="What-is-homology-modelling">What is homology modelling?</h4><blockquote>
+
 
−
<p>Homology modelling is comparative modelling of proteins. It is a comparative protein modelling method designed to find the most probable structure for a sequence given its alignment with related structures. The three-dimensional (3D) model is obtained by optimally satisfying spatial restraints derived from the alignment and expressed as probability density functions (pdfs) for the features restrained. For example, the probabilities for main-chain conformations of a modelled residue may be restrained by its residue type, main-chain conformation of an equivalent residue in a related protein, and the local similarity between the two sequences. Several such pdfs are obtained from the correlations between structural features in 17 families of homologous proteins which have been aligned on the basis of their 3D structures. The pdfs restrain C alpha-C alpha distances, main-chain N-O distances, main-chain and side-chain dihedral angles. A smoothing procedure is used in the derivation of these relationships to minimize the problem of a sparse database. The 3D model of a protein is obtained by optimization of the molecular pdf such that the model violates the input restraints as little as possible. The molecular pdf is derived as a combination of pdfs restraining individual spatial features of the whole molecule. The optimization procedure is a variable target function method that applies the conjugate gradients algorithm topositions of all non-hydrogen atoms. The method is automated. We used modeller to predict all our models.</p>
+
 
−
</blockquote><p><em>Objectives and steps:</em></p>
+
 
 +
<div class='textbody h4'>
 +
<h4>Homology Modelling</h4>
 +
</div>
 +
 
 +
<div class='textbody'>
 +
<p>What is homology modelling?</p>
 +
</div>
 +
 
 +
 
 +
<p>Homology modelling is comparative modelling of proteins. It is a comparative protein modelling method designed to find the most probable structure for a sequence given its alignment with related structures. The three-dimensional (3D) model is obtained by optimally satisfying spatial restraints derived from the alignment and expressed as probability density functions (pdfs) for the features restrained. For example, the probabilities for main-chain conformations of a modelled residue may be restrained by its residue type, main-chain conformation of an equivalent residue in a related protein, and the local similarity between the two sequences. Several such pdfs are obtained from the correlations between structural features in 17 families of homologous proteins which have been aligned on the basis of their 3D structures. The pdfs restrain C alpha-C alpha distances, main-chain N-O distances, main-chain and side-chain dihedral angles. A smoothing procedure is used in the derivation of these relationships to minimize the problem of a sparse database. The 3D model of a protein is obtained by optimization of the molecular pdf such that the model violates the input restraints as little as possible. The molecular pdf is derived as a combination of pdfs restraining individual spatial features of the whole molecule. The optimization procedure is a variable target function method that applies the conjugate gradients algorithm to positions of all non-hydrogen atoms. The method is automated. We used modeller to predict all our models. The steps are: </p>
 +
 
 +
 
 
<table>
 
<table>
 
<thead>
 
<thead>
Line 74: Line 182:
 
</thead>
 
</thead>
 
<tbody></tbody>
 
<tbody></tbody>
−
</table><p><em>Steps</em><br>
+
 
−
For this workflow we will use:</p><ul>
+
 
 +
 
 +
 
 +
<p>For this workflow we will use:</p><ul>
 
<li>UCSF-CHIMERA : download link</li>
 
<li>UCSF-CHIMERA : download link</li>
 
<li>MODELLER: online or installed</li>
 
<li>MODELLER: online or installed</li>
Line 99: Line 210:
 
<p>Position the two structures so that the termini are in a somewhat reasonable place relative to each other to template the fusion protein. In our case, the C-terminal of Defensin monomer attaches with the N-terminal of Ferritin monomer. You can “freeze” one in place by deactivating it and move just the other with the mouse as described here:</p>
 
<p>Position the two structures so that the termini are in a somewhat reasonable place relative to each other to template the fusion protein. In our case, the C-terminal of Defensin monomer attaches with the N-terminal of Ferritin monomer. You can “freeze” one in place by deactivating it and move just the other with the mouse as described here:</p>
 
</li>
 
</li>
−
</ul><p><img src="https://static.igem.org/mediawiki/2018/8/87/T--Paris_Bettencourt--modelling-construct.png" alt="image alt"></p><ul>
+
 
 +
<div class="text3 img">
 +
<img src="https://static.igem.org/mediawiki/2018/8/87/T--Paris_Bettencourt--modelling-construct.png">
 +
</div>
 +
 
 +
 
 
<li>From the sequence alignment window menu choose:  Structure… Modeller (homology) to show the Modeller dialog.  Choose the query as the target and both structures as the template, etc. as in the modeling tutorials.  You may also want to turn on “Use thorough optimization” in the Advanced Options section.</li>
 
<li>From the sequence alignment window menu choose:  Structure… Modeller (homology) to show the Modeller dialog.  Choose the query as the target and both structures as the template, etc. as in the modeling tutorials.  You may also want to turn on “Use thorough optimization” in the Advanced Options section.</li>
−
</ul><p><img src="https://static.igem.org/mediawiki/2018/7/7f/T--Paris_Bettencourt--modelling-optimization.png" alt="image alt"></p><p>NOW YOU CAN SAVE ONE OF THE STRUCTURE AS A PDB FILE AND USE THIS AS A MONOMER FOR FURTHER STEPS.</p><p>Part 2: Make a multimer of the fusion protein using reference PDB structure in Pymol</p><p>Open core assembly ( 4XGS) and fusion monomer (4XGS + 2B68) in pymol.<br>
+
 
 +
<div class="text3 img">
 +
        <img src="https://static.igem.org/mediawiki/2018/7/7f/T--Paris_Bettencourt--modelling-optimization.png">
 +
        </div>
 +
 
 +
 
 +
 
 +
<p>NOW YOU CAN SAVE ONE OF THE STRUCTURE AS A PDB FILE AND USE THIS AS A MONOMER FOR FURTHER STEPS.</p><p>Part 2: Make a multimer of the fusion protein using reference PDB structure in Pymol</p><p>Open core assembly ( 4XGS) and fusion monomer (4XGS + 2B68) in pymol.<br>
 
Use python based code, <a href="http://superpoz.py" target="_blank">superpoz.py</a>. We wrote a small script to generate the assembly structure of our StarCores based on the reference biological assembly of the nude core.</p><h4 id="Result">Result</h4><p>Fusion monomers for all the constructs are developed for molecular dynamic simulation studies and assembly of core with monomer for visualization of scaffold proteins.</p><h2 id="Molecular-dynamic-simulation">Molecular dynamic simulation</h2><h4 id="What-is-Molecular-Dynamic-Simulations-">What is Molecular Dynamic Simulations ?</h4><blockquote>
 
Use python based code, <a href="http://superpoz.py" target="_blank">superpoz.py</a>. We wrote a small script to generate the assembly structure of our StarCores based on the reference biological assembly of the nude core.</p><h4 id="Result">Result</h4><p>Fusion monomers for all the constructs are developed for molecular dynamic simulation studies and assembly of core with monomer for visualization of scaffold proteins.</p><h2 id="Molecular-dynamic-simulation">Molecular dynamic simulation</h2><h4 id="What-is-Molecular-Dynamic-Simulations-">What is Molecular Dynamic Simulations ?</h4><blockquote>
 
<p>MD studies are based on Newton second law of motion :<br>
 
<p>MD studies are based on Newton second law of motion :<br>
Line 128: Line 251:
 
<li>Analyse the result using MD_analyze.mcr</li>
 
<li>Analyse the result using MD_analyze.mcr</li>
 
<li>Analyse each residue  using MD_analyzeres.mcr</li>
 
<li>Analyse each residue  using MD_analyzeres.mcr</li>
−
</ul><h4 id="Results1">Results:</h4><p><img src="https://static.igem.org/mediawiki/2018/4/4f/T--Paris_Bettencourt--modelling-cube.png" alt="image alt"></p><pre><code>       Figure :  Simulation of fusion monomer ferritin + ovispirin (4XGS + 1HU5)
+
 
 +
<div class="text3 img">
 +
        <img src="https://static.igem.org/mediawiki/2018/4/4f/T--Paris_Bettencourt--modelling-cube.png">
 +
        </div>
 +
 
 +
<pre><code>  
 +
Figure :  Simulation of fusion monomer ferritin + ovispirin (4XGS + 1HU5)
 
</code></pre><h4 id="Parameters-">Parameters :</h4><ul>
 
</code></pre><h4 id="Parameters-">Parameters :</h4><ul>
−
<li>pH at which the simulation should be run, by default physiological pH 7.4<br>
+
<li> -pH at which the simulation should be run, by default physiological pH 7.4<br>
 
The ion concentration as a mass fraction, here we use 0.9% NaCl (physiological solution)ions=‘NaCl,0.9’</li>
 
The ion concentration as a mass fraction, here we use 0.9% NaCl (physiological solution)ions=‘NaCl,0.9’</li>
−
<li>Simulation temperature, 298K</li>
+
<li> -Simulation temperature, 298K</li>
−
<li>Water density = 0.997</li>
+
<li> -Water density = 0.997</li>
−
<li>Duration of the simulation = 15ns</li>
+
<li> -Duration of the simulation = 15ns</li>
−
<li>Extension of the cell on each side of the protein ‘10’ means that the cell will be 20 A larger than the protein.<br>
+
<li> -Extension of the cell on each side of the protein ‘10’ means that the cell will be 20 A larger than the protein.<br>
 
Shape of the simulation cell: ‘Cube’.</li>
 
Shape of the simulation cell: ‘Cube’.</li>
−
<li>Forcefield : ForceField AMBER14</li>
+
<li> -Forcefield : ForceField AMBER14</li>
−
<li>Cell boundary : Boundary periodic<br>
+
<li> -Cell boundary : Boundary periodic<br>
 
The simulation speed, ‘fast’ (maximize performance with 2*2.5 fs timestep and constraints)<br>
 
The simulation speed, ‘fast’ (maximize performance with 2*2.5 fs timestep and constraints)<br>
 
The save interval for snapshots. Normally you don’t need more than 500-1000 snapshots of your simulation.<br>
 
The save interval for snapshots. Normally you don’t need more than 500-1000 snapshots of your simulation.<br>
 
Solute from diffusing around and crossing periodic boundaries. Disable that for simulations of crystals.</li>
 
Solute from diffusing around and crossing periodic boundaries. Disable that for simulations of crystals.</li>
−
</ul><p>*Results:$</p><p><img src="https://static.igem.org/mediawiki/2018/b/b5/T--Paris_Bettencourt--modelling-pizza1.png" alt="image alt"><br>
+
 
−
<img src="https://static.igem.org/mediawiki/2018/3/3a/T--Paris_Bettencourt--modelling-core6.png" alt="image alt"></p><p><img src="https://static.igem.org/mediawiki/2018/a/ac/T--Paris_Bettencourt--modelling-core8.png" alt="image alt"><br>
+
<div class="text3 img">
−
<img src="https://static.igem.org/mediawiki/2018/8/82/T--Paris_Bettencourt--modelling-ferritin.png" alt="image alt"><br>
+
        <img src="https://static.igem.org/mediawiki/2018/b/b5/T--Paris_Bettencourt--modelling-pizza1.png">
−
<img src="https://static.igem.org/mediawiki/2018/a/a1/T--Paris_Bettencourt--modelling-pizza2.png" alt="image alt"></p><div dir="ltr" class="resize-sensor" style="position: absolute; left: -10px; top: -10px; right: 0px; bottom: 0px; overflow: hidden; z-index: -1; visibility: hidden;"><div class="resize-sensor-expand" style="position: absolute; left: -10px; top: -10px; right: 0; bottom: 0; overflow: hidden; z-index: -1; visibility: hidden;"><div style="position: absolute; left: 0px; top: 0px; transition: all 0s ease 0s; width: 100000px; height: 100000px;"></div></div><div class="resize-sensor-shrink" style="position: absolute; left: -10px; top: -10px; right: 0; bottom: 0; overflow: hidden; z-index: -1; visibility: hidden;"><div style="position: absolute; left: 0; top: 0; transition: 0s; width: 200%; height: 200%"></div></div></div></div>
+
        </div>
 +
 
 +
 
 +
 
 +
<div class="text3 img">
 +
        <img src="https://static.igem.org/mediawiki/2018/3/3a/T--Paris_Bettencourt--modelling-core6.png">
 +
        </div>
 +
 
 +
 
 +
 
 +
 
 +
<div class="text3 img">
 +
        <img src="https://static.igem.org/mediawiki/2018/a/ac/T--Paris_Bettencourt--modelling-core8.png">
 +
        </div>
 +
 
 +
<div class="text3 img">
 +
        <img src="https://static.igem.org/mediawiki/2018/8/82/T--Paris_Bettencourt--modelling-ferritin.png">
 +
        </div>
 +
 
 +
<div class="text3 img">
 +
        <img src="https://static.igem.org/mediawiki/2018/8/82/T--Paris_Bettencourt--modelling-ferritin.png">
 +
        </div>
 +
 
 +
<div class="text3 img">
 +
        <img src="https://static.igem.org/mediawiki/2018/a/a1/T--Paris_Bettencourt--modelling-pizza2.png">
 +
        </div>
 +
 
 +
 
 +
 
 +
<div dir="ltr" class="resize-sensor" style="position: absolute; left: -10px; top: -10px; right: 0px; bottom: 0px; overflow: hidden; z-index: -1; visibility: hidden;"><div class="resize-sensor-expand" style="position: absolute; left: -10px; top: -10px; right: 0; bottom: 0; overflow: hidden; z-index: -1; visibility: hidden;"><div style="position: absolute; left: 0px; top: 0px; transition: all 0s ease 0s; width: 100000px; height: 100000px;"></div></div><div class="resize-sensor-shrink" style="position: absolute; left: -10px; top: -10px; right: 0; bottom: 0; overflow: hidden; z-index: -1; visibility: hidden;"><div style="position: absolute; left: 0; top: 0; transition: 0s; width: 200%; height: 200%"></div></div></div></div>
  
  

Latest revision as of 03:30, 18 October 2018

Modeling

INTRODUCTION

Our project is based on the consequences of conformational change of antimicrobial peptides. Results generated by the testing group showed that MIC is not a reliable criteria to understand the activity of our StarCores while it has been previously used for species. It was crucial to have some models to:

RESULTS

Our Modelling workflow could be summed up as :

Two hackathons to screen the most promising cores and AMPs based on rationnal criteria

Homology modeling of the 210 constructs with Modeller on Pymol

Molecular Dynamic simulations to test the stability of the predicted structures using Yasara and Inserm Cluster

I. Homology Modelling

Followed by cores and AMPs selection, we were interested in determing the model of the fusion. This led us to the modeling expert, Antoine Tally, who guided us to standardise a protocol for our analysis. Aim: To construct an atomic-resolution model of Star-core monomer via comparative homology modelling

Making fusion protein of 2 known PDB files using CHIMERA

Monomer construct on Chimera using Modeller graphic interface

construct optimization on Chimera

Superpozition on Pymol thanks to handmade script called Superpoz.py



II. Molecular dynamic simulation

OBJECTIVE:

-Study behaviour of fusion molecule in vicinity of cell.

-Assess the stability of Star cores in terms of protein folding, structural characteristics and energy.

-Define constraints and parameters for simulations.

-Use Yasara to analyze MDS data.

CONCLUSION

In this project, we combinatorially fused a set of known AMPs to structurally diverse, self-assembling protein cores to produce star-shaped complexes.

We selected 14 cores that already exist and 15 AMPs, hence over 200 fusions were designed and expressed in a cell-free system, then screened for activity, biocompatibility, and membrane selectivity.

To study the behaviour of in-vitro synthesised molecules we designed the fusion molecules for all constructs using Chimera, Modeller and Pymol. We visualized the assembly and fusion monomer (core + AMP) and studied their behavioural attributes, changes in folding of alpha helix and beta sheath using Yasara as it mimics atom behaviour in real life. Using visualization and MDS studies we confirm that proteins folds well and on an average maintain constant RMSF (Root Mean Square Fluctuations) with all amino acid residues, also expressed well in cell free system.

METHODS

These protocols were defined under certain parameters:

  • AMPs are relatively small, so we assumed that structural changes would be minimal, enabling us to perform Homology modelling.

  • AMPs are only fused on the N/C terminals thus oriented outside the homomultimeric self-assembling protein, the modelisation should keep the geometry of the nude core, displaying the fusionned peptides on the surface of the core. It enabled us to structure superposition.

  • However we verified this hypothesis with our stability MD assay developped with Marc Baaden.

    • Homology Modelling

    What is homology modelling?

    Homology modelling is comparative modelling of proteins. It is a comparative protein modelling method designed to find the most probable structure for a sequence given its alignment with related structures. The three-dimensional (3D) model is obtained by optimally satisfying spatial restraints derived from the alignment and expressed as probability density functions (pdfs) for the features restrained. For example, the probabilities for main-chain conformations of a modelled residue may be restrained by its residue type, main-chain conformation of an equivalent residue in a related protein, and the local similarity between the two sequences. Several such pdfs are obtained from the correlations between structural features in 17 families of homologous proteins which have been aligned on the basis of their 3D structures. The pdfs restrain C alpha-C alpha distances, main-chain N-O distances, main-chain and side-chain dihedral angles. A smoothing procedure is used in the derivation of these relationships to minimize the problem of a sparse database. The 3D model of a protein is obtained by optimization of the molecular pdf such that the model violates the input restraints as little as possible. The molecular pdf is derived as a combination of pdfs restraining individual spatial features of the whole molecule. The optimization procedure is a variable target function method that applies the conjugate gradients algorithm to positions of all non-hydrogen atoms. The method is automated. We used modeller to predict all our models. The steps are:

Determine which constructs would be interesting
Interpret our experiments results

For this workflow we will use:

  • UCSF-CHIMERA : download link
  • MODELLER: online or installed
  • Core: Ferritin: PDB= 4XGS
  • AMP: Cg-Defensin: PDB= 2B68
  • AA seq in FASTA format: download

pDB.001_translation
MGFGCPGNQLKCNNHCKSISCRAGYCDAATLWLRCTCTDCNGKKESSHLKPEMIEKLNEQMNLELYSSLLYQQMSAWCSYHTFEGAAAFLRRHAQEEMTHMQRLFDYLTDTGNLPRINTVESPFAEYSSLDELFQETYKLEQLITQKINELAHAAMTNQDYPTFNFLQWYVSEQHEEEKLFKSIIDKLSLAGKSGEGLYFIDKELSTLDTQN

Part: 1: Making fusion protein of 2 known PDB files using CHIMERA & MODELLER

  • Open the fasta file in Chimera. Then from sequence window menu: Info… Blast Protein to search the PDB for matching structures. In the Blast Protein results, find 4XGS and 2B68 and choose both lines (click, ctrl-click in the results dialog), then click the “Show in MAV” and “Load Structure” buttons at the bottom of the dialog.

  • Now you will have a new sequence alignment window with 3 sequences in it: the fusion protein and the two protein structures, and in the main Chimera window, the two structures.

  • In Chimera, delete all extra protein chains. In other words, if there are extra copies of those structures, just delete them so that you have only one copy to use as the template. Make sure that the remaining copy of each is associated with its sequence in the alignment (sequence alignment window menu: Structures… Associations…)

  • Position the two structures so that the termini are in a somewhat reasonable place relative to each other to template the fusion protein. In our case, the C-terminal of Defensin monomer attaches with the N-terminal of Ferritin monomer. You can “freeze” one in place by deactivating it and move just the other with the mouse as described here:

  • From the sequence alignment window menu choose: Structure… Modeller (homology) to show the Modeller dialog. Choose the query as the target and both structures as the template, etc. as in the modeling tutorials. You may also want to turn on “Use thorough optimization” in the Advanced Options section.

    • NOW YOU CAN SAVE ONE OF THE STRUCTURE AS A PDB FILE AND USE THIS AS A MONOMER FOR FURTHER STEPS.

    Part 2: Make a multimer of the fusion protein using reference PDB structure in Pymol

    Open core assembly ( 4XGS) and fusion monomer (4XGS + 2B68) in pymol.
    Use python based code, superpoz.py. We wrote a small script to generate the assembly structure of our StarCores based on the reference biological assembly of the nude core.

    Result

    Fusion monomers for all the constructs are developed for molecular dynamic simulation studies and assembly of core with monomer for visualization of scaffold proteins.

    Molecular dynamic simulation

    What is Molecular Dynamic Simulations ?

    MD studies are based on Newton second law of motion :
    F = ma
    Where F is the force of an atom, m is the mass of the particle and a is the acceleration.

    As, per universal concept of energy
    F(x) = −∇U(x)

    Where x represents coordinates of all atoms, and U is the potential energy function, Velocity is the derivative of position, and acceleration is the derivative of velocity.

    We can thus write the equations of motion as:
    dx/dt = v
    dv / dt = F(x)/m

    This is a system of ordinary differential equations. For n atoms, we have 3n position coordinates and 3n velocity coordinates. Calculation of “Analytical” (algebraic) solution is impossible but Numerical solution can be given by,
    xi+1 = xi +δ t vi
    vi+1 = vi +δ t F(xi )/m

    This is layman representation of principle used for molecular dynamic simulation studies. Since, in real life and MD also, atoms are in constant motion then motion is used to understand the probability of observing a particular arrangement of atoms as a function of the potential energy.

    Objective :

    To study behaviour of fusion molecule in vicinity of cell.

    Steps :

    • Open YASARA.
    • Go to options > Macro & Movie > Set Target pdb > md_runfast.mcr
    • Wait for the structure to run upto 15ns.
    • Record a video if required.
    • Analyse the result using MD_analyze.mcr
    • Analyse each residue using MD_analyzeres.mcr
      •  
Figure :  Simulation of fusion monomer ferritin + ovispirin (4XGS + 1HU5)

      Parameters :

      • -pH at which the simulation should be run, by default physiological pH 7.4
        The ion concentration as a mass fraction, here we use 0.9% NaCl (physiological solution)ions=‘NaCl,0.9’
      • -Simulation temperature, 298K
      • -Water density = 0.997
      • -Duration of the simulation = 15ns
      • -Extension of the cell on each side of the protein ‘10’ means that the cell will be 20 A larger than the protein.
        Shape of the simulation cell: ‘Cube’.
      • -Forcefield : ForceField AMBER14
      • -Cell boundary : Boundary periodic
        The simulation speed, ‘fast’ (maximize performance with 2*2.5 fs timestep and constraints)
        The save interval for snapshots. Normally you don’t need more than 500-1000 snapshots of your simulation.
        Solute from diffusing around and crossing periodic boundaries. Disable that for simulations of crystals.
        • Centre for Research and Interdisciplinarity (CRI)
        Faculty of Medicine Cochin Port-Royal, South wing, 2nd floor
        Paris Descartes University
        24, rue du Faubourg Saint Jacques
        75014 Paris, France
        paris-bettencourt-2018@cri-paris.org

1. Obtain reference PDB structures representing the core and antimicrobial peptide protein monomers 2. Use MODELLER via CHIMERA interface for homology modelling 3. Choose the best fusion protein model that represents the Star core monomer