Statistical potentials for fold assessment

Heterologous expression, characterisation and 3D-structural insights of GH18 chitinases derived from sago palm (Metroxylon sagu)

Although plants don't have chitins, they produce chitinases to protect themselves from biotic and abiotic stressors. There are two forms of chitinases found in organisms: glycosyl hydrolase 18 (GH18) and 19 (GH19) families. Plant GH19 chitinases are well known for their role in protecting against pathogens, but the roles of GH18 chitinases have not been fully elucidated. This study aimed to produce and characterise two recombinant GH18 chitinases from Metroxylon sagu. Two GH18 chitinase genes, MsChi1 and MsChi2, were identified, with nucleotide sequences of 1009 and 1308 bp, respectively. The proteins encoded by MsChi1 and MsChi2 genes were single polypeptide chains of 310 and 300 amino acids with predicted molecular masses of 31.21 and 30.15 kDa, respectively. Both cDNAs were cloned and expressed in the GS115 strain of Pichia pastoris. Recombinant MsChi1 and MsChi2 exhibited optimal activity at 60 °C with acidic pH 4.0 and 5.0, respectively. Both recombinant enzymes could hydrolyze synthetic and natural substrates (colloidal chitin). rMsChi1 preferred 4-nitrophenol N,N'-diacetyl-β-D chitobioside, while rMsChi2 preferred 4-nitrophenol N,N',N″-triacetyl-β-D chitotriose, suggesting they might function as exochitinase and endochitinase, respectively. They also demonstrated antifungal activities against tested fungi. Homology modeling indicated ASP and GLU as essential residues for proton donation and acceptance.

Lipid turnover through lipophagy in the newly identified extremophilic green microalga Chlamydomonas urium

Autophagy is a central degradative pathway highly conserved among eukaryotes, including microalgae, which remains unexplored in extremophilic organisms. In this study, we described and characterized autophagy in the newly identified extremophilic green microalga Chlamydomonas urium, which was isolated from an acidic environment. The nuclear genome of C. urium was sequenced, assembled and annotated in order to identify autophagy-related genes. Transmission electron microscopy, immunoblotting, metabolomic and photosynthetic analyses were performed to investigate autophagy in this extremophilic microalga. The analysis of the C. urium genome revealed the conservation of core autophagy-related genes. We investigated the role of autophagy in C. urium by blocking autophagic flux with the vacuolar ATPase inhibitor concanamycin A. Our results indicated that inhibition of autophagic flux in this microalga resulted in a pronounced accumulation of triacylglycerols and lipid droplets (LDs). Metabolomic and photosynthetic analyses indicated that C. urium cells with impaired vacuolar function maintained an active metabolism. Such effects were not observed in the neutrophilic microalga Chlamydomonas reinhardtii. Inhibition of autophagic flux in C. urium uncovered an active recycling of LDs through lipophagy, a selective autophagy pathway for lipid turnover. This study provided the metabolic basis by which extremophilic algae are able to catabolize lipids in the vacuole.

GLUT-1DS resistant to ketogenic diet: from clinical feature to in silico analysis. An exemplificative case report with a literature review

Glucose transporter type 1 deficiency syndrome (GLUT-1DS) is characterized by alterations in glucose translocation through the blood-brain barrier (BBB) due to mutation involving the GLUT-1 transporter. The fundamental therapy is ketogenic diet (KD) that provide an alternative energetic substrate - ketone bodies that across the BBB via MCT-1 - for the brain. Symptoms are various and include intractable seizure, acquired microcephalia, abnormal ocular movement, movement disorder, and neurodevelopment delay secondary to an energetic crisis for persistent neuroglycopenia. KD is extremely effective in controlling epileptic seizures and has a positive impact on movement disorders and cognitive impairment. Cases of KD resistance are rare, and only a few of them are reported in the literature, all regarding seizure. Our study describes a peculiar case of GLUT-1DS due to a new deletion involving the first codon of SLC2A1 gene determining a loss of function with a resistance to KD admitted to hospital due to intractable episodes of dystonia. This patient presented a worsening of symptomatology at higher ketonemia values but without hyperketosis and showed a complete resolution of symptomatology while maintaining low ketonemia values. Our study proposes an in-silico genomic and proteomic analysis aimed at explaining the atypical response to KD exhibited by our patient. In this way, we propose a new clinical and research approach based on precision medicine and molecular modelling to be applied to patients with GLUT-1DS resistant to first-line treatment with ketogenic diet by in silico study of genetic and altered protein product.

Biochemical Characterization of a Novel Thermostable 1,4-α-Glucoamylase from Aspergillus brasiliensis Strain Isolated in the Brazilian Atlantic Forest

Glucoamylases are exo-enzymes that cleave the ends of the starch chain, releasing glucose units. In the current work, we described a novel 1,4-α-glucoamylase from an A. brasiliensis strain isolated from an environmental sample. The purified glucoamylase, GlaAb, has a molecular mass of 69 kDa and showed a starch binding domain. GlaAb showed a similar sequence to other fungal glucoamylases, and the molecular 3D model analysis of GlaAb suggests an overall structure as described in the literature, except by elongation in the loop connecting the 4th and 5th α-helices. The enzyme showed activity over a wide range of pH and temperature, with maximum activity at pH 4.5 and 60 °C. GlaAb was stable at 50 °C for 7 h, maintaining 67% residual activity, and it was not inhibited by glucose up to 0.1 M. The glucoamylase was 65% more active in the presence of Mn2+ and showed a Km of 2.21 mg mL-1, Vmax of 155 U mg-1, Kcat 179 s-1, and Kcat/Km 81.06 mg mL-1 s-1 using potato starch as substrate. The results obtained are promising and provide the basis for the development of applications of GlaAb in the industrial process.

Oleic Acid Could Act as a Channel Blocker in the Inhibition of nAChR: Insights from Molecular Dynamics Simulations

The activation of the muscular nicotinic acetylcholine receptor (nAChR) produces the opening of the channel, with the consequent increase in the permeability of cations, triggering an excitatory signal. Free fatty acids (FFA) are known to modulate the activity of the receptor as noncompetitive antagonists, acting at the membrane-AChR interface. We present molecular dynamics simulations of a model of nAChR in a desensitized closed state embedded in a lipid bilayer in which distinct membrane phospholipids were replaced by two different monounsaturated FFA that differ in the position of a double bond. This allowed us to detect and describe that the cis-18:1ω-9 FFA were located at the interface between the transmembrane segments of α2 and γ subunits diffused into the channel lumen with the consequent potential ability to block the channel to the passage of ions.

Atomistic ensemble of active SHP2 phosphatase

SHP2 phosphatase plays an important role in regulating several intracellular signaling pathways. Pathogenic mutations of SHP2 cause developmental disorders and are linked to hematological malignancies and cancer. SHP2 comprises two tandemly-arranged SH2 domains, a catalytic PTP domain, and a disordered C-terminal tail. Under physiological, non-stimulating conditions, the catalytic site of PTP is occluded by the N-SH2 domain, so that the basal activity of SHP2 is low. Whereas the autoinhibited structure of SHP2 has been known for two decades, its active, open structure still represents a conundrum. Since the oncogenic mutant SHP2E76K almost completely populates the active, open state, this mutant has been extensively studied as a model for activated SHP2. By molecular dynamics simulations and accurate explicit-solvent SAXS curve predictions, we present the heterogeneous atomistic ensemble of constitutively active SHP2E76K in solution, encompassing a set of conformational arrangements and radii of gyration in agreement with experimental SAXS data.

Pharmacoinformatic approach to identify potential phytochemicals against SARS-CoV-2 spike receptor-binding domain in native and variants of concern

Severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2) pathogenesis is initiated by the binding of SARS-CoV-2 spike (S) protein with the angiotensin-converting enzyme 2 receptor (ACE2R) on the host cell surface. The receptor-binding domain (RBD) of the S protein mediates the binding and is more prone to mutations resulting in the generation of different variants. Recently, molecules with the potential to inhibit the interaction of S protein with ACE2R have been of interest due to their therapeutic value. In this context, the present work was performed to identify potential RBD binders from the Indian medicinal plant's phytochemical database through virtual screening, molecular docking, and molecular dynamic simulation. Briefly, 1578 compounds filtered from 9596 phytochemicals were chosen for screening against the RBD of the native SARS-CoV-2 S protein. Based on the binding energy, the top 30 compounds were selected and re-docked individually against the native and five variants of concern (VOCs: alpha, beta, gamma, delta, and omicron) of SARS-CoV-2. Four phytochemicals, namely withanolide F, serotobenine, orobanchol, and gibberellin A51, were found to be potential RBD binders in native and all SARS-CoV-2 VOCs. Among the four, withanolide F exhibited lower binding energy (- 10.84 to - 8.56 kcal/mol) and better ligand efficiency (- 0.3 to - 0.25) against all forms of RBD and hence was subjected to a 100 ns MD simulation which confirmed its stringent binding to the RBDs in native and VOCs. The study prioritizes withanolide F as a prospective COVID-19 (Coronavirus disease) therapeutic agent based on the observations. It warrants deeper investigations into the four promising leads for understanding their precise therapeutic value.

