ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Note
Revised

A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling

[version 3; peer review: 3 approved]
PUBLISHED 15 Nov 2017
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

It was demonstrated in recent studies that some rhodopsins can be used in optogenetics as fluorescent indicators of membrane voltage. One of the promising candidates for these applications is archaerhodopsin-3. While it has already shown encouraging results, there is still a large room for improvement. One of possible directions is increasing the intensity of the protein's fluorescent signal. Rational design of mutants with an improved signal is an important task, which requires both experimental and theoretical studies. Herein, we used a homology-based computational approach to predict the three-dimensional structure of archaerhodopsin-3, and a Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid approach with high-level multireference ab initio methodology (SORCI+Q/AMBER) to model optical properties of this protein. We demonstrated that this methodology allows for reliable prediction of structure and spectral properties of archaerhodopsin-3. The results of this study can be utilized for computational molecular design of efficient fluorescent indicators of membrane voltage for modern optogenetics on the basis of archaerhodopsin-3.

Keywords

optogenetics, archaerhodopsin, protein structure prediction, QM/MM, spectral tuning in rhodopsins

Revised Amendments from Version 2

Additional information and data analysis were added in Results section of the article: 1. The comparison of results of different sequence alignment algorithms (new Figure 1 provided); 2. The comparison of the obtained model of archaerhodopsin-3 with the crystallographic structure of archaerhodopsin-1 (new Figure 2, Figure 3); 3. Analysis of the results of absorption maxima calculations for archaerhodopsin-3. Information about analogous results for other rhodopsins is provided. 4. Residues replaced in site-directed mutagenesis study (Mclsaac RS et al,  PNAS, 2014) were visualized (new Figure 4).
 
Modifications made in the Introduction section and Abstract, considering the remarks of our referees. In the Methods section information about sequence identity between archaerhodopsin-1 and archaerhodopsin-3 is provided.

Also, the new variant of the archieve with input and output files, extra grant information (Russian Foundation for Basic Research and Presidium of the RAS) and acknowledgements to two computational research facilities (HPC computing resources at Lomonosov Moscow State University and Computer Center of SPbU) are provided.

The additional grants listed, have aided in the continuation of the project that was supported by previous funding.

See the authors' detailed response to the review by Adam E. Cohen and Samouil Farhi
See the authors' detailed response to the review by Jeremy N. Harvey

Introduction

Precise and quick control of physiological processes using integrated optical and genetic methods is a vast area with a high number of important applications1. One possible approach in this field is using fluorescent voltage-dependent indicators for detecting the activity of mammalian neurons, which allows achievement of the precision of a single neuron without perceptible time delays. It was recently shown that some rhodopsins, especially achaerhodopsin-3, can be potential candidates for such a task25. While the results obtained for archaerhodopsin-3 are already very encouraging, increasing the fluorescent signal of such proteins can lead to a significant progress in this field. It was also shown that the insertion of mutations into these proteins can dramatically improve the signal quality25.

Unfortunately, the fundamental mechanisms underlying the processes that determine fluorescence are not well understood. This lack of knowledge leads to difficulties with design of desired rhodopsin mutants. While rational design is not based on a solid foundation and, for this reason, is not very effective, another experimental approach, random mutagenesis, is very time consuming. Computational studies can provide additional insights into the problem. One of the main obstacles for computational modeling of proteins is the absence of three-dimensional structures of high quality, especially for membrane proteins, which are a challenge for crystallization. On the other hand, computational prediction of three-dimensional structures is not trivial. The goal of this study was to obtain a good-quality structure for achaerhodopsin-3, one of the most used voltage-dependent fluorescent sensors4,5, and based on this structure to predict the optical properties of this protein.

To achieve this goal, we used a homology-based computational approach for structure prediction. As the choice of a structure prediction algorithm is not straightforward, we tested several methods. To evaluate the quality of obtained structures we performed subsequent Quantum Mechanics/Molecular Mechanics (QM/MM) calculations of absorption maxima and compared the results with available experimental data.

Methods

The structure of archaerhodopsin-3 was built using a homology modeling approach. Primary structures for all rhodopsins with crystallographic data are available in the Protein Data Bank library6 (24 structures as of September 2016) were compared with the primary structure of archaerhodopsin-3. Archaerhodopsin-1, which has the highest sequence identity to the target protein, was chosen as a template (RCSB code 1UAZ, sequence identity 84.5%, sequence similarity 91.5%). Three algorithms of homology-based model building were tested: Medeller7, I-TASSER8 and RosettaCM9. All methods of homology modeling heavily rely on externally made target-template alignment of primary sequences, which serves as the main instruction for model building. For this reason, we tested three algorithms of pairwise alignment using their results as an input for each method of model building. Two of the alignment methods are specifically constructed for membrane proteins, MP-T10 and AlignMe11, and the third one, MUSTER12, gains its quality from evolutionary predictions. The latter algorithm is a built-in algorithm of I-TASSER suite and its results were used only for this method of structure prediction.

