Structural Characterization of Heat Shock Protein 90ß and Molecular Interactions with Geldanamycin and Ritonavir: A Computational Study.

Drug repositioning is an important therapeutic strategy for treating breast cancer. Hsp90ß chaperone is an attractive target for inhibiting cell progression. Its structure has a disordered and flexible linker region between the N-terminal and central domains. Geldanamycin was the first Hsp90ß inhibitor to interact specifically at the N-terminal site. Owing to the toxicity of geldanamycin, we investigated the repositioning of ritonavir as an Hsp90ß inhibitor, taking advantage of its proven efficacy against cancer. In this study, we used molecular modeling techniques to analyze the contribution of the Hsp90ß linker region to the flexibility and interaction between the ligands geldanamycin, ritonavir, and Hsp90ß. Our findings indicate that the linker region is responsible for the fluctuation and overall protein motion without disturbing the interaction between the inhibitors and the N-terminus. We also found that ritonavir established similar interactions with the substrate ATP triphosphate, filling the same pharmacophore zone.

Aptamer-Based Recognition of Breast Tumor Cells: A New Era for Breast Cancer Diagnosis.

Breast cancer is one of the leading causes of death among women worldwide and can be classified into four major distinct molecular subtypes based on the expression of specific receptors. Despite significant advances, the lack of biomarkers for detailed diagnosis and prognosis remains a major challenge in the field of oncology. This study aimed to identify short single-stranded oligonucleotides known as aptamers to improve breast cancer diagnosis. The Cell-SELEX technique was used to select aptamers specific to the MDA-MB-231 tumor cell line. After selection, five aptamers demonstrated specific recognition for tumor breast cell lines and no binding to non-tumor breast cells. Validation of aptamer specificity revealed recognition of primary and metastatic tumors of all subtypes. In particular, AptaB4 and AptaB5 showed greater recognition of primary tumors and metastatic tissue, respectively. Finally, a computational biology approach was used to identify potential aptamer targets, which indicated that CSKP could interact with AptaB4. These results suggest that aptamers are promising in breast cancer diagnosis and treatment due to their specificity and selectivity.

Next Generation of Ovarian Cancer Detection Using Aptamers.

Ovarian cancer is among the seven most common types of cancer in women, being the most fatal gynecological tumor, due to the difficulty of detection in early stages. Aptamers are important tools to improve tumor diagnosis through the recognition of specific molecules produced by tumors. Here, aptamers and their potential targets in ovarian cancer cells were analyzed by in silico approaches. Specific aptamers were selected by the Cell-SELEX method using Caov-3 and OvCar-3 cells. The five most frequent aptamers obtained from the last round of selection were computationally modeled. The potential targets for those aptamers in cells were proposed by analyzing proteomic data available for the Caov-3 and OvCar-3 cell lines. Overexpressed proteins for each cell were characterized as to their three-dimensional model, cell location, and electrostatic potential. As a result, four specific aptamers for ovarian tumors were selected: AptaC2, AptaC4, AptaO1, and AptaO2. Potential targets were identified for each aptamer through Molecular Docking, and the best complexes were AptaC2-FXYD3, AptaC4-ALPP, AptaO1-TSPAN15, and AptaO2-TSPAN15. In addition, AptaC2 and AptaO1 could detect different stages and subtypes of ovarian cancer tissue samples. The application of this technology makes it possible to propose new molecular biomarkers for the differential diagnosis of epithelial ovarian cancer.

In silico investigation of riboswitches in fungi: structural and dynamical insights into TPP riboswitches in Aspergillus oryzae.

Riboswitches are RNA sensors affecting post-transcriptional processes through their ability to bind to small molecules. Thiamine pyrophosphate (TPP) riboswitch plays a crucial role in regulating genes involved in synthesizing or transporting thiamine and phosphorylated derivatives in bacteria, archaea, plants, and fungi. Although TPP riboswitch is reasonably well known in bacteria, there is a gap in the knowledge of the fungal TPP riboswitches structure and dynamics, involving mainly sequence variation and TPP interaction with the aptamers. On the other hand, the increase of fungal infections and antifungal resistance raises the need for new antifungal therapies. In this work, we used computational approaches to build three-dimensional models for the three TPP riboswitches identified in Aspergillus oryzae, in which we studied their structure, dynamics, and binding free energy change (ΔGbind) with TPP. Interaction patterns between the TPP and the surrounding nucleotides were conserved among the three models, evidencing high structural conservation. Furthermore, we show that the TPP riboswitch from the A. oryzae NMT1 gene behaves similarly to the E. coli thiA gene concerning the ΔGbind. In contrast, mutations in the fungal TPP riboswitches from THI4 and the nucleoside transporter genes led to structural differences, affecting the binding-site volume, hydrogen bond occupancy, and ΔGbind. Besides, the number of water molecules surrounding TPP influenced the ΔGbind considerably. Notably, our ΔGbind estimation agreed with previous experimental data, reinforcing the relationship between sequence conservation and TPP interaction.

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.

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.

