PACAP key interactions with PAC1, VPAC1, and VPAC2 identified by molecular dynamics simulations.

The neuropeptide pituitary adenylate cyclase-activating polypeptide (PACAP) belongs to the glucagon/secretin family. PACAP interacts with the pituitary adenylate cyclase-activating polypeptide receptor type 1 (PAC1) and vasoactive intestinal peptide receptors 1 and 2 (VPAC1 and VPAC2), exhibiting functions in the immune, endocrine, and nervous systems. This peptide is upregulated in numerous instances of brain injury, acting as a neuroprotective agent. It can also suppress HIV-1 and SARS-CoV-2 viral replication in vitro. This work aimed to identify, in each peptide-receptor system, the most relevant residues for complex stability and interaction energy communication via Molecular Dynamics (MD), Free Energy calculations, and Protein-energy networks, thus revealing in detail the underlying mechanisms of activation of these receptors. Hydrogen bond formation, interaction energies, and computational alanine scanning between PACAP and its receptors showed that His1, Asp3, Arg12, Arg14, and Lys15 are crucial to the peptide's stability. Furthermore, several PACAP interactions with structurally conserved positions deemed necessary in GPCR B1 activation, including Arg2.60, Lys2.67, and Glu7.42, were significant for the peptide's stability within the receptors. According to the protein-energy network, the connection between Asp3 of PACAP and the receptors' conserved Arg2.60 represents a critical energy communication hub in all complexes. Additionally, the ECDs of the receptors were also found to function as energy communication hubs for PACAP. Although the overall binding mode of PACAP in the three receptors was found to be highly conserved, Arg12 and Tyr13 of PACAP were more prominent in complex with PAC1, while Ser2 of PACAP was with VPAC2. The detailed analyses performed in this work pave the way for using PACAP and its receptors as therapeutic targets.Communicated by Ramaswamy H. Sarma.

Thermostabilizing mechanisms of canonical single amino acid substitutions at a GH1 ß-glucosidase probed by multiple MD and computational approaches.

ß-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable ß-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 ß-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized ß-glucosidases for 2G-biofuel production for industry, biotechnology, and science.

RASopathy Cohort of Patients Enrolled in a Brazilian Reference Center for Rare Diseases: A Novel Familial LZTR1 Variant and Recurrent Mutations.

Purpose: Noonan syndrome and related disorders are genetic conditions affecting 1:1000-2000 individuals. Variants causing hyperactivation of the RAS/MAPK pathway lead to phenotypic overlap between syndromes, in addition to an increased risk of pediatric tumors. DNA sequencing methods have been optimized to provide a molecular diagnosis for clinical and genetic heterogeneity conditions. This work aimed to investigate the genetic basis in RASopathy patients through Next Generation Sequencing in a Reference Center for Rare Diseases (IFF/Fiocruz) and implement the precision medicine at a public health institute in Brazil. Patients and Methods: This study comprises 26 cases with clinical suspicion of RASopathies. Sanger sequencing was used to screen variants in exons usually affected in the PTPN11 and HRAS genes for cases with clinical features of Noonan and Costello syndrome, respectively. Posteriorly, negative and new cases with clinical suspicion of RASopathy were analyzed by clinical or whole-exome sequencing. Results: Molecular analysis revealed recurrent variants and a novel LZTR1 missense variant: 24 unrelated individuals with pathogenic variants [PTPN11(11), NF1(2), SOS1(2), SHOC2(2), HRAS(1), BRAF(1), LZTR (1), RAF1(1), KRAS(1), RIT1(1), a patient with co-occurrence of PTPN11 and NF1 mutations (1)]; familial cases carrying a known pathogenic variant in PTPN11 (mother-two children), and a previously undescribed paternally inherited variant in LZTR1. The comparative modeling analysis of the novel LZTR1 variant p.Pro225Leu showed local and global changes in the secondary and tertiary structures, showing a decrease of about 1% in the ß-sheet content. Furthermore, evolutionary conservation indicated that Pro225 is in a highly conserved region, as observed for known dominant pathogenic variants in this protein. Conclusion: Bringing precision medicine through NGS towards congenital syndromes promotes a better understanding of complex clinical and/or undiagnosed cases. The National Policy for Rare Diseases in Brazil emphasizes the importance of incorporating and optimizing diagnostic methodologies in the Unified Brazilian Health System (SUS). Therefore, this work is an important step for the NGS inclusion in diagnostic genetic routine in the public health system.