In silico identification of five binding sites on the SARS-CoV-2 spike protein and selection of seven ligands for such sites

To contribute to the development of products capable of complexing with the SARS-CoV-2 spike protein, and thus preventing the virus from entering the host cell, this work aimed at discovering binding sites in the whole protein structure, as well as selecting substances capable of binding efficiently to such sites. Initially, the three-dimensional structure of the protein, with all receptor binding domains in the closed state, underwent blind docking with 38 substances potentially capable of binding to this protein according to the literature. This allowed the identification of five binding sites. Then, those substances with more affinities for these sites underwent pharmacophoric search in the ZINC15 database. The 14,329 substances selected from ZINC15 were subjected to docking to the five selected sites of the spike protein. The ligands with more affinities for the protein sites, as well as the selected sites themselves, were used in the de novo design of new ligands that were also docked to the binding sites of the protein. The best ligands, regardless of their origins, were used to form complexes with the spike protein, which were subsequently used in molecular dynamics simulations and calculations of ligands affinities to the protein through the molecular mechanics/Poisson-Boltzmann surface area method (MMPBSA). Seven substances with good affinities to the spike protein (-12.9 to -20.6 kcal/mol), satisfactory druggability (Bioavailability score: 0.17 to 0.55), and low acute toxicity to mice (LD50: 751 to 1421 mg/kg) were selected as potentially useful for the future development of new products to manage COVID-19 infections.Communicated by Ramaswamy H. Sarma.

Modeling the Characteristic Residues of Chlorophyll f Synthase (ChlF) from Halomicronema hongdechloris to Determine Its Reaction Mechanism

Photosystem II (PSII) is a quinone-utilizing photosynthetic system that converts light energy into chemical energy and catalyzes water splitting. PsbA (D1) and PsbD (D2) are the core subunits of the reaction center that provide most of the ligands to redox-active cofactors and exhibit photooxidoreductase activities that convert quinone and water into quinol and dioxygen. The performed analysis explored the putative uncoupled electron transfer pathways surrounding P680+ induced by far-red light (FRL) based on photosystem II (PSII) complexes containing substituted D1 subunits in Halomicronema hongdechloris. Chlorophyll f-synthase (ChlF) is a D1 protein paralog. Modeling PSII-ChlF complexes determined several key protein motifs of ChlF. The PSII complexes included a dysfunctional Mn4CaO5 cluster where ChlF replaced the D1 protein. We propose the mechanism of chlorophyll f synthesis from chlorophyll a via free radical chemistry in an oxygenated environment created by over-excited pheophytin a and an inactive water splitting reaction owing to an uncoupled Mn4CaO5 cluster in PSII-ChlF complexes. The role of ChlF in the formation of an inactive PSII reaction center is under debate, and putative mechanisms of chlorophyll f biosynthesis are discussed.

M Protein from Dengue virus oligomerizes to pentameric channel protein: in silico analysis study

The Dengue virus M protein is a 75 amino acid polypeptide with two helical transmembranes (TM). The TM domain oligomerizes to form an ion channel, facilitating viral release from the host cells. The M protein has a critical role in the virus entry and life cycle, making it a potent drug target. The oligomerization of the monomeric protein was studied using ab initio modeling and molecular dynamics (MD) simulation in an implicit membrane environment. The representative structures obtained showed pentamer as the most stable oligomeric state, resembling an ion channel. Glutamic acid, threonine, serine, tryptophan, alanine, isoleucine form the pore-lining residues of the pentameric channel, conferring an overall negative charge to the channel with approximate length of 51.9 Å. Residue interaction analysis (RIN) for M protein shows that Ala94, Leu95, Ser112, Glu124, and Phe155 are the central hub residues representing the physicochemical interactions between domains. The virtual screening with 165 different ion channel inhibitors from the ion channel library shows monovalent ion channel blockers, namely lumacaftor, glipizide, gliquidone, glisoxepide, and azelnidipine to be the inhibitors with high docking scores. Understanding the three-dimensional structure of M protein will help design therapeutics and vaccines for Dengue infection.

Ensemble-based molecular docking and spectrofluorometric analysis of interaction between cytotoxin and tumor necrosis factor receptor 1

Cytotoxin (CTX) is a three-finger toxin presents predominantly in cobra venom. The functional site of the toxin is located at its three hydrophobic loop tips. Its actual mechanism of cytotoxicity remains inconclusive as few conflicting hypotheses have been proposed in addition to direct cytolytic effects. The present work investigated the interaction between CTX and death receptor families via ensemble-based molecular docking and fluorescence titration analysis. Multiple sequence alignments of different CTX isoforms obtained a conserved CTX sequence. The three-dimensional structure of the conserved CTX was later determined using homology modelling, and its quality was validated. Ensemble-based molecular docking of CTX was performed with different death receptors, such as Fas-ligand and tumor necrosis factor receptor families. Our results showed that tumor necrosis factor receptor 1 (TNFR1) was the best receptor interacting with CTX attributed to the interaction of all three functional loops and evinced with low HADDOCK, Z-score and RMSD value. The interaction between CTX and TNFR1 was also supported by a concentration-dependent reduction of fluorescence intensity with increasing binding affinity. The possible intermolecular interactions between CTX and TNFR1 were Van der Waals forces and hydrogen bonding. Our findings suggest a possibility that CTX triggers apoptosis cell death through non-covalent interactions with TNFR1.Communicated by Ramaswamy H. Sarma.

Biochemical Characterization of an Endoglucanase GH7 from Thermophile Thermothielavioides terrestris Expressed on Aspergillus nidulans

Endoglucanases (EC 3.2.1.4) are important enzymes involved in the hydrolysis of cellulose, acting randomly in the β-1,4-glycosidic bonds present in the amorphous regions of the polysaccharide chain. These biocatalysts have been classified into 14 glycosyl hydrolase (GH) families. The GH7 family is of particular interest since it may act on a broad range of substrates, including cellulose, β-glucan, and xylan, an attractive feature for biotechnological applications, especially in the renewable energy field. In the current work, a gene from the thermophilic fungus Thermothielavioides terrestris, encoding an endoglucanase GH7 (TtCel7B), was cloned in the secretion vector pEXPYR and transformed into the high-protein-producing strain Aspergillus nidulans A773. Purified TtCel7B has a molecular weight of approximately 66 kDa, evidenced by SDS-PAGE. Circular dichroism confirmed the high β-strand content consistent with the canonical GH7 family β-jellyroll fold, also observed in the 3D homology model of TtCel7B. Biochemical characterization assays showed that TtCel7B was active over a wide range of pH values (3.5–7.0) and temperatures (45–70 °C), with the highest activity at pH 4.0 and 65 °C. TtCel7B also was stable over a wide range of pH values (3.5–9.0), maintaining more than 80% of its activity after 24 h. The KM and Vmax values in low-viscosity carboxymethylcellulose were 9.3 mg mL−1 and 2.5 × 104 U mg−1, respectively. The results obtained in this work provide a basis for the development of applications of recombinant TtCel7B in the renewable energy field.

Analysis of microbiota-host communication mediated by butyrate in Atlantic Salmon

Butyrate is a microbiota-produced metabolite, sensed by host short-chain fatty acid receptors FFAR2 (Gpr43), FFAR3 (Gpr41), HCAR2 (Gpr109A), and Histone deacetylase (HDAC) that promotes microbiota-host crosstalk. Butyrate influences energy uptake, developmental and immune response in mammals. This microbial metabolite is produced by around 79 anaerobic genera present in the mammalian gut, yet little is known about the role of butyrate in the host-microbiota interaction in salmonid fish. To further our knowledge of this interaction, we analyzed the intestinal microbiota and genome of Atlantic salmon (Salmo salar), searching for butyrate-producing genera and host butyrate receptors. We identified Firmicutes, Proteobacteria, and Actinobacteria as the main butyrate-producing bacteria in the salmon gut microbiota. In the Atlantic salmon genome, we identified an expansion of genes orthologous to FFAR2 and HCAR2 receptors, and class I and IIa HDACs that are sensitive to butyrate. In addition, we determined the expression levels of orthologous of HCAR2 in the gut, spleen, and head-kidney, and FFAR2 in RTgutGC cells. The effect of butyrate on the Atlantic salmon immune response was evaluated by analyzing the pro and anti-inflammatory cytokines response in vitro in SHK-1 cells by RT-qPCR. Butyrate decreased the expression of the pro-inflammatory cytokine IL-1β and increased anti-inflammatory IL-10 and TGF-β cytokines. Butyrate also reduced the expression of interferon-alpha, Mx, and PKR, and decreased the viral load at a higher concentration (4 mM) in cells treated with this molecule before the infection with Infectious Pancreatic Necrosis Virus (IPNV) by mechanisms independent of FFAR2, FFAR3 and HCAR2 expression that probably inhibit HDAC. Moreover, butyrate modified phosphorylation of cytoplasmic proteins in RTgutGC cells. Our data allow us to infer that Atlantic salmon have the ability to sense butyrate produced by their gut microbiota via different specific targets, through which butyrate modulates the immune response of pro and anti-inflammatory cytokines and the antiviral response.

