Estudo estrutural e comparativo entre Sirtuínas 2 de diferentes espécies de Leishmania / Structural and comparative study between Sirtuins 2 of different species of Leishmania by means of Comparative Modeling and Molecular Dynamic

As leishmanioses são doenças negligenciadas que afetam mais de um bilhão e meio de pessoas ao redor do mundo, principalmente nos países em desenvolvimento, provocando grandes impactos socioeconômicos. Os fármacos disponíveis para o tratamento dessas doenças são ineficazes e apresentam graves efeitos adversos. O processo de pesquisa de novos fármacos envolve, entre outras coisas, a seleção de alvos bioquímicos essenciais para a sobrevivência e desenvolvimento do agente causador. Neste sentido, a Sirtuína 2, uma enzima epigenética com atividade hidrolase essencial para a sobrevivência dos parasitas do gênero Leishmania se apresenta como um alvo validado na busca de novos fármacos contra essas parasitoses. O planejamento de fármacos baseado na estrutura do receptor requer o conhecimento da estrutura tridimensional da proteína alvo. Desta forma, a elucidação estrutural e um estudo minucioso das Sirtuínas das várias espécies do gênero Leishmania apresenta-se como uma importante abordagem na aplicação desta estratégia na busca por agentes quimioterápicos. Até o momento, na família Trypanosomatidae, a única estrutura tridimensional resolvida experimentalmente de uma enzima Sirtuína 2 é a da espécie L. infantum. Assim, este trabalho aplicou a abordagem de Modelagem Comparativa utilizando o software Modeller na construção de modelos da Sir2rp1 das espécies L. infantum, L. major e L. braziliensis, cujas sequências de aminoácidos foram extraídas do banco de dados UNIProt. Os modelos construídos foram validados por meio da função de escore DOPE do Modeller e dos servidores PROCHECK, MolProbity e QMEAN, avaliando sua qualidade estereoquímica e seu enovelamento. Os ligantes naturais da enzima foram sobrepostos nos modelos construídos por alinhamento estrutural utilizando o software PyMol e os complexos validados foram submetidos a simulações de Dinâmica Molecular através do pacote GROMACS. Os complexos refinados foram então analisados por meio dos softwares PyMol e LigPlotPlus e dos pacotes GROMACS e gmx_MMPBSA, e foram estudados os sítios de ligação dos substratos e os resíduos de aminoácidos relevantes envolvidos em sua ligação e reconhecimento. A Modelagem Comparativa da Sirtuína 2 humana e seus homólogos das espécies L. infantum, L. major e L. braziliensis, as simulações de Dinâmica Molecular realizadas com os modelos enzimáticos construídos e validados complexados com seus ligantes naturais, os cálculos de energia de interação entre os modelos e seus substratos e o estudo estrutural comparativo realizado entre eles nos fornecem uma base teórica para a busca de novos inibidores da Sirtuína 2 que sejam mais seletivos e potentes contra as enzimas parasitárias, abrindo caminho para o desenvolvimento de candidatos a fármacos leishmanicidas mais seguros e eficazes

Leishmaniasis are neglected diseases that affect more than one and a half billion people around the world, mainly in developing countries, causing major socioeconomic impacts. The drugs available for the treatment of these diseases are ineffective and have serious adverse effects. The process of researching new drugs involves, among other things, the selection of biochemical targets essential for the survival and development of the causative agent. In this sense, Sirtuin 2, an epigenetic enzyme with hydrolase activity essential for the survival of parasites of the Leishmania genus, presents itself as a validated target in the search for new drugs against these parasites. Structure-Based Drug Design requires knowledge of the three-dimensional structure of the target protein. In this way, structural elucidation and a detailed study of Sirtuins from various species of the genus Leishmania presents itself as an important approach in the application of this strategy in the search for chemotherapeutic agents. To date, in the Trypanosomatidae family, the only experimentally resolved three-dimensional structure of a Sirtuin 2 enzyme is that of the species L. infantum. Thus, this work applied the Comparative Modeling approach using the Modeller software in the construction of Sir2rp1 models of the species L. infantum, L. major and L. braziliensis, whose amino acid sequences were retrieved from the UNIProt database. The constructed models were validated using Modeller's DOPE score function and the PROCHECK, MolProbity and QMEAN servers, evaluating their stereochemical quality and folding. The enzyme's natural ligands were superimposed on the built models by structural alignment using the PyMol software and the validated complexes were subjected to Molecular Dynamics simulations using the GROMACS package. The refined complexes were then analyzed using the PyMol and LigPlotPlus softwares and the GROMACS and gmx_MMPBSA packages, and the substrate binding sites and relevant amino acid residues involved in their binding and recognition were studied. The Comparative Modeling of human Sirtuin 2 and its homologues from the species L. infantum, L. major and L. braziliensis, the Molecular Dynamics simulations carried out with the constructed and validated enzymatic models complexed with their natural ligands, the interaction energy calculations between the models and their substrates and the comparative structural study carried out between them provide us with a theoretical basis for the search for new Sirtuin 2 inhibitors that are more selective and potent against the parasitic enzymes, paving the way for the development of safer and more effective leishmanicidal drug candidates