Bacterial 2'-Deoxyguanosine Riboswitch Classes as Potential Targets for Antibiotics: A Structure and Dynamics Study.

The spread of antibiotic-resistant bacteria represents a substantial health threat. Current antibiotics act on a few metabolic pathways, facilitating resistance. Consequently, novel regulatory inhibition mechanisms are necessary. Riboswitches represent promising targets for antibacterial drugs. Purine riboswitches are interesting, since they play essential roles in the genetic regulation of bacterial metabolism. Among these, class I (2'-dG-I) and class II (2'-dG-II) are two different 2'-deoxyguanosine (2'-dG) riboswitches involved in the control of deoxyguanosine metabolism. However, high affinity for nucleosides involves local or distal modifications around the ligand-binding pocket, depending on the class. Therefore, it is crucial to understand these riboswitches' recognition mechanisms as antibiotic targets. In this work, we used a combination of computational biophysics approaches to investigate the structure, dynamics, and energy landscape of both 2'-dG classes bound to the nucleoside ligands, 2'-deoxyguanosine, and riboguanosine. Our results suggest that the stability and increased interactions in the three-way junction of 2'-dG riboswitches were associated with a higher nucleoside ligand affinity. Also, structural changes in the 2'-dG-II aptamers enable enhanced intramolecular communication. Overall, the 2'-dG-II riboswitch might be a promising drug design target due to its ability to recognize both cognate and noncognate ligands.

How does α1Histidine102 affect the binding of modulators to α1ß2γ2 GABAA receptors? molecular insights from in silico experiments.

The activation of GABAA receptors by the neurotransmitter gamma-aminobutyric acid mediates the rapid inhibition response in the central nervous system of mammals. Many neurological and mental health disorders arise from alterations in the structure or function of these pentameric ion channels. GABAA receptors are targets for numerous drugs, including benzodiazepines, which bind to α1ß2γ2 GABAA receptors with high affinity to a site in the extracellular domain, between subunits α1 and γ2. It has been established experimentally that the binding of these drugs depends on the presence of one particular amino acid in the α1 subunit: histidine 102. However, the specific role it plays in the intermolecular interaction has not been elucidated. In this work, we applied in silico methods to understand whether certain protonation and rotamer states of α1His102 facilitate the binding of modulators. We analysed diazepam binding, a benzodiazepine, and the antagonist flumazenil to the GABAA receptor using molecular dynamics simulations and adaptive biasing force simulations. The binding free energy follows changes in the protonation state for both ligands, and rotameric states of α1His102 were specific for the different compounds, suggesting distinct preferences for positive allosteric modulators and antagonists. Moreover, in the presence of diazepam and favoured by a neutral tautomer, we identified a water molecule that links loops A, B, and C and may be relevant to the modulation mechanism.

Shared Binding Mode of Perrottetinene and Tetrahydrocannabinol Diastereomers inside the CB1 Receptor May Incentivize Novel Medicinal Drug Design: Findings from an in Silico Assay.

In recent years, therapeutic compounds derived from phytocannabinoids have brought renewed attention to the benefits they offer to ameliorate chronic disease symptoms. Among cannabinoids, tetrahydrocannabinol (THC) is a well-known component of the Cannabis plant, whose active principles have been studied through the years. Another psychoactive phytocannabinoid, derived from liverworts Radula, perrottetinene (PET), has created interest, especially as a pharmaceutical product and for its legal recreational use. Unfortunately, so far, the interaction mode of these compounds at the type 1 cannabinoid receptors (CB1R) binding site remains unknown, and no experimental three-dimensional structure in complex with THC or PET is available in the Protein Data Bank. Today, many computational methodologies can assist in this crusade and help unveil how these molecules bind, based on the already known pose of a structurally similar compound. In this work, we aim to elucidate the binding mode of THC and PET molecules in both cis and trans conformers, using a combination of several computational methodologies, including molecular docking, molecular dynamics, free energy calculations, and protein-energy network studies. We found that THC and PET interact similarly with the CB1R, in a different conformation depending on the considered diastereomer. We have observed that cis ligands adopted a half-chair conformation of the cycle ring containing the dimethyl group, assuming an axial or equatorial conformation producing a different induced fitting of the surrounding residues compared with trans ligands, with higher interaction energy than the trans conformer. For PET, we have seen that Trp-279 and Trp-356 have a marked influence on the binding. After binding, Trp-279 accommodates its side chain to better interact with the PET's terminal phenyl group, disturbing CB1R residues communication. The interaction with Trp-356 might impair the activation of CB1R and can influence the binding of PET as a partial agonist. Understanding the PET association with CB1R from a molecular perspective can offer a glimpse of preventing potential toxicological or recreational effects since it is an attractive lead for drug development with fewer side effects than trans-THC.

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.

