Comparison of dopamine D2 receptor (homology model and X-ray structure) and virtual screening protocol validation for the antagonism mechanism

The present research work aims to compare the homology model and recent X-ray crystal structure of dopamineD2 receptor (PDB Code: 6CM4) as well as to validate the virtual screening protocol for antagonist compounds. Thecomparison involved the sequence similarity and the capability of both proteins to produce similar risperidone bindingpose with co-crystal structure based on ChemPLP score and Tanimoto Coefficient score generated by PLANTS andPyplif. Homology model failed to give the correct binding pose as the root mean square deviation fell to >2Å evenwith similar sequence and folding. Therefore, 6CM4 should be used for virtual screening instead of the homologymodel. The virtual screening protocol validation of 6CM4 was performed by PLANTS followed by Pyplif filtering.The protocol was able to give EF1% value of 6.238, which was better than the EF1% value of protein dopamine D3receptor that shared >80% similarity with dopamine D2 receptor. Similarity between the docking pose and the actualpose is considered important to obtain better predictivity

In silico analysis of the structural and functional implications of SLC19A1 R27H polymorphism

In view of the documented association of solute carrier family 19 member 1 (SLC19A1) G80A (R27H) polymorphism with the risk for different types of cancers and systemic lupus erythematosus (SLE), we have reanalysed the case–control study on breast cancer to ascertain the conditions in which this polymorphic variant exerts the risk of breast cancer. Association statistics have revealed that this polymorphism exerts the risk for breast cancer under the conditions of low folate intake, and in the absence of well-documented protective polymorphism in cytosolic serine hydroxymethyltransferase. To substantiate this observation, we have developed a homology model of SLC19A1 using glycerol-3-phosphate transporter (d1pw4a) as a template where 73% of the residues were modelled at 90% confidence while 162 residues were modelled ab initio. The wild and mutant proteins shared same topology in S3, S5, S6, S7, S11 and S12 transmembrane domains. The topology varied at S1 (28–43 residue vs 28–44 residue), S2 (66–87 residue vs 69–87 residue), S4 (117–140 residue vs 117–139 residue), S8 (305–325 residue vs 305–324 residue), S9 (336–356 residue vs 336–355residue), and S10 (361–386 residue vs 361–385 residue) transmembrane domains between wild versus mutant proteins. S2 domain is shortened by three amino acid residues in themutantwhile in other domains the difference corresponds to one amino acid residue. The 3DLigandSite analysis revealed that the metallic-ligand-binding sites at 273Trp, 277Asn, 379Leu, 439Phe and 442Leu are although unaffected, there is a loss of active sites corresponding to nonmetallic ligand binding. Tetrahydrofolate and methotrexate have lesseraffinity towards the mutant protein than the wild protein. To conclude, the R27H polymorphism affects the secondary and tertiary structures of SLC19A1 with the significant loss in ligand-binding sites.

Influence of Mg2+ ions on the interaction between 3,5-dicaffeoylquinic acid and HTLV-I integrase / Influencia de los iones Mg+2 sobre la interacción entre el ácido 3,5- dicafeoilquínico y la integrasa del HTLV-I / Influência dos íons Mg2+ sobre a interação entre o ácido 3,5-dicafeoilquínico e a integrase do HTLV-I

Objective. Using molecular simulation, we studied the influence of Mg2+ ions on the binding mode of HTLV-I Integrase (IN) catalytic domain (modeled by homology) with the 3,5- Dicaffeoylquinic Acid (DCQA). HTLV-I Integrase homology model was built using template-like crystallographic data of the IN catalytic domain solved for Avian Sarcoma Virus (VSA, pdb: 1VSD). Materials and methods. In order to analyze the role of Mg2+ in the interaction or coupling between 3,5-DCQA and Integrase, three models were created: i) in the absence of Mg2+ ions, ii) with a Mg2+ ion coordinated at Asp15 and Asp72 and iii) model with two Mg2+ ions coordinated at Asp15-Asp72 and Asp72-Glu108. Coupling force and binding free energy between 3,5-DCQA and HTLV-I IN were assessed in the three models. Results. The lowest docking score and free energy binding were obtained for the second model. Mg2+ ion strongly affected the coupling of the inhibitor 3,5-DCQA with HTLV-I catalytic domain of Integrase, thus revealing a strong interaction in the ligand-protein complex regardless of the ligand-catalytic interaction sites for all three models. Conclusion. Altogether, these results strengthen the hypothesis that the presence of one Mg2+ ion could enhance the interaction in the complex by decreasing free energy, therefore increasing the affinity. Moreover, we propose 3,5-DCQA as an important pharmacophore in the rational design of new antiretroviral drugs.