Development, validation and analysis of a human profurin 3D model using comparative modeling and molecular dynamics simulations.

The emergence of new viruses can lead to the outbreak of pandemics as occurred at the end of 2019 with the coronavirus disease (or COVID-19). The fastest way to effectively control viral infections is to develop broad-spectrum antivirals that can fight at least an entire class of viruses. Profurin, the furin precursor propeptide, is responsible for the autoactivation step which is crucial for the maturation of several viral substrates. This role makes the study of furin and profurin interactions interesting for the development of new potential broad-spectrum antivirals for the treatment against several human viral diseases. Since there is no 3D model of profurin published in the literature or deposited in a database, this work reports the development, validation and analysis of a profurin 3D model using comparative modeling and molecular dynamics. The model is available in ModelArchive at https://www.modelarchive.org/doi/10.5452/ma-ct8l7. The usage of this model will make possible further studies of molecular docking and MD simulations of the profurin-furin system, in the design of new potential broad-spectrum antivirals for the treatment against several human viral diseases.Communicated by Ramaswamy H. Sarma.

In Silico Analysis of a GH3 ß-Glucosidase from Microcystis aeruginosa CACIAM 03.

Cyanobacteria are rich sources of secondary metabolites and have the potential to be excellent industrial enzyme producers. ß-glucosidases are extensively employed in processing biomass degradation as they mediate the most crucial step of bioconversion of cellobiose (CBI), hence controlling the efficiency and global rate of biomass hydrolysis. However, the production and availability of these enzymes derived from cyanobacteria remains limited. In this study, we evaluated the ß-glucosidase from Microcystis aeruginosa CACIAM 03 (MaBgl3) and its potential for bioconversion of cellulosic biomass by analyzing primary/secondary structures, predicting physicochemical properties, homology modeling, molecular docking, and simulations of molecular dynamics (MD). The results showed that MaBgl3 derives from an N-terminal domain folded as a distorted ß-barrel, which contains the conserved His-Asp catalytic dyad often found in glycosylases of the GH3 family. The molecular docking results showed relevant interactions with Asp81, Ala271 and Arg444 residues that contribute to the binding process during MD simulation. Moreover, the MD simulation of the MaBgl3 was stable, shown by analyzing the root mean square deviation (RMSD) values and observing favorable binding free energy in both complexes. In addition, experimental data suggest that MaBgl3 could be a potential enzyme for cellobiose-hydrolyzing degradation.

Differences in Charge Distribution in Leishmania tarentolae Leishmanolysin Result in a Reduced Enzymatic Activity.

Leishmania tarentolae is a non-pathogenic trypanosomatid isolated from lizards widely used for heterologous protein expression and extensively studied to understand the pathogenic mechanisms of leishmaniasis. The repertoire of leishmanolysin genes was reported to be expanded in L. tarentolae genome, but no proteolytic activity was detected. Here, we analyzed L. tarentolae leishmanolysin proteins from the genome to the structural levels and evaluated the enzymatic activity of the wild-type and overexpressing mutants of leishmanolysin. A total of 61 leishmanolysin sequences were retrieved from the L. tarentolae genome. Five of them were selected for phylogenetic analysis, and for three of them, we built 3D models based on the crystallographic structure of L. major ortholog. Molecular dynamics simulations of these models disclosed a less negative electrostatic potential compared to the template. Subsequently, L. major LmjF.10.0460 and L. tarentolae LtaP10.0650 leishmanolysins were cloned in a pLEXSY expression system into L. tarentolae. Proteins from the wild-type and the overexpressing parasites were submitted to enzymatic analysis. Our results revealed that L. tarentolae leishmanolysins harbor a weak enzymatic activity about three times less abundant than L. major leishmanolysin. Our findings strongly suggest that the less negative electrostatic potential of L. tarentolae leishmanolysin can be the reason for the reduced proteolytic activity detected in this parasite.

