Comparative Protein Structure Modeling Using MODELLER

B. Oliveira C, Fadhloun R, Torabi E, Shahbazzadeh D, Pooshang Bagheri K, Salgado Ferreira R, Borchani L. Antigenic and Substrate Preference Differences between Scorpion and Spider Dermonecrotic Toxins, a Comparative Investigation

The Hemiscorpius lepturus scorpion and brown spider Loxosceles intermedia represent a public health problem in Asia and America, respectively. Although distinct, these organisms contain similar toxins responsible for the principal clinical signs of envenomation. To better understand the properties of these toxins, we designed a study to compare recombinant Heminecrolysin (rHNC) and rLiD1, the major phospholipase D toxins of scorpion and spider venom, respectively. Using a competitive ELISA and a hemolytic inhibition test, we come to spot a cross reaction between scorpion and spider venoms along with an epitopic similarity between rHNC and rLiD1 associated with neutralizing antibodies. Results show that the ability of the rHNC to hydrolyze lysophosphatidylcholine (LPC) is equivalent to that of rLiD1 to hydrolyze sphingomyelin and vice-versa. rHNC exclusively catalyze transphosphatidylation of LPC producing cyclic phosphatidic acid (cPA). The in-silico analysis of hydrogen bonds between LPC and toxins provides a possible explanation for the higher transphosphatidylase activity of rHNC. Interestingly, for the first time, we reveal that lysophosphatidic acid (LPA) can be a substrate for both enzymes using cellular and enzymatic assays. The finding of the usage of LPA as a substrate as well as the formation of cPA as an end product could shed more light on the molecular basis of Hemiscorpius lepturus envenomation as well as on loxoscelism.

DSS1 allosterically regulates the conformation of the tower domain of BRCA2 that has dsDNA binding specificity for homologous recombination

DSS1 is an evolutionary conserved, small intrinsically disordered protein that regulates various cellular functions. Although several studies have elucidated the role of DSS1 in stabilizing BRCA2 and its importance in homologous recombination repair (HRR), yet the structural mechanism behind the stability and HRR remains elusive. In this study, using molecular dynamics simulation we show that DSS1 stabilizes linearly arranged DNA/DSS1 binding domains of BRCA2 with many native contacts. These contacts are absent in the complexes with two missense DSS1 mutants associated with germline breast cancer and somatic mouth carcinoma. Most importantly, our protein energy-based network models show DSS1 allosterically regulates the conformation of the distant tower domain of BRCA2 that has dsDNA binding specificity for HRR. We further postulate that the unique conformation of the tower domain with kinked-helices might be responsible for DNA strand invasion and initiation of HRR. Induced conformation of the tower domain by the kinked-helices is absent in the unbound BRCA2, as well as in the two mutant DSS1-BRCA2 complexes. This suggests that DSS1 allosterically regulates the tower domain conformations of BRCA2 that affects dsDNA binding, essential for HRR. Our results add a new dimension to the function of DSS1 and its role in regulating HRR.

Experimental and computational analysis of biased agonism on full-length and a C-terminally truncated adenosine A2A receptor

Biased agonism, the ability of agonists to differentially activate downstream signaling pathways by stabilizing specific receptor conformations, is a key issue for G protein-coupled receptor (GPCR) signaling. The C-terminal domain might influence this functional selectivity of GPCRs as it engages G proteins, GPCR kinases, β-arrestins, and several other proteins. Thus, the aim of this paper is to compare the agonist-dependent selectivity for intracellular pathways in a heterologous system expressing the full-length (A_2AR) and a C-tail truncated (A_2A^Δ40R lacking the last 40 amino acids) adenosine A_2A receptor, a GPCR that is already targeted in Parkinson’s disease using a first-in-class drug. Experimental data such as ligand binding, cAMP production, β-arrestin recruitment, ERK1/2 phosphorylation and dynamic mass redistribution assays, which correspond to different aspects of signal transduction, were measured upon the action of structurally diverse compounds (the agonists adenosine, NECA, CGS-21680, PSB-0777 and LUF-5834 and the SCH-58261 antagonist) in cells expressing A_2AR and A_2A^Δ40R. The results show that taking cAMP levels and the endogenous adenosine agonist as references, the main difference in bias was obtained with PSB-0777 and LUF-5834. The C-terminus is dispensable for both G-protein and β-arrestin recruitment and also for MAPK activation. Unrestrained molecular dynamics simulations, at the μs timescale, were used to understand the structural arrangements of the binding cavity, triggered by these chemically different agonists, facilitating G protein binding with different efficacy.

Mycobacterium tuberculosis Calcium Pump CtpF Modulates the Autophagosome in an mTOR-Dependent Manner

Calcium is a very important second messenger, whose concentration in various cellular compartments is under tight regulation. A disturbance in the levels of calcium in these compartments can play havoc in the cell, as it regulates various cellular processes by direct or indirect mechanisms. Here, we have investigated the functional importance of a calcium transporting P2A ATPase, CtpF of Mycobacterium tuberculosis (Mtb) in the pathogen's interaction with the host. Among its uncanny ways of dealing with the host with umpteen strategies for survival and persistence in humans, CtpF is identified as a new player. The levels of ctpF are upregulated in macrophage stresses like hypoxia, high nitric oxide levels and acidic pH. Using confocal microscopy and fluorimetry, we show that CtpF effluxes calcium in macrophages in early stages of Mtb infection. Downregulation of ctpF expression by conditional knockdown resulted in perturbation of host calcium levels and consequent decreased activation of mTOR. We present a mechanism how calcium efflux by the pathogen inhibits mTOR-dependent autophagy and enhances bacterial survival. Our work highlights how Mtb engages its metal efflux pumps to exploit host autophagic process for its proliferation.

IMPROvER: the Integral Membrane Protein Stability Selector

Identifying stabilising variants of membrane protein targets is often required for structure determination. Our new computational pipeline, the Integral Membrane Protein Stability Selector (IMPROvER) provides a rational approach to variant selection by employing three independent approaches: deep-sequence, model-based and data-driven. In silico tests using known stability data, and in vitro tests using three membrane protein targets with 7, 11 and 16 transmembrane helices provided measures of success. In vitro, individual approaches alone all identified stabilising variants at a rate better than expected by random selection. Low numbers of overlapping predictions between approaches meant a greater success rate was achieved (fourfold better than random) when approaches were combined and selections restricted to the highest ranked sites. The mix of information IMPROvER uses can be extracted for any helical membrane protein. We have developed the first general-purpose tool for selecting stabilising variants of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\upalpha$$\end{document}α-helical membrane proteins, increasing efficiency and reducing workload. IMPROvER can be accessed at http://improver.ddns.net/IMPROvER/.

Modeling, Mutagenesis and In-silico Structural Stability Assay of the Model of Superoxide Dismutase of Lactococcus Lactis Subsp. Cremoris MG1363

Background:

Characterizing the structure and function of superoxide dismutase (SOD), as an antioxidant enzyme providing opportunities for its application in food supplements.

Objectives:

In this study, the features of the Manganese-SOD of Lactococcus Lactis (SDLL), subsp. Cremoris MG1363, as probiotic bacteria, were determined on the basis of the computational methods.

Materials and Methods:

The protein’s physicochemical properties and the prediction of its secondary structure were determined via the ProtParam server and the GOR program respectively. Moreover, the 3D structures of the proteins were constructed via the MODELLER on the basis of the homology method and the threading algorithm MUSTER. On the other hand, the structural stability of the models was assayed under the quasi-physiological conditions by the GROMACS program via the GROMOS96 43a1 force field in Linux system. Finally, using the Molecular Docking Studies, the functionality features of the models were predicted through their affinity with the corresponding substrates.

Results:

The results revealed the physicochemical properties of the SDLL and a 3D model of a chain of the enzyme being similar to the SOD from the Bacillus Subtilis (SDBS). The model of the SDLL was checked for quality control purposes including the Ramachandran plot, the ERRAT and the Verifiy3D. The model was suggestive of the structural stability in quasi-physiological conditions; yet, less than that of the SDBS. Assessing the cause of the instability in the SDLL model was indicative of two unstable regions in the area far from the enzyme’s active position, they were considered suitable for mutagenesis. Accordingly, the loop substitution for the corresponding region of SDBS and the deletion of the loop positioned at the C terminal of SDLL resulted in a mutant of SDLL with more stability and appropriate affinity with the corresponding substrate.

