Polymorphism in F pocket affects peptide selection and stability of type 1 diabetes-associated HLA-B39 allotypes.

HLA-B*39:06, HLA-B*39:01, and HLA-B*38:01 are closely related HLA allotypes differentially associated with type 1 diabetes (T1D) risk and progression. B*39:06 is highly predisposing, while B*39:01 and B*38:01 are weakly predisposing and protective allotypes, respectively. Here, we aimed to decipher molecular mechanisms underlying the differential association of these allotypes with T1D pathogenesis. We addressed peptide binding and conformational stability of HLA-B allotypes using computational and experimental approaches. Computationally, we found that B*39:06 and B*39:01 allotypes had more rigid peptide-binding grooves and were more promiscuous in binding peptides than B*38:01. Peptidomes of B*39:06 and B*39:01 contained fewer strong binders and were of lower affinity than that of B*38:01. Experimentally, we demonstrated that B*39:06 and B*39:01 had a higher capacity to bind peptides and exit to the cell surface but lower surface levels and were degraded faster than B*38:01. In summary, we propose that promiscuous B*39:06 and B*39:01 may bind suboptimal peptides and transport them the cell surface, where such unstable complexes may contribute to the pathogenesis of T1D.

In Vitro, In Vivo, and In Silico Analysis of Pyraclostrobin and Cyprodinil and Their Mixture Reveal New Targets and Signaling Mechanisms.

Pyraclostrobin and cyprodinil are broad-spectrum fungicides that are used in crops to control diseases. However, they are excessively used and, as a result, end up in the environment and threaten human health and ecosystems. Hence, knowledge of their mechanisms of action is critical to revealing their environmental fate and negative effects and regulating their use. In the present study, we conducted a comprehensive study to show the adverse effects of pyraclostrobin, cyprodinil, and their mixture using zebrafish larvae and different cell lines. Several end points were investigated, including mortality, development, gene expression, reporter assays, and molecular docking simulations. We found that both compounds and their mixture caused developmental delays and mortality in zebrafish, with a higher effect displayed by pyraclostrobin. Both compounds altered the expression of genes involved in several signaling pathways, including oxidative stress and mitochondrial function, lipid and drug metabolisms, the cell cycle, DNA damage, apoptosis, and inflammation. A noteworthy result of this study is that cyprodinil and the mixture group acted as NFκB activators, while pyraclostrobin demonstrated antagonist activity. The AHR activity was also upregulated by cyprodinil and the mixture group; however, pyraclostrobin did not show any effect. For the first time, we also demonstrated that pyraclostrobin had androgen receptor antagonist activity.

Unraveling the Allosteric Communication Mechanisms in T-Cell Receptor-Peptide-Loaded Major Histocompatibility Complex Dynamics Using Molecular Dynamics Simulations: An Approach Based on Dynamic Cross Correlation Maps and Residue Interaction Energy Calculations.

Antigen presentation by major histocompatibility complex (MHC) proteins to T-cell receptors (TCRs) plays a crucial role in triggering the adaptive immune response. Most of our knowledge on TCR-peptide-loaded major histocompatibility complex (pMHC) interaction stemmed from experiments yielding static structures, yet the dynamic aspects of this molecular interaction are equally important to understand the underlying molecular mechanisms and to develop treatment strategies against diseases such as cancer and autoimmune diseases. To this end, computational biophysics studies including all-atom molecular dynamics simulations have provided useful insights; however, we still lack a basic understanding of an overall allosteric mechanism that results in conformational changes in the TCR and subsequent T-cell activation. Previous hydrogen-deuterium exchange and nuclear magnetic resonance studies provided clues regarding these molecular mechanisms, including global rigidification and allosteric effects on the constant domain of TCRs away from the pMHC interaction site. Here, we show that molecular dynamics simulations can be used to identify how this overall rigidification may be related to the allosteric communication within TCRs upon pMHC interaction via essential dynamics and nonbonded residue-residue interaction energy analyses. The residues taking part in the rigidification effect are highlighted with an intricate analysis on residue interaction changes, which lead to a detailed outline of the complex formation event. Our results indicate that residues of the Cß domain of TCRs show significant differences in their nonbonded interactions upon complex formation. Moreover, the dynamic cross correlations between these residues are also increased, in line with their nonbonded interaction energy changes. Altogether, our approach may be valuable for elucidating intramolecular allosteric changes in the TCR structure upon pMHC interaction in molecular dynamics simulations.

ProSNEx: a web-based application for exploration and analysis of protein structures using network formalism.

