Hydropathic AF-2 variants in the androgen receptor gene among androgen insensitivity patients.

BACKGROUND: Androgen insensitivity syndrome (AIS) is a common condition among individuals with differences of sexual development (DSD) and results from germline allelic variants in the androgen receptor (AR) gene. Understanding the phenotypic consequences of AR allelic variants that disrupt the activation function 2 (AF2) region is essential to grasping its clinical significance. OBJECTIVES: This study aims to provide insights into the phenotypic characteristics and clinical impact of AR mutations affecting the AF2 region in AIS patients. We achieve this by reviewing reported AR variants in the AF2 region among individuals with AIS, including identifying a new phenotype associated with the c.2138T>C variant (p.Leu713Pro) in the AR gene. MATERIALS AND METHODS: We comprehensively reviewed AR variants within the AF2 region reported in AIS and applied molecular dynamics simulations to assess the impact of the p.Leu713Pro variant on protein dynamics. RESULTS: Our review of reported AR variants in the AF2 region revealed a spectrum of phenotypic outcomes in AIS patients. Molecular dynamics simulations indicated that the p.Leu713Pro variant significantly alters the local dynamics of the AR protein and disrupts the correlation and covariance between variables. DISCUSSION: The diverse phenotypic presentations observed among individuals with AR variants in the AF2 region highlight the complexity of AIS. The altered protein dynamics resulting from the p.Leu713Pro variant further emphasize the importance of the AF2 region in AR function. CONCLUSION: Our study provides valuable insights into AR mutations' phenotypic characteristics and clinical impact on the AF2 region in AIS. Moreover, the disruption of protein dynamics underscores the significance of the AF2 region in AR function and its role in the pathogenesis of AIS.

Epigenetic regulation exerted by Caliphruria subedentata and galantamine: an in vitro and in silico approach for mimic Alzheimer's disease.

At the interface between genes and environment, epigenetic mechanisms, including DNA methylation and histone modification, regulate neurogenic processes such as differentiation, proliferation, and maturation of neural stem cells. However, these mechanisms are altered in Alzheimer's disease (AD), a neurodegenerative condition that mainly affects older adults. Since epigenetic mechanisms are known to be reversible, a number of molecules from natural sources are being studied as epigenetic regulators in AD. Recently, in vitro and in silico studies have shown that C. subedentata and its alkaloids modulated neurotoxicity. However, studies exploring the epigenetic activity of these alkaloids are limited. We conducted a set of bioassays to evaluate neuronal differentiation and the sensitivity of undifferentiated SH-SY5 cells against a neurotoxic stimulus. In addition, we analyzed the methylation profiles in genes such as APP, PSI, and BACE1 due to their role in amyloid processing. Docking and molecular dynamic analysis were used to explore the effect exerted by C. subedentata alkaloids on the regulation of histone deacetylases (HDAC2, HDAC3 and HDAC7). The results demonstrated that C. subedentata and galantamine induce neuronal differentiation and protect the undifferentiated SH-SY5Y cells against Aß(1-42)-induced neurotoxicity. The methylation profiles of the studied genes show no statistically significant differences between C. subedentata, galantamine. However, these findings should be interpreted with caution, since small changes in methylation promoters in the brain could not be easily detected. Results from in silico approaches describe for the first time the potential promissing epigenetic effects of galantamine by regulating HDAC3 and HDAC7 modification.Communicated by Ramaswamy H. Sarma.

Molecular Docking of Lac_CB10: Highlighting the Great Potential for Bioremediation of Recalcitrant Chemical Compounds by One Predicted Bacteroidetes CopA-Laccase.