Conclusion:

In general, the study provides a new model of SDLL with certain thermostable features, and a new mutant with suitable stability and functionality on the basis of the direct mutagenesis being used in different applications.

Intrinsically Disordered Regions of the DNA-Binding Domain of Human FoxP1 Facilitate Domain Swapping

Forkhead box P (FoxP) proteins are unique transcription factors that spatiotemporally regulate gene expression by tethering two chromosome loci together via functional domain-swapped dimers formed through their DNA-binding domains. Further, the differential kinetics on this dimerization mechanism underlie an intricate gene regulation network at physiological conditions. Nonetheless, poor understanding of the structural dynamics and steps of the association process impedes to link the functional domain swapping to human-associated diseases. Here, we have characterized the DNA-binding domain of human FoxP1 by integrating single-molecule Förster resonance energy transfer and hydrogen-deuterium exchange mass spectrometry data with molecular dynamics simulations. Our results confirm the formation of a previously postulated domain-swapped (DS) FoxP1 dimer in solution and reveal the presence of highly populated, heterogeneous, and locally disordered dimeric intermediates along the dimer dissociation pathway. The unique features of FoxP1 provide a glimpse of how intrinsically disordered regions can facilitate domain swapping oligomerization and other tightly regulated association mechanisms relevant in biological processes.

Potent novel inhibitors against hepatitis C virus NS3 (HCV NS3 GT-3a) protease domain

HCV NS3, a non-structural hepatitis C viral protein is used as one of the potential targets for inhibition by direct-acting antivirals. It is known that the success rate for HCV genotype-1 treatment remained very high, however, treatment of genotype-3a (GT-3a), is still quite challenging. In the current study, the HCV GT-3a full-length NS3 gene was amplified and sequenced. The complete nucleotide sequence was translated into the amino acid sequence and homology models of HCV-NS3 GT-3a were generated by HCV-NS3 genotype-1b as a template. The objective of the study was to screen novel therapeutic hits from large databases. For this aim, various small molecule databases including, BindingDB (∼45.000 compounds), NCI (∼265.000 compounds), and Specs-SC (∼212.000 compounds) were used. Firstly, all of the compounds were screened using binary-QSAR models from the MetaCore/MetaDrug server, and compounds were filtered based on therapeutic activity predictions by the anti-viral QSAR model. Filtered molecules were used in 26 different toxicity QSAR models and active non-toxic compounds were identified. These selected molecules were then used in docking and molecular dynamics (MD) simulations studies at the binding cavities of the NS3 protease domain of the GT-3a. Results were compared with known inhibitors and novel molecules are proposed against HCV-NS3 GT-3a. These molecules have high ligand efficiencies as compared to the reference molecules suggesting a better alternate to the existing suite of inhibitors. Thus, this study will be a step ahead in the development of new potential compounds as antiviral drugs for the GT-3a target.

Sodium-coupled glucose transport, the SLC5 family, and therapeutically relevant inhibitors: from molecular discovery to clinical application

Sodium glucose transporters (SGLTs) belong to the mammalian solute carrier family SLC5. This family includes 12 different members in human that mediate the transport of sugars, vitamins, amino acids, or smaller organic ions such as choline. The SLC5 family belongs to the sodium symporter family (SSS), which encompasses transporters from all kingdoms of life. It furthermore shares similarity to the structural fold of the APC (amino acid-polyamine-organocation) transporter family. Three decades after the first molecular identification of the intestinal Na+-glucose cotransporter SGLT1 by expression cloning, many new discoveries have evolved, from mechanistic analysis to molecular genetics, structural biology, drug discovery, and clinical applications. All of these advances have greatly influenced physiology and medicine. While SGLT1 is essential for fast absorption of glucose and galactose in the intestine, the expression of SGLT2 is largely confined to the early part of the kidney proximal tubules, where it reabsorbs the bulk part of filtered glucose. SGLT2 has been successfully exploited by the pharmaceutical industry to develop effective new drugs for the treatment of diabetic patients. These SGLT2 inhibitors, termed gliflozins, also exhibit favorable nephroprotective effects and likely also cardioprotective effects. In addition, given the recent finding that SGLT2 is also expressed in tumors of pancreas and prostate and in glioblastoma, this opens the door to potential new therapeutic strategies for cancer treatment by specifically targeting SGLT2. Likewise, further discoveries related to the functional association of other SGLTs of the SLC5 family to human pathologies will open the door to potential new therapeutic strategies. We furthermore hope that the herein summarized information about the physiological roles of SGLTs and the therapeutic benefits of the gliflozins will be useful for our readers to better understand the molecular basis of the beneficial effects of these inhibitors, also in the context of the tubuloglomerular feedback (TGF), and the renin-angiotensin system (RAS). The detailed mechanisms underlying the clinical benefits of SGLT2 inhibition by gliflozins still warrant further investigation that may serve as a basis for future drug development.

The Human TSHβ Subunit Proteins and Their Binding Sites on the TSH Receptor Using Molecular Dynamics Simulation

To gain further insight into the binding of the normal and variant human TSHβ subunits (TSHβ and TSHβv), we modeled the 2 monomeric proteins and studied their interaction with the TSH receptor ectodomain (TSHR-ECD) using molecular dynamics simulation Furthermore, analyzed their bioactivity in vitro using recombinant proteins to confirm that such binding was physiologically relevant. Examining the interaction of TSHβ and TSHβv with the TSHR-ECD model using molecular dynamic simulation revealed strong binding of these proteins to the receptor ECD. The specificity of TSHβ and TSHβv binding to the TSHR-ECD was examined by analyzing the hydrogen-bonding residues of these subunits to the FSH receptor ECD, indicating the inability of these molecules to bind to the FSH receptors. Furthermore, the modelling suggests that TSHβ and TSHβv proteins clasped the concave surface of the leucine rich region of the TSHR ECD in a similar way to the native TSH using dynamic hydrogen bonding. These mutually exclusive stable interactions between the subunits and ECD residues included some high-affinity contact sites corresponding to binding models of native TSH. Furthermore, we cloned TSHβ and TSHβv proteins using the entire coding ORF and purified the flag-tagged proteins. The expressed TSHβ subunit proteins retained bioactivity both in a coculture system as well as with immune-purified proteins. In summary, we showed that such interactions can result in a functional outcome and may exert physiological or pathophysiological effects in immune cells.

Computer-aided Discovery of a New Nav1.7 Inhibitor for Treatment of Pain and Itch

BACKGROUND

Voltage-gated sodium channel Nav1.7 has been validated as a perspective target for selective inhibitors with analgesic and anti-itch activity. The objective of this study was to discover new candidate compounds with Nav1.7 inhibitor properties. The authors hypothesized that their approach would yield at least one new compound that inhibits sodium currents in vitro and exerts analgesic and anti-itch effects in mice.

METHODS

In silico structure-based similarity search of 1.5 million compounds followed by docking to the Nav1.7 voltage sensor of Domain 4 and molecular dynamics simulation was performed. Patch clamp experiments in Nav1.7-expressing human embryonic kidney 293 cells and in mouse and human dorsal root ganglion neurons were conducted to test sodium current inhibition. Formalin-induced inflammatory pain model, paclitaxel-induced neuropathic pain model, histamine-induced itch model, and mouse lymphoma model of chronic itch were used to confirm in vivo activity of the selected compound.

RESULTS

After in silico screening, nine compounds were selected for experimental assessment in vitro. Of those, four compounds inhibited sodium currents in Nav1.7-expressing human embryonic kidney 293 cells by 29% or greater (P < 0.05). Compound 9 (3-(1-benzyl-1H-indol-3-yl)-3-(3-phenoxyphenyl)-N-(2-(pyrrolidin-1-yl)ethyl)propanamide, referred to as DA-0218) reduced sodium current by 80% with a 50% inhibition concentration of 0.74 μM (95% CI, 0.35 to 1.56 μM), but had no effects on Nav1.5-expressing human embryonic kidney 293 cells. In mouse and human dorsal root ganglion neurons, DA-0218 reduced sodium currents by 17% (95% CI, 6 to 28%) and 22% (95% CI, 9 to 35%), respectively. The inhibition was greatly potentiated in paclitaxel-treated mouse neurons. Intraperitoneal and intrathecal administration of the compound reduced formalin-induced phase II inflammatory pain behavior in mice by 76% (95% CI, 48 to 100%) and 80% (95% CI, 68 to 92%), respectively. Intrathecal administration of DA-0218 produced acute reduction in paclitaxel-induced mechanical allodynia, and inhibited histamine-induced acute itch and lymphoma-induced chronic itch.