Objetivo: Usando simulación molecular, estudiamos la influencia de los iones Mg2+ en la interacción del dominio catalítico de la integrasa HTLV-I (IN) (modelado por homología) con el ácido 3,5-Dicafeoilquínico (DCQA). Materiales y métodos. El modelo por homología de la HTLV-I IN fue construido usando como molde la estructura cristalina de la IN del virus del sarcoma aviar (VSA, pdb: 1VSD). Para analizar el rol de los iones Mg2+ en la interacción con el DCQA y la integrasa, tres modelos fueron creados: i) en ausencia de iones Mg2+, ii) con un ion Mg2+ coordinado con Asp15 y Asp72 y iii) con dos iones Mg2+ coordinados con Asp15-Asp72 y Asp72-Glu108. Las fuerzas de interacción y la energía libre de unión entre el DCQA y la HTLV-I IN fueron calculadas en los tres modelos. Resultados. El puntaje más bajo en el docking y la menor energía libre fueron obtenidos para el modelo con un solo ion de Mg2+. El Mg2+ afecta fuertemente el acoplamiento del inhibidor DCQA con el dominio catalítico de la HTLV-I IN, revelando una fuerte interacción entre el complejo ligando-proteína que es independiente del sitio catalítico de los tres modelos usados. Conclusión. Los resultados sugieren que la presencia de un ion de Mg² + podría incrementar la interacción en el complejo debido a la disminución de la energía libre, intensificando así la afinidad. Por lo que, proponemos al DCQA como farmacóforo para el diseño de drogas antiretrovirales.

Objetivo. Usando simulação molecular, estudamos a influência dos íons Mg2+ na interação do domínio catalítico da integrase HTLV-I (IN) (modelado por homologia) com o ácido 3,5-dicafeoilquínico (3,5-diCQA). Materiais e métodos. O modelo de homologia da HTLV-I IN foi construído utilizando como fôrma a estrutura cristalina da IN de vírus do sarcoma aviário (VSA, pdb: 1VSD). Para analisar o papel dos íons Mg2+ na interação com o 3,5-diCQA e a integrase, três modelos foram criados: i) na ausência de íons Mg2+, ii) com um íon Mg2+ coordenado com Asp15 e Asp72 e, iii) com dois íons Mg2+ coordenados com Asp15-Asp72 e Asp72-Glu108. As forças de interação e a energia livre de ligação entre o 3,5-diCQA e a HTLV-I IN foram calculadas nos três modelos. Resultados. A pontuação mais baixa no docking e a menor energia livre foram obtidas para o modelo com um único íon de Mg2+. O Mg2+ afeta fortemente o acoplamento do inibidor 3,5-diCQA com o domínio catalítico da HTLV-I IN, revelando uma forte interação entre o complexo ligando-proteína que é independente do sítio catalítico dos três modelos utilizados. Conclusão. Os resultados sugerem que a presença de um íon Mg2+ poderia aumentar a interação no complexo devido à diminuição da energia livre, aumentando assim a afinidade. Assim, propomos ao 3,5-diCQA como farmacóforo para a fabricação de medicamentos antirretrovirais.

In silico analysis identifies a C3HC4-RING finger domain of a putative E3 ubiquitin-protein ligase located at the C-terminus of a polyglutamine-containing protein

Almost identical polyglutamine-containing proteins with unknown structures have been found in human, mouse and rat genomes (GenBank AJ277365, AF525300, AY879229). We infer that an identical new gene (RING) finger domain of real interest is located in each C-terminal segment. A three-dimensional (3-D) model was generated by remote homology modeling and the functional implications are discussed. The model consists of 65 residues from terminal position 707 to 772 of the human protein with a total length of 796 residues. The 3-D model predicts a ubiquitin-protein ligase (E3) as a binding site for ubiquitin-conjugating enzyme (E2). Both enzymes are part of the ubiquitin pathway to label unwanted proteins for subsequent enzymatic degradation. The molecular contact specificities are suggested for both the substrate recognition and the residues at the possible E2-binding surface. The predicted structure, of a ubiquitin-protein ligase (E3, enzyme class number 6.3.2.19, CATH code 3.30.40.10.4) may contribute to explain the process of ubiquitination. The 3-D model supports the idea of a C3HC4-RING finger with a partially new pattern. The putative E2-binding site is formed by a shallow hydrophobic groove on the surface adjacent to the helix and one zinc finger (L722, C739, P740, P741, R744). Solvent-exposed hydrophobic amino acids lie around both zinc fingers (I717, L722, F738, or P765, L766, V767, V733, P734). The 3-D structure was deposited in the protein databank theoretical model repository (2B9G, RCSB Protein Data Bank, NJ).