Laccases are multicopper oxidases (MCOs) with a broad application spectrum, particularly in second-generation ethanol biotechnology and the bioremediation of xenobiotics and other highly recalcitrant compounds. Synthetic pesticides are xenobiotics with long environmental persistence, and the search for their effective bioremediation has mobilized the scientific community. Antibiotics, in turn, can pose severe risks for the emergence of multidrug-resistant microorganisms, as their frequent use for medical and veterinary purposes can generate constant selective pressure on the microbiota of urban and agricultural effluents. In the search for more efficient industrial processes, some bacterial laccases stand out for their tolerance to extreme physicochemical conditions and their fast generation cycles. Accordingly, to expand the range of effective approaches for the bioremediation of environmentally important compounds, the prospection of bacterial laccases was carried out from a custom genomic database. The best hit found in the genome of Chitinophaga sp. CB10, a Bacteroidetes isolate obtained from a biomass-degrading bacterial consortium, was subjected to in silico prediction, molecular docking, and molecular dynamics simulation analyses. The putative laccase CB10_180.4889 (Lac_CB10), composed of 728 amino acids, with theoretical molecular mass values of approximately 84 kDa and a pI of 6.51, was predicted to be a new CopA with three cupredoxin domains and four conserved motifs linking MCOs to copper sites that assist in catalytic reactions. Molecular docking studies revealed that Lac_CB10 had a high affinity for the molecules evaluated, and the affinity profiles with multiple catalytic pockets predicted the following order of decreasing thermodynamically favorable values: tetracycline (-8 kcal/mol) > ABTS (-6.9 kcal/mol) > sulfisoxazole (-6.7 kcal/mol) > benzidine (-6.4 kcal/mol) > trimethoprim (-6.1 kcal/mol) > 2,4-dichlorophenol (-5.9 kcal/mol) mol. Finally, the molecular dynamics analysis suggests that Lac_CB10 is more likely to be effective against sulfisoxazole-like compounds, as the sulfisoxazole-Lac_CB10 complex exhibited RMSD values lower than 0.2 nm, and sulfisoxazole remained bound to the binding site for the entire 100 ns evaluation period. These findings corroborate that LacCB10 has a high potential for the bioremediation of this molecule.

Evidence for Epistatic Interaction between HLA-G and LILRB1 in the Pathogenesis of Nonsegmental Vitiligo.

Vitiligo is the most frequent cause of depigmentation worldwide. Genetic association studies have discovered about 50 loci associated with disease, many with immunological functions. Among them is HLA-G, which modulates immunity by interacting with specific inhibitory receptors, mainly LILRB1 and LILRB2. Here we investigated the LILRB1 and LILRB2 association with vitiligo risk and evaluated the possible role of interactions between HLA-G and its receptors in this pathogenesis. We tested the association of the polymorphisms of HLA-G, LILRB1, and LILRB2 with vitiligo using logistic regression along with adjustment by ancestry. Further, methods based on the multifactor dimensionality reduction (MDR) approach (MDR v.3.0.2, GMDR v.0.9, and MB-MDR) were used to detect potential epistatic interactions between polymorphisms from the three genes. An interaction involving rs9380142 and rs2114511 polymorphisms was identified by all methods used. The polymorphism rs9380142 is an HLA-G 3'UTR variant (+3187) with a well-established role in mRNA stability. The polymorphism rs2114511 is located in the exonic region of LILRB1. Although no association involving this SNP has been reported, ChIP-Seq experiments have identified this position as an EBF1 binding site. These results highlight the role of an epistatic interaction between HLA-G and LILRB1 in vitiligo pathogenesis.

Computational and atomistic studies applied to the understanding of the structural and behavioral features of the immune checkpoint HLA-G molecule and gene.

We took advantage of the increasingly evolving approaches for in silico studies concerning protein structures, protein molecular dynamics (MD), protein-protein and protein-DNA docking to evaluate: (i) the structure and MD characteristics of the HLA-G well-recognized isoforms, (ii) the impact of missense mutations at HLA-G receptor genes (LILRB1/2), and (iii) the differential binding of the hypoxia-inducible factor 1 (HIF1) to hypoxia-responsive elements (HRE) at the HLA-G gene. Besides reviewing these topics, they were revisited including the following novel results: (i) the HLA-G6 isoforms were unstable docked or not with ß2-microglobulin or peptide, (ii) missense mutations at LILRB1/2 genes, exchanging amino acids at the intracellular domain, particularly those located within and around the ITIM motifs, may impact the HLA-G binding strength, and (iii) HREs motifs at the HLA-G promoter or exon 2 regions exhibiting a guanine at their third position present a higher affinity for HIF1 when compared to an adenine at the same position. These data shed some light into the functional aspects of HLA-G, particularly how polymorphisms may influence the role of the molecule. Computational and atomistic studies have provided alternative tools for experimental physical methodologies, which are time-consuming, expensive, demanding large quantities of purified proteins, and exhibit low output.