Mycobacterium tuberculosis and M. bovis BCG Moreau Fumarate Reductase Operons Produce Different Polypeptides That May Be Related to Non-canonical Functions.

Tuberculosis is a world widespread disease, caused by Mycobacterium tuberculosis (M.tb). Although considered an obligate aerobe, this organism can resist life-limiting conditions such as microaerophily mainly due to its set of enzymes responsible for energy production and coenzyme restoration under these conditions. One of these enzymes is fumarate reductase, an heterotetrameric complex composed of a catalytic (FrdA), an iron-sulfur cluster (FrdB) and two transmembrane (FrdC and FrdD) subunits involved in anaerobic respiration and important for the maintenance of membrane potential. In this work, aiming to further characterize this enzyme function in mycobacteria, we analyzed the expression of FrdB-containing proteins in M.tb and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) Moreau, the Brazilian vaccine strain against tuberculosis. We identified three isoforms in both mycobacteria, two of them corresponding to the predicted encoded polypeptides of M.tb (27 kDa) and BCG Moreau (40 kDa) frd sequences, as due to an insertion on the latter's operon a fused FrdBC protein is expected. The third 52 kDa band can be explained by a transcriptional slippage event, typically occurring when mutation arises in a repetitive region within a coding sequence, thought to reduce its impact allowing the production of both native and variant forms. Comparative modeling of the M.tb and BCG Moreau predicted protein complexes allowed the detection of subtle overall differences, showing a high degree of structure and maybe functional resemblance among them. Axenic growth and macrophage infection assays show that the frd locus is important for proper bacterial development in both scenarios, and that both M.tb's and BCG Moreau's alleles can partially revert the hampered phenotype of the knockout strain. Altogether, our results show that the frdABCD operon of Mycobacteria may have evolved to possess other yet non-described functions, such as those necessary during aerobic logarithmic growth and early stage steps of infection.

Effects of the Q80K Polymorphism on the Physicochemical Properties of Hepatitis C Virus Subtype 1a NS3 Protease.

Hepatitis C virus genotype 1a (HCV-1a) comprises clades I and II. The Q80K polymorphism is found predominantly in clade I but rarely in clade II. Here, we investigated whether natural polymorphisms in HCV-1a clade II entailed structural protein changes when occurrence of the Q80K variant was simulated. Based on HCV-1a clade I and II protein sequences, the structure of the HCV-1a Q80K mutant NS3-4A was obtained by comparative modeling. Its physicochemical properties were studied by molecular dynamics simulations and network analysis. Results demonstrate that, in the presence of the K80 variant, clade II protease polymorphisms A91 and S/G174 led to variations in hydrogen bond occupancies. Structural analyses revealed differences in (i) flexibility of the H57 catalytic residue on the NS3 protease and (ii) correlations between amino acids on the NS3 protease and the NS4A cofactor. The latter indicated possible destabilization of interactions, resulting in increased separation of these proteins. The present findings describe how the relationships between different HCV-1a NS3 protease amino acid residues could affect the appearance of viral variants and the existence of distinct genetic barriers to HCV-1a isolates.

In silico analysis of the V66M variant of human BDNF in psychiatric disorders: An approach to precision medicine.