ProSNEx (Protein Structure Network Explorer) is a web service for construction and analysis of Protein Structure Networks (PSNs) alongside amino acid flexibility, sequence conservation and annotation features. ProSNEx constructs a PSN by adding nodes to represent residues and edges between these nodes using user-specified interaction distance cutoffs for either carbon-alpha, carbon-beta or atom-pair contact networks. Different types of weighted networks can also be constructed by using either (i) the residue-residue interaction energies in the format returned by gRINN, resulting in a Protein Energy Network (PEN); (ii) the dynamical cross correlations from a coarse-grained Normal Mode Analysis (NMA) of the protein structure; (iii) interaction strength. Upon construction of the network, common network metrics (such as node centralities) as well as shortest paths between nodes and k-cliques are calculated. Moreover, additional features of each residue in the form of conservation scores and mutation/natural variant information are included in the analysis. By this way, tool offers an enhanced and direct comparison of network-based residue metrics with other types of biological information. ProSNEx is free and open to all users without login requirement at http://prosnex-tool.com.

Identification of novel inhibitors of the ABC transporter BmrA.

The resistance of microbes to commonly used antibiotics has become a worldwide health problem. A major underlying mechanism of microbial antibiotic resistance is the export of drugs from bacterial cells. Drug efflux is mediated through the action of multidrug resistance efflux pumps located in the bacterial cell membranes. The critical role of bacterial efflux pumps in antibiotic resistance has directed research efforts to the identification of novel efflux pump inhibitors that can be used alongside antibiotics in clinical settings. Here, we aimed to find potential inhibitors of the archetypical ATP-binding cassette (ABC) efflux pump BmrA of Bacillus subtilis via virtual screening of the Mu.Ta.Lig. Chemotheca small molecule library. Molecular docking calculations targeting the nucleotide-binding domain of BmrA were performed using AutoDock Vina. Following a further drug-likeness filtering step based on Lipinski's Rule of Five, top 25 scorers were identified. These ligands were then clustered into separate groups based on their contact patterns with the BmrA nucleotide-binding domain. Six ligands with distinct contact patterns were used for further in vitro inhibition assays based on intracellular ethidium bromide accumulation. Using this methodology, we identified two novel inhibitors of BmrA from the Chemotheca small molecule library.

gRINN: a tool for calculation of residue interaction energies and protein energy network analysis of molecular dynamics simulations.

Atomistic molecular dynamics (MD) simulations generate a wealth of information related to the dynamics of proteins. If properly analyzed, this information can lead to new insights regarding protein function and assist wet-lab experiments. Aiming to identify interactions between individual amino acid residues and the role played by each in the context of MD simulations, we present a stand-alone software called gRINN (get Residue Interaction eNergies and Networks). gRINN features graphical user interfaces (GUIs) and a command-line interface for generating and analyzing pairwise residue interaction energies and energy correlations from protein MD simulation trajectories. gRINN utilizes the features of NAMD or GROMACS MD simulation packages and automatizes the steps necessary to extract residue-residue interaction energies from user-supplied simulation trajectories, greatly simplifying the analysis for the end-user. A GUI, including an embedded molecular viewer, is provided for visualization of interaction energy time-series, distributions, an interaction energy matrix, interaction energy correlations and a residue correlation matrix. gRINN additionally offers construction and analysis of Protein Energy Networks, providing residue-based metrics such as degrees, betweenness-centralities, closeness centralities as well as shortest path analysis. gRINN is free and open to all users without login requirement at http://grinn.readthedocs.io.

A computational docking study on the pH dependence of peptide binding to HLA-B27 sub-types differentially associated with ankylosing spondylitis.

A single amino acid difference (Asp116His), having a key role in a pathogenesis pathway, distinguishes HLA-B*27:05 and HLA-B*27:09 sub-types as associated and non-associated with ankylosing spondylitis, respectively. In this study, molecular docking simulations were carried out with the aim of comprehending the differences in the binding behavior of both alleles at varying pH conditions. A library of modeled peptides was formed upon single point mutations aiming to address the effect of 20 naturally occurring amino acids at the binding core peptide positions. For both alleles, computational docking was applied using Autodock 4.2. Obtained free energies of binding (FEB) were compared within the peptide library and between the alleles at varying pH conditions. The amino acid preferences of each position were studied enlightening the role of each on binding. The preferred amino acids for each position of pVIPR were found to be harmonious with experimental studies. Our results indicate that, as the pH is lowered, the capacity of HLA-B*27:05 to bind peptides in the library is largely lost. Hydrogen bonding analysis suggests that the interaction between the main anchor positions of pVIPR and their respective binding pocket residues are affected from the pH the most, causing an overall shift in the FEB profiles.