Genome-wide analysis of the chromatin sites targeted by HEX 70a storage protein in the honeybee brain and fat body.

Hexamerins, the proteins massively stored in the larval haemolymph of insects, are gradually used throughout metamorphosis as a source of raw material and energy for the development of adult tissues. Such behaviour defined hexamerins as storage proteins. Immunofluorescence experiments coupled with confocal microscopy show a hexamerin, HEX 70a, in the nucleus of the brain and fat body cells from honeybee workers, an unexpected localization for a storage protein. HEX 70a colocalizes with fibrillarin, a nucleolar-specific protein and H3 histone, thus suggesting a potential role as a chromatin-binding protein. This was investigated through chromatin immunoprecipitation and high-throughput DNA sequencing (ChIP-seq). The significant HEX 70a-DNA binding sites were mainly localized at the intergenic, promoter and intronic regions. HEX 70a targeted DNA stretches mapped to the genomic regions encompassing genes with relevant functional attributes. Several HEX 70a targeted genes were associated with H3K27ac or/and H3K27me3, known as active and repressive histone marks. Brain and fat body tissues shared a fraction of the HEX 70 targeted genes, and tissue-specific targets were also detected. The presence of overrepresented DNA motifs in the binding sites is consistent with specific HEX 70a-chromatin association. In addition, a search for HEX 70a targets in RNA-seq public libraries of fat bodies from nurses and foragers revealed differentially expressed targets displaying hex 70a-correlated developmental expression, thus supporting a regulatory activity for HEX 70a. Our results support the premise that HEX 70a is a moonlighting protein that binds chromatin and has roles in the brain and fat body cell nuclei, apart from its canonical role as a storage protein.

Virtual screening and molecular dynamics study of natural products against Rab10 for the treatment of Alzheimer's disease.

Alzheimer's disease (AD) is a neurodegenerative disorder associated with aging. Various enzymatic targets have been and are still being studied in an attempt to discover new drugs for the treatment of AD; however, Rab GTPases are still relatively unexplored. These enzymes regulate cellular processes by alternating of GDP and GTP nucleotides. In vitro studies have shown that the knockdown of Rab10 reduces the production of Aß40 and Aß42 peptides, making it a promising target for the treatment of AD. In order to identify potential Rab10 inhibitors, the structure-based virtual screening (SBVS) was used considering a subset of 80763 natural products obtained from ZINC15 database. Tertiary structure of Rab10 was obtained from the Protein Data Bank and the Autodock Vina program was used in the SBVS to filter potential bioactive substances against this enzyme. The SBVS protocol was validated by redocking the co-crystallized GNP and the binding energies of the GDP and GTP were used as controls in the pharmacodynamic analysis. Thus, it was possible to select 45 compounds with binding energy less or equal -11 kcal.mol-1. ADME/T properties of these compounds were evaluated by the SwissADME program, where it was possible to identify 6 promising molecules. The resulting complexes were subjected to molecular dynamics simulations to analyze the pharmacodynamics over time. The results suggest that the compound ZINC4090657 (derived from quinolizidine) and the compounds ZINC4000106 and ZINC0630250 (derived from coumarin) have favorable pharmacological characteristics for the inhibition of Rab10, with ZINC4090657 being the most promising one.Communicated by Ramaswamy H. Sarma.

Genetic diversity of the LILRB1 and LILRB2 coding regions in an admixed Brazilian population sample.