Brain-derived neurotrophic factor (BDNF) plays an important role in neurogenesis and synapse formation. The V66M is the most prevalent BDNF mutation in humans and impairs the function and distribution of BDNF. This mutation is related to several psychiatric disorders. The pro-region of BDNF, particularly position 66 and its adjacent residues, are determinant for the intracellular sorting and activity-dependent secretion of BDNF. However, it has not yet been fully elucidated. The present study aims to analyze the effects of the V66M mutation on BDNF structure and function. Here, we applied nine algorithms, including SIFT and PolyPhen-2, for functional and stability prediction of the V66M mutation. The complete theoretical model of BNDF was generated by Rosetta and validated by PROCHECK, RAMPAGE, ProSa, QMEAN and Verify-3D algorithms. Structural alignment was performed using TM-align. Phylogenetic analysis was performed using the ConSurf server. Molecular dynamics (MD) simulations were performed and analyzed using the GROMACS 2018.2 package. The V66M mutation was predicted as deleterious by PolyPhen-2 and SIFT in addition to being predicted as destabilizing by I-Mutant. According to SNPeffect, the V66M mutation does not affect protein aggregation, amyloid propensity, and chaperone binding. The complete theoretical structure of BDNF proved to be a reliable model. Phylogenetic analysis indicated that the V66M mutation of BDNF occurs at a non-conserved position of the protein. MD analyses indicated that the V66M mutation does not affect the BDNF flexibility and surface-to-volume ratio, but affects the BDNF essential motions, hydrogen-bonding and secondary structure particularly at its pre and pro-domain, which are crucial for its activity and distribution. Thus, considering that these parameters are determinant for protein interactions and, consequently, protein function; the alterations observed throughout the MD analyses may be related to the functional impairment of BDNF upon V66M mutation, as well as its involvement in psychiatric disorders.

Unraveling RNA dynamical behavior of TPP riboswitches: a comparison between Escherichia coli and Arabidopsis thaliana.

Riboswitches are RNA sensors that affect post-transcriptional processes through their ability to bind to small molecules. Thiamine pyrophosphate (TPP) riboswitch class is the most widespread riboswitch occurring in all three domains of life. Even though it controls different genes involved in the synthesis or transport of thiamine and its phosphorylated derivatives in bacteria, archaea, fungi, and plants, the TPP aptamer has a conserved structure. In this study, we aimed at understanding differences in the structural dynamics of TPP riboswitches from Escherichia coli and Arabidopsis thaliana, based on their crystallographic structures (TPPswec and TPPswat, respectively) and dynamics in aqueous solution, both in apo and holo states. A combination of Molecular Dynamics Simulations and Network Analysis empowered to find out slight differences in the dynamical behavior of TPP riboswitches, although relevant for their dynamics in bacteria and plants species. Our results suggest that distinct interactions in the microenvironment surrounding nucleotide U36 of TPPswec (and U35 in TPPswat) are related to different responses to TPP. The network analysis showed that minor structural differences in the aptamer enable enhanced intramolecular communication in the presence of TPP in TPPswec, but not in TPPswat. TPP riboswitches of plants present subtler and slower regulation mechanisms than bacteria do.

Genomic and structural features of the yellow fever virus from the 2016-2017 Brazilian outbreak.

Southeastern Brazil has been suffering a rapid expansion of a severe sylvatic yellow fever virus (YFV) outbreak since late 2016, which has reached one of the most populated zones in Brazil and South America, heretofore a yellow fever-free zone for more than 70 years. In the current study, we describe the complete genome of 12 YFV samples from mosquitoes, humans and non-human primates from the Brazilian 2017 epidemic. All of the YFV sequences belong to the modern lineage (sub-lineage 1E) of South American genotype I, having been circulating for several months prior to the December 2016 detection. Our data confirm that viral strains associated with the most severe YF epidemic in South America in the last 70 years display unique amino acid substitutions that are mainly located in highly conserved positions in non-structural proteins. Our data also corroborate that YFV has spread southward into Rio de Janeiro state following two main sylvatic dispersion routes that converged at the border of the great metropolitan area comprising nearly 12 million unvaccinated inhabitants. Our original results can help public health authorities to guide the surveillance, prophylaxis and control measures required to face such a severe epidemiological problem. Finally, it will also inspire other workers to further investigate the epidemiological and biological significance of the amino acid polymorphisms detected in the Brazilian 2017 YFV strains.

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.

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.

