Difference between revisions of "Team:Madrid-OLM/Computationaimprovement"

Line 1: Line 1:
 
{{Madrid-OLM/1}}
 
{{Madrid-OLM/1}}
 
<html lang="en">
 
<html lang="en">
 +
   
 
     <style>
 
     <style>
        .tittle-secc{
+
    .lesspad{
             padding-top: 12em !important;
+
             padding-right: 1em !important;
 +
            padding-left: 1em !important;
 
         }
 
         }
        .figureimage{
+
   
            margin-bottom: 0.5em
+
            }
+
 
     </style>
 
     </style>
 
     <head>
 
     <head>
 
         <meta charset="utf-8">
 
         <meta charset="utf-8">
         <title>Computational improvement of the aptamer</title>
+
         <title>Aptamers</title>
 
         <meta name="viewport" content="width=device-width, initial-scale=1.0">
 
         <meta name="viewport" content="width=device-width, initial-scale=1.0">
         <meta name="description" content="Site Description Here">
+
         <meta name="description" content="Page where the team expose how to optains the aptamers">
 
     </head>
 
     </head>
 
     <body class=" ">
 
     <body class=" ">
 
         <a id="start"></a>
 
         <a id="start"></a>
 
         <div class="main-container">
 
         <div class="main-container">
            <section class="page-navigator">
 
                <ul>
 
                    <li>
 
                        <a href="#home" class="inner-link" data-title="Home"></a>
 
                    </li>
 
                    <li>
 
                        <a href="#protein" class="inner-link" data-title="Protein structure"></a>
 
                    </li>
 
                    <li>
 
                        <a href="#aptamer" class="inner-link" data-title="Aptamer sctucture"></a>
 
                    </li>
 
                    <li>
 
                        <a href="#dockin" class="inner-link" data-title="Docking and mutation proposal"></a>
 
                    </li>
 
                    <li>
 
                        <a href="#links" class="inner-link" data-title="Links of interest"></a>
 
                    </li>
 
                </ul>
 
            </section>
 
 
              
 
              
             <section id="home"class="tittle-secc text-center switchable feature-large">
+
           
                <div class="container">
+
             <section class="text-center imagebg space--lg" data-overlay="3">
                    <div class="row justify-content-around">
+
                <div class="background-image-holder">
                        <div class="col-md-10 col-lg-8 ">
+
                    <img alt="background" src="https://static.igem.org/mediawiki/2018/4/4a/T--Madrid-OLM--Aptamer--Mainimage.png" />
                            <h1 id="Teamtittle">Computational improvement of the aptamer</h1>
+
                            <p class="lead">
+
                                The target of this protocol is to improve the affinity of the union between the protein and the aptamer, obtained by the process that had been explained previously, by bioinformatics methods.
+
                            </p>
+
                            <p class="lead">
+
                                To carry out this affinity improvement, you can start from only 2 elements: the DNA aptamer sequence and the name or preferably the amino acid protein sequence. However, any extra information that one could find in previous studies could be so helpful as to save time and reduce the margin of cumulative error. The complete protocol has four different sections: 
+
                            </p>
+
                            <ul>
+
                                <li>
+
                                    <p class="lead">• Obtaining the 3D protein structure.
+
                                    </p>     
+
                                </li>
+
                                <li>
+
                                    <p class="lead">• Obtaining the 3D aptamer structure (most critical point).
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">• A Search of the binding site between protein and aptamer through a docking process.  
+
                                    </p>     
+
                                </li>
+
                                <li>
+
                                    <p class="lead">• Study of the union and proposal of mutation in the sequence.
+
                                    </p>     
+
                                </li>
+
                            </ul>
+
                           
+
                            <p class="lead">