Leukocyte immunoglobulin (Ig)-like receptors (LILR) LILRB1 and LILRB2 may play a pivotal role in maintaining self-tolerance and modulating the immune response through interaction with classical and nonclassical HLA molecules. Although both diversity and natural selection patterns over HLA genes have been extensively evaluated, little information is available concerning the genetic diversity and selection signatures on the LILRB1/2 regions. Therefore, we identified the LILRB1/2 genetic diversity using next-generation sequencing in a population sample from São Paulo State, Brazil. We identified 58 LILRB1 Single Nucleotide Variants (SNVs), which gave rise to 13 haplotypes, and 41 LILRB2 SNVs arranged into 11 haplotypes. Although we may not exclude as a possible effect of population structure, we found evidence of either positive or purifying selection on LILRB1/2 coding regions. Some residues in both proteins showed to be under the effect of positive selection, suggesting that amino acid replacements in these proteins resulted in beneficial functional changes. Finally, we have revealed that allelic variation (six and five amino acid exchanges in LILRB1 and LILRB2, respectively) affects the structure and/or stability of both molecules. Nonetheless, LILRB2 has shown higher average stability, with no D1/D2 residue affecting protein structure. Overall, our findings demonstrate that LILRB1 and LILRB2 are as polymorphic as HLA class Ib genes and provide strong evidence supporting the directional selection regime hypothesis.

Identification of novel potential cyclooxygenase-2 inhibitors using ligand- and structure-based virtual screening approaches.

Cyclooxygenase 2 (COX-2) is a well-established target for the design of anti-inflammatory intermediates. Celecoxib was selected as a template molecule to perform ligand-based virtual screening, i.e. to search for structures with similarity in shape and electrostatic potential, with a gradual increase in accuracy through the combined fitting of several steps using eight commercial databases. The molecules ZINC408709 and ZINC2090319 reproduced values within the limits established in an initial study of absorption and distribution in the body. No alert was fired for possible toxic groups when these molecules were subjected to toxicity prediction. Molecular docking results with these compounds showed a higher binding affinity in comparison to rofecoxib for the COX-2 target. Additionally, ZINC408709 and ZINC2090319 were predicted to be potentially biologically active. In in silico prediction of endocrine disruption potential, it was established that the molecule ZINC2090319 binds strongly to the target related to cardiovascular risk in a desirable way as a non-steroidal antagonist and ZINC408709 binds strongly to the target that is associated with the treatment of inflammatory pathologies and similar to celecoxib. Metabolites generated from these compounds are less likely to have side effects. Simulations were used to evaluate the interaction of compounds with COX-1 and COX-2 during 200 ns. Despite the differences, ZINC408709 molecule showed better stability for COX-2 during molecular dynamics simulation. In the calculations of free energy MM/PBSA, the molecule ZINC408709 ΔGbind value has a higher affinity to celecoxib and rofecoxib COX-2. This demonstrates that the selected substances can be considered as promising COX-2 inhibitors. Communicated by Ramaswamy H. Sarma.

The Role of Celecoxib as a Potential Inhibitor in the Treatment of Inflammatory Diseases - A Review.

This article aims at reviewing celecoxib as a potential inhibitor in the treatment of inflammatory diseases. The enzyme cyclooxygenase (COX) predominantly has two isoforms called cyclooxygenase 1 (COX-1) and cyclooxygenase 2 (COX-2). The former plays a constitutive role related to homeostatic effects in renal and platelets, while the latter is mainly responsible for the induction of inflammatory effects. Since COX-2 plays an important role in the pathogenesis of inflammatory diseases, it has been signaled as a target for the planning of anti-inflammatory intermediates. Many inhibitors developed and planned for COX-2 inhibition have presented side effects to humans, mainly in the gastrointestinal and/or cardiovascular tract. Therefore, it is necessary to design new potential COX-2 inhibitors, which are relatively safe and have no side effects. In this sense, celecoxib is the only potent, selective COX-2 inhibitor that is still commercially available (within the "coxib" family). Thus, celecoxib became a commercial prototype inhibitor for the development of anti-inflammatory agents for the COX-2 enzyme. This review provides inhibition highlights that should provide a structural basis for the design of promising new non-steroidal anti-inflammatory drugs (NSAIDs), which act as COX-2 inhibitors with lesser side effects on the human body.