CONCLUSIONS

This study's computer-aided drug discovery approach yielded a new Nav1.7 inhibitor that shows analgesic and anti-pruritic activity in mouse models.

Continuous Assembly of β-Roll Structures Is Implicated in the Type I-Dependent Secretion of Large Repeat-in-Toxins (RTX) Proteins

Repeats-in-Toxin (RTX) proteins of Gram-negative bacteria are excreted through the type I secretion system (T1SS) that recognizes non-cleavable C-terminal secretion signals. These are preceded by arrays of glycine and aspartate-rich nonapeptide repeats grouped by four to eight β strands into blocks that fold into calcium-binding parallel β-roll structures. The β-rolls are interspersed by linkers of variable length and sequence and the organization of multiple RTX repeat blocks within large RTX domains remains unknown. Here we examined the structure and function of the RTX domain of Bordetella pertussis adenylate cyclase toxin (CyaA) that is composed of five β-roll RTX blocks. We show that the non-folded RTX repeats maintain the stability of the CyaA polypeptide in the Ca²⁺-depleted bacterial cytosol and thereby enable its efficient translocation through the T1SS apparatus. The efficacy of secretion of truncated CyaA constructs was dictated by the number of retained RTX repeat blocks and depended on the presence of extracellular Ca²⁺ ions. We further describe the crystal structure of the RTX blocks IV-V of CyaA (CyaA_1372-1681) that consists of a contiguous assembly of two β-rolls that differs substantially from the arrangement of the RTX blocks observed in RTX lipases or other RTX proteins. These results provide a novel structural insight into the architecture of the RTX domains of large RTX proteins and support the "push-ratchet" mechanism of the T1SS-mediated secretion of very large RTX proteins.

The Molecular Basis of JAZ-MYC Coupling, a Protein-Protein Interface Essential for Plant Response to Stressors

The jasmonic acid (JA) signaling pathway is one of the primary mechanisms that allow plants to respond to a variety of biotic and abiotic stressors. Within this pathway, the JAZ repressor proteins and the basic helix-loop-helix (bHLH) transcription factor MYC3 play a critical role. JA is a volatile organic compound with an essential role in plant immunity. The increase in the concentration of JA leads to the decoupling of the JAZ repressor proteins and the bHLH transcription factor MYC3 causing the induction of genes of interest. The primary goal of this study was to identify the molecular basis of JAZ-MYC coupling. For this purpose, we modeled and validated 12 JAZ-MYC3 3D in silico structures and developed a molecular dynamics/machine learning pipeline to obtain two outcomes. First, we calculated the average free binding energy of JAZ-MYC3 complexes, which was predicted to be -10.94 +/-2.67 kJ/mol. Second, we predicted which ones should be the interface residues that make the predominant contribution to the free energy of binding (molecular hotspots). The predicted protein hotspots matched a conserved linear motif SL••FL•••R, which may have a crucial role during MYC3 recognition of JAZ proteins. As a proof of concept, we tested, both in silico and in vitro, the importance of this motif on PEAPOD (PPD) proteins, which also belong to the TIFY protein family, like the JAZ proteins, but cannot bind to MYC3. By mutating these proteins to match the SL••FL•••R motif, we could force PPDs to bind the MYC3 transcription factor. Taken together, modeling protein-protein interactions and using machine learning will help to find essential motifs and molecular mechanisms in the JA pathway.

Does SARS-CoV-2 Bind to Human ACE2 More Strongly Than Does SARS-CoV?

The 2019 novel coronavirus (SARS-CoV-2) epidemic, which was first reported in December 2019 in Wuhan, China, was declared a pandemic by the World Health Organization in March 2020. Genetically, SARS-CoV-2 is closely related to SARS-CoV, which caused a global epidemic with 8096 confirmed cases in more than 25 countries from 2002 to 2003. Given the significant morbidity and mortality rate, the current pandemic poses a danger to all of humanity, prompting us to understand the activity of SARS-CoV-2 at the atomic level. Experimental studies have revealed that spike proteins of both SARS-CoV-2 and SARS-CoV bind to angiotensin-converting enzyme 2 (ACE2) before entering the cell for replication. However, the binding affinities reported by different groups seem to contradict each other. Wrapp et al. (Science2020, 367, 1260–1263) showed that the spike protein of SARS-CoV-2 binds to the ACE2 peptidase domain (ACE2-PD) more strongly than does SARS-CoV, and this fact may be associated with a greater severity of the new virus. However, Walls et al. (Cell2020, 181, 281–292) reported that SARS-CoV-2 exhibits a higher binding affinity, but the difference between the two variants is relatively small. To understand the binding mechnism and experimental results, we investigated how the receptor binding domain (RBD) of SARS-CoV (SARS-CoV-RBD) and SARS-CoV-2 (SARS-CoV-2-RBD) interacts with a human ACE2-PD using molecular modeling. We applied a coarse-grained model to calculate the dissociation constant and found that SARS-CoV-2 displays a 2-fold higher binding affinity. Using steered all-atom molecular dynamics simulations, we demonstrate that, like a coarse-grained simulation, SARS-CoV-2-RBD was associated with ACE2-PD more strongly than was SARS-CoV-RBD, as evidenced by a higher rupture force and larger pulling work. We show that the binding affinity of both viruses to ACE2 is driven by electrostatic interactions.

In vitro and in silico analysis of L. donovani enoyl acyl carrier protein reductase - A possible drug target

The emergence of increased resistance to the available drugs has created a situation that demands to find out more specific molecular drug targets for Leishmaniasis. The enoyl acyl carrier protein reductase (ENR), a regulatory enzyme in type II fatty acid synthesis, was confirmed as a novel drug target and triclosan as its specific inhibitor in many microorganisms. In this study, the triclosan was tested for the leishmanicidal property against Leishmania donovani (L. donovani) and the results of in vitro and ex vivo drug assays on promastigotes and amastigotes showed that triclosan possessed antileishmanial activity with a half minimal inhibitory concentration (IC₅₀) of 30 µM. Consequently, adopting in silico approach, we have tested the triclosan's ability to bind with the L. donovani enoyl acyl carrier protein reductase (LdENR). The 3D structure of LdENR was modelled, triclosan and cofactors were docked in LdENR model and molecular dynamic simulations were performed to observe the protein-ligands interactions, stability, compactness and binding energy calculation of the ligands-LdENR complexes. The observation showed that triclosan stably interacted with LdENR in presence of both the cofactors (NADPH and NADH), however, simulation results favor NADH as a preferred co-factor for LdENR. These results support that the reduction of L. donovani growth in the in vitro and ex vivo drug assays may be due to the interaction of triclosan with LdENR, which should be confirmed through enzymatic assays. The results of this study suggest that LdENR could be a potential drug target and triclosan as a lead for Leishmaniasis.Communicated by Ramaswamy H. Sarma.

Insight into natural inhibitors and bridging docking to dynamic simulation against sugar Isomerase (SIS) domain protein

The pathogen Legionella longbeachae is a causative agent of legionellosis. The antibiotic resistance is the major problem of this modern world. Thus, selective pressure warrants the need for identification of newer drug target. In current study, subtractive proteomics approach screen out SIS (sugar isomerase) domain protein as an attractive receptor molecule for rational drug design. This protein is involved in lipopolysaccharide biosynthesis and catalyzes the isomerization of sedoheptulose 7-phosphate in D-glycero-D-manno-heptose 7-phosphate. Molecular docking revealed compound 1 (2-(6-(N,N-dimethyl sulfamoyl)pipridin-4-yl)pyrazin-2-yl)imidazol-3-ium-1-ide) as the potent inhibitor having GOLD fitness score of 69. The complex is affirmed by half-site effect via simulation analysis. Complex stability was investigated via several approaches that follows dynamic simulation and binding energies. Trajectory analysis revealed slight change in ring positioning of inhibitor inside the active pocket during 130 ns (nanosecond). Interestingly, it was affirmed via binding interactions' density distribution. Hence, radial distribution function (RDF) inferred that SER55 and SER83 are the major residues that take part in hydrogen bonding and complex stability. Furthermore, an indigenously developed method axial frequency distribution (AFD) has revealed that ligand moved closer to the active site with both the residues SER55 and SER83 binding to the ligand. The phenomena was observed via rotating motion with respect to receptor center cavity. Thus, inhibitor movement towards allosteric site was observed at the end of simulations. Finally, binding free energy calculations by MMPB/GBSA predicts high compound affinity for the complex. Hence, findings from the current study will aid in the novel drug discovery and future experimental studies. Graphical abstract.