Functional annotation and comparative modeling of ligninolytic enzymes from Trametes villosa (SW.) Kreisel for biotechnological applications.

Functional annotation of Trametes villosa genome was performed to search Class II peroxidase proteins in this white-rot fungus, which can be valuable for several biotechnological processes. After sequence identification and manual curation, five proteins were selected to build 3 D models by comparative modeling. Analysis of sequential and structural sequences from selected targets revealed the presence of two putative Lignin Peroxidase and three putative Manganese Peroxidase on this fungal genome. All 3 D models had a similar folding pattern from selected 3 D structure templates. After minimization and validation steps, the best 3 D models were subjected to docking studies and molecular dynamics to identify structural requirements and the interactions required for molecular recognition. Two reliable 3 D models of Class II peroxidases, with typical catalytic site and architecture, and its protein sequences are indicated to recombinant production in biotechnological applications, such as bioenergy.Communicated by Ramaswamy H. Sarma.

Structural insights into NS5B protein of novel equine hepaciviruses and pegiviruses complexed with polymerase inhibitors.

Infections produced by hepaciviruses have been associated with liver disease in horses. Currently, at least three viruses belonging to the Flaviviridae family are capable of producing a chronic infection in equines: non-primate hepacivirus (NPHV), Theiler's disease-associated virus (TDAV), and equine pegivirus (EPgV). The RNA-dependent RNA polymerases of viruses (RdRp) (NS5 protein), from the flavivirus family, use de novo RNA synthesis to initiate synthesis. The two antiviral drugs currently used to treat hepatitis C (HCV), sofosbuvir and dasabuvir, act on the viral NS5B polymerase as nucleoside and non-nucleoside inhibitors, respectively. Both drugs have shown significant clinical inhibition of viral response. In this work, we aimed to model the NS5B polymerase of the equine hepacivirus (EHCV) subtypes 1 and 2, TDAV and EPgV, to assess whether current direct-acting antiviral drugs against HCV interact with these proteins. Crystal structures of HCV-NS5B were used as templates for modeling target sequences in both conformations (open and closed). Also, molecular docking of sofosbuvir and dasabuvir were performed to predict their possible binding modes at the modeled NS5B polymerase binding sites. We observed that the NS5B models of the EHCV and EPgV shared well-conserved 3D structures to HCV-NS5B and other RdRps, suggesting functional conservation. Interactions of EHCV subtypes 1, 2 and TDAV polymerases with sofosbuvir showed a similar molecular interaction pattern compared to HCV-NS5B, while interactions with dasabuvir were less conserved. In silico studies of molecular interactions between these modeled structures and sofosbuvir suggest that this compound could be efficient in combating equine pathogens, thus contributing to animal welfare.

Experimental Approaches and Computational Modeling of Rat Serum Albumin and Its Interaction with Piperine.

The bioactive piperine (1-piperoyl piperidine) compound found in some pepper species (Piper nigrum linn and Piper sarmentosum Roxb) has been shown to have therapeutic properties and to be useful for well-being. The tests used to validate these properties were performed in vitro or with small rats. However, in all these assays, the molecular approach was absent. Although the first therapeutic trials relied on the use of rats, no proposal was mentioned either experimentally or computationally at the molecular level regarding the interaction between piperine and rat serum albumin (RSA). In the present study, several spectroscopic techniques were employed to characterize rat serum albumin and, aided by computational techniques, the protein modeling was proposed. From the spectroscopic results, it was possible to estimate the binding constant (3.9 × 104 M-1 at 288 K) using the Stern-Volmer model and the number of ligands (three) associated with the protein applying interaction density function model. The Gibbs free energy, an important thermodynamic parameter, was determined (-25 kJ/mol), indicating that the interaction was spontaneous. This important set of experimental results served to parameterize the computational simulations. The results of molecular docking and molecular dynamics matched appropriately made it possible to have detailed microenvironments of RSA accessed by piperine.

GABAA receptor family: overview on structural characterization.

The pentameric γ-aminobutyric acid type A receptors are ion channels activated by ligands, which intervene in the rapid inhibitory transmission in the mammalian CNS. Due to their rich pharmacology and therapeutic potential, it is essential to understand their structure and function thoroughly. This deep characterization was hampered by the lack of experimental structural information for many years. Thus, computational techniques have been extensively combined with experimental data, in order to undertake the study of γ-aminobutyric acid type A receptors and their interaction with drugs. Here, we review the exciting journey made to assess the structures of these receptors and outline major outcomes. Finally, we discuss the brand new structure of the α1ß2γ2 subtype and the amazing advances it brings to the field.