Structural Characterization of the Interaction of Hypoxia Inducible Factor-1 with Its Hypoxia Responsive Element at the -964G > A Variation Site of the HLA-G Promoter Region.

Human Antigen Leukocyte-G (HLA-G) gene encodes an immune checkpoint molecule that has restricted tissue expression in physiological conditions; however, the gene may be induced in hypoxic conditions by the interaction with the hypoxia inducible factor-1 (HIF1). Hypoxia regulatory elements (HRE) located at the HLA-G promoter region and at exon 2 are the major HIF1 target sites. Since the G allele of the -964G > A transversion induces higher HLA-G expression when compared to the A allele in hypoxic conditions, here we analyzed HIF1-HRE complex interaction at the pair-atom level considering both -964G > A polymorphism alleles. Mouse HIF2 dimer crystal (Protein Data Bank ID: 4ZPK) was used as template to perform homology modelling of human HIF1 quaternary structure using MODELLER v9.14. Two 3D DNA structures were built from 5'GCRTG'3 HRE sequence containing the -964G/A alleles using x3DNA. Protein-DNA docking was performed using the HADDOCK v2.4 server, and non-covalent bonds were computed by DNAproDB server. Molecular dynamic simulation was carried out per 200 ns, using Gromacs v.2019. HIF1 binding in the HRE containing -964G allele results in more hydrogen bonds and van der Waals contact formation than HRE with -964A allele. Protein-DNA complex trajectory analysis revealed that HIF1-HRE-964G complex is more stable. In conclusion, HIF1 binds in a more stable and specific manner at the HRE with G allele.

GDAP1 mutations are frequent among Brazilian patients with autosomal recessive axonal Charcot-Marie-Tooth disease.

Mutations in ganglioside-induced differentiation-associated-protein 1 (GDAP1) are associated with several subtypes of Charcot-Marie-Tooth (CMT) disease, including autosomal recessive and demyelinating (CMT4A); autosomal recessive and axonal (AR-CMT2K); autosomal dominant and axonal (CMT2K); and an intermediate and recessive form (CMTRIA). To date, at least 103 mutations in this gene have been described, but the relative frequency of GDAP1 mutations in the Brazilian CMT population is unknown. In this study, we investigated the frequency of GDAP1 mutations in a cohort of 100 unrelated Brazilian CMT patients. We identified five variants in unrelated axonal CMT patients, among which two were novel and probably pathogenic (N64S, P119T) one was novel and was classified as VUS (K207L) and two were known pathogenic variants (R125* and Q163*). The prevalence rate of GDAP1 among the axonal CMT cases was 7,14% (5/70), all of them of recessive inheritance, thus suggesting that the prevalence was higher than what is observed in most countries. All patients exhibited severe early-onset CMT that was rapidly progressive. Additionally, this study widens the mutational spectrum of GDAP1-related CMT through identification of two novel likely pathogenic variants.

Warming and elevated CO2 induces changes in the reproductive dynamics of a tropical plant species.