Mechanisms of pathogenesis of missense mutations on the KDM6A-H3 interaction in type 2 Kabuki Syndrome

Mutations in genes encoding for histone methylation proteins are associated with several developmental disorders. Among them, KDM6A is the disease causative gene of type 2 Kabuki Syndrome, a rare multisystem disease. While nonsense mutations and short insertions/deletions are known to trigger pathogenic mechanisms, the functional effects of missense mutations are still uncharacterized. In this study, we demonstrate that a selected set of missense mutations significantly hamper the interaction between KDM6A and the histone H3, by modifying the dynamics of the linker domain, and then causing a loss of function effect.

Genome-wide and structural analyses of pseudokinases encoded in the genome of Arabidopsis thaliana provide functional insights

Protein Kinase-Like Non-Kinases (PKLNKs), commonly known as "pseudokinases", are homologous to eukaryotic Ser/Thr/Tyr protein kinases (PKs) but lack the crucial aspartate residue in the catalytic loop, indispensable for phosphotransferase activity. Therefore, they are predicted to be "catalytically inactive" enzyme homologs. Analysis of protein-kinase like sequences from Arabidopsis thaliana led to the identification of more than 120 pseudokinases lacking catalytic aspartate, majority of which are closely related to the plant-specific receptor-like kinase family. These pseudokinases engage in different biological processes, enabled by their diverse domain architectures and specific subcellular localizations. Structural comparison of pseudokinases with active and inactive conformations of canonical PKs, belonging to both plant and animal origin, revealed unique structural differences. The currently available crystal structures of pseudokinases show that the loop topologically equivalent to activation segment of PKs adopts a distinct-folded conformation, packing against the pseudoenzyme core, in contrast to the extended and inhibitory geometries observed for active and inactive states, respectively, of catalytic PKs. Salt-bridge between ATP-binding Lys and DFG-Asp as well as hydrophobic interactions between the conserved nonpolar residue C-terminal to the equivalent DFG motif and nonpolar residues in C-helix mediate such a conformation in pseudokinases. This results in enhanced solvent accessibility of the pseudocatalytic loop in pseudokinases that can possibly serve as an interacting surface while associating with other proteins. Specifically, our analysis identified several residues that may be involved in pseudokinase regulation and hints at the repurposing of pseudocatalytic residues to achieve mechanistic control over noncatalytic functions of pseudoenzymes.

Serpins in Fasciola hepatica: insights into host-parasite interactions

Protease inhibitors play crucial roles in parasite development and survival, modulating the immune responses of their vertebrate hosts. Members of the serpin family are irreversible inhibitors of serine proteases and regulate systems related to defence against parasites. Limited information is currently available on protease inhibitors from the liver fluke Fasciola hepatica. In this study, we characterised four serpins from F. hepatica (FhS-1-FhS-4). Biochemical characterisation revealed that recombinant FhS-2 (rFhS) inhibits the activity of human neutrophil cathepsin G, while rFhS-4 inhibits the activity of bovine pancreatic chymotrypsin and cathepsin G. Consistent with inhibitor function profiling data, rFhS-4 inhibited cathepsin G-activated platelet aggregation in a dose-responsive manner.Similar to other serpins, rFhS2 and rFhS-4 bind to heparin with high affinity. Tissue localisation demonstrated that these serpins have different spatial distributions. FhS-2 is localised in the ovary, while FhS-4 was found in gut cells. Both of them co-localised in the spines within the tegument. These findings provide the basis for study of functional roles of these proteins as part of an immune evasion mechanism in the adult fluke, and in protection of eggs to ensure parasite life cycle continuity. Further understanding of serpins from the liver fluke may lead to the discovery of novel anti-parasitic interventions.

The development of the dog heartworm is highly sensitive to sterols which activate the orthologue of the nuclear receptor DAF-12

Prevention therapy against Dirofilaria immitis in companion animals is currently threatened by the emergence of isolates resistant to macrocyclic lactone anthelmintics. Understanding the control over developmental processes in D. immitis is important for elucidating new approaches to heartworm control. The nuclear receptor DAF-12 plays a role in the entry and exit of dauer stage in Caenorhabditis elegans and in the development of free-living infective third-stage larvae (iL3) of some Clade IV and V parasitic nematodes. We identified a DAF-12 ortholog in the clade III nematode D. immitis and found that it exhibited a much higher affinity for dafachronic acids than described with other nematode DAF-12 investigated so far. We also modelled the DimDAF-12 structure and characterized the residues involved with DA binding. Moreover, we showed that cholesterol derivatives impacted the molting process from the iL3 to the fourth-stage larvae. Since D. immitis is unable to synthesize cholesterol and only completes its development upon host infection, we hypothesize that host environment contributes to its further molting inside the host vertebrate. Our discovery contributes to a better understanding of the developmental checkpoints of D. immitis and offers new perspectives for the development of novel therapies against filarial infections.

In-silico analysis of Calcium Dependent Protein Kinase 6 of Cryptosporidium parvum through molecular modeling, docking, and dynamics simulation study

Calcium Dependent Protein Kinases are found in the Apicomplexan, algae, and plants; however, they are not reported in vertebrates and are regarded as excellent drug targets for pharmaceutical interventions. Calcium Dependent Protein Kinases of Cryptosporidium are probably involved in the regulation of invasion and egress process during the infection of the host cells. The previous study reported that after the Calcium Dependent Protein Kinase 1 gene, Calcium Dependent Protein Kinase 6 of Cryptosporidium parvum is expressed in all stages of the parasite (merozoites/schizonts as well as sexual stages) at a comparable level and makes it as a valid drug target. In this study, an attempt is made to address the similarity in sequences and phylogenetic study of Calcium Dependent Protein Kinase 6 (CDPK6) among Calcium Dependent Protein Kinases of Apicomplexans. Further, the three-dimensional structure determination of CDPK6 of C. parvum was performed through a molecular modeling approach followed by virtual screening of small-molecule inhibitors from different datasets. The best inhibitor from Tres Cantos Antimalarial Set with ID 11730 reported a binding affinity of -8.2 kcal/mol against CDPK6 of C. parvum. Furthermore, the reliability of the binding mode of the inhibitor is validated through a complex molecular dynamics simulation study for a time interval of 100 ns. The simulation study advocates that the inhibitor Tres Cantos Antimalarial Set_11730 formed a stable interaction with the predicted active site residues and can be considered for industrial pharmaceutical research in future.Communicated by Ramaswamy H. Sarma.

MDockPeP: A Web Server for Blind Prediction of Protein-Peptide Complex Structures

Protein-peptide interactions mediate a wide range of important cellular tasks. In silico prediction of protein-peptide complex structure is highly desirable for mechanistic investigation of these processes and for therapeutic design. Recently, we developed a docking-based method for predicting protein-peptide complex structures, which starts with the peptide sequence and globally docks the all-atom, flexible peptide onto the protein structure. The produced modes are then evaluated with a statistical potential-based scoring function. The method has been implemented into an online server, MDockPeP server, which is freely available at http://zougrouptoolkit.missouri.edu/mdockpep . The server can be used for protein-peptide complex structure prediction. The server can also be used for initial-stage sampling of the protein-peptide binding modes for computational-demanding simulation or docking methods.

Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins

