Journal of Chemistry

Journal of Chemistry / 2018 / Article

Research Article | Open Access

Volume 2018 |Article ID 4173479 | 10 pages |

Identification of Lutzomyia longipalpis Odorant Binding Protein Modulators by Comparative Modeling, Hierarchical Virtual Screening, and Molecular Dynamics

Academic Editor: Fabio Polticelli
Received19 May 2018
Accepted22 Jul 2018
Published24 Sep 2018


Visceral leishmaniasis (VL) is the second most important vector-borne disease in the world. It is transmitted by Lutzomyia longipalpis in America; therefore, controlling the vector is essential to prevent the disease, especially using traps with chemical attractants. It is known that odorant binding proteins (OBPs) act at the first odor selection level, so in silico methodology was used to identify putative vector chemical modulators based on OBPs on known ligand structures. Therefore, 3D structures of L. longipalpis OBP were predicted through different comparative modeling methods. The best model was subjected to molecular dynamics studies. Then, a hierarchical virtual screening approach filtered OBP modulator-like compounds from ZINC12 biogenic database based in global chemical space, using principal components from ChemGPS-NP server. Such compounds then were evaluated and ranked according to their affinity with the OBP orthosteric site by molecular docking in DOCK 6.7. The compounds were scored by Grid Score function and top five ranked poses had their intermolecular complex interactions analyzed in PLIP server. Most ligands in the top of the rank were lysophospholipids, which could potentially interact with the OBP hydrophobic pocket through Phe72, Tyr76, Ile79, Ala87, Lys88, Asp92, Phe61, Leu75, Trp113, His120, and Phe122 residues and H-bonding with His120 and Phe122. Next, compounds in the top of the rank were evaluated by 50 ns MD and the results showed that the phosphate group of these compounds could set a salt bridge with His110. Additionally, Tyr76, Ala87, Met91, Trp113, and Phe122 were important to hydrophobic interactions with the ligand. These results highlight the importance of accurate assessments such as MD studies in order to analyze the docking results in the identification of new odorant modulators.

1. Introduction

Vector-borne diseases are one of the major and most concerning human epidemic issues to be solved. They affect more than half of the world population and cause more than a million deaths per year [1]. Visceral leishmaniasis (VL) is an important vector-borne disease manifested through systemic infections caused by parasites belonging to genus Leishmania. It is transmitted by the bites of female phlebotomine Lutzomyia longipalpis, only [2]. Estimates assume that 500,000 new people get VL in a yearly basis and that 60,000 of them die due to the disease [3, 4]. The main vector-control strategy applied to prevent and combat vector-borne diseases rely on using insecticides such as dichlorodiphenyltrichloroethane (DDT) [5]. However, damages to the environment and to human health are caused by DDT accumulation restraint due to the application of these insecticides [68].

Chemical ecology is an alternative approach to insecticide application, besides being a nonpolluting procedure. This approach is based on volatile chemical compounds used to disturb and manage the chemical odorant sensing of insect vectors in order to modulate their behavior [911]. The odorant binding proteins (OBPs) are the first molecular filters favoring insects’ sensitivity because they act in odor selection [1113]. These proteins are responsible for the selective captation and transport of small hydrophobic compounds (odors) into the sensilar lymph, for an aqueous solution that fills the sense organs into membrane olfactive receptors in the neuron cell and for carrying pheromones into insects glands [14, 15]. Therefore, OBPs structures can guide the identification of putative-behavior modulator compounds according to their affinity with the transporter binding site [16, 17].

Furthermore, there are some advantages on dealing with OBPs because they are extracellular soluble proteins, and their structure is lesser complex than that of olfactory receptors with seven transmembrane domains (7TM). So, these molecular targets have been used in a new method known as chemical reverse ecology since it uses the protein-target structure rather than the structure of known similar compounds to develop chemicals that are attractive and/or repellent to insects by using computational tools [16, 17].

Computational tools allow the fast and cheaper assessment to big virtual compound libraries and help identifying bioactive/modulative molecules. Overall, this virtual screening (VS) process may be guided by the structure of the protein binding site or by chemical similarity to the structures of known modulators [18]. Hierarchical virtual screening (HVS) methods use all the information about the structure (from protein and ligands) to strategically filter compound libraries.