Comparative Surface Electrostatics and Normal Mode Analysis of High and Low Pathogenic H7N7 Avian Influenza Viruses

Influenza A viruses are rarely symptomatic in wild birds, while representing a higher threat to poultry and mammals, where they can cause a variety of symptoms, including death. H5 and H7 subtypes of influenza viruses are of particular interest because of their pathogenic potential and reported capacity to spread from poultry to mammals, including humans. The identification of molecular fingerprints for pathogenicity can help surveillance and early warning systems, which are crucial to prevention and protection from such potentially pandemic agents. In the past decade, comparative analysis of the surface features of hemagglutinin, the main protein antigen in influenza viruses, identified electrostatic fingerprints in the evolution and spreading of H5 and H9 subtypes. Electrostatic variation among viruses from avian or mammalian hosts was also associated with host jump. Recent findings of fingerprints associated with low and highly pathogenic H5N1 viruses, obtained by means of comparative electrostatics and normal modes analysis, prompted us to check whether such fingerprints can also be found in the H7 subtype. Indeed, evidence presented in this work showed that also in H7N7, hemagglutinin proteins from low and highly pathogenic strains present differences in surface electrostatics, while no meaningful variation was found in normal modes.

Design of Novel IRAK4 Inhibitors Using Molecular Docking, Dynamics Simulation and 3D-QSAR Studies

Treatment of several autoimmune diseases and types of cancer has been an intense area of research over the past two decades. Many signaling pathways that regulate innate and/or adaptive immunity, as well as those that induce overexpression or mutation of protein kinases, have been targeted for drug discovery. One of the serine/threonine kinases, Interleukin-1 Receptor Associated Kinase 4 (IRAK4) regulates signaling through various Toll-like receptors (TLRs) and interleukin-1 receptor (IL1R). It controls diverse cellular processes including inflammation, apoptosis, and cellular differentiation. MyD88 gain-of-function mutations or overexpression of IRAK4 has been implicated in various types of malignancies such as Waldenström macroglobulinemia, B cell lymphoma, colorectal cancer, pancreatic ductal adenocarcinoma, breast cancer, etc. Moreover, over activation of IRAK4 is also associated with several autoimmune diseases. The significant role of IRAK4 makes it an interesting target for the discovery and development of potent small molecule inhibitors. A few potent IRAK4 inhibitors such as PF-06650833, RA9 and BAY1834845 have recently entered phase I/II clinical trial studies. Nevertheless, there is still a need of selective inhibitors for the treatment of cancer and various autoimmune diseases. A great need for the same intrigued us to perform molecular modeling studies on 4,6-diaminonicotinamide derivatives as IRAK4 inhibitors. We performed molecular docking and dynamics simulation of 50 ns for one of the most active compounds of the dataset. We also carried out MM-PBSA binding free energy calculation to identify the active site residues, interactions of which are contributing to the total binding energy. The final 50 ns conformation of the most active compound was selected to perform dataset alignment in a 3D-QSAR study. Generated RF-CoMFA (q² = 0.751, ONC = 4, r² = 0.911) model revealed reasonable statistical results. Overall results of molecular dynamics simulation, MM-PBSA binding free energy calculation and RF-CoMFA model revealed important active site residues of IRAK4 and necessary structural properties of ligand to design more potent IRAK4 inhibitors. We designed few IRAK4 inhibitors based on these results, which possessed higher activity (predicted pIC₅₀) than the most active compounds of the dataset selected for this study. Moreover, ADMET properties of these inhibitors revealed promising results and need to be validated using experimental studies.

The Flexible, Extended Coil of the PDZ-Binding Motif of the Three Deadly Human Coronavirus E Proteins Plays a Role in Pathogenicity

The less virulent human (h) coronaviruses (CoVs) 229E, NL63, OC43, and HKU1 cause mild, self-limiting respiratory tract infections, while the more virulent SARS-CoV-1, MERS-CoV, and SARS-CoV-2 have caused severe outbreaks. The CoV envelope (E) protein, an important contributor to the pathogenesis of severe hCoV infections, may provide insight into this disparate severity of the disease. We, therefore, generated full-length E protein models for SARS-CoV-1 and -2, MERS-CoV, HCoV-229E, and HCoV-NL63 and docked C-terminal peptides of each model to the PDZ domain of the human PALS1 protein. The PDZ-binding motif (PBM) of the SARS-CoV-1 and -2 and MERS-CoV models adopted a more flexible, extended coil, while the HCoV-229E and HCoV-NL63 models adopted a less flexible alpha helix. All the E peptides docked to PALS1 occupied the same binding site and the more virulent hCoV E peptides generally interacted more stably with PALS1 than the less virulent ones. We hypothesize that the increased flexibility of the PBM in the more virulent hCoVs facilitates more stable binding to various host proteins, thereby contributing to more severe disease. This is the first paper to model full-length 3D structures for both the more virulent and less virulent hCoV E proteins, providing novel insights for possible drug and/or vaccine development.

Structural and functional insights into ABHD5, a ligand-regulated lipase co-activator

Alpha/beta hydrolase domain-containing protein 5 (ABHD5) is a highly conserved protein that regulates various lipid metabolic pathways via interactions with members of the perilipin (PLIN) and Patatin-like phospholipase domain-containing protein (PNPLA) protein families. Loss of function mutations in ABHD5 result in Chanarin-Dorfman Syndrome (CDS), characterized by ectopic lipid accumulation in numerous cell types and severe ichthyosis. Recent data demonstrates that ABHD5 is the target of synthetic and endogenous ligands that might be therapeutic beneficial for treating metabolic diseases and cancers. However, the structural basis of ABHD5 functional activities, such as protein-protein interactions and ligand binding is presently unknown. To address this gap, we constructed theoretical structural models of ABHD5 by comparative modeling and topological shape analysis to assess the spatial patterns of ABHD5 conformations computed in protein dynamics. We identified functionally important residues on ABHD5 surface for lipolysis activation by PNPLA2, lipid droplet targeting and PLIN-binding. We validated the computational model by examining the effects of mutating key residues in ABHD5 on an array of functional assays. Our integrated computational and experimental findings provide new insights into the structural basis of the diverse functions of ABHD5 as well as pathological mutations that result in CDS.

Structure‐conditioned amino‐acid couplings: how contact geometry affects pairwise sequence preferences

Relating a protein's sequence to its conformation is a central challenge for both structure prediction and sequence design. Statistical contact potentials, as well as their more descriptive versions that account for side‐chain orientation and other geometric descriptors, have served as simplistic but useful means of representing second‐order contributions in sequence–structure relationships. Here we ask what happens when a pairwise potential is conditioned on the fully defined geometry of interacting backbones fragments. We show that the resulting structure‐conditioned coupling energies more accurately reflect pair preferences as a function of structural contexts. These structure‐conditioned energies more reliably encode native sequence information and more highly correlate with experimentally determined coupling energies. Clustering a database of interaction motifs by structure results in ensembles of similar energies and clustering them by energy results in ensembles of similar structures. By comparing many pairs of interaction motifs and showing that structural similarity and energetic similarity go hand‐in‐hand, we provide a tangible link between modular sequence and structure elements. This link is applicable to structural modeling, and we show that scoring CASP models with structured‐conditioned energies results in substantially higher correlation with structural quality than scoring the same models with a contact potential. We conclude that structure‐conditioned coupling energies are a good way to model the impact of interaction geometry on second‐order sequence preferences.

Structures of the TMC-1 complex illuminate mechanosensory transduction

The initial step in the sensory transduction pathway underpinning hearing and balance in mammals involves the conversion of force into the gating of a mechanosensory transduction channel¹. Despite the profound socioeconomic impacts of hearing disorders and the fundamental biological significance of understanding mechanosensory transduction, the composition, structure and mechanism of the mechanosensory transduction complex have remained poorly characterized. Here we report the single-particle cryo-electron microscopy structure of the native transmembrane channel-like protein 1 (TMC-1) mechanosensory transduction complex isolated from Caenorhabditis elegans. The two-fold symmetric complex is composed of two copies each of the pore-forming TMC-1 subunit, the calcium-binding protein CALM-1 and the transmembrane inner ear protein TMIE. CALM-1 makes extensive contacts with the cytoplasmic face of the TMC-1 subunits, whereas the single-pass TMIE subunits reside on the periphery of the complex, poised like the handles of an accordion. A subset of complexes additionally includes a single arrestin-like protein, arrestin domain protein (ARRD-6), bound to a CALM-1 subunit. Single-particle reconstructions and molecular dynamics simulations show how the mechanosensory transduction complex deforms the membrane bilayer and suggest crucial roles for lipid-protein interactions in the mechanism by which mechanical force is transduced to ion channel gating.