The current outbreak of coronavirus disease-2019 (COVID-19) poses unprecedented challenges to global health¹. The new coronavirus responsible for this outbreak-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-shares high sequence identity to SARS-CoV and a bat coronavirus, RaTG13². Although bats may be the reservoir host for a variety of coronaviruses^3,4, it remains unknown whether SARS-CoV-2 has additional host species. Here we show that a coronavirus, which we name pangolin-CoV, isolated from a Malayan pangolin has 100%, 98.6%, 97.8% and 90.7% amino acid identity with SARS-CoV-2 in the E, M, N and S proteins, respectively. In particular, the receptor-binding domain of the S protein of pangolin-CoV is almost identical to that of SARS-CoV-2, with one difference in a noncritical amino acid. Our comparative genomic analysis suggests that SARS-CoV-2 may have originated in the recombination of a virus similar to pangolin-CoV with one similar to RaTG13. Pangolin-CoV was detected in 17 out of the 25 Malayan pangolins that we analysed. Infected pangolins showed clinical signs and histological changes, and circulating antibodies against pangolin-CoV reacted with the S protein of SARS-CoV-2. The isolation of a coronavirus from pangolins that is closely related to SARS-CoV-2 suggests that these animals have the potential to act as an intermediate host of SARS-CoV-2. This newly identified coronavirus from pangolins-the most-trafficked mammal in the illegal wildlife trade-could represent a future threat to public health if wildlife trade is not effectively controlled.

Bioactivities of synthetic salicylate-substituted carboxyl (E)-β-Farnesene derivatives as ecofriendly agrochemicals and their binding mechanism with potential targets in aphid olfactory system

BACKGROUND

The aphid alarm pheromone, (E)-β-farnesene (EβF), is a natural product secreted from the aphid cornicle as a signal to warn companions of danger. Odorant binding proteins (OBPs) are the vital targets in insect signal transduction pathways. To improve bioactivity of EβF as more economic and stable aphid control agents, EβF derivatives containing an active substructure, salicylic acid moiety, were designed, synthesized, and evaluated for their bioactivities in a structure-function study under laboratory conditions.

RESULTS

EβF derivatives, (E)-3,7-dimethylocta-2,6-dien-1-yl-2-hydroxy-3-methylbenzoate and (E)-3,7-dimethylocta-2,6-dien-1-yl-2-hydroxy-3-methoxybenzoate showed outstanding aphid-repellent activity at a dose of 5 μg against Acyrthosiphon pisum (repellency proportions of 67.3% and 71.2%, respectively) and Myzus persicae (repellency proportions of 80.0% and 74.4%, respectively) in laboratory. EβF and most of its derivatives bound strongly to ApisOBP9 with a higher affinity than those of the reported potential targets AphisOBP3 and ApisOBP7. The binding affinities to these three ApisOBPs were generally consistent with the in vivo aphid-repellent activity. A molecular docking study suggested that the hydrophobic effect was crucial for the interactions between the derivatives and the OBPs.

CONCLUSION

New EβF derivatives containing salicylic acid moiety and their repellent activity, binding mechanism with three potential OBPs are presented. A new OBP, ApisOBP9, was characterized as a potential EβF and EβF derivatives binding protein for the first time. The hydrophobic nature of these analogues is responsible for their activity. Two analogues 3b and 3e with outstanding aphid-repellent activity could be new leads for aphid control agents.

Site of Action of Brain Neurosteroid Pregnenolone Sulfate at the N-Methyl-D-Aspartate Receptor

N-methyl-D-aspartate receptor (NMDAR) hypofunction has been implicated in several neurodevelopmental disorders. NMDAR function can be augmented by positive allosteric modulators, including endogenous compounds, such as cholesterol and neurosteroid pregnenolone sulfate (PES). Here we report that PES accesses the receptor via the membrane, and its binding site is different from that of cholesterol. Alanine mutagenesis has identified residues that disrupt the steroid potentiating effect at the rat GluN1 (G638; I642) and GluN2B (W559; M562; Y823; M824) subunit. Molecular dynamics simulation indicates that, in the absence of PES, the GluN2B M1 helix residue W559 interacts with the M4 helix residue M824. In the presence of PES, the M1 and M4 helices of agonist-activated receptor rearrange, forming a tighter interaction with the GluN1 M3 helix residues G638 and I642. This stabilizes the open-state position of the GluN1 M3 helices. Together, our data identify a likely binding site for the NMDAR-positive allosteric modulator PES and describe a novel molecular mechanism by which NMDAR activity can be augmented.SIGNIFICANCE STATEMENT There is considerable interest in drugs that enhance NMDAR function and could compensate for receptor hypofunction associated with certain neuropsychiatric disorders. Positive allosteric modulators of NMDARs include an endogenous neurosteroid pregnenolone sulfate (PES), but the binding site of PES on the NMDAR and the molecular mechanism of potentiation are unknown. We use patch-clamp electrophysiology in combination with mutagenesis and in silico modeling to describe the interaction of PES with the NMDAR. Our data indicate that PES binds to the transmembrane domain of the receptor at a discrete group of residues at the GluN2B membrane helices M1 and M4 and the GluN1 helix M3, and that PES potentiates NMDAR function by stabilizing the open-state position of the GluN1 M3 helices.

Glutantβase: a database for improving the rational design of glucose-tolerant β-glucosidases

Β-glucosidases are key enzymes used in second-generation biofuel production. They act in the last step of the lignocellulose saccharification, converting cellobiose in glucose. However, most of the β-glucosidases are inhibited by high glucose concentrations, which turns it a limiting step for industrial production. Thus, β-glucosidases have been targeted by several studies aiming to understand the mechanism of glucose tolerance, pH and thermal resistance for constructing more efficient enzymes. In this paper, we present a database of β-glucosidase structures, called Glutantβase. Our database includes 3842 GH1 β-glucosidase sequences collected from UniProt. We modeled the sequences by comparison and predicted important features in the 3D-structure of each enzyme. Glutantβase provides information about catalytic and conserved amino acids, residues of the coevolution network, protein secondary structure, and residues located in the channel that guides to the active site. We also analyzed the impact of beneficial mutations reported in the literature, predicted in analogous positions, for similar enzymes. We suggested these mutations based on six previously described mutants that showed high catalytic activity, glucose tolerance, or thermostability (A404V, E96K, H184F, H228T, L441F, and V174C). Then, we used molecular docking to verify the impact of the suggested mutations in the affinity of protein and ligands (substrate and product). Our results suggest that only mutations based on the H228T mutant can reduce the affinity for glucose (product) and increase affinity for cellobiose (substrate), which indicates an increment in the resistance to product inhibition and agrees with computational and experimental results previously reported in the literature. More resistant β-glucosidases are essential to saccharification in industrial applications. However, thermostable and glucose-tolerant β-glucosidases are rare, and their glucose tolerance mechanisms appear to be related to multiple and complex factors. We gather here, a set of information, and made predictions aiming to provide a tool for supporting the rational design of more efficient β-glucosidases. We hope that Glutantβase can help improve second-generation biofuel production. Glutantβase is available at http://bioinfo.dcc.ufmg.br/glutantbase .

Molecular Topology of RNA Polymerase I Upstream Activation Factor

Upstream activation factor (UAF) is a multifunctional transcription factor in Saccharomyces cerevisiae that plays dual roles in activating RNA polymerase I (Pol I) transcription and repression of Pol II. For Pol I, UAF binds to a specific upstream element in the ribosomal DNA (rDNA) promoter and interacts with two other Pol I initiation factors, the TATA-binding protein (TBP) and core factor (CF). We used an integrated combination of chemical cross-linking mass spectrometry (CXMS), molecular genetics, protein biochemistry, and structural modeling to understand the topological framework responsible for UAF complex formation. Here, we report the molecular topology of the UAF complex, describe new structural and functional domains that play roles in UAF complex integrity, assembly, and biological function, and provide roles for previously identified UAF domains that include the Rrn5 SANT and histone fold domains. We highlight the role of new domains in Uaf30 that include an N-terminal winged helix domain and a disordered tethering domain as well as a BORCS6-like domain found in Rrn9. Together, our results reveal a unique network of topological features that coalesce around a histone tetramer-like core to form the dual-function UAF complex.

Structure-based drug design, synthesis and screening of MmaA1 inhibitors as novel anti-TB agents