In silico characterization of a cyanobacterial plant-type isoaspartyl aminopeptidase/asparaginase.

Asparaginases are found in a range of organisms, although those found in cyanobacteria have been little studied, in spite of their great potential for biotechnological application. This study therefore sought to characterize the molecular structure of an L-asparaginase from the cyanobacterium Limnothrix sp. CACIAM 69d, which was isolated from a freshwater Amazonian environment. After homology modeling, model validation was performed using a Ramachandran plot, VERIFY3D, and the RMSD. We also performed molecular docking and dynamics simulations based on binding free-energy analysis. Structural alignment revealed homology with the isoaspartyl peptidase/asparaginase (EcAIII) from Escherichia coli. When compared to the template, our model showed full conservation of the catalytic site. In silico simulations confirmed the interaction of cyanobacterial isoaspartyl peptidase/asparaginase with its substrate, ß-Asp-Leu dipeptide. We also observed that the residues Thr154, Thr187, Gly207, Asp218, and Gly237 were fundamental to protein-ligand complexation. Overall, our results suggest that L-asparaginase from Limnothrix sp. CACIAM 669d has similar properties to E. coli EcAIII asparaginase. Our study opens up new perspectives for the biotechnological exploitation of cyanobacterial asparaginases.

Structural insights into leishmanolysins encoded on chromosome 10 of Leishmania (Viannia) braziliensis

BACKGROUND Leishmanolysins have been described as important parasite virulence factors because of their roles in the infection of promastigotes and resistance to host's defenses. Leishmania (Viannia) braziliensis contains several leishmanolysin genes in its genome, especially in chromosome 10. However, the functional impact of such diversity is not understood, but may be attributed partially to the lack of structural data for proteins from this parasite. OBJECTIVES This works aims to compare leishmanolysin sequences from L. (V.) braziliensis and to understand how the diversity impacts in their structural and dynamic features. METHODS Leishmanolysin sequences were retrieved from GeneDB. Subsequently, 3D models were built using comparative modeling methods and their dynamical behavior was studied using molecular dynamic simulations. FINDINGS We identified three subgroups of leishmanolysins according to sequence variations. These differences directly affect the electrostatic properties of leishmanolysins and the geometry of their active sites. We identified two levels of structural heterogeneity that might be related to the ability of promastigotes to interact with a broad range of substrates. MAIN CONCLUSION Altogether, the structural plasticity of leishmanolysins may constitute an important evolutionary adaptation rarely explored when considering the virulence of L. (V.) braziliensis parasites.

Comparative modeling and molecular dynamics suggest high carboxylase activity of the Cyanobium sp. CACIAM14 RbcL protein.

Rubisco catalyzes the first step reaction in the carbon fixation pathway, bonding atmospheric CO2/O2 to ribulose 1,5-bisphosphate; it is therefore considered one of the most important enzymes in the biosphere. Genetic modifications to increase the carboxylase activity of rubisco are a subject of great interest to agronomy and biotechnology, since this could increase the productivity of biomass in plants, algae and cyanobacteria and give better yields in crops and biofuel production. Thus, the aim of this study was to characterize in silico the catalytic domain of the rubisco large subunit (rbcL gene) of Cyanobium sp. CACIAM14, and identify target sites to improve enzyme affinity for ribulose 1,5-bisphosphate. A three-dimensional model was built using MODELLER 9.14, molecular dynamics was used to generate a 100 ns trajectory by AMBER12, and the binding free energy was calculated using MM-PBSA, MM-GBSA and SIE methods with alanine scanning. The model obtained showed characteristics of form-I rubisco, with 15 beta sheets and 19 alpha helices, and maintained the highly conserved catalytic site encompassing residues Lys175, Lys177, Lys201, Asp203, and Glu204. The binding free energy of the enzyme-substrate complexation of Cyanobium sp. CACIAM14 showed values around -10 kcal mol(-1) using the SIE method. The most important residues for the interaction with ribulose 1,5-bisphosphate were Arg295 followed by Lys334. The generated model was successfully validated, remaining stable during the whole simulation, and demonstrated characteristics of enzymes with high carboxylase activity. The binding analysis revealed candidates for directed mutagenesis sites to improve rubisco's affinity.

Computational Studies of Benzoxazinone Derivatives as Antiviral Agents against Herpes Virus Type 1 Protease.