Structurally Guided Reprogramming of Valerenadiene Synthase

Valerena-1,10-diene synthase (VDS) catalyzes the conversion of the universal precursor farnesyl diphosphate into the unusual sesquiterpene valerena-1,10-diene (VLD), which possesses a unique isobutenyl substituent group. In planta, one of VLD's isobutenyl terminal methyl groups becomes oxidized to a carboxylic acid forming valerenic acid (VA), an allosteric modulator of the GABA_A receptor. Because a structure-activity relationship study of VA for its modulatory activity is desired, we sought to manipulate the VDS enzyme for the biosynthesis of structurally diverse scaffolds that could ultimately lead to the generation of VA analogues. Using three-dimensional structural homology models, phylogenetic sequence comparisons to well-characterized sesquiterpene synthases, and a substrate-active site contact mapping approach, the contributions of specific amino acid residues within or near the VDS active site to possible catalytic cascades for VLD and other sesquiterpene products were assessed. An essential role of Tyr535 in a germacrenyl route to VLD was demonstrated, while its contribution to a family of other sesquiterpenes derived from a humulyl route was not. No role for Cys415 or Cys452 serving as a proton donor to reaction intermediates in VLD biosynthesis was observed. However, a gatekeeper role for Asn455 in directing farnesyl carbocations down all-trans catalytic cascades (humulyl and germacrenyl routes) versus a cisoid cascade (nerolidyl route) was demonstrated. Altogether, these results have mapped residues that establish a context for the catalytic cascades operating in VDS and future manipulations for generating more structurally constrained scaffolds.

Second-generation non-hematopoietic erythropoietin-derived peptide for neuroprotection

Erythropoietin (EPO) is a well-known erythropoietic cytokine having a tissue-protective effect in various tissues against hypoxic stress, including the brain. Thus, its recombinants may function as neuroprotective compounds. However, despite considerable neuroprotective effects, the EPO-based therapeutic approach has side effects, including hyper-erythropoietic and tumorigenic effects. Therefore, some modified forms and derivatives of EPO have been proposed to minimize the side effects. In this study, we generated divergently modified new peptide analogs derived from helix C of EPO, with several amino acid replacements that interact with erythropoietin receptors (EPORs). This modification resulted in unique binding potency to EPOR. Unlike recombinant EPO, among the peptides, ML1-h3 exhibited a potent neuroprotective effect against oxidative stress without additional induction of cell-proliferation, owing to a differential activating mode of EPOR signaling. Furthermore, it inhibited neuronal death and brain injury under hypoxic stress in vitro and in an in vivo ischemic brain injury model. Therefore, the divergent modification of EPO-derivatives for affinity to EPOR could provide a basis for a more advanced and optimal neuroprotective strategy.

•

Short peptides derived from helix C of EPO have a neuroprotective effect.

•

Divergent modification of EPO-derived peptides has a differential affinity to EPOR.

•

ML1 and its analogs have differential cell protective and proliferative effects.

•

ML1-h3 protects neurons by suppressing in vitro oxidative stress.

•

ML1-h3 mitigates brain injury in the in vivo mouse ischemic model without hematopoietic effect.

Molecular Docking Simulation Studies Identifies Potential Natural Product Derived-Antiwolbachial Compounds as Filaricides against Onchocerciasis

Onchocerciasis is the leading cause of blindness and severe skin lesions which remain a major public health problem, especially in tropical areas. The widespread use of antibiotics and the long duration required for effective treatment continues to add to the increasing global menace of multi-resistant pathogens. Onchocerca volvulus harbors the endosymbiont bacteria Wolbachia, essential for the normal development of embryos, larvae and long-term survival of the adult worm, O. volvulus. We report here results of using structure-based drug design (SBDD) approach aimed at identifying potential novel Wolbachia inhibitors from natural products against the Wolbachia surface protein (WSP). The protein sequence of the WSP with UniProtKB identifier Q0RAI4 was used to model the three-dimensional (3D) structure via homology modelling techniques using three different structure-building algorithms implemented in Modeller, I-TASSER and Robetta. Out of the 15 generated models of WSP, one was selected as the most reasonable quality model which had 82, 15.5, 1.9 and 0.5% of the amino acid residues in the most favored regions, additionally allowed regions, generously allowed regions and disallowed regions, respectively, based on the Ramachandran plot. High throughput virtual screening was performed via Autodock Vina with a library comprising 42,883 natural products from African and Chinese databases, including 23 identified anti-Onchocerca inhibitors. The top six compounds comprising ZINC000095913861, ZINC000095486235, ZINC000035941652, NANPDB4566, acetylaleuritolic acid and rhemannic acid had binding energies of −12.7, −11.1, −11.0, −11, −10.3 and −9.5 kcal/mol, respectively. Molecular dynamics simulations including molecular mechanics Poisson-Boltzmann (MMPBSA) calculations reinforced the stability of the ligand-WSP complexes and plausible binding mechanisms. The residues Arg45, Tyr135, Tyr148 and Phe195 were predicted as potential novel critical residues required for ligand binding in pocket 1. Acetylaleuritolic acid and rhemannic acid (lantedene A) have previously been shown to possess anti-onchocercal activity. This warrants the need to evaluate the anti-WSP activity of the identified molecules. The study suggests the exploitation of compounds which target both pockets 1 and 2, by investigating their potential for effective depletion of Wolbachia. These compounds were predicted to possess reasonably good pharmacological profiles with insignificant toxicity and as drug-like. The compounds were computed to possess biological activity including antibacterial, antiparasitic, anthelmintic and anti-rickettsials. The six natural products are potential novel antiwolbachial agents with insignificant toxicities which can be explored further as filaricides for onchocerciasis.

A genetic approach reveals different modes of action of prefoldins

The prefoldin complex (PFDc) was identified in humans as a co-chaperone of the cytosolic chaperonin T-COMPLEX PROTEIN RING COMPLEX (TRiC)/CHAPERONIN CONTAINING TCP-1 (CCT). PFDc is conserved in eukaryotes and is composed of subunits PFD1-6, and PFDc-TRiC/CCT folds actin and tubulins. PFDs also participate in a wide range of cellular processes, both in the cytoplasm and in the nucleus, and their malfunction causes developmental alterations and disease in animals and altered growth and environmental responses in yeast and plants. Genetic analyses in yeast indicate that not all of their functions require the canonical complex. The lack of systematic genetic analyses in plants and animals, however, makes it difficult to discern whether PFDs participate in a process as the canonical complex or in alternative configurations, which is necessary to understand their mode of action. To tackle this question, and on the premise that the canonical complex cannot be formed if one subunit is missing, we generated an Arabidopsis (Arabidopsis thaliana) mutant deficient in the six PFDs and compared various growth and environmental responses with those of the individual mutants. In this way, we demonstrate that the PFDc is required for seed germination, to delay flowering, or to respond to high salt stress or low temperature, whereas at least two PFDs redundantly attenuate the response to osmotic stress. A coexpression analysis of differentially expressed genes in the sextuple mutant identified several transcription factors, including ABA INSENSITIVE 5 (ABI5) and PHYTOCHROME-INTERACTING FACTOR 4, acting downstream of PFDs. Furthermore, the transcriptomic analysis allowed assigning additional roles for PFDs, for instance, in response to higher temperature.

An Analysis of Structure-function Co-relation between GLI Oncoprotein and HLA Immune-gene Transcriptional Regulation through Molecular Docking

Background:

GLI proteins play a significant role in the transduction of the Hedgehog (Hh) signaling pathway. A variety of human cancers, including the brain, gastrointestinal, lung, breast, and prostate cancers, demonstrate inappropriate activation of this pathway. GLI helps in proliferation and has an inhibitory role in the differentiation of hematopoietic stem cells. Malignancies may have a defect in differentiation. Different types of malignancies and undifferentiated cells have a low level of HLA expression on their cell surface.

Objective:

Human Leukocytic Antigen (HLA) downregulation is frequently observed in cancer cells. This work is aimed to hypothesize whether this downregulation of HLA molecules is GLI oncoprotein mediated or not. To understand the roles of different types of GLI oncoproteins on different classes of HLA transcriptional machinery was carried out through structure-based modeling and molecular docking studies.

Methods:

To investigate the role of GLI in HLA expression /downregulation is Hh-GLI mediated or not, molecular docking based computational interaction studies were performed between different GLI proteins (GLI1, GLI2, and GLI3) with TATA box binding protein (TBP) and compare the binding efficiencies of different HLA gene (both HLA class I and –II) regulating transcription factors (RelA, RFX5, RFXAP, RFXANK, CIITA, CREB1, and their combinations) with TBP. Due to unavailability of 3D protein structures of GLI2 and cyclin D2 (a natural ligand of GLI1) were modelled followed by structural validation by Ramachandran plot analysis.