Tuberculosis is one of the leading causes of death across the world. The treatment regimens for tuberculosis are well established, but still the control of the disease faces many challenges such as lengthy treatment protocols, drug resistance and toxicity. In the present work, mycolic acid methyl transferase (MmaA1), a protein involved in the maturation of mycolic acids in the biochemical pathway of the Mycobacterium, was studied for novel drug discovery. The homology model of the MmaA1 protein was built and validated by using computational techniques. The MmaA1 protein has 286 amino acid residues consisting of 10 α-helices and 7 β-sheets. The active site of the MmaA1 protein was identified using CASTp, SiteMap and PatchDock. Virtual screening studies were performed with two small molecule ligand databases: Asinex synergy and Diverse_Elite_Gold_Platinum databases having a total of 43,446 molecules and generated 1,30,814 conformers against the predicted and validated active site of the MmaA1 protein. Binding analysis showed that the residues ASP 19, PHE 22, TRP 30, TYR 32, TRP 74 and ALA 77 of MmaA1 protein have consistent interactions with the ligands. The hit ligands were further filtered by in silico ADME properties to eliminate potentially toxic molecules. Of the top 10 molecules, 3-(2-morpholinoacetamido)-N-(1,4-dihydro-4-oxoquinazolin-6-yl) benzamide was synthesised and screened for in vitro anti-TB activity against Mtb H37Rv using MABA assay. The compound and its intermediates exhibited good in vitro anti-TB activity which can be taken up for future lead optimisation studies. Structure based virtual screening study was performed using a validated homology model against small molecules from two virtual compound libraries. Synthesised the lead compound 3-(2-morpholinoacetamido)-N-(1,4-dihydro-4-oxoquinazolin-6-yl)benzamide obtained from virtual screening. In vitro activity against Mtb H37Rv has given a promising result.

Computer-Aided Ligand Discovery for Estrogen Receptor Alpha

Breast cancer (BCa) is one of the most predominantly diagnosed cancers in women. Notably, 70% of BCa diagnoses are Estrogen Receptor α positive (ERα+) making it a critical therapeutic target. With that, the two subtypes of ER, ERα and ERβ, have contrasting effects on BCa cells. While ERα promotes cancerous activities, ERβ isoform exhibits inhibitory effects on the same. ER-directed small molecule drug discovery for BCa has provided the FDA approved drugs tamoxifen, toremifene, raloxifene and fulvestrant that all bind to the estrogen binding site of the receptor. These ER-directed inhibitors are non-selective in nature and may eventually induce resistance in BCa cells as well as increase the risk of endometrial cancer development. Thus, there is an urgent need to develop novel drugs with alternative ERα targeting mechanisms that can overcome the limitations of conventional anti-ERα therapies. Several functional sites on ERα, such as Activation Function-2 (AF2), DNA binding domain (DBD), and F-domain, have been recently considered as potential targets in the context of drug research and discovery. In this review, we summarize methods of computer-aided drug design (CADD) that have been employed to analyze and explore potential targetable sites on ERα, discuss recent advancement of ERα inhibitor development, and highlight the potential opportunities and challenges of future ERα-directed drug discovery.

A variant in IL6ST with a selective IL-11 signaling defect in human and mouse

The GP130 cytokine receptor subunit encoded by IL6ST is the shared receptor for ten cytokines of the IL-6 family. We describe a homozygous non-synonymous variant in IL6ST (p.R281Q) in a patient with craniosynostosis and retained deciduous teeth. We characterize the impact of the variant on cytokine signaling in vitro using transfected cell lines as well as primary patient-derived cells and support these findings using a mouse model with the corresponding genome-edited variant Il6st p.R279Q. We show that human GP130 p.R281Q is associated with selective loss of IL-11 signaling without affecting IL-6, IL-27, OSM, LIF, CT1, CLC, and CNTF signaling. In mice Il6st p.R279Q lowers litter size and causes facial synostosis and teeth abnormalities. The effect on IL-11 signaling caused by the GP130 variant shows incomplete penetrance but phenocopies aspects of IL11RA deficiency in humans and mice. Our data show that a genetic variant in a pleiotropic cytokine receptor can have remarkably selective defects.

Linking Biochemical and Structural States of SERCA: Achievements, Challenges, and New Opportunities

Sarcoendoplasmic reticulum calcium ATPase (SERCA), a member of the P-type ATPase family of ion and lipid pumps, is responsible for the active transport of Ca²⁺ from the cytoplasm into the sarcoplasmic reticulum lumen of muscle cells, into the endoplasmic reticulum (ER) of non-muscle cells. X-ray crystallography has proven to be an invaluable tool in understanding the structural changes of SERCA, and more than 70 SERCA crystal structures representing major biochemical states (defined by bound ligand) have been deposited in the Protein Data Bank. Consequently, SERCA is one of the best characterized components of the calcium transport machinery in the cell. Emerging approaches in the field, including spectroscopy and molecular simulation, now help integrate and interpret this rich structural information to understand the conformational transitions of SERCA that occur during activation, inhibition, and regulation. In this review, we provide an overview of the crystal structures of SERCA, focusing on identifying metrics that facilitate structure-based categorization of major steps along the catalytic cycle. We examine the integration of crystallographic data with different biophysical approaches and computational methods to link biochemical and structural states of SERCA that are populated in the cell. Finally, we discuss the challenges and new opportunities in the field, including structural elucidation of functionally important and novel regulatory complexes of SERCA, understanding the structural basis of functional divergence among homologous SERCA regulators, and bridging the gap between basic and translational research directed toward therapeutic modulation of SERCA.

On the Conformational Dynamics of β-Amyloid Forming Peptides: A Computational Perspective

Understanding the conformational dynamics of proteins and peptides involved in important functions is still a difficult task in computational structural biology. Because such conformational transitions in β-amyloid (Aβ) forming peptides play a crucial role in many neurological disorders, researchers from different scientific fields have been trying to address issues related to the folding of Aβ forming peptides together. Many theoretical models have been proposed in the recent years for studying Aβ peptides using mathematical, physicochemical, and molecular dynamics simulation, and machine learning approaches. In this article, we have comprehensively reviewed the developmental advances in the theoretical models for Aβ peptide folding and interactions, particularly in the context of neurological disorders. Furthermore, we have extensively reviewed the advances in molecular dynamics simulation as a tool used for studying the conversions between polymorphic amyloid forms and applications of using machine learning approaches in predicting Aβ peptides and aggregation-prone regions in proteins. We have also provided details on the theoretical advances in the study of Aβ peptides, which would enhance our understanding of these peptides at the molecular level and eventually lead to the development of targeted therapies for certain acute neurological disorders such as Alzheimer's disease in the future.

CRiSP: accurate structure prediction of disulfide-rich peptides with cystine-specific sequence alignment and machine learning

MOTIVATION

High-throughput sequencing discovers many naturally occurring disulfide-rich peptides or cystine-rich peptides (CRPs) with diversified bioactivities. However, their structure information, which is very important to peptide drug discovery, is still very limited.

RESULTS

We have developed a CRP-specific structure prediction method called Cystine-Rich peptide Structure Prediction (CRiSP), based on a customized template database with cystine-specific sequence alignment and three machine-learning predictors. The modeling accuracy is significantly better than several popular general-purpose structure modeling methods, and our CRiSP can provide useful model quality estimations.

AVAILABILITY AND IMPLEMENTATION

The CRiSP server is freely available on the website at http://wulab.com.cn/CRISP.

CONTACT

wuyd@pkusz.edu.cn or jiangfan@pku.edu.cn.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A structurally minimized yet fully active insulin based on cone-snail venom insulin principles

Human insulin and its current therapeutic analogs all show propensity, albeit varyingly, to self-associate into dimers and hexamers, which delays their onset of action and makes blood glucose management difficult for people with diabetes. Recently, we described a monomeric, insulin-like peptide in cone-snail venom with moderate human insulin-like bioactivity. Here, with insights from structural biology studies, we report the development of mini-Ins-a human des-octapeptide insulin analog-as a structurally minimal, full-potency insulin. Mini-Ins is monomeric and, despite the lack of the canonical B-chain C-terminal octapeptide, has similar receptor binding affinity to human insulin. Four mutations compensate for the lack of contacts normally made by the octapeptide. Mini-Ins also has similar in vitro insulin signaling and in vivo bioactivities to human insulin. The full bioactivity of mini-Ins demonstrates the dispensability of the PheB24-PheB25-TyrB26 aromatic triplet and opens a new direction for therapeutic insulin development.

Effect of the Met148Leu mutation on the structure and dynamics of the rusticyanin protein from Acidithiobacillus sp. FJ2