+
                                <br/>It is very important to understand that this process does not guarantee a better result. However it gives the opportunity to improve the aptamer obtained. Although there are many works with proteins structures and the interaction between them, there are almost no studies that work with nucleic acids in simulation or structure prediction of the structures beyond double strand helix. Therefore any final result obtained should be checked in the laboratory, either to confirm or discard the new sequence obtained. important to understand that this process not guarantee a better result. However it gives the opportunity of improve the aptamer obtained. Although there are so many works with proteins structures and the interaction between each other’s, there are almost no studies that work with nucleic acids in simulation or structure prediction terms beyond double strand helix. Therefore any final result obtained should be checked in the laboratory, either to confirm or discard the new sequence obtained.
+
                            </p>
+
                            <p class="lead">
+
                                <br/>You could see all the server that we will mention in the link section:
+
                            </p>
+
                            <a class="btn btn--primary type--uppercase " href="#links" class="inner-link">
+
                                <span class="btn__text">
+
                                    Links of interest
+
                                </span>
+
                            </a>
+
                        </div>
+
                    </div>
+
                    <!--end of row-->
+
 
                 </div>
 
                 </div>
                 <!--end of container-->
+
                 <div class="container pos-vertical-center">
            </section>
+
           
+
            <section id="protein" class="text-center">
+
                <div class="container">
+
 
                     <div class="row">
 
                     <div class="row">
                         <div class="col-md-10 col-lg-8 boxed boxed--border bg--secondary boxed--lg box-shadow">
+
                         <div class="col-md-12">
                             <h2>Obtaining the 3D protein structure</h2>
+
                             <h1 class="h1--large">
                            <p class="lead">
+
                                 Aptamers
                                The final result of this section should be a file in .pdb format where we can find the coordinates of every single atom of the protein. That file could be interpreted by bioinformatics programs like Phymol or Chimera and represent the protein in different graphically. To obtain it, there are two different paths depending on one fact: has the protein been previously crystallographed?
+
                             </h1>
                            </p>
+
                            <h4>A structure obtained experimentally exists</h4>
+
                            <p class="lead">
+
                                We will simply use this proven structure. This will be the most reliable method due to it being the result of a group's investigation. To check if there is any previous structure we recommend the following simple steps:
+
                            </p>
+
                            <ul>
+
                                <li>
+
                                    <p class="lead">1. Access the UniProt server and search for the protein of interest.
+
                                    </p>     
+
                                </li>
+
                                <li>
+
                                    <p class="lead">2. Search for the structure description area and check if it has an entry to the Protein Data Bank (PDB) database.
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">3. In the positive case, access that PDB entry and download the file that is needed.
+
                                    </p>     
+
                                </li>
+
                                <li>
+
                                    <p class="lead">4. In the case of no entry, the steps explained in the next section must be performed.
+
                                </li>
+
                            </ul>
+
                            <p class="lead">
+
                                <br/>Clarification, the PDB is a universal database where are the 3D structure of all proteins which have been obtained its structure experimentally. That is, any structure that we find there is quite reliable and proven. However, predictions of structures are not included, owing to this procedure is not experimental.
+
                            </p>
+
                            <h4>A structure obtained experimentally don’t exists.</h4>
+
                            <p class="lead">
+
                                In this case is necessary to start with a structure prediction process. Although this is not the optimal method because we could start to accumulate error, we can obtain very reliable results. Several different methods can be used for this purpose. We recommend the following two and in this order:
+
                            </p>
+
                            <ul>
+
                                <li>
+
                                    <p class="lead">1. <b>Homology methods:</b> This method looks for sequences of proteins that are included in the PDB that have a high similarity with our protein and adopt its 3D structure. If the result shows a similarity of more than 30% in the entire sequence, we could be certain the structure. For this step we recommend the Swiss-Model server. Otherwise, if it gives an inferior value, we should move on to the next step.
+
                                    </p>     
+
                                </li>
+
                                <li>
+
                                    <p class="lead">2. <b>I-Tasser (Iterative Threading ASSEmbly Refinement): </b> This server uses several methods, as they define themselves: “a hierarchical approach to protein structure and function prediction. It first identifies structural templates from the PDB by multiple threading approach LOMETS, with full-length atomic models constructed by iterative template fragment assembly simulations.” This method is currently the most reliable and they have been the winners in predicting structure in the last 6 editions of the CAP tests.
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">3. As an extra option we mention <b>Robetta: </b> It is an ab-initio simulation method, where it tries to create the protein from 0, using physical simulation and gross computational power. It works best for short sequences and can take a long time to get reliable results.
+
                                    </p>     