Before QM/MM calculations, several preparation steps were performed: hydrogen atoms were added using pdb2pqr package13 version 2.1.1 using CHARMM force field version 2714, pH=7; hydrogen atoms were equilibrated by energy minimization in NAMD package15, version 2.11. The retinal chromophore was bound to the lysine residue Lys226, whole lysine + retinal system was parameterized in CHARMM force field16. The protein was inserted in the POPC membrane17, the whole system was inserted in a water solvent box, the TIP3P water model18 was used, the size of water box was selected so that there were at least 10 Å from any atom of protein to the edge of the system, and the system was neutralized by addition of Na+ and Cl- ions.

Relaxation of the system was performed in several steps: relaxation of retinal + lysine complex with all other atoms fixed, relaxation of all atoms that were within 6 Å of the chromophore system, relaxation of whole protein and water box. During all these steps the following parameters were used: a 10 Å cutoff with switching starting at 8.5 Å was applied to the electrostatics and van der Waals interactions; Particle Mesh Ewald method19 was used for dealing with electrostatics interactions, grid spacing 1Å. Equilibrated protein structure was extracted; internal waters were added into protein cavities using WaterDock program20.

To calculate absorption maxima, we used the methodology that has been proven as efficient in a number of our previous studies for different kind of rhodopsins and rhodopsin mimics2127. The structures of the archaerhodopsin-3 obtained at the previous step were optimized using two-layer ONIOM (QM:MM-EE) scheme. (QM=B3LYP/6-31G*; MM= AMBER for aminoacids and TIP3P for water, EE=electronic embedding) was implemented in the Gaussian09 package28. To calculate the spectral properties of the chromophore in the presence of the protein environment (described as AMBER point charges) SORCI+Q/6-31G* level of the theory was used, as it is implemented in ORCA6.0 package29. For details of previously performed QM/MM methodology see Altun et al.30,31.

Results

On the first step of structure prediction of archaerhodopsin-3 we performed the alignment of its amino acid sequence with the sequence of archaerhodopsin-1, comparing three different algorithms: AlignMe, MUSTER and MP-T. The results of AlignMe and MUSTER algorithms were identical; they differ from the MP-T alignment result in the flexible N-terminal part (Figure 1, a–c).

On the next step we built 7 three-dimensional models of archaerhodopsin-3 using I-TASSER, Medeller and RosettaCM algorithms of template-based structure prediction. Comparison of crystallographic structure of archaerhodopsin-1 with the predicted structure of archaerhodopsin-3 showed that these two proteins differ only in loop regions and on the edges of the alpha-helices (Figure 2, Figure 3). These regions are located on relatively large distance from the retinal chromophore. We also highlighted the residues, which were the sites for the site-directed mutagenesis of the archaerhodopsin-35 (Figure 4). The computational characterization of these mutants is an interesting problem for further investigations.

For all seven predicted models we calculated the absorption maximum wavelengths. All obtained models provide reasonable results -- the deviation range is only 31 nm from the experiment (Table 1). This deviation range is sensible considering the approximations used in our methodologies. In our previous studies, we obtained a 29 nm shift for the absorption maximum of halorhodopsin from N.pharaonis23, a 25 nm shift for the absorption maximum of archaerhodopsin-2 and a 39 nm shift for the channelrhodopsin-2. The model with the smallest deviation was predicted by I-TASSER algorithm with AlignMe alignment (Figure 2).

62feff92-316f-4d5b-80d7-4fbf21221b0b_figure1.gif

Figure 1.

The results of the pairwise sequence alignment algorithms: a) AlignMe; b) MP-T; c) MUSTER.

62feff92-316f-4d5b-80d7-4fbf21221b0b_figure2.gif

Figure 2. The predicted structure of archaerhodopsin-3 (green) superimposed on crystallographic structure of archaerhodopsin-1 (blue).

62feff92-316f-4d5b-80d7-4fbf21221b0b_figure3.gif

Figure 3. The illustration of the differences between structures of archaerhodopsin-1 and archaerhodopsin-3.

62feff92-316f-4d5b-80d7-4fbf21221b0b_figure4.gif

Figure 4. The visualization of residues, which were replaced in the site-directed mutagenesis study of archaerhodopsin-35.

Table 1. Absorption spectrum maximum of achaerhodopsin-3 for different models.

Alignment
method
Model building methodλmax, nm
--Experimental wild-type structure556
AlignMeI-TASSER578
MP-TI-TASSER581
MUSTERI-TASSER581
AlignMeRosettaCM587
MP-TRosettaCM585
AlignMeMedeller580
MP-TMedeller586

Conclusions

In this study, we predicted the structure of fluorescent voltage-dependent sensor achaerhodopsin-3 and evaluated its quality with subsequent QM/MM high level ab initio calculations of spectral properties. The calculated absorption maximum is within 31 nm from the experimental value. Several methods of model building were tested and spectral characteristics were calculated for all resulting models. We showed that our methodology allowed for reliable prediction of optical properties of archaerhodopsin-3. The results of this study can be utilized for high-level QM/MM investigation of different aspects of photochemistry of this voltage-dependent fluorescent sensor and, therefore, to contribute in development of the efficient molecular tools for modern optogenetics.

Data availability