The rare lncRNA GOLLD is widespread and structurally conserved among Mycobacterium tRNA arrays.

Noncoding RNA (ncRNA) genes produce transcripts involved in a wide range of functions, including catalytic and regulatory functions. Besides, some transcripts have highly complex structures that may impact their activities. Among the largest bacterial ncRNAs, there is the rare GOLLD RNA, which is associated with tRNA genes and supposed to be chromosome- and phage-encoded in specialized groups of bacteria, including those from Lactobacillales and Actinomycetales orders. The only GOLLD structure was inferred from a variety of sequences, including many marine metagenomes. To explore GOLLD RNA in bacterial genomes, we mined the GOLLD gene in thousands of Mycobacterium and virus genomes using Infernal software. We identified this gene in 350 mycobacteria, including megaplasmids, and 39 bacteriophages, mainly in the genomic context of tRNA arrays. Mycobacterium GOLLD genes presented a high diversity and were distributed in three phylogenetic groups: (i) Mycobacterium exclusive; (ii) Mycobacterium and mycobacteriophages; and (iii) mycobacteriophage exclusive. We also determined the GOLLD secondary structure of each group using R2 R software based on GOLLD alignments generated by Infernal software. All GOLLD groups displayed a 3' half conserved structure, including utter E-loops pseudoknots substructures, also shared by non-Mycobacterium GOLLD while the 5' half motif was different among the groups. Here, we showed that the lncRNA GOLLD is widespread in Mycobacterium within tRNA arrays and corroborated the previously predicted GOLLD secondary structure.

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.

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.

Exploring Novel Antimalarial Compounds Targeting Plasmodium falciparum Enoyl-ACP Reductase: Computational and Experimental Insights.

Malaria, caused by Plasmodium protozoa with Plasmodium falciparum as the most virulent species, continues to pose significant health challenges. Despite the availability of effective antimalarial drugs, the emergence of resistance has heightened the urgency for developing novel therapeutic compounds. In this study, we investigated the enoyl-ACP reductase enzyme of P. falciparum (PfENR) as a promising target for antimalarial drug discovery. Through a comprehensive analysis, we conducted a comparative evaluation of two lead compounds, LD1 (CID: 44405336, lead compounds 1) and LD2 (CID: 72703246, lead compounds 2), obtained from the PubChem/NCBI ligand database, to serve as reference molecules in the identification of potential derivatives using virtual screening assays. Among the newly identified candidates, Ligand 1 (LG1) and Ligand 2 (LG2) exhibited intriguing characteristics and underwent further investigation through docking and molecular dynamics simulations. Ligand 1 (LG1) demonstrated interactions similar to LD1, including hydrogen bonding with Asp218, while Ligand 2 (LG2) displayed superior binding energy comparable to LD1 and LD2, despite lacking hydrogen bonding interactions observed in the control compounds triclosan and its derivative 7-(4-chloro-2-hydroxyphenoxy)-4-methyl-2H-chromen-2-one (CHJ). Following computational validation using the MM/GBSA method to estimate binding free energy, commercially acquired LG1 and LG2 ligands were subjected to in vitro testing. Inhibition assays were performed to evaluate their potential as PfENR inhibitors alongside triclosan as a control compound. LG1 exhibited no inhibitory effects, while LG2 demonstrated inhibitory effects like triclosan. In conclusion, this study contributes valuable insights into developing novel antimalarial drugs by identifying LG2 as a potential ligand and employing a comprehensive approach integrating computational and experimental methodologies.

Characterization of Mycobacterium smegmatis Glutaminase-Free Asparaginase (MSMEG_3173).

l-asparaginase is an enzyme catalyzing the hydrolysis of l-asparagine into l-aspartate and ammonia, which is of great therapeutic importance in tumor treatment. However, commercially available enzymes are associated with adverse effects, and searching for a new l-asparaginase with better pharmaceutical properties was the aim of this work. The coding sequence for Mycobacterium smegmatisl-asparaginase (MsA) was cloned and expressed. The recombinant protein showed high activity toward l-asparagine, whereas none was detected for l-glutamine. The enzymatic properties (K m = 1.403 ± 0.24 mM and k cat = 708.1 ± 25.05 s-1) indicate that the enzyme would be functional within the expected blood l-asparagine concentration, with good activity, as shown by k cat. The pH and temperature profiles suggest its use as a biopharmaceutical in humans. Molecular dynamics analysis of the MsA model reveals the formation of a hydrogen bond network involving catalytic residues with l-asparagine. However, the same is not observed with l-glutamine, mainly due to steric hindrance. Additionally, the structural feature of residue 119 being a serine rather than a proline has significant implications. These findings help explain the low glutaminase activity observed in MsA, like what is described for the Wolinella succinogenes enzyme. This establishes mycobacterial asparaginases as key scaffolds to develop biopharmaceuticals against acute lymphocytic leukemia.

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.