Herpes simplex virus infections have been described in the medical literature for centuries, yet the the drugs available nowadays for therapy are largely ineffective and low oral bioavailability plays an important role on the inefficacy of the treatments. Additionally, the details of the inhibition of Herpes Virus type 1 are still not fully understood. Studies have shown that several viruses encode one or more proteases required for the production new infectious virions. This study presents an analysis of the interactions between HSV-1 protease and benzoxazinone derivatives through a combination of structure-activity relationships, comparative modeling and molecular docking studies. The structure activity relationship results showed an important contribution of hydrophobic and polarizable groups and limitations for bulky groups in specific positions. Two Herpes Virus type 1 protease models were constructed and compared to achieve the best model which was obtained by MODELLER. Molecular docking results pointed to an important interaction between the most potent benzoxazinone derivative and Ser129, consistent with previous mechanistic data. Moreover, we also observed hydrophobic interactions that may play an important role in the stabilization of inhibitors in the active site. Finally, we performed druglikeness and drugscore studies of the most potent derivatives and the drugs currently used against Herpes virus.

3-D modeling and molecular dynamics simulation of interleukin-22 from the So-iny mullet, Liza haematocheila

Background: Interleukin-22 (IL-22) plays an important role in the regulation of immune responses. However, little is known about its function or structure in fish. Results: The IL-22 gene was first cloned from So-iny mullet (Liza haematocheila), one of commercially important fish species in China. Then, 3-D structure model of the mullet IL-22 was constructed by comparative modeling method using human IL-22 (1M4R) as template, and a 5 ns molecular dynamics (MD) was studied. The open reading frame (ORF) of mullet IL-22 cDNA was 555 bp, encoding 184 amino acids. The mullet IL-22 shared higher identities with the other fish IL-22 homologs and possessed a conserved IL-10 signature motif at its C-terminal. The mullet IL-22 model possessed six conserved helix structure. PROCHECK, SAVES and Molprobity server analysis confirmed that this model threaded well with human IL-22. Strikingly, analysis with CastP, cons-PPISP server suggested that the cysteines in mullet IL-22 might not be involved in the forming of disulfide bond for structural stabilization, but related to protein-protein interactions. Conclusions: The structure of IL-22 in So-iny mullet (Liza haematocheila) was constructed using comparative modeling method which provide more information for studying the function of fish IL-22.

Construção do receptor CYSLTR1 por modelagem comparativa / Building of CYSLTR1 receptor by comparative modeling

A asma é caracterizada como um distúrbio inflamatório crônico das vias aéreas, provocada pela contração da musculatura lisa dos bronquíolos, ocasionando obstrução parcial dos mesmos e dificultando a respiração. A Organização Mundial de Saúde (OMS) estima que 300 milhões de pessoas atualmente sofram de asma, sendo que as crianças estão entre as mais acometidas. Somente em 2005, 255.000 pessoas morreram de asma. Dentre os mediadores envolvidos no processo asmático, os leucotrienos cisteínicos, derivados do ácido araquidônico, são considerados os mais potentes entre aqueles envolvidos no processo asmático, indicados como principais mediadores da inflamação reversível das vias aéreas. Nos últimos 20 anos, grandes esforços têm sido realizados para identificar e desenvolver antagonistas dos receptores de leucotrienos na busca de melhorar o tratamento da asma, limitar a sua morbidade, e reduzir os efeitos dos medicamentos atuais. Portanto, o presente estudo propôs a construção de um modelo teórico do receptor do leucotrienos cisteínicos denominado CysLT1 por Modelagem Comparativa. Conclui-se que o modelo obtido através das metodologias computacionais e apresentado no presente estudo pode auxiliar em futuros testes, principalmente em metodologias que empregam ancoragem molecular e de novo design testando, in silico, ligantes de diversas fontes contra o receptor CysLT1.

Asthma is a chronic inflammatory disorder of the airways characterized by contraction of the smooth muscle of the bronchioles, causing their partial obstruction and making it difficult to breathe. The World Health Organization (WHO) estimates that 300 million people currently suffer from asthma, which is more common among children. In 2005 alone, 255,000 people died of asthma. The cysteinyl leukotrienes, derived from arachidonic acid, are considered the most potent mediators of the asthmatic process, being indicated as key mediators of reversible airway inflammation. In the past 20 years, great efforts have been made to identify and develop leukotriene receptor antagonists, in the quest to improve the treatment of asthma, limit its morbidity and reduce the side-effects of current drugs. Therefore, the objective of this study was to build a theoretical model of the cysteinyl leukotriene receptor, CysLT1, by Comparative Modeling. We conclude that the model generated by computational methods and presented in this paper may help in future studies, especially where docking and de novo design are involved, in which new ligands from diverse sources are tested in silico against the CysLT1 receptor.