Tropical plant species are vulnerable to climate change and global warming. Since flowering is a critical factor for plant reproduction and seed-set, warming and elevated atmospheric carbon dioxide concentrations (eCO2) are crucial climate change factors that can affect plant reproductive dynamics and flowering related events in the tropics. Using a combined free-air CO2 enrichment and a free-air temperature-controlled enhancement system, we investigate how warming (+2 °C above ambient, eT) and elevated [CO2] (~600 ppm, eCO2) affect the phenological pattern, plant-insect interactions, and outcrossing rates in the tropical legume forage species Stylosanthes capitata Vogel (Fabaceae). In comparison to the control, a significantly greater number of flowers (NF) per plot (+62%) were observed in eT. Furthermore, in warmed plots flowers began opening approximately 1 h earlier (~09:05), with a canopy temperature of ~23 °C, than the control (~09:59) and eCO2 (~09:55) treatments. Flower closure occurred about 3 h later in eT (~11:57) and control (~13:13), with a canopy temperature of ~27 °C. These changes in flower phenology increased the availability of floral resources and attractiveness for pollinators such as Apis mellifera L. and visitors such as Paratrigona lineata L., with significant interactions between eT treatments and insect visitation per hour/day, especially between 09:00-10:40. In comparison to the control, the additive effects of combined eCO2 + eT enhanced the NF by 137%, while the number of A. mellifera floral visits per plot/week increased by 83% during the period of greatest flower production. Although we found no significant effect of treatments on mating system parameters, the overall mean multilocus outcrossing rate (tm = 0.53 ± 0.03) did confirm that S. capitata has a mixed mating system. The effects of elevated [CO2] and warming on plant-pollinator relationships observed here may have important implications for seed production of tropical forage species in future climate scenarios.

Precursor-directed Biosynthesis in Tabernaemontana catharinensis as a New Avenue for Alzheimer's Disease-modifying Agents.

Plants produce a high diversity of metabolites that can act as regulators of cholinergic dysfunction. Among plants, the potential of species of the genus Tabernaemontana to treat neurological disorders has been linked to iboga-type alkaloids that are biosynthesized by those species. In this context, precursor-directed biosynthesis approaches were carried out using T. catharinensis plantlets to achieve new-to-nature molecules as promising agents against Alzheimer's disease. Aerial parts of T. catharinensis, cultured in vitro, produced 7 unnatural alkaloids (5-fluoro-ibogamine, 5-fluoro-voachalotine, 5-fluoro-12-methoxy-Nb-methyl-voachalotine, 5-fluoro-isovoacangine, 5-fluoro-catharanthine, 5-fluoro-19-(S)-hydroxy-ibogamine, and 5-fluoro-coronaridine), while root extracts showed the presence of the same unnatural iboga-type alkaloids and 2 additional ones: 5-fluoro-voafinine and 5-fluoro-affinisine. Moreover, molecular docking approaches were carried out to evaluate the potential inhibition activity of T. catharinensis' natural and unnatural alkaloids against AChE and BChE enzymes. Fluorinated iboga alkaloids (5-fluoro-catharanthine, 5-fluoro-voachalotine, 5-fluoro-affinisine, 5-fluoro-isovoacangine, 5-fluoro-corinaridine) were more active than natural ones and controls against AchE, while 5-fluoro-19-(S)-hydroxy-ibogamine, 5-fluoro-catharanthine, 5-fluoro-isovoacangine, and 5-fluoro-corinaridine showed better activity than natural ones and controls against BChE. Our findings showed that precursor-directed biosynthesis strategies generated "new-to-nature" alkaloids that are promising Alzheimer's disease drug candidates. Furthermore, the isotopic experiments also allowed us to elucidate the initial steps of the biosynthetic pathway for iboga-type alkaloids, which are derived from the MEP and shikimate pathways.

Assessment of the hypoglycemic effect of Bixin in alloxan-induced diabetic rats: in vivo and in silico studies.

The objectives of this study were to extract and purify Bixin from the seeds of Bixa orellana and to evaluate its hypoglycemic activity in vivo, as well as, to conduct an in silico study of selectivity on peroxisome proliferator-activated receptors via molecular docking and molecular dynamics simulations. Oral administration of Bixin (10 mg/kg) significantly reduced their glucose level that was alloxan-induced diabetic rats. Bixin showed in silico selectivity on peroxisome proliferator-activated receptors (PPARs), particularly by the peroxisome proliferator-activated receptor gamma (PPARγ), which supports the hypoglycemic activity of Bixin. From the results obtained, it can be inferred that Bixin presents hypoglycemic characteristics, which was confirmed by the results obtained from the in vivo and in silico tests. Bixin may act by other pathways to control blood glucose and thus it is possible that it presents a different toxicity profile than troglitazone, rosiglitazone and pioglitazone. However, more studies on the activity and toxicity of Bixin are needed to evaluate for further clinical use. Communicated by Ramaswamy H. Sarma.