Unraveling the ligand specificity and promiscuity of the Staphylococcus aureus NorA efflux pump: a computational study.

Staphylococcus aureus, a gram-positive bacterial pathogen, develops antibiotic resistance partly through enhanced activity of transmembrane multi-drug efflux pump proteins like NorA. Being a prominent member of the Major Facilitator Superfamily (MFS), NorA transports various small molecules including hydrophilic fluoroquinolone antibiotics across the cell membrane. Intriguingly, NorA is inhibited by a structurally diverse set of small molecule inhibitors as well, indicating a highly promiscuous ligand/inhibitor recognition. Our study aims to elucidate the structural facets of this promiscuity. Known NorA inhibitors were grouped into five clusters based on chemical class and docked into ligand binding pockets on NorA conformations generated via molecular dynamics simulations. We discovered that several key residues, such as I23, E222, and F303, are involved in inhibitor binding. Additionally, residues I244, T223, F303, and F140 were identified as prominent in interactions with specific ligand clusters. Our findings suggest that NorA's substrate binding site, encompassing residues aiding ligand recognition based on chemical nature, facilitates the recognition of chemically diverse ligands. This insight into NorA's structural promiscuity in ligand recognition not only enhances understanding of antibiotic resistance mechanisms in S. aureus but also sets the stage for the development of more effective efflux pump inhibitors, vital for combating multidrug resistance.Communicated by Ramaswamy H. Sarma.

Non-apoptotic cell death induction via sapogenin based supramolecular particles.

The discovery of novel chemotherapeutics that act through different mechanisms is critical for dealing with tumor heterogeneity and therapeutic resistance. We previously reported a saponin analog (AG-08) that induces non-canonical necrotic cell death and is auspicious for cancer therapy. Here, we describe that the key element in triggering this unique cell death mechanism of AG-08 is its ability to form supramolecular particles. These self-assembled particles are internalized via a different endocytosis pathway than those previously described. Microarray analysis suggested that AG-08 supramolecular structures affect several cell signaling pathways, including unfolded protein response, immune response, and oxidative stress. Finally, through investigation of its 18 analogs, we further determined the structural features required for the formation of particulate structures and the stimulation of the unprecedented cell death mechanism of AG-08. The unique results of AG-08 indicated that supramolecular assemblies of small molecules are promising for the field of anticancer drug development, although they have widely been accepted as nuisance in drug discovery studies.

In silico and in vitro assessment of androgen receptor antagonists.

There is a growing concern for male reproductive health as studies suggest that there is a sharp increase in prostate cancer and other fertility related problems. Apart from lifestyle, pollutants are also known to negatively affect the reproductive system. In addition to many other compounds that have been shown to alter androgen signaling, several environmental pollutants are known to disrupt androgen signaling via binding to androgen receptor (AR) or indirectly affecting the androgen synthesis. We analyzed here the molecular mechanism of the interaction between the human AR Ligand Binding Domain (hAR-LBD) and two environmental pollutants, linuron (a herbicide) and procymidone (a pesticide), and compared with the steroid agonist dihydrotestosterone (DHT) and well-known hAR antagonists bicalutamide and enzalutamide. Using molecular docking and dynamics simulations, we showed that the co-activator interaction site of the hAR-LBD is disrupted in different ways by different ligands. Binding free energies of the ligands were also ordered in increasing order as follows: linuron, procymidone, DHT, bicalutamide, and enzalutamide. These data were confirmed by in vitro assays. Reporter assay with MDA-kb2 cells showed that linuron, procymidone, bicalutamide and enzalutamide can inhibit androgen mediated activation of luciferase activity. Gene expression analysis further showed that these compounds can inhibit the expression of prostate specific antigen (PSA) and microseminoprotein beta (MSMB) in prostate cell line LNCaP. Comparative analysis showed that procymidone is more potent than linuron in inhibiting AR activity. Furthermore, procymidone at 10 µM dose showed equivalent and higher activity to AR inhibitor enzalutamide and bicalutamide respectively.

Sequence-structure-function relationships in class I MHC: A local frustration perspective.