Results:

GLI proteins especially, GLI1 and GLI2, have almost similar binding energy of RFX5-RFXANK- RFXAP and CIITA multi-protein complex to TBP but has lower binding energy between RelA to TBP.

Conclusion:

This study suggests that HLA class I may not be downregulated by GLI; however, over-expression of GLI1 is may be responsible for HLA class II downregulation. Thus this protein may be responsible for the maintenance of the undifferentiated state of malignant cells. This study also suggests the implicative role of GLI1 in the early definitive stage of hematopoiesis.

Variants of human CLDN9 cause mild to profound hearing loss

Hereditary deafness is clinically and genetically heterogeneous. We investigated deafness segregating as a recessive trait in two families. Audiological examinations revealed an asymmetric mild to profound hearing loss with childhood or adolescent onset. Exome sequencing of probands identified a homozygous c.475G>A;p.(Glu159Lys) variant of CLDN9 (NM_020982.4) in one family and a homozygous c.370_372dupATC;p.(Ile124dup) CLDN9 variant in an affected individual of a second family. Claudin 9 (CLDN9) is an integral membrane protein and constituent of epithelial bicellular tight junctions (TJs) that form semipermeable, paracellular barriers between inner ear perilymphatic and endolymphatic compartments. Computational structural modeling predicts that substitution of a lysine for glutamic acid p.(Glu159Lys) alters one of two cis-interactions between CLDN9 protomers. The p.(Ile124dup) variant is predicted to locally misfold CLDN9 and mCherry tagged p.(Ile124dup) CLDN9 is not targeted to the HeLa cell membrane. In situ hybridization shows that mouse Cldn9 expression increases from embryonic to postnatal development and persists in adult inner ears coinciding with prominent CLDN9 immunoreactivity in TJs of epithelia outlining the scala media. Together with the Cldn9 deaf mouse and a homozygous frameshift of CLDN9 previously associated with deafness, the two bi-allelic variants of CLDN9 described here point to CLDN9 as a bona fide human deafness gene.

Structural model and functional properties of an exo-polygalacturonase from Neosartorya glabra

A purified exo-polygalacturonase of Neosartorya glabra (EplNg) was successfully characterized. EplNg native presented 68.2 kDa, with 32% carbohydrate content. The deglycosylated form showed 46.3 kDa and isoelectric point of 5.4. The identity of EplNg was confirmed as an exo-polygalacturonase class I (EC 3.2.1.67) using mass spectrometry and Western-Blotting. Capillary electrophoresis indicated that only galacturonic acid was released by the action of EplNg on sodium polypectate, confirming an exoenzyme character. The structural model confers that EplNg has a core formed by twisted parallel β-sheets structure. Among twelve putative cysteines, ten were predicted to form disulfide bridges. The catalytic triad predicted is composed of Asp²²³, Asp²⁴⁵, and Asp²⁴⁶ aligned along with a distance in 4-5 Å, suggesting that EplNg probably does not perform the standard inverting catalytic mechanism described for the GH28 family. EplNg was active from 30 to 90 °C, with maximum activity at 65 °C, pH 5.0. The Km and Vmax determined using sodium polypectate were 6.9 mg·mL^-1 and Vmax 690 μmol·min^-1.mg^-1, respectively. EplNg was active and stable over a wide range of pH values and temperatures, confirming the interesting properties EplNg and provide a basis for the development of the enzyme in different biotechnological processes.

Similarities and differences in the structures and proteoform profiles of the complement proteins C6 and C7

The human complement system provides a first line of defence against pathogens. It requires a well-orchestrated sequential assembly of an array of terminal complement components (C5, C6, C7, C8, and C9), ultimately forming the membrane attack complex (MAC). Although much information about MAC assembly is available, the structure of the soluble C7 has remained elusive. The complement proteins C7 and C6 share very high sequence homology and exhibit several conserved domains, disulphide bridges, and C-mannosylation sites. Here, we used an integrative structural MS-based approach combining native MS, glycopeptide-centric MS, in-gel cross-linking MS (IGX-MS) and structural modelling to describe structural features, including glycosylation, of human serum soluble C7. We compare this data with structural and glycosylation data for human serum C6. The new structural model for C7 shows that it adopts a compact conformation in solution. Although C6 and C7 share many similarities, our data reveals distinct O-, and N-linked glycosylation patterns in terms of location and glycan composition. Cumulatively, our data provide valuable new insight into the structure and proteoforms of C7, solving an essential piece of the puzzle in our understanding of MAC assembly.

Identifying Primate ACE2 Variants That Confer Resistance to SARS-CoV-2

SARS-CoV-2 infects humans through the binding of viral S-protein (spike protein) to human angiotensin I converting enzyme 2 (ACE2). The structure of the ACE2-S-protein complex has been deciphered and we focused on the 27 ACE2 residues that bind to S-protein. From human sequence databases, we identified nine ACE2 variants at ACE2-S-protein binding sites. We used both experimental assays and protein structure analysis to evaluate the effect of each variant on the binding affinity of ACE2 to S-protein. We found one variant causing complete binding disruption, two and three variants, respectively, strongly and mildly reducing the binding affinity, and two variants strongly enhancing the binding affinity. We then collected the ACE2 gene sequences from 57 nonhuman primates. Among the 6 apes and 20 Old World monkeys (OWMs) studied, we found no new variants. In contrast, all 11 New World monkeys (NWMs) studied share four variants each causing a strong reduction in binding affinity, the Philippine tarsier also possesses three such variants, and 18 of the 19 prosimian species studied share one variant causing a strong reduction in binding affinity. Moreover, one OWM and three prosimian variants increased binding affinity by >50%. Based on these findings, we proposed that the common ancestor of primates was strongly resistant to and that of NWMs was completely resistant to SARS-CoV-2 and so is the Philippine tarsier, whereas apes and OWMs, like most humans, are susceptible. This study increases our understanding of the differences in susceptibility to SARS-CoV-2 infection among primates.

Protein model accuracy estimation empowered by deep learning and inter-residue distance prediction in CASP14

The inter-residue contact prediction and deep learning showed the promise to improve the estimation of protein model accuracy (EMA) in the 13th Critical Assessment of Protein Structure Prediction (CASP13). To further leverage the improved inter-residue distance predictions to enhance EMA, during the 2020 CASP14 experiment, we integrated several new inter-residue distance features with the existing model quality assessment features in several deep learning methods to predict the quality of protein structural models. According to the evaluation of performance in selecting the best model from the models of CASP14 targets, our three multi-model predictors of estimating model accuracy (MULTICOM-CONSTRUCT, MULTICOM-AI, and MULTICOM-CLUSTER) achieve the averaged loss of 0.073, 0.079, and 0.081, respectively, in terms of the global distance test score (GDT-TS). The three methods are ranked first, second, and third out of all 68 CASP14 predictors. MULTICOM-DEEP, the single-model predictor of estimating model accuracy (EMA), is ranked within top 10 among all the single-model EMA methods according to GDT-TS score loss. The results demonstrate that inter-residue distance features are valuable inputs for deep learning to predict the quality of protein structural models. However, larger training datasets and better ways of leveraging inter-residue distance information are needed to fully explore its potentials.

Tryptophan Operon Diversity Reveals Evolutionary Trends among Geographically Disparate Chlamydia trachomatis Ocular and Urogenital Strains Affecting Tryptophan Repressor and Synthase Function