An Animal Able To Tolerate D2 O.

It is possible to gain a deeper insight into the role of water in biology by using physicochemical variant molecules, such as deuterium oxide (D2 O); however, D2 O is toxic to multicellular organisms in high concentrations. By using a unique desiccation-rehydration process, we demonstrate that the anhydrobiotic nematode Panagrolaimus superbus is able to tolerate and proliferate in 99 % D2 O. Moreover, we analysed P. superbus' water-channel protein (aquaporin; AQP), which is associated with dehydration/rehydration, by comparing its primary structure and modelling its tertiary structure in silico. Our data evidence that P. superbus' AQP is an aquaglyceroporin, a class of water channel known to display a wider pore; this helps to explain the rapid and successful organismal influx of D2 O into this species. This is the first demonstration of an animal able to withstand high D2 O levels, thus paving a way for the investigation of the effects D2 O on higher levels of biological organization.

Structural Modeling and Molecular Dynamics of the Immune Checkpoint Molecule HLA-G.

HLA-G is considered to be an immune checkpoint molecule, a function that is closely linked to the structure and dynamics of the different HLA-G isoforms. Unfortunately, little is known about the structure and dynamics of these isoforms. For instance, there are only seven crystal structures of HLA-G molecules, being all related to a single isoform, and in some cases lacking important residues associated to the interaction with leukocyte receptors. In addition, they lack information on the dynamics of both membrane-bound HLA-G forms, and soluble forms. We took advantage of in silico strategies to disclose the dynamic behavior of selected HLA-G forms, including the membrane-bound HLA-G1 molecule, soluble HLA-G1 dimer, and HLA-G5 isoform. Both the membrane-bound HLA-G1 molecule and the soluble HLA-G1 dimer were quite stable. Residues involved in the interaction with ILT2 and ILT4 receptors (α3 domain) were very close to the lipid bilayer in the complete HLA-G1 molecule, which might limit accessibility. On the other hand, these residues can be completely exposed in the soluble HLA-G1 dimer, due to the free rotation of the disulfide bridge (Cys42/Cys42). In fact, we speculate that this free rotation of each protomer (i.e., the chains composing the dimer) could enable alternative binding modes for ILT2/ILT4 receptors, which in turn could be associated with greater affinity of the soluble HLA-G1 dimer. Structural analysis of the HLA-G5 isoform demonstrated higher stability for the complex containing the peptide and coupled ß2-microglobulin, while structures lacking such domains were significantly unstable. This study reports for the first time structural conformations for the HLA-G5 isoform and the dynamic behavior of HLA-G1 molecules under simulated biological conditions. All modeled structures were made available through GitHub (https://github.com/KavrakiLab/), enabling their use as templates for modeling other alleles and isoforms, as well as for other computational analyses to investigate key molecular interactions.

Genistein and Galantamine Combinations Decrease ß-Amyloid Peptide (1-42)-Induced Genotoxicity and Cell Death in SH-SY5Y Cell Line: an In Vitro and In Silico Approach for Mimic of Alzheimer's Disease.

Alzheimer's disease (AD) is the primary dementia-causing disease worldwide, involving a multifactorial combination of environmental, genetic, and epigenetic factors, with essential participation of age and sex. Biochemically, AD is characterized by the presence of abnormal deposition of beta amyloid peptide (Aß(1-42)), which in the brain is strongly correlated with oxidative stress, inflammation, DNA damage, and cholinergic impairment. The multiple mechanisms involved in its etiology create significant difficulty in producing an effective treatment. Neuroprotective properties of genistein and galantamine have been widely demonstrated through different mechanisms; however, it is unknown a possible synergistic neuroprotective effect against Aß(1-42). In order to understand how genistein and galantamine combinations regulate the mechanisms of neuroprotection, we conducted a set of bioassays in vitro to evaluate cell viability, clonogenic survival, cell death, and anti-genotoxicity. Through molecular docking and therapeutic viability assays, we analyzed the inhibitory activity exerted by genistein on three major protein targets (AChE, BChE, and NMDA) involved in AD. The results showed that genistein and galantamine afforded significant protection at higher concentrations; however, combinations of sub-effective concentrations of both compounds provided marked neuroprotection when they were combined. In silico approaches showed that genistein has higher scores than the positive controls and low toxicity levels; nevertheless, the therapeutic viability indicated that unlike galantamine, genistein cannot undergo the action by P glycoprotein (PGP) and probably may be unable to cross the blood-brain barrier. In conclusion, our results show that genistein and galantamine exert neuroprotective by decreasing genotoxicity and cell death. In silico analysis, suggest that genistein modulates positively the expression of AChE, BChE, and NMDA. In this context, a combination of two or more drugs could inspire an attractive therapeutic strategy.

Computational design of new protein kinase 2 inhibitors for the treatment of inflammatory diseases using QSAR, pharmacophore-structure-based virtual screening, and molecular dynamics.

Receptor-interacting protein kinase 2 (RIPK2) plays an essential role in autoimmune response and is suggested as a target for inflammatory diseases. A pharmacophore model was built from a dataset with ponatinib (template) and 18 RIPK2 inhibitors selected from BindingDB database. The pharmacophore model validation was performed by multiple linear regression (MLR). The statistical quality of the model was evaluated by the correlation coefficient (R), squared correlation coefficient (R2), explanatory variance (adjusted R2), standard error of estimate (SEE), and variance ratio (F). The best pharmacophore model has one aromatic group (LEU24 residue interaction) and two hydrogen bonding acceptor groups (MET98 and TYR97 residues interaction), having a score of 24.739 with 14 aligned inhibitors, which were used in virtual screening via ZincPharmer server and the ZINC database (selected in function of the RMSD value). We determined theoretical values of biological activity (logRA) by MLR, pharmacokinetic and toxicology properties, and made molecular docking studies comparing binding affinity (kcal/mol) results with the most active compound of the study (ponatinib) and WEHI-345. Nine compounds from the ZINC database show satisfactory results, yielding among those selected, the compound ZINC01540228, as the most promising RIPK2 inhibitor. After binding free energy calculations, the following molecular dynamics simulations showed that the receptor protein's backbone remained stable after the introduction of ligands.

Aire Disruption Influences the Medullary Thymic Epithelial Cell Transcriptome and Interaction With Thymocytes.

The function of medullary thymic epithelial cells (mTECs) is associated with thymocyte adhesion, which is crucial for the negative selection of autoreactive thymocytes in the thymus. This process represents the root of central tolerance of self-components and prevents the onset of autoimmune diseases. Since thymic epithelia correspond to an important target of donor T cells during the onset of chronic graft-vs-host-disease, mTEC-thymocyte adhesion may have implications for alloimmunity. The Aire and Fezf2 genes function as transcriptome controllers in mTECs. The central question of this study is whether there is a mutual relationship between mTEC-thymocyte adhesion and the control of the mTEC transcriptome and whether Aire is involved in this process. Here, we show that in vitro mTEC-thymocyte adhesion causes transcriptome changes in mTECs and upregulates the transcriptional expression of Aire and Fezf2, as well as cell adhesion-related genes such as Cd80 or Tcf7, among others. Crispr-Cas9-mediated Aire gene disruption demonstrated that this gene plays a role in the process of mTEC-thymocyte adhesion. Consistent with the nuclear localization signal (NLS) encoded by Aire exon 3, which was targeted, we demonstrate that Aire KO-/- mTECs impair AIRE protein localization in the nucleus. Consequently, the loss of function of Aire reduced the ability of these cells to adhere to thymocytes. Their transcriptomes differed from their wild-type Aire+/+ counterparts, even during thymocyte adhesion. A set of mRNA isoforms that encode proteins involved in cell adhesion were also modulated during this process. This demonstrates that both thymocyte interactions and Aire influence transcriptome profiling of mTEC cells.