Thereby, the main aim of the present study was to search for putative chemical modulators of Lutzomyia longipalpis behavior by using the structure of an OBP found in the transcriptome of the insect’s male sexual organs (RAAPBAR014D07) and also in the cDNA of its antennae (BAM011G02) [19, 20]. This work rose from previous computational studies focused in L. longipalpis OBPs [21] in which 3D OBP comparative models were generated through different methods, and the best ones were chosen based on their overall quality (validation). Then, they were optimized through minimization procedures and subjected to molecular dynamics (MD) simulations. Finally, these comparative models were used to screen a virtual library of natural compounds, which was previously filtered by ligand-based virtual screening, i.e., a scheme composed of different hierarchical steps.

2. Materials and Methods

2.1. Comparative Modeling: Template Selection

The peptide sequence of Lutzomyia longipalpis OBP4 (RAAPBAR014D07) was subjected to the SignalIP v. 4.1 server [22] in order to identify and further remove the signal peptide. Template selection was achieved by using a sequence of similar searches at PSI-BLAST [23] against sequences of 3D structures experimentally solved and deposited in the Protein Data Bank (PDB) repository [24]. Only proteins presenting aligned sequences with at least 45% identity, structure with 2 Å resolution, and 20% R-in complex with ligand (holo structure) were taken into consideration [2527].

2.2. Comparative Modeling: Construction of L. longipalpis OBP 3D Model