The obligate intracellular pathogen Chlamydia trachomatis (Ct) is the leading cause of bacterial sexually transmitted infections and blindness globally. To date, Ct urogenital strains are considered tryptophan prototrophs, utilizing indole for tryptophan synthesis within a closed-conformation tetramer comprised of two α (TrpA)- and two β (TrpB)-subunits. In contrast, ocular strains are auxotrophs due to mutations in TrpA, relying on host tryptophan pools for survival. It has been speculated that there is strong selective pressure for urogenital strains to maintain a functional operon. Here, we performed genetic, phylogenetic, and novel functional modeling analyses of 595 geographically diverse Ct ocular, urethral, vaginal, and rectal strains with complete operon sequences. We found that ocular and urogenital, but not lymphogranuloma venereum, TrpA-coding sequences were under positive selection. However, vaginal and urethral strains exhibited greater nucleotide diversity and a higher ratio of nonsynonymous to synonymous substitutions [Pi(a)/Pi(s)] than ocular strains, suggesting a more rapid evolution of beneficial mutations. We also identified nonsynonymous amino acid changes for an ocular isolate with a urogenital backbone in the intergenic region between TrpR and TrpB at the exact binding site for YtgR-the only known iron-dependent transcription factor in Chlamydia-indicating that selective pressure has disabled the response to fluctuating iron levels. In silico effects on protein stability, ligand-binding affinity, and tryptophan repressor (TrpR) affinity for single-stranded DNA (ssDNA) measured by calculating free energy changes (ΔΔG) between Ct reference and mutant tryptophan operon proteins were also analyzed. We found that tryptophan synthase function was likely suboptimal compared to other bacterial tryptophan prototrophs and that a diversity of urogenital strain mutations rendered the synthase nonfunctional or inefficient. The novel mutations identified here affected active sites in an orthosteric manner but also hindered α- and β-subunit allosteric interactions from distant sites, reducing efficiency of the tryptophan synthase. Importantly, strains with mutant proteins were inclined toward energy conservation by exhibiting an altered affinity for their respective ligands compared to reference strains, indicating greater fitness. This is not surprising as l-tryptophan is one of the most energetically costly amino acids to synthesize. Mutations in the tryptophan repressor gene (trpR) among urogenital strains were similarly detrimental to function. Our findings indicate that urogenital strains are evolving more rapidly than previously recognized with mutations that impact tryptophan operon function in a manner that is energetically beneficial, providing a novel host-pathogen evolutionary mechanism for intracellular survival.IMPORTANCE Chlamydia trachomatis (Ct) is a major global public health concern causing sexually transmitted and ocular infections affecting over 130 million and 260 million people, respectively. Sequelae include infertility, preterm birth, ectopic pregnancy, and blindness. Ct relies on available host tryptophan pools and/or substrates to synthesize tryptophan to survive. Urogenital strains synthesize tryptophan from indole using their intact tryptophan synthase (TS). Ocular strains contain a trpA frameshift mutation that encodes a truncated TrpA with loss of TS function. We found that TS function is likely suboptimal compared to other tryptophan prototrophs and that urogenital stains contain diverse mutations that render TS nonfunctional/inefficient, evolve more rapidly than previously recognized, and impact operon function in a manner that is energetically beneficial, providing an alternative host-pathogen evolutionary mechanism for intracellular survival. Our research has broad scientific appeal since our approach can be applied to other bacteria that may explain evolution/survival in host-pathogen interactions.

A blueprint for gene function analysis through Base Editing in the model plant Physcomitrium (Physcomitrella) patens

CRISPR-Cas9 has proven to be highly valuable for genome editing in plants, including the model plant Physcomitrium patens. However, the fact that most of the editing events produced using the native Cas9 nuclease correspond to small insertions and deletions is a limitation. CRISPR-Cas9 base editors enable targeted mutation of single nucleotides in eukaryotic genomes and therefore overcome this limitation. Here, we report two programmable base-editing systems to induce precise cytosine or adenine conversions in P. patens. Using cytosine or adenine base editors, site-specific single-base mutations can be achieved with an efficiency up to 55%, without off-target mutations. Using the APT gene as a reporter of editing, we could show that both base editors can be used in simplex or multiplex, allowing for the production of protein variants with multiple amino-acid changes. Finally, we set up a co-editing selection system, named selecting modification of APRT to report gene targeting (SMART), allowing up to 90% efficiency site-specific base editing in P. patens. These two base editors will facilitate gene functional analysis in P. patens, allowing for site-specific editing of a given base through single sgRNA base editing or for in planta evolution of a given gene through the production of randomly mutagenised variants using multiple sgRNA base editing.

Representations of protein structure for exploring the conformational space: A speed-accuracy trade-off

•

We compare ten structural representations, either atomistic or coarse-grained.

•

Thus, ten distance-dependent statistical potentials of mean force (PMF) were built.

•

The Cβ-only and Cα + Cβ representations provide the best speed–accuracy trade-off.

•

Including glycines through Cα, in a Cβ-only representation, yields a higher accuracy.

•

We generalize the conclusions to the total information gain (TIG) scoring function.

The recent breakthrough in the field of protein structure prediction shows the relevance of using knowledge-based based scoring functions in combination with a low-resolution 3D representation of protein macromolecules. The choice of not using all atoms is barely supported by any data in the literature, and is mostly motivated by empirical and practical reasons, such as the computational cost of assessing the numerous folds of the protein conformational space. Here, we present a comprehensive study, carried on a large and balanced benchmark of predicted protein structures, to see how different types of structural representations rank in either accuracy or calculation speed, and which ones offer the best compromise between these two criteria. We tested ten representations, including low-resolution, high-resolution, and coarse-grained approaches. We also investigated the generalization of the findings to other formalisms than the widely-used “potential of mean force” (PMF) method. Thus, we observed that representing protein structures by their β carbons—combined or not with Cα—provides the best speed–accuracy trade-off, when using a “total information gain” scoring function. For statistical PMFs, using MARTINI backbone and side-chains beads is the best option. Finally, we also demonstrated the necessity of training the reference state on all atom types, and of including the Cα atoms of glycine residues, in a Cβ-based representation.

Insights into the Molecular Inhibition of the Oncogenic Channel KV10.1 by Globular Toxins

Inhibition of the expression of the human ether-à-go-go (hEAG1 or hK_V10.1) channel is associated with a dramatic reduction in the growth of several cancerous tumors. The modulation of this channel's activity is a promising target for the development of new anticancer drugs. Although some small molecules have shown inhibitory activity against K_V10.1, their lack of specificity has prevented their use in humans. In vitro studies have recently identified a limited number of peptide toxins with proven specificity in their hK_V10.1 channel inhibitory effect. These peptide toxins have become desirable candidates to use as lead compounds to design more potent and specific hK_V10.1 inhibitors. However, the currently available studies lack the atomic resolution needed to characterize the molecular features that favor their binding to hK_V10.1. In this work, we present the first attempt to locate the possible hK_V10.1 binding sites of the animal peptide toxins APETx4, Aa1a, Ap1a, and k-hefutoxin 1, all of which described as hK_V10.1 inhibitors. Our studies incorporated homology modeling to construct a robust three-dimensional (3D) model of hK_V10.1, applied protein docking, and multiscale molecular dynamics techniques to reveal in atomic resolution the toxin-channel interactions. Our approach suggests that some peptide toxins bind in the outer vestibule surrounding the pore of hK_V10.1; it also identified the channel residues Met397 and Asp398 as possible anchors that stabilize the binding of the evaluated toxins. Finally, a description of the possible mechanism for inhibition and gating is presented.

A Molecular Modeling Approach to Identify Potential Antileishmanial Compounds Against the Cell Division Cycle (cdc)-2-Related Kinase 12 (CRK12) Receptor of Leishmania donovani

The huge burden of leishmaniasis caused by the trypanosomatid protozoan parasite Leishmania is well known. This illness was included in the list of neglected tropical diseases targeted for elimination by the World Health Organization. However, the increasing evidence of resistance to existing antimonial drugs has made the eradication of the disease difficult to achieve, thus warranting the search for new drug targets. We report here studies that used computational methods to identify inhibitors of receptors from natural products. The cell division cycle-2-related kinase 12 (CRK12) receptor is a plausible drug target against Leishmania donovani. This study modelled the 3D molecular structure of the L. donovani CRK12 (LdCRK12) and screened for small molecules with potential inhibitory activity from African flora. An integrated library of 7722 African natural product-derived compounds and known inhibitors were screened against the LdCRK12 using AutoDock Vina after performing energy minimization with GROMACS 2018. Four natural products, namely sesamin (NANPDB1649), methyl ellagic acid (NANPDB1406), stylopine (NANPDB2581), and sennecicannabine (NANPDB6446) were found to be potential LdCRK12 inhibitory molecules. The molecular docking studies revealed two compounds NANPDB1406 and NANPDB2581 with binding affinities of -9.5 and -9.2 kcal/mol, respectively, against LdCRK12 which were higher than those of the known inhibitors and drugs, including GSK3186899, amphotericin B, miltefosine, and paromomycin. All the four compounds were predicted to have inhibitory constant (Ki) values ranging from 0.108 to 0.587 μM. NANPDB2581, NANPDB1649 and NANPDB1406 were also predicted as antileishmanial with Pa and Pi values of 0.415 and 0.043, 0.391 and 0.052, and 0.351 and 0.071, respectively. Molecular dynamics simulations coupled with molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) computations reinforced their good binding mechanisms. Most compounds were observed to bind in the ATP binding pocket of the kinase domain. Lys488 was predicted as a key residue critical for ligand binding in the ATP binding pocket of the LdCRK12. The molecules were pharmacologically profiled as druglike with inconsequential toxicity. The identified molecules have scaffolds that could form the backbone for fragment-based drug design of novel leishmanicides but warrant further studies to evaluate their therapeutic potential.

Biological impact of mutually exclusive exon switching

Alternative splicing can expand the diversity of proteomes. Homologous mutually exclusive exons (MXEs) originate from the same ancestral exon and result in polypeptides with similar structural properties but altered sequence. Why would some genes switch homologous exons and what are their biological impact? Here, we analyse the extent of sequence, structural and functional variability in MXEs and report the first large scale, structure-based analysis of the biological impact of MXE events from different genomes. MXE-specific residues tend to map to single domains, are highly enriched in surface exposed residues and cluster at or near protein functional sites. Thus, MXE events are likely to maintain the protein fold, but alter specificity and selectivity of protein function. This comprehensive resource of MXE events and their annotations is available at: http://gene3d.biochem.ucl.ac.uk/mxemod/. These findings highlight how small, but significant changes at critical positions on a protein surface are exploited in evolution to alter function.