+
                                 </li>
+
                             </ul>
+
                            <br/>
+
                            <h4>Energy minimization.</h4>
+
                            <p class="lead">
+
                                This last step is completely recommendable whatever the origin of the file with the structure be. It is necessary to have the structure in a state, as stable as possible. To minimize the structure we can employ 2 strategies:
+
                            </p>
+
                            <ul>
+
                                <li>
+
                                    <p class="lead">1. <b>Quick vacuum minimization</b> This is the simplest method. For this purpose we will use the Chimera program and within its tools we will choose the "Energy minimization". Then we will adjust the parameters to what is most convenient for us
+
                                    </p>     
+
                                </li>
+
                                <li>
+
                                    <p class="lead">2. <b>Minimization in aqueous medium: </b> Complex method, for which we recommend looking for a tutorial to understand the management of the programs used and the necessary concepts. In this, more complex case, we will use the VMD program (Visual Molecular Dynamics). It will be necessary to build a structure around the protein with the "Automatic PSF builder" tool, which includes water molecules and ions. Subsequently a simulation of energy minimization will be carried out with the tool "AutoIMD" which in turn calls an external program (NAND) that must be installed on the PC. Finally, the resulting file should be saved and opened in Chimera to eliminate all water molecules and leave only the structure of the protein.
+
                                    </p>
+
                                </li>
+
                            </ul>
+
                            <br/>
+
                            <img class= "figureimage" alt="Figure1" src="https://static.igem.org/mediawiki/2018/b/bc/T--Madrid-OLM--Aptamer--Modelization--OleE1Comparition.png" style="width:60%;"/>
+
                            <p class="lead" style="margin-left:10%; margin-right:10%;">Figure 1. Comparison of the structure of Ole E1 before (purple) and after (blue) an energy minimization process.
+
                            </p>
+
                            <img class= "figureimage" alt="Figure2" src="https://static.igem.org/mediawiki/2018/2/23/T--Madrid-OLM--Aptamer--Modelization--OleE1ribbsurf.png" style="width:75%;"/>
+
                            <p class="lead" style="margin-left:10%; margin-right:10%;">Figure 2. Ole E1 after an energy minimization process represented as ribbon and as ASA (Accessible Surface Area). Made by Chimera tool..
+
                            </p>
+
                           
+
 
                         </div>
 
                         </div>
 
                     </div>
 
                     </div>
Line 166: Line 37:
 
             </section>
 
             </section>
 
              
 
              
             <section id="aptamer" class="text-center">
+
             <section id="story" class="text-center">
 
                 <div class="container">
 
                 <div class="container">
 
                     <div class="row">
 
                     <div class="row">
                         <div class="col-md-10 col-lg-8 boxed boxed--border bg--secondary boxed--lg box-shadow">
+
                         <div class="col-md-10 col-lg-8">
                             <h2>Obtaining the 3D aptamer structure</h2>
+
                             <h2>What is an Aptamer?</h2>
 
                             <p class="lead">
 
                             <p class="lead">
                                 Previously there were two options, that the protein had a known structure or not. In this case it is a DNA chain that we have personally obtained in the laboratory, so the possibility of finding a structure published for that exact sequence it's almost zero. In this way the only case to consider is to try to predict its structure. This action presents a series of drawbacks nowadays. While for proteins there are numerous methods and servers with high reliability, for nucleic acids there are hardly any servers and prediction methods. This is due to the structure of proteins has always been very important and varied for its function, but from the genetic material, the important thing has always been its code, whether in the form of DNA or RNA. For the case of RNAs, there is more than one function, reason why in certain occasions if it has been attended to the need to know their structures but it has not been developed correctly. Even so, we propose a series of steps to predict the structure of the aptamer. Although this predicted structure will not be reliable and it is important to understand this, it can serve us for our objective of proposing an improvement in the initial sequence.