The rusticyanin protein, a blue monomeric copper protein type-1, is one of the main components in the iron-electron transfer chain of the Acidithiobacillus ferrooxidans, and is the product of the rus gene expression. Herein, first the bacterial DNA of Acidithiobacillus sp. FJ2 was extracted. Then, the rus gene sequence and the sequence amino acid rusticyanin protein were determined. The Met148Leu mutation increased the oxidase activity of the rusticyanin protein, thereby enhancing the efficiency of the bioleaching process by bacteria Acidithiobacillus ferroxidans. Met148Leu mutation was created in the rusticyanin protein, then molecular dynamics (MD) simulations and structural analysis were performed. The MD analysis of the wild-type and mutant protein demonstrated a slight instability in the mutant protein and significant instability in the active site of the mutant protein. The usefulness of this study is the genetic manipulation of the native Acidithiobacillus sp. FJ2 bacterium, which can boost the bioleaching efficiency of the bacterium to some extent, and investigating its effects on the structure of a mutant protein using computational methods.

Knowledge-based structural models of SARS-CoV-2 proteins and their complexes with potential drugs

The World Health Organization (WHO) has declared the coronavirus disease 2019 (COVID‐19) caused by the novel coronavirus SARS‐CoV‐2 a pandemic. There is, however, no confirmed anti‐COVID‐19 therapeutic currently. In order to assist structure‐based discovery efforts for repurposing drugs against this disease, we constructed knowledge‐based models of SARS‐CoV‐2 proteins and compared the ligand molecules in the template structures with approved/experimental drugs and components of natural medicines. Our theoretical models suggest several drugs, such as carfilzomib, sinefungin, tecadenoson, and trabodenoson, that could be further investigated for their potential for treating COVID‐19.

In Silico Study and Optimization of Bacillus megaterium alpha-Amylases Production Obtained from Honey Sources

This study aimed to screen alpha-amylase producing microorganisms from honey as a low water activity medium, a suitable source for selecting stable and cost-beneficial bacterial enzyme production systems. Plackett-Burman method was used to select twelve effective factors including pH, inoculum size, temperature, time, corn starch, KH₂PO₄, peptone, MgSO₄, CaCl₂, NaCl, glycerin, and yeast extract concentrations on bacterial alpha-amylases production yield. The Box-Behnken method was utilized to optimize the level of selected significant factors. The stability of bacterial alpha-amylases was also determined in low pH and high-temperature conditions. In addition, in silico study was used to create the alpha-amylase structure and study the stability in high-temperature and low water available condition. Among all isolated and characterized microorganisms, Bacillus megaterium produced the highest amount of alpha-amylases. The in silico data showed the enzyme 3D structure similarity to alpha-amylase from Halothermothrix orenii and highly negative charge amino acids on its surface caused the enzyme activity and stability in low water conditions. Based on Box-Behnken results, the temperature 35 °C, pH 6 and starch 40 g/l were determined as the optimum level of significant factors to achieve the highest alpha-amylases unit (101.44 U/ml). This bacterial alpha-amylases enzyme showed stability at pH 5 and a range of temperatures from 40 to 60 °C that indicates this enzyme may possess the potential for using in industrial processes.

An Evolutionary Marker of the Ribokinase Superfamily Is Responsible for Zinc-Mediated Regulation of Human Pyridoxal Kinase

The ribokinase superfamily catalyzes the phosphorylation of a vast diversity of substrates, and its members are characterized by the conservation of a common structural fold along with highly conserved sequence motifs responsible for phosphoryl transfer (GXGD) and stabilization of the metal-nucleotide complex (NXXE). Recently, a third motif (HXE) exclusive from ADP-dependent enzymes was identified, with its glutamic acid participating in water-mediated interactions with the metal-nucleotide complex and in stabilization of the ternary complex during catalysis. In this work, we bioinformatically determine that the aspartic acid of another motif (DPV), exclusively found in hydroxyethyl thiazole (THZK), hydroxymethyl pyrimidine (HMPK) and pyridoxal kinases (PLK), is structurally equivalent to the acidic residue in the HXE motif. Moreover, this residue is highly conserved among all ribokinase superfamily members. To determine whether the functional role of the DPV motif is similar to the HXE motif, we employed molecular dynamics simulations using crystal structures of phosphoryl donor substrate-complexed THZK and PLK, showing that its aspartic acid participated in water-mediated or direct interactions with the divalent metal of the metal-nucleotide complex. Lastly, enzyme kinetic assays on human PLK, an enzyme that utilizes zinc, showed that site-directed mutagenesis of the aspartic acid from the DPV motif abolishes the inhibition of this enzyme by increasing free zinc concentrations. Altogether, our results highlight that the DPV and HXE motifs are evolutionary markers of the functional and structural divergence of the ribokinase superfamily and evidence the role of the DPV motif in the interaction with both free and nucleotide-complexed divalent metals in the binding site of these enzymes.

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.

Salt Enhances the Thermostability of Enteroviruses by Stabilizing Capsid Protein Interfaces

Enteroviruses are common agents of infectious disease that are spread by the fecal-oral route. They are readily inactivated by mild heat, which causes the viral capsid to disintegrate or undergo conformational change. While beneficial for the thermal treatment of food or water, this heat sensitivity poses challenges for the stability of enterovirus vaccines. The thermostability of an enterovirus can be modulated by the composition of the suspending matrix, though the effects of the matrix on virus stability are not understood. Here, we determined the thermostability of four enterovirus strains in solutions with various concentrations of NaCl and different pH values. The experimental findings were combined with molecular modeling of the protein interaction forces at the pentamer and the protomer interfaces of the viral capsids. While pH only had a modest effect on thermostability, increasing NaCl concentrations raised the breakpoint temperatures of all viruses tested by up to 20°C. This breakpoint shift could be explained by an enhancement of the van der Waals attraction forces at the two protein interfaces. In comparison, the (net repulsive) electrostatic interactions were less affected by NaCl. Depending on the interface considered, the breakpoint temperature shifted by 7.5 or 5.6°C per 100-kcal/(mol·Å) increase in protein interaction force.IMPORTANCE The genus Enterovirus encompasses important contaminants of water and food (e.g., coxsackieviruses), as well as viruses of acute public health concern (e.g., poliovirus). Depending on the properties of the surrounding matrix, enteroviruses exhibit different sensitivities to heat, which in turn influences their persistence in the environment, during food treatment, and during vaccine storage. Here, we determined the effect of NaCl and pH on the heat stability of different enteroviruses and related the observed effects to changes in protein interaction forces in the viral capsid. We demonstrate that NaCl renders enteroviruses thermotolerant and that this effect stems from an increase in van der Waals forces at different protein subunits in the viral capsid. This work sheds light on the mechanism by which salt enhances virus stability.

Temporal Regulation of Dendritic Spines Through NrCAM-Semaphorin3F Receptor Signaling in Developing Cortical Pyramidal Neurons

Neuron-glial related cell adhesion molecule NrCAM is a newly identified negative regulator of spine density that genetically interacts with Semaphorin3F (Sema3F), and is implicated in autism spectrum disorders (ASD). To investigate a role for NrCAM in spine pruning during the critical adolescent period when networks are established, we generated novel conditional, inducible NrCAM mutant mice (Nex1Cre-ERT2: NrCAMflox/flox). We demonstrate that NrCAM functions cell autonomously during adolescence in pyramidal neurons to restrict spine density in the visual (V1) and medial frontal cortex (MFC). Guided by molecular modeling, we found that NrCAM promoted clustering of the Sema3F holoreceptor complex by interfacing with Neuropilin-2 (Npn2) and PDZ scaffold protein SAP102. NrCAM-induced receptor clustering stimulated the Rap-GAP activity of PlexinA3 (PlexA3) within the holoreceptor complex, which in turn, inhibited Rap1-GTPase and inactivated adhesive β1 integrins, essential for Sema3F-induced spine pruning. These results define a developmental function for NrCAM in Sema3F receptor signaling that limits dendritic spine density on cortical pyramidal neurons during adolescence.

Isolation of SARS-CoV-2-related coronavirus from Malayan pangolins