Computational Modulation of the V3 Region of Glycoprotein gp125 of HIV-2

HIV-2 infection is frequently neglected in HIV/AIDS campaigns. However, a special emphasis must be given to HIV-2 as an untreated infection that also leads to AIDS and death, and for which the efficacy of most available drugs is limited against HIV-2. HIV envelope glycoproteins mediate binding to the receptor CD4 and co-receptors at the surface of the target cell, enabling fusion with the cell membrane and viral entry. Here, we developed and optimized a computer-assisted drug design approach of an important HIV-2 glycoprotein that allows us to explore and gain further insights at the molecular level into protein structures and interactions crucial for the inhibition of HIV-2 cell entry. The 3D structure of a key HIV-2ROD gp125 region was generated by a homology modeling campaign. To disclose the importance of the main structural features and compare them with experimental results, 3D-models of six mutants were also generated. These mutations revealed the selective impact on the behavior of the protein. Furthermore, molecular dynamics simulations were performed to optimize the models, and the dynamic behavior was tackled to account for structure flexibility and interactions network formation. Structurally, the mutations studied lead to a loss of aromatic features, which is very important for the establishment of π-π interactions and could induce a structural preference by a specific coreceptor. These new insights into the structure-function relationship of HIV-2 gp125 V3 and surrounding regions will help in the design of better models and the design of new small molecules capable to inhibit the attachment and binding of HIV with host cells.

Molecular Modelling Studies on Pyrazole Derivatives for the Design of Potent Rearranged during Transfection Kinase Inhibitors

RET (rearranged during transfection) kinase, one of the receptor tyrosine kinases, plays a crucial role in the development of the human nervous system. It is also involved in various cell signaling networks responsible for the normal cell division, growth, migration, and survival. Previously reported clinical studies revealed that deregulation or aberrant activation of RET signaling can cause several types of human cancer. For example, medullary thyroid carcinoma (MTC) and multiple endocrine neoplasia (MEN2A, MEN2B) occur due to sporadic mutation or germline RET mutation. A number of RET kinase inhibitors have been approved by the FDA for the treatment of cancer, such as cabozantinib, vandetanib, lenvatinib, and sorafenib. However, each of these drugs is a multikinase inhibitor. Hence, RET is an important therapeutic target for cancer drug design. In this work, we have performed various molecular modelling studies, such as molecular docking and dynamics simulation for the most active compound of the pyrazole series as RET kinase inhibitors. Furthermore, molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) free energy calculation and 3-dimensional quantitative structure-activity relationship (3D-QSAR) were performed using g_mmpbsa and SYBYL-X 2.1 package. The results of this study revealed the crucial binding site residues at the active site of RET kinase and contour map analysis showed important structural characteristics for the design of new highly active inhibitors. Therefore, we have designed ten RET kinase inhibitors, which showed higher inhibitory activity than the most active compound of the series. The results of our study provide insights to design more potent and selective RET kinase inhibitors.

In silico design of novel FAK inhibitors using integrated molecular docking, 3D-QSAR and molecular dynamics simulation studies

Focal adhesion kinase (FAK) is a cytoplasmic tyrosine kinase that plays a crucial role in integrin signaling that regulates essential cellular functions including growth, motility, proliferation and survival in different types of cells. Interestingly, it has also shown to be up-regulated in various types of tumors, hence it has emerged as a significant therapeutic target for the development of selective inhibitors. In present work, with the aim of achieving further insight into the structural characteristics required for the FAK inhibitory activity, a combined approach of molecular modeling studies including molecular docking, three-dimensional quantitative structure activity relationship (3D-QSAR) and molecular dynamics (MD) simulation were carried out on a series of 7H-pyrrolo[2,3-d]pyrimidine and thieno[3,2-d]pyrimidine FAK inhibitors. The probable binding modes and interactions of inhibitors into the FAK active site were predicted by molecular docking. The 3D-QSAR models were developed using the comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) methods, with three ligand-based, docking-based and receptor-based alignment techniques. Both CoMFA and CoMSIA models obtained from receptor-based alignment were superior to the ones obtained by other alignment methods. However, the CoMSIA model (q² = 0.679, r² = 0.954 and r²_pred = 0.888) depicted almost better predictive ability than the CoMFA model (q² = 0.617, r² = 0.932 and r²_pred = 0.856). The contour map analysis revealed the relationship between the structural features and inhibitory activity. The docking results and CoMFA and CoMSIA contour maps were in good accordance. Based on the information obtained from the molecular docking and contour map analysis, a series of novel FAK inhibitors were designed that showed better predicted inhibitory activity than the most potent compound 31 in the data set. Finally, the stability of the reference molecule 31 and the designed compounds D15 and D27 were evaluated through a 30 ns of MD simulation and their binding free energies were calculated using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The result of MD simulation and binding free energy decomposition demonstrated the important role of van der Waals interactions alongside H-bond ones that were in consistent with the docking and contour maps analysis results. In sum, the results from this study may provide a significant insight for developing more effective novel FAK inhibitors.Communicated by Ramaswamy H. Sarma.

Selective cross-linking of coinciding protein assemblies by in-gel cross-linking mass spectrometry

Cross-linking mass spectrometry has developed into an important method to study protein structures and interactions. The in-solution cross-linking workflows involve time and sample consuming steps and do not provide sensible solutions for differentiating cross-links obtained from co-occurring protein oligomers, complexes, or conformers. Here we developed a cross-linking workflow combining blue native PAGE with in-gel cross-linking mass spectrometry (IGX-MS). This workflow circumvents steps, such as buffer exchange and cross-linker concentration optimization. Additionally, IGX-MS enables the parallel analysis of co-occurring protein complexes using only small amounts of sample. Another benefit of IGX-MS, demonstrated by experiments on GroEL and purified bovine heart mitochondria, is the substantial reduction of undesired over-length cross-links compared to in-solution cross-linking. We next used IGX-MS to investigate the complement components C5, C6, and their hetero-dimeric C5b6 complex. The obtained cross-links were used to generate a refined structural model of the complement component C6, resembling C6 in its inactivated state. This finding shows that IGX-MS can provide new insights into the initial stages of the terminal complement pathway.

Identification of a SARS-like bat coronavirus that shares structural features with the spike glycoprotein receptor-binding domain of SARS-CoV-2

SARS-CoV-2 is a recently emerged coronavirus that binds angiotensin-converting enzyme 2 (ACE2) for cell entry via its receptor-binding domain (RBD) on a surface-expressed spike glycoprotein. Studies show that despite its similarities to severe acute respiratory syndrome (SARS) coronavirus, there are critical differences in key RBD residues when compared to SARS-CoV-2. Here we present a short in silico study, showing that SARS-like bat coronavirus Rs3367 shares a high conservation with SARS-CoV-2 in important RBD residues for ACE2 binding: SARS-CoV-2's Phe486, Thr500, Asn501 and Tyr505; implicated in receptor-binding strength and host-range determination. These features were not shared with other studied bat coronaviruses belonging to the betacoronavirus genus, including RaTG13, the closest reported bat coronavirus to SARS-CoV-2's spike protein. Sequence and phylogeny analyses were followed by the computation of a reliable model of the RBD of SARS-like bat coronavirus Rs3367, which allowed structural insight of the conserved residues. Superimposition of this model on the SARS-CoV-2 ACE2-RBD complex revealed critical ACE2 contacts are also maintained. In addition, residue Asn488Rs3367 interacted with a previously defined pocket on ACE2 composed of Tyr41, Lys353 and Asp355. When compared to available SARS-CoV-2 crystal structure data, Asn501SARS-CoV-2 showed a different interaction with the ACE2 pocket. Taken together, this study offers molecular insights on RBD-receptor interactions with implications for vaccine design.

An information gain-based approach for evaluating protein structure models

For three decades now, knowledge-based scoring functions that operate through the "potential of mean force" (PMF) approach have continuously proven useful for studying protein structures. Although these statistical potentials are not to be confused with their physics-based counterparts of the same name-i.e. PMFs obtained by molecular dynamics simulations-their particular success in assessing the native-like character of protein structure predictions has lead authors to consider the computed scores as approximations of the free energy. However, this physical justification is a matter of controversy since the beginning. Alternative interpretations based on Bayes' theorem have been proposed, but the misleading formalism that invokes the inverse Boltzmann law remains recurrent in the literature. In this article, we present a conceptually new method for ranking protein structure models by quality, which is (i) independent of any physics-based explanation and (ii) relevant to statistics and to a general definition of information gain. The theoretical development described in this study provides new insights into how statistical PMFs work, in comparison with our approach. To prove the concept, we have built interatomic distance-dependent scoring functions, based on the former and new equations, and compared their performance on an independent benchmark of 60,000 protein structures. The results demonstrate that our new formalism outperforms statistical PMFs in evaluating the quality of protein structural decoys. Therefore, this original type of score offers a possibility to improve the success of statistical PMFs in the various fields of structural biology where they are applied. The open-source code is available for download at https://gitlab.rpbs.univ-paris-diderot.fr/src/ig-score.