+
                                 Aptamers are single-stranded nucleic acids with defined and stable tertiary structures that bind with high affinity and specificity to target molecules. Aptamers adopt a wide variety of structural motifs, such as stems, hairpins, bulges, pseudoknots, and G-quadruplexes, recognizing with their target via diverse interactions such as hydrogen bonds, hydrophobic interactions, van der Waals interactions, aromatic stacking, and so on.
 
                             </p>
 
                             </p>
                            <h4>Step 1: Obtaining the secondary structure</h4>
 
 
                             <p class="lead">
 
                             <p class="lead">
                                 As there is no single server as in proteins to obtain a final result giving only the sequence, we have to divide the process into different stages. In this first stage we need to obtain a secondary structure of the single strand of DNA. The server we recommend is the MFold server of the University of Albany. It provides both a graphic map of the links between bases, and the secondary structure in text format (in the Vienna format that we need later), as we show in this example:
+
                                 Their affinity can be compared to antibodies, but the cost of their production is significantly less expensive than the antibodies.
 
                             </p>
 
                             </p>
                            <p class="lead" style="text-align: center;">
 
                                GTGACGTAGGTTGGTGTGGTTGGGGCGTCAC <br/>
 
                                (((((((.................)))))))
 
                            </p>
 
                            <img class= "figureimage" alt="Figure3" src="https://static.igem.org/mediawiki/2018/5/59/T--Madrid-OLM--Aptamer--Modelization--MapAptamer2Structure.png" style="width:40%;"/>
 
                            <p class="lead" style="margin-left:20%; margin-right:10%;">Figure 3. Example image of the links between base pairs that the server provides.
 
                            </p>
 
                            <h4>Step 2: Prediction of 3D structure</h4>
 
 
                             <p class="lead">
 
                             <p class="lead">
                                 The next step is to predict the tertiary structure. The only server capable of doing that nowadays is ROSIE and it requires the secondary structure that we have previously obtained in Vienna format in addition to the sequence. This server predicts RNA structures (not DNA), so it will be necessary change the T for the U in the sequence.
+
                                 The steps we have been gone through to obtain aptamers are:
                            </p>
+
                            <p class="lead">
+
                                The result obtained will be a series of three-dimensional RNA structure proposals in .pdb format, so a posteriori it will be necessary to make file modifications to eliminate the excess atoms and modify the links so that it can be identified as a DNA chain.
+
                            </p>
+
                            <p class="lead">
+
                                We must take into account that the algorithm used by this server is not as tight as in proteins and usually commit some deviations in the ab initio simulation. It is particularly very important to know if our aptamer has incorporated some ion because this method is not taken into account by aptamers that incorporate external molecules and in those cases, it usually deviates a lot in the final result.
+
                            </p>
+
                            <p class="lead">
+
                                Subsequently it will be necessary to perform a process of energy minimization as explained in the previous section, choosing the method that is preferred to have the most stable molecule.
+
                            </p>
+
                            <img class= "figureimage" alt="Figure4" src="https://static.igem.org/mediawiki/2018/e/e9/T--Madrid-OLM--Aptamer--Modelization--Aptamer3Dcompare.png" style="width:80%;"/>
+
                            <p class="lead" style="margin-left:10%; margin-right:10%;">Figure 4. Trombine aptamer structure comparison between the crystallographed in a lab (blue) and the prediction one (grey). We can see that the blue one have a ion that conditions the structure.
+
                            </p>
+
                            <p class="lead">
+
                                <b>Fact: </b>We tried to fine-tune this step with an aptamer already characterized by thrombin. Although theoretically, the protocol was logical, the results obtained with the set-up led us to conclude that we did not have it prepared yet due to the high error rate. We tried to fine-tune this step with an aptamer already characterized by thrombin. Although theoretically, the protocol was logical, the results obtained with the set-up led us to conclude that we did not have it ready yet due to the high error rate.