Evaluation of 7-arylaminopyrazolo[1,5-a]pyrimidines as anti-Plasmodium falciparum, antimalarial, and Pf-dihydroorotate dehydrogenase inhibitors.

Malaria remains one of the most serious global infectious diseases. An important target for antimalarial chemotherapy is the enzyme dihydroorotate dehydrogenase from Plasmodium falciparum (PfDHODH), which is responsible for the conversion of dihydroorotate to orotate in the de novo pyrimidine biosynthetic pathway. In this study, we have designed and synthesized fifteen 7-arylpyrazolo[1,5-a]pyrimidine derivatives using ring bioisosteric replacement and molecular hybridization of functional groups based on the highly active 5-methyl-N-(naphthalen-2-yl)-2-(trifluoromethyl)- [1,2,4]triazolo[1,5-a]pyrimidin-7-amine. The compounds were tested against Plasmodium falciparum, as antimalarials in mice with P. berghei, and as inhibitors of PfDHODH. Thirteen compounds were found to be active against P. falciparum, with IC50 values ranging from 1.2 ± 0.3 to 92 ± 26 µM in the anti-HRP2 and hypoxanthine assays. Four compounds showed the highest selective index (SI), which is a ratio between cytotoxicity and activity in vitro. The inhibition of PfDHODH showed that compound 30 (R2 = CH3; R5 = CF3; Ar = 7-ß-naphthyl) displayed higher and selective inhibitory activity, with IC50 = 0.16 ± 0.01 µM, followed by 25 (R2 = CH3; R5 = CH3; Ar = 7-ß-Naphthyl) and 19 (R2 = CF3; R5 = CF3; Ar = 7-ß-naphthyl), with IC50 = 4 ± 1 µM and 6 ± 1 µM, respectively. The trifluoromethyl group at the 2- or 5-positions of the pyrazolo[1,5-a]pyrimidine ring led to increased drug activity. The docking results agreed with the values obtained from enzymatic assays.

Boosting Docking-Based Virtual Screening with Deep Learning.

In this work, we propose a deep learning approach to improve docking-based virtual screening. The deep neural network that is introduced, DeepVS, uses the output of a docking program and learns how to extract relevant features from basic data such as atom and residues types obtained from protein-ligand complexes. Our approach introduces the use of atom and amino acid embeddings and implements an effective way of creating distributed vector representations of protein-ligand complexes by modeling the compound as a set of atom contexts that is further processed by a convolutional layer. One of the main advantages of the proposed method is that it does not require feature engineering. We evaluate DeepVS on the Directory of Useful Decoys (DUD), using the output of two docking programs: Autodock Vina1.1.2 and Dock 6.6. Using a strict evaluation with leave-one-out cross-validation, DeepVS outperforms the docking programs, with regard to both AUC ROC and enrichment factor. Moreover, using the output of Autodock Vina1.1.2, DeepVS achieves an AUC ROC of 0.81, which, to the best of our knowledge, is the best AUC reported so far for virtual screening using the 40 receptors from the DUD.

Structural and dynamic insights into the C-terminal extension of cysteine proteinase B from Leishmania amazonensis.

Cysteine proteinase B (CPB) is a significant virulence factor for Leishmania infections. Upon processing from its zymogen form, it happens a release of the immunomodulatory CPB C-terminal extension (cyspep) into the cytoplasm of the macrophage. Epitopes derived from this fragment were shown to influence the proportion of lymphocytes CD8+ upon infection, favoring the parasite escaping from the host́s immune system. At present, there is no available structural data of cyspep, which impairs a proper understanding of its biological functions. Here, we attempted to build molecular models for this fragment and subsequently evaluate their stabilities in aqueous solution from molecular dynamics simulations analysis. Characterization of our models obtained with distinct techniques (comparative modeling, threading, and ab initio) indicates a prevalence of ß-sheets in agreement with consensus secondary structure predictions. Simulation data supported this finding since the formation of new strands, along with a rapid disruption of helical content, were observed. Overall, this study provides a rationalization of epitope mapping data and an improved understanding of cyspep antigenicity.

Caracterização do envelhecimento populacional no município do Rio de Janeiro: contribuições para políticas públicas sustentáveis / Demographic aging characterization in the city of Rio de Janeiro: contributions to sustainable public policies