The current outbreak of coronavirus disease-2019 (COVID-19) poses unprecedented challenges to global health¹. The new coronavirus responsible for this outbreak-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-shares high sequence identity to SARS-CoV and a bat coronavirus, RaTG13². Although bats may be the reservoir host for a variety of coronaviruses^3,4, it remains unknown whether SARS-CoV-2 has additional host species. Here we show that a coronavirus, which we name pangolin-CoV, isolated from a Malayan pangolin has 100%, 98.6%, 97.8% and 90.7% amino acid identity with SARS-CoV-2 in the E, M, N and S proteins, respectively. In particular, the receptor-binding domain of the S protein of pangolin-CoV is almost identical to that of SARS-CoV-2, with one difference in a noncritical amino acid. Our comparative genomic analysis suggests that SARS-CoV-2 may have originated in the recombination of a virus similar to pangolin-CoV with one similar to RaTG13. Pangolin-CoV was detected in 17 out of the 25 Malayan pangolins that we analysed. Infected pangolins showed clinical signs and histological changes, and circulating antibodies against pangolin-CoV reacted with the S protein of SARS-CoV-2. The isolation of a coronavirus from pangolins that is closely related to SARS-CoV-2 suggests that these animals have the potential to act as an intermediate host of SARS-CoV-2. This newly identified coronavirus from pangolins-the most-trafficked mammal in the illegal wildlife trade-could represent a future threat to public health if wildlife trade is not effectively controlled.

Secondary Structure and Contact Guided Differential Evolution for Protein Structure Prediction

Ab initio protein tertiary structure prediction is one of the long-standing problems in structural bioinformatics. With the help of residue-residue contact and secondary structure prediction information, the accuracy of ab initio structure prediction can be enhanced. In this study, an improved differential evolution with secondary structure and residue-residue contact information referred to as SCDE is proposed for protein structure prediction. In SCDE, two score models based on secondary structure and contact information are proposed, and two selection strategies, namely, secondary structure-based selection strategy and contact-based selection strategy, are designed to guide conformation space search. A probability distribution function is designed to balance these two selection strategies. Experimental results on a benchmark dataset with 28 proteins and four free model targets in CASP12 demonstrate that the proposed SCDE is effective and efficient.

The SARS-CoV-2 Exerts a Distinctive Strategy for Interacting with the ACE2 Human Receptor

The COVID-19 disease has plagued over 200 countries with over three million cases and has resulted in over 200,000 deaths within 3 months. To gain insight into the high infection rate of the SARS-CoV-2 virus, we compare the interaction between the human ACE2 receptor and the SARS-CoV-2 spike protein with that of other pathogenic coronaviruses using molecular dynamics simulations. SARS-CoV, SARS-CoV-2, and HCoV-NL63 recognize ACE2 as the natural receptor but present a distinct binding interface to ACE2 and a different network of residue–residue contacts. SARS-CoV and SARS-CoV-2 have comparable binding affinities achieved by balancing energetics and dynamics. The SARS-CoV-2–ACE2 complex contains a higher number of contacts, a larger interface area, and decreased interface residue fluctuations relative to the SARS-CoV–ACE2 complex. These findings expose an exceptional evolutionary exploration exerted by coronaviruses toward host recognition. We postulate that the versatility of cell receptor binding strategies has immediate implications for therapeutic strategies.

Innovative IgG Biomarkers Based on Phage Display Microbial Amyloid Mimotope for State and Stage Diagnosis in Alzheimer's Disease

An innovative approach to identify new conformational antigens of Aβ_1–42 recognized by IgG autoantibodies as biomarkers of state and stage in Alzheimer’s disease (AD) patients is described. In particular, through the use of bioinformatics modeling, conformational similarities between several Aβ_1–42 forms and other amyloid-like proteins with F1 capsular antigen (Caf1) of Yersinia pestis were first found. pVIII M13 phage display libraries were then screened against YPF19, anti-Caf1 monoclonal antibody, and IgGs of AD patients, in alternate biopanning cycles of a so-called “double binding” selection. From the selected phage clones, one, termed 12III1, was found to be able to prevent in vitro Aβ_1–42-induced cytotoxicity in SH-SY5Y cells, as well as to promote disaggregation of preformed fibrils, to a greater extent with respect to wild-type phage (pC89). IgG levels detected by 12III1 provided a significant level of discrimination between diseased and nondemented subjects, as well as a good correlation with the state progression of the disease. These results give significant impact in AD state and stage diagnosis, paving the way for the development not only for an innovative blood diagnostic assay for AD precise diagnosis, progressive clinical assessment, and screening but also for new effective treatments.

Specifying the role of BAK1-interacting receptor-like kinase 3 in brassinosteroid signaling

Brassinosteroids (BR) are involved in the control of several developmental processes ranging from root elongation to senescence and adaptation to environmental cues. Thus, BR perception and signaling have to be precisely regulated. One regulator is BRI1-associated kinase 1 (BAK1)-interacting receptor-like kinase 3 (BIR3). In the absence of BR, BIR3 forms complexes with BR insensitive 1 (BRI1) and BAK1. However, the biophysical and energetic requirements for complex formation in the absence of the ligand have yet to be determined. Using computational modeling, we simulated the potential complexes between the cytoplasmic domains of BAK1, BRI1 and BIR3. Our calculations and experimental data confirm the interaction of BIR3 with BAK1 and BRI1, with the BAK1 BIR3 interaction clearly favored. Furthermore, we demonstrate that BIR3 and BRI1 share the same interaction site with BAK1. This suggests a competition between BIR3 and BRI1 for binding to BAK1, which results in preferential binding of BIR3 to BAK1 in the absence of the ligand thereby preventing the active participation of BAK1 in BR signaling. Our model also suggests that BAK1 and BRI1 can interact even while BAK1 is in complex with BIR3 at an additional binding site of BAK1 that does not allow active BR signaling.

Carbohydrate-Binding Capability and Functional Conformational Changes of AbnE, an Arabino-oligosaccharide Binding Protein

ABC importers are membrane proteins responsible for the transport of nutrients into the cells of prokaryotes. Although the structures of ABC importers vary, all contain four conserved domains: two nucleotide-binding domains (NBDs), which bind and hydrolyze ATP, and two transmembrane domains (TMDs), which help translocate the substrate. ABC importers are also dependent on an additional protein component, a high-affinity substrate-binding protein (SBP) that specifically binds the target ligand for delivery to the appropriate ABC transporter. AbnE is a SBP belonging to the ABC importer for arabino-oligosaccharides in the Gram-positive thermophilic bacterium Geobacillus stearothermophilus. Using isothermal titration calorimetry (ITC), purified AbnE was shown to bind medium-sized arabino-oligosaccharides, in the range of arabino-triose (A3) to arabino-octaose (A8), all with K_d values in the nanomolar range. We describe herein the 3D structure of AbnE in its closed conformation in complex with a wide range of arabino-oligosaccharide substrates (A2-A8). These structures provide the basis for the detailed structural analysis of the AbnE-sugar complexes, and together with complementary quantum chemical calculations, site-specific mutagenesis, and isothermal titration calorimetry (ITC) experiments, provide detailed insights into the AbnE-substrate interactions involved. Small-angle X-ray scattering (SAXS) experiments and normal mode analysis (NMA) are used to study the conformational changes of AbnE, and these data, taken together, suggest clues regarding its binding mode to the full ABC importer.

Peptide cargo tunes a network of correlated motions in human leucocyte antigens

Most biomolecular interactions are typically thought to increase the (local) rigidity of a complex, for example, in drug‐target binding. However, detailed analysis of specific biomolecular complexes can reveal a more subtle interplay between binding and rigidity. Here, we focussed on the human leucocyte antigen (HLA), which plays a crucial role in the adaptive immune system by presenting peptides for recognition by the αβ T‐cell receptor (TCR). The role that the peptide plays in tuning HLA flexibility during TCR recognition is potentially crucial in determining the functional outcome of an immune response, with obvious relevance to the growing list of immunotherapies that target the T‐cell compartment. We have applied high‐pressure/temperature perturbation experiments, combined with molecular dynamics simulations, to explore the drivers that affect molecular flexibility for a series of different peptide–HLA complexes. We find that different peptide sequences affect peptide–HLA flexibility in different ways, with the peptide cargo tuning a network of correlated motions throughout the pHLA complex, including in areas remote from the peptide‐binding interface, in a manner that could influence T‐cell antigen discrimination.