+
 
                             </p>
 
                             </p>
 
                         </div>
 
                         </div>
 
                     </div>
 
                     </div>
 
                     <!--end of row-->
 
                     <!--end of row-->
                </div>
+
                   
                <!--end of container-->
+
                   
            </section>
+
                    <div class="row">
 +
                        <div class="col-md-1">
 +
                        </div>   
 +
                        <div class="col-md-3 ">
 +
                            <div class="boxed boxed--border bg--secondary boxed--lg box-shadow lesspad">
 
              
 
              
            <section id="dockin" class="text-center">
+
                                <a class="btn btn--primary-2 btn--sm type--uppercase" href="https://2018.igem.org/Team:Madrid-OLM/AptDiscovery">
                <div class="container">
+
                                    <span class="btn__text">
                    <div class="row">
+
                                        Aptamer discovery
                        <div class="col-md-10 col-lg-8 boxed boxed--border bg--secondary boxed--lg box-shadow">
+
                                    </span>
                            <h2>Docking procces and mutation proposal</h2>
+
                                </a>
                            <p class="lead">
+
                                <p class="lead">Find a binding aptamer with our target protein.</p>
                                After obtaining the structures of the protein and the aptamer, we need to know in which area and with which orientation the interaction between them takes place; Studying which pairs of bases and which amino acids are involved in this interaction and propose a change of 1 or 2 bases in order to make this interaction stronger.
+
                                           
                            </p>
+
                             </div>
                            <h4>Docking</h4>
+
                            <p class="lead">
+
                                The docking is the computational simulation method in which the ideal orientation of an organic molecule is calculated when interacting with another one forming a stable complex (if it exists). Generally this process is designed for protein-protein interactions, the reason why most servers are optimized for these cases. However, we recommend using the NPDock server. This server allows the simulation of the interaction between a protein and a single-stranded nucleic acid molecule, having to identify whether it is DNA or RNA before sending the work.
+
                            </p>
+
                            <img class= "figureimage" alt="Figure5" src="https://static.igem.org/mediawiki/2018/a/a6/T--Madrid-OLM--Aptamer--Modelization--NPDockResult.png" style="width:80%;"/>
+
                            <p class="lead" style="margin-left:10%; margin-right:10%;">Figure 5. Interaction between thrombine and an aptamer simulated by the server NPDock.
+
                            </p>
+
                            <h4>Proposal of mutation</h4>
+
                            <p class="lead">
+
                                With the .pdb file obtained from the previous step, we only have the most subjective step. In this case we have to perform a study of the interaction and see which atoms interact with each other. For this, the most useful thing is to use a visualizer like the ones we mentioned at the beginning (Chimera or phymol). With the evaluation of the different options will have to interpret if making any change in the aptamer, these interactions could improve
+
                            </p>
+
                            <p class="lead">
+
                                The proposed mutation cannot be accepted blindly. It will be necessary to synthesize the proposed aptamer and check in the laboratory if its affinity is superior to that of the original aptamer.
+
                             </p>
+
 
                         </div>
 
                         </div>
                    </div>
+
                        <div class="col-md-3 ">
                    <!--end of row-->
+
                            <div class="boxed boxed--border bg--secondary boxed--lg box-shadow lesspad">
                </div>
+
                <!--end of container-->
+
            </section>
+
 
              
 
              
 +
                                <a class="btn btn--primary-2 btn--sm type--uppercase" href="https://2018.igem.org/Team:Madrid-OLM/AptCharacterization">
 +
                                    <span class="btn__text">
 +
                                        Aptamer Characterization
 +
                                    </span>
 +
                                </a>
 +
                                <p class="lead">Test the affinity of the aptamers we have selected.</p>
 +
                                           
 +
                            </div>
 +
                        </div>
 +
                        <div class="col-md-3 ">
 +
                            <div class="boxed boxed--border bg--secondary boxed--lg box-shadow lesspad">
 
              
 
              
            <section id="links" class="text-center">