Ligand Pose Predictions for Human G Protein-Coupled Receptors: Insights from the Amber-Based Hybrid Molecular Mechanics/Coarse-Grained Approach

Human G protein-coupled receptors (hGPCRs) are the most frequent targets of Food and Drug Administration (FDA)-approved drugs. Structural bioinformatics, along with molecular simulation, can support structure-based drug design targeting hGPCRs. In this context, several years ago, we developed a hybrid molecular mechanics (MM)/coarse-grained (CG) approach to predict ligand poses in low-resolution hGPCR models. The approach was based on the GROMOS96 43A1 and PRODRG united-atom force fields for the MM part. Here, we present a new MM/CG implementation using, instead, the Amber 14SB and GAFF all-atom potentials for proteins and ligands, respectively. The new implementation outperforms the previous one, as shown by a variety of applications on models of hGPCR/ligand complexes at different resolutions, and it is also more user-friendly. Thus, it emerges as a useful tool to predict poses in low-resolution models and provides insights into ligand binding similarly to all-atom molecular dynamics, albeit at a lower computational cost.

MASS: predict the global qualities of individual protein models using random forests and novel statistical potentials

BACKGROUND

Protein model quality assessment (QA) is an essential procedure in protein structure prediction. QA methods can predict the qualities of protein models and identify good models from decoys. Clustering-based methods need a certain number of models as input. However, if a pool of models are not available, methods that only need a single model as input are indispensable.

RESULTS

We developed MASS, a QA method to predict the global qualities of individual protein models using random forests and various novel energy functions. We designed six novel energy functions or statistical potentials that can capture the structural characteristics of a protein model, which can also be used in other protein-related bioinformatics research. MASS potentials demonstrated higher importance than the energy functions of RWplus, GOAP, DFIRE and Rosetta when the scores they generated are used as machine learning features. MASS outperforms almost all of the four CASP11 top-performing single-model methods for global quality assessment in terms of all of the four evaluation criteria officially used by CASP, which measure the abilities to assign relative and absolute scores, identify the best model from decoys, and distinguish between good and bad models. MASS has also achieved comparable performances with the leading QA methods in CASP12 and CASP13.

CONCLUSIONS

MASS and the source code for all MASS potentials are publicly available at http://dna.cs.miami.edu/MASS/ .

Putative circumsporozoite protein (CSP) of Plasmodium vivax is considerably distinct from the well-known CSP and plays a role in the protein ubiquitination pathway

Amidst technical challenges which limit successful culture and genetic manipulation of P. vivax parasites, we used a computational approach to identify a critical target with evolutionary significance. The putative circumsporozoite protein on chromosome 13 of P. vivax (PvpuCSP)is distinct from the well-known vaccine candidate PfCSP. The aim of this study was to understand the role of PvpuCSP and its relatedness to the well-known CSP. The study revealed PvpuCSP as a membrane bound E3 ubiquitin ligase involved in ubiquitination. It has a species-specific tetra-peptide unit which is differentially repeated in various P. vivax strains. The PvpuCSP is different from CSP in terms of stage-specific expression and function. Since E3 ubiquitin ligases are known antimalarial drug targets targeting the proteasome pathway, PvpuCSP, with evolutionary connotation and a key role in orchestrating protein degradation in P. vivax, can be explored as a potential drug target.

•

PvpuCSP is predicted as E3 ubiquitin ligase, a part of ubiquitination pathway.

•

Tetra-peptide tandem repeat at C terminal of PvpuCSP is exclusive to P. vivax.

•

Moderately expressed during all parasitic stages in host and vector

•

Partially disordered protein with both structured domains and two distinct IDRs

•

A transmembrane protein with highly conserved functional domain across Apicomplexa

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.

The Mechanism of Cholesterol Modification of Hedgehog Ligand

Hedgehog (Hh) proteins are important components of signal transduction pathways involved in animal development, and their defects are implicated in carcinogenesis. Their N-terminal domain (HhN) acts as a signaling ligand, and their C-terminal domain (HhC) performs an autocatalytic function of cleaving itself away, while adding a cholesterol moiety to HhN. HhC has two sub-domains: a hedgehog/intein (hint) domain that primarily performs the autocatalytic activity, and a sterol-recognition region (SRR) that binds to cholesterol and properly positions it with respect to HhN. The three-dimensional details of this autocatalytic mechanism remain unknown, as does the structure of the precursor Hh protein. In this study, a complete cholesterol-bound precursor form of the drosophila Hh precursor is modeled using known crystal structures of HhN and the hint domain, and a hypothesized similarity of SRR to an unrelated but similar-sized cholesterol binding protein. The restrained geometries and topology switching (RGATS) strategy is then used to predict atomic-detail pathways for the full autocatalytic reaction starting from the precursor and ending in a cholesterol-linked HhN domain and a cleaved HhC domain. The RGATS explicit solvent simulations indicate the roles of individual HhC residues in facilitating the reaction, which can be confirmed through mutational experiments. These simulations also provide plausible structural models for the N/S acyl transfer intermediate and the product states of this reaction. This study thus provides a good framework for future computational and experimental studies to develop a full structural and dynamic understanding of Hh autoprocessing. © 2019 Wiley Periodicals, Inc.

3D-QSAR Assisted Design of Novel 7-Deazapurine Derivatives as TNNI3K Kinase Inhibitors Using Molecular Docking and Molecular Dynamics Simulation

Background:

Cardiac troponin I-interacting kinase (TNNI3K) is a cardiac-specific kinase that belongs to MAPKKK family. It is a dual-function kinase with tyrosine and serine/threonine kinase activity. Over-expression of TNNI3K results in various cardiovascular diseases such as cardiomyopathy, ischemia/reperfusion injury, heart failure, etc. Since, it is a cardiac-specific kinase and expressed only in heart tissue, it is an ideal molecular target to treat cardiac diseases. The main objective of the work is to study and understand the structure-activity relationship of the reported deazapurine derivatives and to use the 3D-QSAR and docking results to design potent and novel TNNI3K inhibitors of this series.

Methods:

In the present study, we have used molecular docking 3D QSAR, and molecular dynamics simulation to understand the structure-activity correlation of reported TNNI3K inhibitors and to design novel compounds of deazapurine derivatives with increased activity.

Results:

Both CoMFA (q2=0.669, NOC=5, r2=0.944) and CoMSIA (q2=0.783, NOC=5, r2=0.965) have resulted in satisfactory models. The models were validated using external test set, Leave-out- Five, bootstrapping, progressive scrambling, and rm2 metrics calculations. The validation procedures showed the developed models were robust and reliable. The docking results and the contour maps analysis helped in the better understanding of the structure-activity relationship.

Conclusion:

This is the first report on 3D-QSAR modeling studies of TNNI3K inhibitors. Both docking and MD results were consistent and showed good correlation with the previous experimental data. Based on the information obtained from contour maps, 31 novel TNNI3K inhibitors were designed. These designed compounds showed higher activity than the existing dataset compounds.

Characterisation of a catalytic triad and reaction selectivity in the dual mechanism of the catalyse hydride transfer in xylitol phosphate dehydrogenase

Xylitol is a high-value low-calorie sweetener used as sugar substitute in food and pharmaceutical industry. Xylitol phosphate dehydrogenase (XPDH) catalyses the conversion of d-xylulose 5-phosphate (XU5P) and d-ribulose 5-phosphate (RU5P) to xylitol and ribitol respectively in the presence of nicotinamide adenine dinucleotide hydride (NADH). Although these enzymes have been shown to produce xylitol and ribitol, there is an incomplete understanding of the mechanism of the catalytic events of these reactions and the detailed mechanism has yet to be elucidated. The main goal of this work is to analyse the conformational changes of XPDH-bound ligands such as zinc^, NADH, XU5P, and RU5P to elucidate the key amino acids involved in the substrate binding. In silico modelling, comparative molecular dynamics simulations, interaction analysis and conformational study were carried out on three XPDH enzymes of the Medium-chain dehydrogenase (MDR) family in order to elucidate the atomistic details of conformational transition, especially on the open and closed state of XPDH. The analysis also revealed the possible mechanism of substrate specificity that are responsible in the catalyse hydride transfer are the residues His58 and Ser39 which would act as the proton donor for reduction of XU5P and RU5P respectively. The structural comparison and MD simulations displayed a significant difference in the conformational dynamics of the catalytic and coenzyme loops between Apo and XPDH-complexes and highlight the contribution of newly found triad residues. This study would assist future mutagenesis study and enzyme modification work to increase the catalysis efficiency of xylitol production in the industry.