The target-template sequence alignment was subjected to the SWISS-MODEL server ( [28], to the MODELLER software [29, 30], and to the I-TASSER server [31] ( for structure prediction methods applied to rigid body assembly, spatial restraints, and threading, respectively.

2.3. Comparative Modeling: L. longipalpis OBP Model Validation

The OBP4 models early developed from different prediction methods were visually examined in the UCSF Chimera 1.9.8 [30] and Discovery Studio Visualizer 4.1 software [32]. They had their steric quality estimated through Ramachandran plots generated in the PROCHECK v.3.5.4 software [33]. Energy quality examinations were conducted through ANOLEA in order to corroborate the overall quality of the models [34] (; global quality estimations were performed through QMEAN6 [35] and Z-SCORE [36]. The best-quality model was chosen to optimize the procedures through minimization cycles and molecular dynamics (MD) simulation of the model in the complex (with template bound ligand in holo crystal structure). It was done to reproduce solvation effects such as ligand-protein intermolecular interactions [37].

2.4. Molecular Dynamics Simulations

The calculations described below were performed in GROMACS v.5 with GROMOS 53a6 force field, in a NPT ensemble [38]. Firstly, the protonation states of early-model ionizable residues were calculated in PROPKA in order to mimic the pH (pH 7.5) of the sensilar lymph [39]. Next, the OBP model was complexed with the ligand found in the template structure, and the system was solvated in 7496 SPC-E water models [40] (without the addition of counterions) in a dodecahedron crystal cell of 71.92 Å lengths (A, B, and C) with each angle set at 60° (alpha, beta, and gamma). The protein was placed in the center of the box (13 Å in relation to box vectors), and the energy in the system was minimized through two steepest descent algorithm (SD) routines: first, all protein heavy atoms were restrained to their initial positions using 200 kJ/mol force and after that, without restraints. LINCS algorithm restricted all the atom bond movements in the protein [41]. Then, the system was further minimized through conjugated gradient (CG) algorithm, without any restraint. Verlet integrator [42] was used to force the system to reach the equilibrium when the system was heated from 0 to 300 K at a constant pressure through short MD simulation (1 ns). MD simulation procedures were performed within periodic boundary conditions assuring that the atoms would remain inside the simulation box. Long-range bonds (at maximum 13 Å) of each atom were taken into account, and the electrostatic interactions were treated through the particle mesh ewald (PME) method [43]. An extensive MD simulation (100 ns) was conducted under the aforementioned conditions without any restraint.

The choice was made for a representative structure to be used in the structure-based virtual screening procedures (docking). It was done by analyzing frames from every 10 ps of simulation after the system reached stability and by subsequently selecting the central frame of the most populated cluster of system coordinates according to the RMSD values at cutoff 2 Å of RMSD for cluster differentiation based on the GROMOS method.

2.5. Ligand-Based Virtual Screening

Ligands required for ligand-based screening were selected from crystallographic holo structures (protein + ligand) of mosquito OBPs, because of their evolutionary and structural link to phlebotominae OBPs, since there is no phlebotomine OBP structure solved through experimental methods such as crystallography and NMR (nuclear magnetic resonance) deposited in the protein data bank. The biogenic (natural) chemical compound data bank of ZINC12 [44] was filtered according to its structural similarity to previously selected ligands. Therefore, similarity was based on the spatial chemical coordinates (from axes PC1, PC2, and PC3) collected in the CHEMGPS-NP Web [45] through Euclidean distance calculations. The spatial location of biogenic chemical structures was compared to the location of mosquito OBP ligands, but only the closest cluster to the OBP ligand structures was selected for the second screening filter known as molecular docking (receptor-based virtual screening).

2.6. Structure-Based Virtual Screening

Lutzomyia longipalpis OBP4 was the protein structure used in docking procedures; it was obtained through comparative modeling and optimized through MD calculations. The adopted chemical structure database was the one derived from a previously ligand-based filtered biogenic set of compounds. The protein orthosteric site was predicted through the superpositioning of template and homologue modeled structures in the Discovery Studio 4.5 [32]. The position of the backbone ligand complexed in the template structure was used to define a central point. Thus, the docking site was composed of amino acid residues located 8 Å from the central point. The protein structure was prepared in the UFSC CHIMERA 1.9.8 [30]. First, hydrogen atoms and the Gasteiger charge were added to the PDB protein file. A copy of this charged file (without the hydrogens) was used to find the molecular surface (dms file), which was used for sphere generation in all the protein cavities. The spheres within 8.0 Å of the ligand complexed structure were selected, and the docking calculation regions were delimited. Subsequently, grid generation was conducted by using the charged protein file and the selected spheres. The docking calculations were performed in the DOCK 6.7 package [46].

The binding energy of all docked compounds was scored through the Grid Score. The top virtual hit compounds were these compounds complex interactions calculated in the PLIP server (, which is a fully automated protein-ligand interaction profiler [47]. PLIP generated 3D interaction maps depicting the OBP orthosteric site and the docked ligand structures presenting the best energy values.

2.7. MD Simulation for Complex

The top virtual hit compound classified after hierarchical virtual screening procedures was subjected to assess molecular dynamics studies with the Lutzomyia longipalpis OBP model in silico under the same conditions adopted to previous MD studies implemented for OBP optimization. All the complex interactions were analyzed by using 3D interaction maps generated in PLIP [47].

3. Results and Discussion

The local alignment (BLAST) of target sequences without signal peptide amino acids allowed choosing the Culex quinquefasciatus OBP1 holo structure complexed with an oviposition pheromone crystalized at 8.2 pH, (PDB ID 3OGN-resolution: 1,3 Å; R-Value Free: 0.174) as the best template for homology modeling. The alignment had 61% of overall sequence identity and 57% orthosteric site conservation at 98% alignment coverage and 3e − 49 e-value (Figure 1). These identity values (overall and binding site) corroborate the extremely well conserved secondary structure folding the OBP family; on the other hand, they highlight the specificity of the olfactory chemical sense modulation. Sequence identity and similarity levels in the target-template alignment were satisfactory for homology prediction, since it only presented a six-residue region poorly conserved in the alignment (in red) and one noncorresponding additional residue at the beginning of the target sequence (Figure 1). The OBP target had the cysteine signature of the protein family, as well as presented the same typical C-terminal ionizable residues (HIS110, ASP117, and HIS120) found in OBPs that capture and release ligands according to pH medium variations [48] (Figure 1).

The different strategies used to generate comparative models provided three 3D OBP models. Energetic parameters of the I-Tasser model presented the best overall value depending on the quality of the taken measurements; however, it was the model recording less residues in favorable Ramachandran regions. The MODELLER model with the lowest % of high energy had the best Ramachandran result at one steric clash in a flexible zone of the structure away from the orthosteric site. The model deriving from the SWISS-MODEL had energy and steric values raging between the values recorded for the other two mentioned strategies; therefore, the model resulting from MODELLER due to the adoption of the restriction method allowed the best steric evaluation results: 96.4% of the residues in favorable regions of the Ramachandran plot. The global MODELLER measurement was also satisfactory, this model recorded normalized discrete optimized protein energy (zDOPE) score (−1.94), and it indicated that the atom pair distances in the model comply the distances recorded for a large sample of known protein structures; therefore, it was the model of choice (Table 1).

Comparative modeling software% aminoacides in Ramachandran plot regionsET (E/kT)% aa (high energy)Qmean6Z-score
FavorableAllowedGen. allowedNonallowed


The structural alignment between the homology model and the template (3OGN) OBPs showed C-α root mean square deviation (RMSD) 0.179 Å, which was within the range of experimental models deriving from the same protein (0.5 Å), and it indicated the high quality of the early-built homology model [49]. The built model presented all the common characteristics of insect odorant binding proteins, namely, globular shape formed by secondary alpha-helices structures presenting a hydrophobic cavity and a tail at the c-terminal chain. This tail helps to enclosure the ligand in the cavity (Figure 2).

The protonation prediction of OBP4 amino acid residues indicates that all acid residues and histidines were deprotonated at pH = 7.5 and that the basic residues were protonated. It is worth noticing that both tail termini histidines (HIS120 and HIS110) could be protonated in more acidic environments (at pH 6.5), and it reinforced the suggested mechanism of releasing odor through pH dropping close to the olfactory neuron membrane.

3.1. Model Refinement

The system solvated at the bound state (holo) was subjected to minimization runs and to 150 ns of MD simulation, but it did not require any charge to be neutralized. The RMSD trajectory of all 150 ns recorded standard deviation (4.9). The last 50 ns reached stability path at standard deviation 1.52 and mean value 3.7 Å (Figure 3).

The overall structure had good packing at the mean radius of gyration value 13.9 Å (standard deviation 1.11) when the stabilized period was taken into account, and it met the values found for others OBPs [50] (Figure 4). The overall energy had inconspicuous change, fact that confirmed a more stable state for the model. Furthermore, the residues presenting prominent RMS fluctuation (RMSF larger than 2 Å) were mainly located in flexible regions such as the coils and alpha-helices extremities, however, only in residues close to the cavity of the binding site (alpha-helices 5 and 6) (Figure 5).

Structure clustering within the trajectory was achieved according to the RMSD variation in the last 50 ns in 11 clusters. The largest cluster comprised 2332 structure frames from which the representative structure at 139.09 ns was retrieved. The optimized protein model selected for virtual screening showed satisfactory quality, it recorded 99.1% of residues in favored regions of the Ramachandran plot and only 0.9% (one residue) of residues in disallowed regions (Figure 6).

3.2. Ligand-Based Virtual Screening

A tiny fraction of all 120 million compounds available for virtual Screening [51] was selected, namely, 119,624 molecules of natural compounds from the biogenic zinc database [44]. Roughly 60,000 of these compounds were very close to the chemical space of typical OBP ligands (decanal, indol, PEG6, mosquito pheromone oviposition (MOP), and picaridin); therefore, they composed a database for the next step: the receptor-based virtual screening (Figure 7).

3.3. Structure-Based Virtual Screening

The docking results of ligands in the orthosteric OBP site were scored through the Grid Score and ranked by number; the top five compounds are shown in Table 2. All molecules were lysophospholipids and phospholipids missing one of their two O-acyl chains; they are commonly found in human metabolome since they are used for plant communication.

Top-ranked 1 (ZINC ID: 000040165350 PubChem annotation ID: 53480942 ) grid score energy: −98.73 KJ/mol
Top-ranked 2 (ZINC ID: 000040165373 PubChem annotation ID: 52925143 ) grid score energy: −97.82 KJ/mol
Top-ranked 3 (ZINC ID: 000040164463 PubChem annotation ID: 5283562; 11454491; 194606 ) grid score energy: −97.42 KJ/mol
Top-ranked 4 (ZINC ID: 000040165351 PubChem annotation ID: 53480943 ) grid score energy: −96.64 KJ/mol
Top-ranked 5 (ZINC ID: 000040165341 PubChem annotation ID: 53480932) grid score energy: −96.21 KJ/mol

Z40165350 (Table 2) showed hydrophobic interactions with the amino acids PHE72, TYR76, ILE79, ALA87, LYS88, and ASP92 due to its best scored energy: two hydrogen bonds, a donor h-bond between the HIS110 residue and the N(3) amino terminal of the ligand, and another h-bond donor between the amine of the PHE122 core chain and the ligand phosphate O(3). Additionally, two donor pi-type interactions with residues PHE61 and TRP113 (visualized in Discovery Studio) were observed. Despite the aforementioned residues, PHE61, LEU75, TRP113, HIS120, and PHE122 completed the hydrophobic interactions in other top-ranked ligand complexes and in an unusual complex. The ILE18, LEU14, and GLU13 residues set hydrophobic interactions with a longer fatty acid chain ligand (z40165351). The same residues (HIS120 and PHE120) could interact with h-donor ligand tails in h-bonds. Moreover, the terminal loop residue ALA124 was subtly facing the hydrophobic cavity in contrast to the 3OGN template, and it could allow greater anchoring for longer chain ligands.

A stabilized MD trajectory was reached when the Lutzomyia longipalpis OBP complex presenting top virtual hit compound (Z40165350) was evaluated (∼1.37 Å) (Figure 8(a)). A representative structure, which evidenced that the hydrogen bonds previously identified through molecular docking did not stayed in place, was chosen from this trajectory. However, the h-bond between the nitrogen of PHE122 and the oxygen of ligand appeared to be more consistent (Figure 8(b)). A salt bridge was also detected between HIS110 and the ligand phosphate group in addition to hydrophobic interactions between residues TYR76, ALA87, MET91, TRP113, and PHE122 and the aliphatic chain of ligand.

4. Conclusions

Understanding the molecular mechanisms responsible for the olfaction modulation in insects is a complex and long path capable of revolutionizing the manner we deal with pollinators, pests, and disease vectors in general. The effect of neglect and technological backwardness on the control of disease vectors has been devastating in recent years. Vector control is extremely difficult to perform when it comes to leishmaniasis; therefore, strategies focused on the olfactory modulation of the vector’s behavior could help developing new tools to prevent this disease.

The use of odor binding proteins (OBPs) as important molecular target in studies focused on vectors’ olfactory modulation has been gaining more space in chemical ecology. Some of these proteins have already been identified for Lutzomyia longipalpis; however, the three-dimensional structure of these molecules has never been investigated, and there is no information about their expression or interaction with ligands. Hence, this is the first study showing detailed aspects of the structural biology of L. longipalpis OBPs through comparative modeling techniques. The validation of the quality of the built models allowed inferring that MODELLER is very useful to predict 3D models of high-quality OBPs and made it possible prioritizing the best refinement model (OBP4), including the use of molecular dynamics simulations. The final structure, which is the most stable one and corresponds to the production of the MD trajectory, is expected to be the closest to the molecular bioactive conformation of OBP.

The use of the hierarchical approach applied to virtual screening has allowed systematically, efficiently, and roughly filtering 120,000 natural compounds by selecting compounds in the chemical space similar to the known OBP modulators, as well as facilitated the selection of relevant ligands (top-ranked in receptor-based VS) for future evaluation in vitro or in vivo.

Data Availability

Data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


The authors thank for academic support of Graduate Program in Biotechnology, Universidade Estadual de Feira de Santana, which allowed the development of “Comparative Modeling and Hierarchical Virtual Screening for Identification of OBPs Modulators of Lutzomyia Longipalpis” (dissertation) that includes the data shown on this manuscript.


  1. WHO, A Global Brief on Vector-Borne Diseases, World Health Organization, Geneva, Switzerland, 2014.
  2. P. D. Ready, “Biology of phlebotomine sand flies as vectors of disease agents,” Annual Review of Entomology, vol. 58, no. 1, pp. 227–250, 2013. View at: Publisher Site | Google Scholar
  3. P. D. Ready, “Epidemiology of visceral leishmaniasis,” Clinical Epidemiology, vol. 6, pp. 147–154, 2014. View at: Publisher Site | Google Scholar
  4. WHO, Control of The Leishmaniases, vol. 949, World Health Organization, Geneva, Switzerland, 2010, Technical Report series.
  5. WHO, The Use of DDT in Malaria Vector Control, World Health Organization, Geneva, Switzerland, 2011.
  6. S. M. Snedeker, “Pesticides and breast cancer risk: a review of DDT, DDE and dieldrin,” Environmental Health Perspectives, vol. 109, no. 1, pp. 35–47, 2001. View at: Publisher Site | Google Scholar
  7. W. Stuetz, “Global surveillance of DDT and DDE levels in human tissues,” International Journal of Occupational Medicine and Environmental Health, vol. 19, no. 1, p. 83, 2006. View at: Google Scholar
  8. H. Van Den Berg, “DDT and malaria prevention: Van Den Berg responds,” Environmental Health Perspectives, vol. 118, pp. A15–A16, 2010. View at: Publisher Site | Google Scholar
  9. G. M. Leal and W. S. Leal, “Binding of a fluorescence reporter and a ligand to an odorant-binding protein of the yellow fever mosquito, Aedes aegypti,” F1000Research, vol. 3, p. 305, 2015. View at: Publisher Site | Google Scholar
  10. W. S. Leal, “Pheromone reception,” Topics in Current Chemistry, vol. 240, pp. 1–36, 2005. View at: Google Scholar
  11. W. S. Leal, “Odorant reception in insects: roles of receptors, binding proteins, and degrading enzymes,” Annual Review of Entomology, vol. 58, no. 1, pp. 373–391, 2013. View at: Publisher Site | Google Scholar
  12. R.G. Vogt and L.M. Riddiford, “Pheromone binding and inactivation by moth antennae,” Nature, vol. 293, no. 5828, pp. 161–163, 1981. View at: Publisher Site | Google Scholar
  13. P. Xu, R. Atkinson, D. N. M. Jones, and D. P. Smith, “Drosophila OBP LUSH is required for activity of pheromone-sensitive neurons,” Neuron, vol. 45, no. 2, pp. 193–200, 2005. View at: Publisher Site | Google Scholar
  14. P. Pelosi, I. Iovinella, A. Felicioli, and F. R. Dani, “Soluble proteins of chemical communication: an overview across arthropods,” Frontiers in Physiology, vol. 5, p. 87, 2014. View at: Publisher Site | Google Scholar
  15. P. Pelosi, M. Calvello, and L. Ban, “Diversity of odorant-binding proteins and chemosensory proteins in insects,” Chemical Senses, vol. 30, no. 1, pp. i291–i292, 2005. View at: Publisher Site | Google Scholar
  16. P. D. Kamala, V. Kempraj, R. M. Aurade, T. K. Roy, K. S. Hivashankara, and A. Verghese, “Computational reverse chemical ecology: virtual screening and predicting behaviorally active semiochemicals for Bactrocera dorsalis,” BMC Genomics, vol. 15, no. 1, p. 209, 2014. View at: Publisher Site | Google Scholar
  17. W. S. Leal, R. M. R. Barbosa, W. Xu et al., “Reverse and conventional chemical ecology approaches for the development of oviposition attractants for Culex mosquitoes,” PLoS One, vol. 3, no. 8, Article ID e3045, 2008. View at: Publisher Site | Google Scholar
  18. R. P. Rodrigues, S. P. Mantoani, J. R. de Almeida et al., “Estratégias de triagem virtual no planejamento de fármacos,” Revista Virtual de Quimica, vol. 4, no. 6, pp. 739–776, 2012. View at: Google Scholar
  19. R. V. Azevedo, D. B. Dias, J. A. Bretãs et al., “The transcriptome of Lutzomyia longipalpis (Diptera: Psychodidae) male reproductive organs,” PLoS One, vol. 7, no. 4, Article ID e34495, 2012. View at: Publisher Site | Google Scholar
  20. D. B. S. Dias, Análise de Genes Expressos Nas Antenas de Lutzomyia Longipalpis (Lutz, Neiva, 1912) (Diptera: Psychodidae: Phlebotominae), Instituto Oswaldo Cruz, Rio de Janeiro, Brazil, 2008.
  21. I. B. Santana, “Modelagem comparativa e triagem virtual hierárquica para identificação de moduladores das OBPs de Lutzomyia longipalpis,” Universidade Estadual de Feira de Santana, Feira de Santana, Brazil, 2016, Dissertation. View at: Google Scholar
  22. T. N. Petersen, S. Brunak, G. von Heijne, and H. Nielsen, “SignalP 4.0: discriminating signal peptides from transmembrane regions,” Nature Methods, vol. 8, no. 10, pp. 785-786, 2011. View at: Publisher Site | Google Scholar
  23. S. Altschul, “Gapped BLAST and PSI-BLAST: a new generation of protein database search programs,” Nucleic Acids Research, vol. 25, no. 17, pp. 3389–3402, 1997. View at: Publisher Site | Google Scholar
  24. H. M. Berman, “The Protein Data Bank,” Nucleic Acids Research, vol. 28, no. 1, pp. 235–242, 2000. View at: Publisher Site | Google Scholar
  25. A. Fiser, “Template-based protein structure modeling,” Methods in Molecular Biology, vol. 673, pp. 73–94, 2010. View at: Publisher Site | Google Scholar
  26. S. Goldsmith-Fischman and B. Honig, “Structural genomics: computational methods for structure analysis,” Protein Science, vol. 12, no. 9, pp. 1813–1821, 2003. View at: Publisher Site | Google Scholar
  27. Z. Xiang, “Advances in homology protein structure modeling,” Current Protein and Peptide Science, vol. 7, no. 3, pp. 217–227, 2006. View at: Publisher Site | Google Scholar
  28. M. Biasini, S. Bienert, A. Waterhouse et al., “SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information,” Nucleic Acids Research, vol. 42, no. W1, pp. W252–W258, 2014. View at: Publisher Site | Google Scholar
  29. N. Eswar, B. Webb, M. A. Marti-Renom et al., “Comparative protein structure modeling using MODELLER,” Current Protocols in Bioinformatics, vol. 8, no. 5-6, 2006. View at: Publisher Site | Google Scholar
  30. E. F. Pettersen, T. D. Goddard, C. C. Huang et al., “UCSF Chimera? A visualization system for exploratory research and analysis,” Journal of Computational Chemistry, vol. 25, no. 13, pp. 1605–1612, 2004. View at: Publisher Site | Google Scholar
  31. Y. Zhang, “I-TASSER server for protein 3D structure prediction,” BMC Bioinformatics, vol. 9, no. 1, p. 40, 2008. View at: Publisher Site | Google Scholar
  32. D. S. Biovia, Discovery Studio Modeling Environment, San Diego Dassault Systèmes, San Diego, CA, USA, 2015.
  33. R. A. Laskowski, M. W. MacArthur, D. S. Moss, and J. M. Thornton, “PROCHECK: a program to check the stereochemical quality of protein structures,” Journal of Applied Crystallography, vol. 26, no. 2, pp. 283–291, 1993. View at: Publisher Site | Google Scholar
  34. F. Melo, D. Devos, E. Depiereux, and E. Feytmans, “ANOLEA: a www server to assess protein structures,” International Conference on Intelligent Systems for Molecular Biology, vol. 5, pp. 187–190, 1997. View at: Google Scholar
  35. P. Benkert, M. Künzli, and T. Schwede, “QMEAN server for protein model quality estimation,” Nucleic Acids Research, vol. 37, no. 2, pp. W510–W514, 2009. View at: Publisher Site | Google Scholar
  36. L. Zhang and J. Skolnick, “What should the Z-score of native protein structures be?” Protein Science, vol. 7, no. 5, pp. 1201–1207, 1998. View at: Publisher Site | Google Scholar
  37. D. Pal and P. Chakrabarti, “On residues in the disallowed region of the Ramachandran map,” Biopolymers, vol. 63, no. 3, pp. 195–206, 2002. View at: Publisher Site | Google Scholar
  38. C. Oostenbrink, A. Villa, A. E. Mark, and W. F. Van Gunsteren, “A biomolecular force field based on the free enthalpy of hydration and solvation: the GROMOS force-field parameter sets 53A5 and 53A6,” Journal of Computational Chemistry, vol. 25, no. 13, pp. 1656–1676, 2004. View at: Publisher Site | Google Scholar
  39. M. Rostkowski, M. H. M. Olsson, C. R. Søndergaard, and J. H. Jensen, “Graphical analysis of pH-dependent properties of proteins predicted using PROPKA,” BMC Structural Biology, vol. 11, no. 1, p. 6, 2011. View at: Publisher Site | Google Scholar
  40. H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, “The missing term in effective pair potentials,” Journal of Physical Chemistry, vol. 91, no. 24, pp. 6269–6271, 1987. View at: Publisher Site | Google Scholar
  41. B. Hess, H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije, “LINCS: A linear constraint solver for molecular simulations,” Journal of Computational Chemistry, vol. 18, no. 12, pp. 1463–1472, 1997. View at: Publisher Site | Google Scholar
  42. L. Verlet, “Computer “experiments” on classical fluids. I. Thermodynamical properties of Lennard-Jones molecules,” Physical Review, vol. 159, no. 1, pp. 98–103, 1967. View at: Publisher Site | Google Scholar
  43. T. Darden, D. York, and L. Pedersen, “Particle mesh Ewald: an N·log(N) method for Ewald sums in large systems,” Journal of Chemical Physics, vol. 98, no. 12, Article ID 10089, 1993. View at: Publisher Site | Google Scholar
  44. J. J. Irwin, T. Sterling, M. M. Mysinger, E. S. Bolstad, and R. G. Coleman, “ZINC: a free tool to discover chemistry for biology,” Journal of Chemical Information and Modeling, vol. 52, no. 7, pp. 1757–1768, 2012. View at: Publisher Site | Google Scholar
  45. J. Larsson, J. Gottfries, S. Muresan, and A. Backlund, “ChemGPS-NP: tuned for navigation in biologically relevant chemical space,” Journal of Natural Products, vol. 70, no. 5, pp. 789–794, 2007. View at: Publisher Site | Google Scholar
  46. W. J. Allen, T. E. Balius, S. Mukherjee et al., “DOCK 6: impact of new features and current docking performance,” Journal of Computational Chemistry, vol. 36, no. 15, pp. 1132–1156, 2015. View at: Publisher Site | Google Scholar
  47. S. Salentin, S. Schreiber, V. J. Haupt et al., “PLIP: fully automated protein-ligand interaction profiler,” Nucleic Acids Research, vol. 43, pp. W443–W447, 2015. View at: Publisher Site | Google Scholar
  48. F. F. Damberger, E. Michel, Y. Ishida, W. S. Leal, and K. Wuthrich, “Pheromone discrimination by a pH-tuned polymorphism of the Bombyx mori pheromone-binding protein,” Proceedings of the National Academy of Sciences, vol. 110, no. 46, pp. 18680–18685, 2013. View at: Publisher Site | Google Scholar
  49. C. Chothia and A. M. Lesk, “The relation between the divergence of sequence and structure in proteins,” EMBO Journal, vol. 5, no. 4, pp. 823–826, 1986. View at: Google Scholar
  50. L. Han, Y. J. Zhang, L. Zhang et al., “Operating mechanism and molecular dynamics of pheromone-binding protein ASP1 as influenced by pH,” PLoS One, vol. 9, no. 10, Article ID e110565, 2014. View at: Publisher Site | Google Scholar
  51. T. Sterling and J. J. Irwin, “ZINC 15-ligand discovery for everyone,” Journal of Chemical Information and Modeling, vol. 55, no. 11, pp. 2324–2337, 2015. View at: Publisher Site | Google Scholar

Copyright © 2018 Isis Bugia Santana et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

705 Views | 288 Downloads | 1 Citation
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.