Revisiting Jatropha curcas Monomeric Esterase: A Dienelactone Hydrolase Compatible with the Electrostatic Catapult Model.

Jatropha curcas contains seeds with a high oil content, suitable for biodiesel production. After oil extraction, the remaining mass can be a rich source of enzymes. However, data from the literature describing physicochemical characteristics for a monomeric esterase from the J. curcas seed did not fit the electrostatic catapult model for esterases/lipases. We decided to reevaluate this J. curcas esterase and extend its characterization to check this apparent discrepancy and gain insights into the enzyme's potential as a biocatalyst. After anion exchange chromatography and two-dimensional gel electrophoresis, we identified the enzyme as belonging to the dienelactone hydrolase family, characterized by a cysteine as the nucleophile in the catalytic triad. The enzyme displayed a basic optimum hydrolysis pH of 9.0 and an acidic pI range, in contrast to literature data, making it well in line with the electrostatic catapult model. Furthermore, the enzyme showed low hydrolysis activity in an organic solvent-containing medium (isopropanol, acetonitrile, and ethanol), which reverted when recovering in an aqueous reaction mixture. This enzyme can be a valuable tool for hydrolysis reactions of short-chain esters, useful for pharmaceutical intermediates synthesis, due to both its high hydrolytic rate in basic pH and its stability in an organic solvent.

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.

Studying effects of different protonation states of His11 and His102 in ribose-5-phosphate isomerase of Trypanosoma cruzi: an example of cooperative behavior.

The Trypanosoma cruzi ribose-5-phosphate isomerase B (TcRpiB) is a crucial piece in the pentose phosphate pathway and thus is a potential drug target for treatment of Chagas' disease. TcRpiB residues, such as Cys69, Asp45, Glu149 and Pro47, have confirmed their roles in substrate recognition, catalytic reaction and binding site conformation. However, the joint performance of His11 and His102, in the D-ribose-5-phosphate (R5P) in the catalysis is not well understood. In this work, we probed the influence of different protonation states of His11 and His102 on the behavior of the ligand R5P using molecular dynamics simulations, network analysis and thermodynamic integration. Simulations revealed that a protonated His11 combined with a neutral His102 (His11+‒His102) was able to stabilize the ligand R5P in the binding site. Moreover, calculated relative free energy differences showed that when protonated His11 was coupled to a neutral His102 an exergonic process takes place. On the other hand, neutral His11 combined with a protonated His102 (His11‒His102+), sampled conformations that resembled the catalyzed product D-ribulose-5-phosphate (Ru5P). Network analysis also demonstrated some peculiarities for these systems with some negatively correlated nodes in the binding site for His11‒His102+, and exclusive suboptimal paths for His11+‒His102. Therefore, the combined approach presented in this paper proposes two suitable protonation states for the TcRpiB catalytic mechanism, where an extra proton in either histidines might favor R5P binding or influence isomerization reaction to Ru5P. Our results may guide further in silico drug discovery studies. Communicated by Ramaswamy H. Sarma.

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.

Identification of Novel and Recurrent RMRP Variants in a Series of Brazilian Patients with Cartilage-Hair Hypoplasia: McKusick Syndrome.

Cartilage-hair hypoplasia syndrome (CHH) is an autosomal recessive disorder caused by pathogenic variants of the RMRP gene and characterized by metaphyseal bone dysplasia associated with hypotrichosis, immunodeficiency, and predisposition to malignancy. However, the genotype-phenotype correlation in CHH is not well understood. Here, we report a single country cohort of 23 Brazilian patients with clinical and radiological features consistent with CHH. We found 23 different pathogenic variants in the RMRP gene - 12 novel and 11 previously described in the literature. Interestingly, the most frequent Finnish pathogenic variant related to CHH (g.71A>G) was not found in our cohort. In contrast, more than 50% of the patients carried the rare g.196C>T variant suggesting a possible founder effect in the Brazilian population. In silico analysis showed that pathogenic variants occurred either in the regions conserved in mammalian species or within essential domains for the ribonucleoprotein complex. Pathogenicity prediction studies can improve the understanding of how these variants affect RNA.

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.