The sequence of archaerhodopsin-3 was taken from Uniprot database: http://www.uniprot.org/uniprot/P96787

The template for homology modeling was taken from PDB database (rcsb code 1UAZ): http://www.rcsb.org/pdb/explore/explore.do?structureId=1UAZ

The input files of I-TASSER suite, RosettaCM, Medeller algorithms with corresponding README files, zipped output of I-TASSER suite, scripts for processing structure after homology modeling stage (with instructions in README file), input files for spectra calculations are available: doi, https://doi.org/10.5281/zenodo.83002532.

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 11 Jan 2017
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Nikolaev DM, Emelyanov A, Boitsov VM et al. A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling [version 3; peer review: 3 approved]. F1000Research 2017, 6:33 (https://doi.org/10.12688/f1000research.10541.3)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 3
VERSION 3
PUBLISHED 15 Nov 2017
Revised
Views
5
Cite
Reviewer Report 22 Nov 2017
Jeremy N. Harvey, Division of Quantum Chemistry and Physical Chemistry, Department of Chemistry, KU Leuven, Leuven, Belgium 
Approved
VIEWS 5
I am satisfied ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Harvey JN. Reviewer Report For: A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling [version 3; peer review: 3 approved]. F1000Research 2017, 6:33 (https://doi.org/10.5256/f1000research.14321.r28011)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
8
Cite
Reviewer Report 22 Nov 2017
Adam E. Cohen, Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA 
Approved
VIEWS 8
The authors have improved the manuscript in response ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Cohen AE. Reviewer Report For: A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling [version 3; peer review: 3 approved]. F1000Research 2017, 6:33 (https://doi.org/10.5256/f1000research.14321.r28010)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Version 2
VERSION 2
PUBLISHED 17 Jan 2017
Revised
Views
28
Cite
Reviewer Report 31 Mar 2017
Jeremy N. Harvey, Division of Quantum Chemistry and Physical Chemistry, Department of Chemistry, KU Leuven, Leuven, Belgium 
Approved with Reservations
VIEWS 28
The very short manuscript reports homology modelling of the membrane protein, coupled with 'validation' by using vertical excitation energies computed with a QM/MM protocol. I think the work is interesting, particularly since the authors have made their data available, making ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Harvey JN. Reviewer Report For: A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling [version 3; peer review: 3 approved]. F1000Research 2017, 6:33 (https://doi.org/10.5256/f1000research.11514.r20540)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 15 Nov 2017
    Mikhail N. Ryazantsev, Saint-Petersburg State University, St. Petersburg, Russian Federation
    15 Nov 2017
    Author Response
    Dear Prof. Harvey,
     
    Thank you for your review of our article “A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling”. We have carefully read it and made ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 15 Nov 2017
    Mikhail N. Ryazantsev, Saint-Petersburg State University, St. Petersburg, Russian Federation
    15 Nov 2017
    Author Response
    Dear Prof. Harvey,
     
    Thank you for your review of our article “A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling”. We have carefully read it and made ... Continue reading
Views
36
Cite
Reviewer Report 20 Mar 2017
Adam E. Cohen, Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA, USA 
Samouil Farhi, Harvard University, Cambridge, MA, USA 
Approved with Reservations
VIEWS 36
Nikolaev et al. attempt to predict archaerhodopsin-3 structural and optical properties using in silico approaches, layering QM/MM calculations of absorption spectra over homology based structural modeling. As the authors note, such work would be invaluable for guiding the design of ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Cohen AE and Farhi S. Reviewer Report For: A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling [version 3; peer review: 3 approved]. F1000Research 2017, 6:33 (https://doi.org/10.5256/f1000research.11514.r20866)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 15 Nov 2017
    Mikhail N. Ryazantsev, Saint-Petersburg State University, St. Petersburg, Russian Federation
    15 Nov 2017
    Author Response
    Dear Prof. Cohen and Samouil Farhi.
     
    Thank you for your review of our article “A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling”. We have carefully read ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 15 Nov 2017
    Mikhail N. Ryazantsev, Saint-Petersburg State University, St. Petersburg, Russian Federation
    15 Nov 2017
    Author Response
    Dear Prof. Cohen and Samouil Farhi.
     
    Thank you for your review of our article “A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling”. We have carefully read ... Continue reading
Views
26
Cite
Reviewer Report 31 Jan 2017
Kirill A. Afonin, Nanoscale Science Program, Department of Chemistry, The Center for Biomedical Engineering and Science (CBES), The University of North Carolina at Charlotte, Charlotte, NC, USA 
Approved
VIEWS 26
Overall the paper by Nikolaev et al., is well written and an interesting part of the photo- and biochemical realm. The paper provides very promising results and a unique way of analyzing structure functionality. I am interested to see what ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Afonin KA. Reviewer Report For: A voltage-dependent fluorescent indicator for optogenetic applications, archaerhodopsin-3: Structure and optical properties from in silico modeling [version 3; peer review: 3 approved]. F1000Research 2017, 6:33 (https://doi.org/10.5256/f1000research.11514.r19223)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.

Comments on this article Comments (0)

Version 3
VERSION 3 PUBLISHED 11 Jan 2017
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.