Class I Major Histocompatibility Complex (MHC) binds short antigenic peptides with the help of Peptide Loading Complex (PLC), and presents them to T-cell Receptors (TCRs) of cytotoxic T-cells and Killer-cell Immunglobulin-like Receptors (KIRs) of Natural Killer (NK) cells. With more than 10000 alleles, human MHC (Human Leukocyte Antigen, HLA) is the most polymorphic protein in humans. This allelic diversity provides a wide coverage of peptide sequence space, yet does not affect the three-dimensional structure of the complex. Moreover, TCRs mostly interact with HLA in a common diagonal binding mode, and KIR-HLA interaction is allele-dependent. With the aim of establishing a framework for understanding the relationships between polymorphism (sequence), structure (conserved fold) and function (protein interactions) of the human MHC, we performed here a local frustration analysis on pMHC homology models covering 1436 HLA I alleles. An analysis of local frustration profiles indicated that (1) variations in MHC fold are unlikely due to minimally-frustrated and relatively conserved residues within the HLA peptide-binding groove, (2) high frustration patches on HLA helices are either involved in or near interaction sites of MHC with the TCR, KIR, or tapasin of the PLC, and (3) peptide ligands mainly stabilize the F-pocket of HLA binding groove.

Computational investigation of peptide binding stabilities of HLA-B*27 and HLA-B*44 alleles.

Major Histocompatibility Complex (MHC) is a cell surface glycoprotein that binds to foreign antigens and presents them to T lymphocyte cells on the surface of Antigen Presenting Cells (APCs) for appropriate immune recognition. Recently, studies focusing on peptide-based vaccine design have allowed a better understanding of peptide immunogenicity mechanisms, which is defined as the ability of a peptide to stimulate CTL-mediated immune response. Peptide immunogenicity is also known to be related to the stability of peptide-loaded MHC (pMHC) complex. In this study, ENCoM server was used for structure-based estimation of the impact of single point mutations on pMHC complex stabilities. For this purpose, two human MHC molecules from the HLA-B*27 group (HLA-B*27:05 and HLA-B*27:09) in complex with four different peptides (GRFAAAIAK, RRKWRRWHL, RRRWRRLTV and IRAAPPPLF) and three HLA-B*44 molecules (HLA-B*44:02, HLA-B*44:03 and HLA-B*44:05) in complex with two different peptides (EEYLQAFTY and EEYLKAWTF) were analyzed. We found that the stability of pMHC complexes is dependent on both peptide sequence and MHC allele. Furthermore, we demonstrate that allele-specific peptide-binding preferences can be accurately revealed using structure-based computational methods predicting the effect of mutations on protein stability.

How do mutations and allosteric inhibitors modulate caspase-7 activity? A molecular dynamics study.

Caspases are members of a highly regulated aspartate-cysteine protease family which have important roles in apoptosis. Pharmaceutical studies focused on these molecules since they are involved in diseases such as cancer and neurodegenerative disorders. A small molecule which binds to the dimeric interface away from the binding site induces a conformational change that resembles the pro-caspase form of the molecule by shifting loop positions. The fluctuation mechanisms caused by mutations or binding of a ligand can explain the key mechanism for the function of that molecule. In this study, we performed molecular dynamics simulations on wild-type and mutated structures (C290N, R187M, Y223A, G188L and G188P) as well as allosterically inhibited structure (DICA-bound caspase-7) to observe the effects of the single mutations on intrinsic dynamics. The results show that previously known changes in catalytic activity upon mutations or allosteric ligand binding are reflected in corresponding changes in the global dynamics of caspase-7. Communicated by Ramaswamy H. Sarma.

Computational characterization of residue couplings and micropolymorphism-induced changes in the dynamics of two differentially disease-associated human MHC class-I alleles.

Human major histocompatibility complex class I (MHC I) - or human leukocyte antigen (HLA) - proteins present intracellularly processed peptides to cytotoxic T lymphocytes in the adaptive immune response to pathogens. A high level of polymorphism in human MHC I proteins defines the peptide-binding specificity of thousands of different MHC alleles. However, polymorphism as well as the peptide ligand can also affect the global dynamics of the complex. In this study, we conducted classical molecular dynamics simulations of two HLA alleles, the ankylosing spondylitis (AS) associated/tapasin-dependent HLA-B*27:05 and nondisease-associated/tapasin-independent HLA-B*27:09, both in peptide-free forms as well as complex with four different peptides ligands. Our results indicate that in peptide-free form, the single amino acid substitution distinguishing the two alleles (D116H), leads to a weaker dynamic coupling of residues in the tapasin-dependent HLA-B*27:05. In peptide-bound form, several residues of the binding-groove, mostly in A and B pockets, show hinge-like behavior in the global motion of the MHC. Moreover, allele-dependent changes are shown in residue interactions, affecting the B-pocket as well as the beta-2-microglobulin (ß2m)-facing residues of the HLA chain.