+
                                <a class="btn btn--primary-2 btn--sm type--uppercase" href="https://2018.igem.org/Team:Madrid-OLM/Computationaimprovement">
                <div class="container">
+
                                     <span class="btn__text">
                    <div class="row">
+
                                         Computational improvement
                        <div class="col-md-10 col-lg-8 boxed boxed--border bg--secondary boxed--lg box-shadow">
+
                                     </span>
                            <h2>Links of interest</h2>
+
                                 </a>
                            <ul>
+
                                 <p class="lead">Try to improve the affinity between the aptamers and the target.</p>
                                <li>
+
                                           
                                    <p class="lead">
+
                             </div>
                                        <a href="https://www.uniprot.org/" >UniProt:</a>Database with all proteins and the related information.
+
                                    </p> 
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="https://www.rcsb.org/" >PDB:</a>Database with all 3D shapes of protein, nucleic acids, and complex assemblies.
+
                                     </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                         <a href="https://pymol.org/2/" >Phymol:</a>Is a user-sponsored molecular visualization system on an open-source foundation, maintained and distributed by Schrödinger.
+
                                     </p>
+
                                 </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="https://www.cgl.ucsf.edu/chimera/" >Chimera:</a>A highly extensible program for interactive visualization and analysis of molecular structures and related data, including density maps, supramolecular assemblies, sequence alignments, docking results, trajectories, and conformational ensembles.
+
                                    </p>
+
                                 </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="https://www.ks.uiuc.edu/Research/vmd/" >VMD:</a>A molecular visualization program for displaying, animating, and analyzing large biomolecular systems using 3-D graphics and built-in scripting.
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="https://swissmodel.expasy.org/" >Swiss-Model:</a>A fully automated protein structure homology-modelling server, accessible via the ExPASy web server, or from the program DeepView (Swiss Pdb-Viewer).
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="https://zhanglab.ccmb.med.umich.edu/I-TASSER/" >I-tasser:</a>A hierarchical approach to protein structure and function prediction.
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="http://robetta.bakerlab.org/" >Robbeta:</a>    Ab-inition full chain protein structures prediction server.
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="http://unafold.rna.albany.edu/?q=mfold/DNA-Folding-Form" >Albany University MFold:</a>    Web server for the prediction of secundary structures in DNA single-strand chains.
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="http://rosie.rosettacommons.org/rna_denovo/submit" >Rossie:</a>    Web server for ab-initio RNA single-stand structure prediction.
+
                                    </p>
+
                                </li>
+
                                <li>
+
                                    <p class="lead">
+
                                        <a href="http://genesilico.pl/NPDock" >NPDock:</a>  Web server for modeling of RNA-protein and DNA-protein complex structures.
+
                                    </p>
+
                                </li>
+
                             </ul>
+
 
                         </div>
 
                         </div>
                     </div>
+
                     </div>    
                    <!--end of row-->
+
 
                 </div>
 
                 </div>
 
                 <!--end of container-->
 
                 <!--end of container-->
 
             </section>
 
             </section>
 
+
           
                     
+
           
 +
           
 +
           
 +
           
 +
                   
 
         </div>
 
         </div>
 
         <!--<div class="loader"></div>-->
 
         <!--<div class="loader"></div>-->

Revision as of 21:20, 13 October 2018

Madrid-OLM

Aptamers
background

Aptamers

What is an Aptamer?

Aptamers are single-stranded nucleic acids with defined and stable tertiary structures that bind with high affinity and specificity to target molecules. Aptamers adopt a wide variety of structural motifs, such as stems, hairpins, bulges, pseudoknots, and G-quadruplexes, recognizing with their target via diverse interactions such as hydrogen bonds, hydrophobic interactions, van der Waals interactions, aromatic stacking, and so on.

Their affinity can be compared to antibodies, but the cost of their production is significantly less expensive than the antibodies.

The steps we have been gone through to obtain aptamers are:

Aptamer discovery

Find a binding aptamer with our target protein.

Aptamer Characterization

Test the affinity of the aptamers we have selected.

Computational improvement

Try to improve the affinity between the aptamers and the target.