Resumo Objetivo Descrever o estágio do envelhecimento populacional no município do Rio de Janeiro. Métodos Estudo ecológico tendo como unidades de observação os 160 bairros que compõem o município, utilizando indicadores sociais e demográficos construídos a partir de informações do Censo 2010. Realizou-se análise exploratória por meio de mapas temáticos e determinou-se a dependência espacial pelo Índice de Moran Global. Para agrupar bairros em estágios semelhantes do envelhecimento foi realizada uma análise de agrupamento a partir do método K-means. Resultados Encontraram-se três grupos de bairros em estágios diferentes de envelhecimento populacional, identificando-se uma tendência espacial no sentido oeste-leste com os bairros da “Zona Sul” se encontrando no estágio mais avançado de envelhecimento. Conclusão O estudo identificou as diferenças no processo de envelhecimento populacional e na composição etária dos bairros, apontando para a necessidade de políticas de saúde pública específicas que contemplem as particularidades desse processo em cada localidade, visando garantir um envelhecimento sustentável.

Abstract Objective Describe the stage of population aging in the city of Rio de Janeiro. Methods Ecological study with the 160 city neighborhoods as observational units, using social and demographic indicators built with information from the 2010 census. The exploratory analysis was undertaken with thematic maps, and the spatial dependence was measured with the Global Moran’s Index. K-means clustering was used for grouping neighborhoods with similar aging stages. Results Three neighborhood clusters in different stages of population aging were found and a spatial trend in the west-east direction was identified, with neighborhoods in the ‘South Zone’ in a more advanced stage of population aging. Conclusion The study identified differences in the population aging process and in the age composition of neighborhoods, indicating the need for specific public health policies that allow for the particularities of this process in each location, aiming a sustainable population aging.

Exploring the unbinding of Leishmania (L.) amazonensis CPB derived-epitopes from H2 MHC class I proteins.

New strategies to control Leishmania disease demand an extensive knowledge about several aspects of infection including the understanding of its molecular events. In murine models, cysteine proteinase B from Leishmania amazonensis promotes regulation of immune response, and fragments from its C-terminus extension (cyspep) can play a decisive role in the host-parasite interaction. The interaction between cyspep-derived peptides and major histocompatibility complex (MHC) proteins is a crucial factor in Leishmania infections. Seven cyspep-derived peptides, previously identified as capable of interacting with H-2 (murine) MHC class I proteins, were studied in this work. We established a protocol to simulate the unbinding of these peptides from the cleft of H-2 receptors. From the simulations, we estimated the corresponding free energy of dissociation (ΔGd ) and described the molecular events that occur during the exit of peptides from the cleft. To test the reliability of this method, we first applied it to a calibration set of four crystallographic MHC/peptide complexes. Next, we explored the unbinding of the seven complexes mentioned above. Results were consistent with ΔGd values obtained from surface plasmon resonance (SPR) experiments. We also identified some of the primary interactions between peptides and H-2 receptors, and we detected three regions of influence for the interaction. This pattern was systematically observed for the peptides and helped determine a minimum distance for the real interaction between peptides and H-2 proteins occurring at ∼ 25 Å.

Computational drug design strategies applied to the modelling of human immunodeficiency virus-1 reverse transcriptase inhibitors.

Reverse transcriptase (RT) is a multifunctional enzyme in the human immunodeficiency virus (HIV)-1 life cycle and represents a primary target for drug discovery efforts against HIV-1 infection. Two classes of RT inhibitors, the nucleoside RT inhibitors (NRTIs) and the nonnucleoside transcriptase inhibitors are prominently used in the highly active antiretroviral therapy in combination with other anti-HIV drugs. However, the rapid emergence of drug-resistant viral strains has limited the successful rate of the anti-HIV agents. Computational methods are a significant part of the drug design process and indispensable to study drug resistance. In this review, recent advances in computer-aided drug design for the rational design of new compounds against HIV-1 RT using methods such as molecular docking, molecular dynamics, free energy calculations, quantitative structure-activity relationships, pharmacophore modelling and absorption, distribution, metabolism, excretion and toxicity prediction are discussed. Successful applications of these methodologies are also highlighted.