How well can the accuracy of comparative protein structure models be predicted?

Design, synthesis, and biological characterization of proteolysis targeting chimera (PROTACs) for the ataxia telangiectasia and RAD3-related (ATR) kinase

The Ataxia telangiectasia and RAD3-related (ATR) kinase is a key regulator of DNA replication stress responses and DNA-damage checkpoints. Several potent and selective ATR inhibitors are reported and four of them are currently in clinical trials in combination with radio- or chemotherapy. Based on the idea of degrading target proteins rather than inhibiting them, we designed, synthesized and biologically characterized a library of ATR-targeted proteolysis targeting chimera (PROTACs). Among the synthesized compounds, the lenalidomide-based PROTAC 42i was the most promising. In pancreatic and cervix cancer cells cancer cells, it reduced ATR to 40 % of the levels in untreated cells. 42i selectively degraded ATR through the proteasome, dependent on the E3 ubiquitin ligase component cereblon, and without affecting the associated kinases ATM and DNA-PKcs. 42i may be a promising candidate for further optimization and biological characterization in various cancer cells.

Structural coverage of the human interactome

Complex biological processes in cells are embedded in the interactome, representing the complete set of protein-protein interactions. Mapping and analyzing the protein structures are essential to fully comprehending these processes' molecular details. Therefore, knowing the structural coverage of the interactome is important to show the current limitations. Structural modeling of protein-protein interactions requires accurate protein structures. In this study, we mapped all experimental structures to the reference human proteome. Later, we found the enrichment in structural coverage when complementary methods such as homology modeling and deep learning (AlphaFold) were included. We then collected the interactions from the literature and databases to form the reference human interactome, resulting in 117 897 non-redundant interactions. When we analyzed the structural coverage of the interactome, we found that the number of experimentally determined protein complex structures is scarce, corresponding to 3.95% of all binary interactions. We also analyzed known and modeled structures to potentially construct the structural interactome with a docking method. Our analysis showed that 12.97% of the interactions from HuRI and 73.62% and 32.94% from the filtered versions of STRING and HIPPIE could potentially be modeled with high structural coverage or accuracy, respectively. Overall, this paper provides an overview of the current state of structural coverage of the human proteome and interactome.

Antibiofilm potential of gallic acid against Klebsiella pneumoniae and Enterobacter hormaechei: in-vitro and in-silico analysis

Biofilm refers to a community of microorganisms that adhere to a substrate and play a crucial role in microbial pathogenesis and developing infections associated with medical devices. Enterobacter hormaechei and Klebsiella pneumoniae are classified as significant nosocomial pathogens within the ESKAPE category and cause diverse infections. In addition to their reputation as prolific biofilm formers, these pathogens are increasingly becoming drug-resistant and pose a substantial threat to the healthcare setting. Due to the inherent resistance of biofilms to conventional therapies, novel strategies are imperative for effectively controlling E. hormaechei and K. pneumoniae biofilms. This study aimed to assess the anti-biofilm activity of gallic acid (GA) against E. hormaechei and K. pneumoniae. The results of biofilm quantification assays demonstrated that GA exhibited significant antibiofilm activity against E. hormaechei and K. pneumoniae at concentrations of 4 mg mL-1, 2 mg mL-1, 1 mg mL-1, and 0.5 mg mL-1. Similarly, GA exhibited a dose-dependent reduction in violacein production, a QS-regulated purple pigment, indicating its ability to suppress violacein production and disrupt QS mechanisms in Chromobacterium violaceum. Additionally, computational tools were utilized to identify the potential target involved in the biofilm formation pathway. The computational analysis further indicated the strong binding affinity of GA to essential biofilm regulators, MrkH and LuxS, suggesting its potential in targeting the c-di-GMP and quorum sensing (QS) pathways to hinder biofilm formation in K. pneumoniae. These compelling findings strongly advocate GA as a promising drug candidate against biofilm-associated infections caused by E. hormaechei and K. pneumoniae.

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.

The novel anti-cancer feature of Brazzein through activating of hTLR5 by integration of biological evaluation: molecular docking and molecular dynamics simulation

Many of plant proteins exhibit the properties similar to the antitumor proteins although the anticancer activity of Brazzein on modulating the autophagy signaling pathway has not been determined so far. The present study aimed to develop a simplified system to enable the rational design of the activating extracellular domain of human Toll-like receptor 5 (hTLR5). To identify the anticancer effect of Brazzein, HADDOCK program and molecular dynamics (MD) simulation were applied to examine the binding of the wild type (WT) and p.A19K mutant of Brazzein to the TLR5. The expression of MAP1S and TNF-α genes was estimated based on real-time PCR. The results clearly confirmed that the WT of Brazzein activated hTLR5 in the MCF-7 cell line since the genes were more and significantly less expressed in the cells treated with the WT and p.A19K mutant than the control, respectively. The snapshots of MD simulation exhibit the consistent close interactions of hTLR5 with the two helices of Brazzein on its lateral side. The results of per residue-free energy decomposition analysis substantiate those of intermolecular contact analysis perfectly one. We propose that the WT of Brazzein can act as an antitumor drug candidate.

The evolution of the HIV-1 protease folding stability

The evolution of structural proteins is generally constrained by the folding stability. However, little is known about the particular capacity of viral proteins to accommodate mutations that can potentially affect the protein stability and, in general, the evolution of the protein stability over time. As an illustrative model case, here, we investigated the evolution of the stability of the human immunodeficiency virus (HIV-1) protease (PR), which is a common HIV-1 drug target, under diverse evolutionary scenarios that include (1) intra-host virus evolution in a cohort of seventy-five patients sampled over time, (2) intra-host virus evolution sampled before and after specific PR-based treatments, and (3) inter-host evolution considering extant and ancestral (reconstructed) PR sequences from diverse HIV-1 subtypes. We also investigated the specific influence of currently known HIV-1 PR resistance mutations on the PR folding stability. We found that the HIV-1 PR stability fluctuated over time within a constant and wide range in any studied evolutionary scenario, accommodating multiple mutations that partially affected the stability while maintaining activity. We did not identify relationships between change of PR stability and diverse clinical parameters such as viral load, CD4⁺ T-cell counts, and a surrogate of time from infection. Counterintuitively, we predicted that nearly half of the studied HIV-1 PR resistance mutations do not significantly decrease stability, which, together with compensatory mutations, would allow the protein to adapt without requiring dramatic stability changes. We conclude that the HIV-1 PR presents a wide structural plasticity to acquire molecular adaptations without affecting the overall evolution of stability.

Definition of the Acceptor Substrate Binding Specificity in Plant Xyloglucan Endotransglycosylases Using Computational Chemistry

Xyloglucan endotransglycosylases (XETs) play key roles in the remodelling and reconstruction of plant cell walls. These enzymes catalyse homo-transglycosylation reactions with xyloglucan-derived donor and acceptor substrates and hetero-transglycosylation reactions with a variety of structurally diverse polysaccharides. In this work, we describe the basis of acceptor substrate binding specificity in non-specific Tropaeolum majus (TmXET6.3) and specific Populus tremula x tremuloides (PttXET16A) XETs, using molecular docking and molecular dynamics (MD) simulations combined with binding free energy calculations. The data indicate that the enzyme-donor (xyloglucan heptaoligosaccharide or XG-OS7)/acceptor complexes with the linear acceptors, where a backbone consisted of glucose (Glc) moieties linked via (1,4)- or (1,3)-β-glycosidic linkages, were bound stably in the active sites of TmXET6.3 and PttXET16A. Conversely, the acceptors with the (1,6)-β-linked Glc moieties were bound stably in TmXET6.3 but not in PttXET16A. When in the (1,4)-β-linked Glc containing acceptors, the saccharide moieties were replaced with mannose or xylose, they bound stably in TmXET6.3 but lacked stability in PttXET16A. MD simulations of the XET-donor/acceptor complexes with acceptors derived from (1,4;1,3)-β-glucans highlighted the importance of (1,3)-β-glycosidic linkages and side chain positions in the acceptor substrates. Our findings explain the differences in acceptor binding specificity between non-specific and specific XETs and associate theoretical to experimental data.

All-Atom Simulations Uncover Structural and Dynamical Properties of STING Proteins in the Membrane System

Recent studies have shown that the stimulator of interferon gene (STING) protein plays a central role in the immune system by facilitating the production of type I interferons in cells. The STING signaling pathway is also a prominent activator of cancer-killing T cells that initiate a powerful adaptive immune response. Since biomolecular signaling pathways are complicated and not easily identified through traditional experiments, molecular dynamics (MD) has often been used to study structural and dynamical responses of biological pathways. Here, we carried out MD simulations for full-length chicken and human STING (chSTING and hSTING) proteins. Specifically, we investigated ligand-bound closed (holo) and ligand-unbound open (apo) forms of STING in the membrane system by comparing their conformational and dynamical differences. Our research provides clues for understanding the mechanism of the STING signaling pathway by uncovering detailed insights for the examined systems: the residues from each chain in the binding pocket are strongly correlated to one another in the open STING structure compared with those in the closed STING structure. Ligand-bound closed STING displays ∼174° rotation of the ligand-binding domain (LBD) relative to the open STING structure. The dynamical analysis of residue Cys148 located in the linker region of hSTING does not support the earlier hypothesis that Cys148 can form disulfide bonds between adjacent STING dimers. We also demonstrate that using the full-length proteins is critical, since the MD simulations of the LBD portion alone cannot properly describe the global conformational properties of STING.

The tomato cytochrome P450 CYP712G1 catalyses the double oxidation of orobanchol en route to the rhizosphere signalling strigolactone, solanacol

Strigolactones (SLs) are rhizosphere signalling molecules and phytohormones. The biosynthetic pathway of SLs in tomato has been partially elucidated, but the structural diversity in tomato SLs predicts that additional biosynthetic steps are required. Here, root RNA-seq data and co-expression analysis were used for SL biosynthetic gene discovery. This strategy resulted in a candidate gene list containing several cytochrome P450s. Heterologous expression in Nicotiana benthamiana and yeast showed that one of these, CYP712G1, can catalyse the double oxidation of orobanchol, resulting in the formation of three didehydro-orobanchol (DDH) isomers. Virus-induced gene silencing and heterologous expression in yeast showed that one of these DDH isomers is converted to solanacol, one of the most abundant SLs in tomato root exudate. Protein modelling and substrate docking analysis suggest that hydroxy-orbanchol is the likely intermediate in the conversion from orobanchol to the DDH isomers. Phylogenetic analysis demonstrated the occurrence of CYP712G1 homologues in the Eudicots only, which fits with the reports on DDH isomers in that clade. Protein modelling and orobanchol docking of the putative tobacco CYP712G1 homologue suggest that it can convert orobanchol to similar DDH isomers as tomato.

Enzyme mediated transformation of CO2 into calcium carbonate using purified microbial carbonic anhydrase

In this study, a bacterial carbonic anhydrase (CA) was purified from Corynebacterium flavescens for the CO₂ conversion into CaCO₃. The synthesized CaCO₃ can be utilized in the papermaking industry as filler material, construction material and in steel industry. Herein, the CA was purified by using a Sephadex G-100 column chromatography having 29.00 kDa molecular mass in SDS-PAGE analysis. The purified CA showed an optimal temperature of 35 °C and pH 7.5. In addition, a kinetic study of CA using p-NPA as substrate showed V_max (166.66 μmoL/mL/min), K_m (5.12 mM), and K_cat (80.56 sec^-1) using Lineweaver Burk plot. The major inhibitors of CA activity were Na²⁺, K⁺, Mn²⁺, and Al³⁺, whereas Zn²⁺ and Fe²⁺ slightly enhanced it. The purified CA showed a good efficacy to convert the CO₂ into CaCO₃ with a total conversion rate of 65.05 mg CaCO₃/mg of protein. In silico analysis suggested that the purified CA has conserved Zn²⁺ coordinating residues such as His 111, His 113, and His 130 in the active site center. Further analysis of the CO₂ binding site showed conserved residues such as Val 132, Val 142, Leu 196, Thr 197, and Val 205. However, a substitution has been observed where Trp 208 of its closest structural homolog T. ammonificans CA is replaced with Arg 207 of C. flavescens. The presence of a hydrophilic mutation in the CO₂ binding hydrophobic region is a further subject of investigation.

Improvement of FGF7 Thermal Stability by Introduction of Mutations in Close Vicinity to Disulfide Bond and Surface Salt Bridge

Fibroblast Growth Factor 7 (FGF7), a growth factor specific to epithelial cells, has attracted attention as a therapeutic protein. However, FGF7 has a limitation in its use due to low protein stability. Here, the mutations were designed to increase the stability of FGF7 by analyzing its 3D structure and sequence of other FGFs. Palifermin, N-terminal truncated FGF7 is known to have improved stability and was used as control protein in our study. The K126 and K178 were substituted into glutamate to form salt bridge with the neighboring residue R175 respectively and A120C mutation was introduced in close vicinity to disulfide bond between C133 and C137. The data of Circular Dichroism (CD) showed that all mutant proteins tested had higher Tm value than Palifermin and Tm of A120C/K126E/K178E FGF7 mutant protein was 15.24 °C higher than that of Palifermin. The results of cell proliferation activity and soluble protein analyzed by sodium dodecyl sulfate–polyacrylamide gel electrophoresis (SDS-PAGE) after 37 °C or 45 °C incubation exhibited that the stability of A120C mutant protein and A120C-including mutant proteins was improved. These results suggest that the mutation of amino acid in close vicinity to disulfide bond and the salt bridge at the surface of FGF7 enhanced thermal stability and make FGF7 more useful for pharmaceutical and cosmetical application.

A Benchmark Dataset for Evaluating Practical Performance of Model Quality Assessment of Homology Models

Protein structure prediction is an important issue in structural bioinformatics. In this process, model quality assessment (MQA), which estimates the accuracy of the predicted structure, is also practically important. Currently, the most commonly used dataset to evaluate the performance of MQA is the critical assessment of the protein structure prediction (CASP) dataset. However, the CASP dataset does not contain enough targets with high-quality models, and thus cannot sufficiently evaluate the MQA performance in practical use. Additionally, most application studies employ homology modeling because of its reliability. However, the CASP dataset includes models generated by de novo methods, which may lead to the mis-estimation of MQA performance. In this study, we created new benchmark datasets, named a homology models dataset for model quality assessment (HMDM), that contain targets with high-quality models derived using homology modeling. We then benchmarked the performance of the MQA methods using the new datasets and compared their performance to that of the classical selection based on the sequence identity of the template proteins. The results showed that model selection by the latest MQA methods using deep learning is better than selection by template sequence identity and classical statistical potentials. Using HMDM, it is possible to verify the MQA performance for high-accuracy homology models.

Broad Specific Xyloglucan:Xyloglucosyl Transferases Are Formidable Players in the Re-Modelling of Plant Cell Wall Structures

Plant xyloglucan:xyloglucosyl transferases, known as xyloglucan endo-transglycosylases (XETs) are the key players that underlie plant cell wall dynamics and mechanics. These fundamental roles are central for the assembly and modifications of cell walls during embryogenesis, vegetative and reproductive growth, and adaptations to living environments under biotic and abiotic (environmental) stresses. XET enzymes (EC 2.4.1.207) have the β-sandwich architecture and the β-jelly-roll topology, and are classified in the glycoside hydrolase family 16 based on their evolutionary history. XET enzymes catalyse transglycosylation reactions with xyloglucan (XG)-derived and other than XG-derived donors and acceptors, and this poly-specificity originates from the structural plasticity and evolutionary diversification that has evolved through expansion and duplication. In phyletic groups, XETs form the gene families that are differentially expressed in organs and tissues in time- and space-dependent manners, and in response to environmental conditions. Here, we examine higher plant XET enzymes and dissect how their exclusively carbohydrate-linked transglycosylation catalytic function inter-connects complex plant cell wall components. Further, we discuss progress in technologies that advance the knowledge of plant cell walls and how this knowledge defines the roles of XETs. We construe that the broad specificity of the plant XETs underscores their roles in continuous cell wall restructuring and re-modelling.

Atypical Hemolytic Uremic Syndrome-Associated FHR1 Isoform FHR1*B Enhances Complement Activation and Inflammation

Atypical hemolytic uremic syndrome (aHUS) is a rare but severe type of thrombotic microangiopathy that is triggered by the abnormal activation of the alternative complement pathway. Previous studies have reported that three completely linked coding variants of CFHR1 form two haplotypes, namely, CFHR1*A (c.469C, c.475C, c.523G) and CFHR1*B (c.469T, c.475G, c.523C). CFHR1*B is associated with susceptibility to aHUS. To explore the genetic mechanism by which CFHR1 isoforms contribute to aHUS, we compared the structures of FHR1*A and FHR1*B by homology modeling and found differences in the angles between SCR3 and SCR4-SCR5, as FHR1*B had a larger angle than FHR1*A. Then, we expressed FHR1*A and FHR1*B recombinant proteins and compared their functions in complement system regulation and inflammation. We found that FHR1*B presented a significantly higher capacity for binding C3b and necrotic cells than FHR1*A. In a cofactor assay, the FHR-1*B showed stronger influence on FH mediated cofactor function than the FHR-1*A, resulted in fewer C3b cleavage products. In the C3 convertase assays, FHR1*B showed more powerful effect compared with FHR1*A regarding to de-regulate FH function of inhibition the assembling of C3bBb. Additionally, we also found that FHR1*B triggered monocytes to secrete higher levels of IL-1β and IL-6 than FHR1*A. In the present study, we showed that variants of CFHR1 might differently affect complement activation and sterile inflammation. Our findings provide a possible mechanism underlying the predisposition to aHUS caused by CFHR1 isoform CFHR1*B.

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.

Identification of the key functional genes in salt-stress tolerance of Cyanobacterium Phormidium tenue using in silico analysis

The development of artificial biocrust using cyanobacterium Phormidium tenue has been suggested as an effective strategy to prevent soil degradation. Here, a combination of in silico approaches with growth rate, photosynthetic pigment, morphology, and transcript analysis was used to identify specific genes and their protein products in response to 500 mM NaCl in P. tenue. The results show that 500 mM NaCl induces the expression of genes encoding glycerol-3-phosphate dehydrogenase (glpD) as a Flavoprotein, ribosomal protein S12 methylthiotransferase (rimO), and a hypothetical protein (sll0939). The constructed co-expression network revealed a group of abiotic stress-responsive genes. Using the Basic Local Alignment Search Tool (BLAST), the homologous proteins of rimO, glpD, and sll0939 were identified in the P. tenue genome. Encoded proteins of glpD, rimO, and DUF1622 genes, respectively, contain (DAO and DAO C), (UPF0004, Radical SAM and TRAM 2), and (DUF1622) domains. The predicted ligand included 22B and MG for DUF1622, FS5 for rimO, and FAD for glpD protein. There was no direct disruption in ligand-binding sites of these proteins by Na+, Cl-, or NaCl. The growth rate, photosynthetic pigment, and morphology of P. tenue were investigated, and the result showed an acceptable tolerance rate of this microorganism under salt stress. The quantitative real-time polymerase chain reaction (qRT-PCR) results revealed the up-regulation of glpD, rimO, and DUF1622 genes under salt stress. This is the first report on computational and experimental analyses of the glpD, rimO, and DUF1622 genes in P. tenue under salt stress to the best of our knowledge.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-03050-w.

Homology Modeling and Molecular Dynamics Simulations of Trypanosoma cruzi Phosphodiesterase b1

Cyclic nucleotide phosphodiesterases have been implicated in the proliferation, differentiation and osmotic regulation of trypanosomatids; in some trypanosomatid species, they have been validated as molecular targets for the development of new therapeutic agents. Because the experimental structure of Trypanosoma cruzi PDEb1 (TcrPDEb1) has not been solved so far, an homology model of the target was created using the structure of Trypanosoma brucei PDEb1 (TbrPDEb1) as a template. The model was refined by extensive enhanced sampling molecular dynamics simulations, and representative snapshots were extracted from the trajectory by combined clustering analysis. This structural ensemble was used to develop a structure-based docking model of the target. The docking accuracy of the model was validated by redocking and cross-docking experiments using all available crystal structures of TbrPDEb1, whereas the scoring accuracy was validated through a retrospective screen, using a carefully curated dataset of compounds assayed against TbrPDEb1 and/or TcrPDEb1. Considering the results from in silico validations, the model may be applied in prospective virtual screening campaigns to identify novel hits, as well as to guide the rational design of potent and selective inhibitors targeting this enzyme.

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.

Evidence that specific interactions play a role in the cholesterol sensitivity of G protein-coupled receptors

G protein-coupled receptors (GPCRs) are known to be modulated by membrane cholesterol levels, but whether or not the effects are caused by specific receptor-cholesterol interactions or cholesterol's general effects on the membrane is not well-understood. We performed coarse-grained molecular dynamics (CGMD) simulations coupled with structural bioinformatics approaches on the β2-adrenergic receptor (β2AR) and the cholecystokinin (CCK) receptor subfamily. The β2AR has been shown to be sensitive to membrane cholesterol and cholesterol molecules have been clearly resolved in numerous β2AR crystal structures. The two CCK receptors are highly homologous and preserve similar cholesterol recognition motifs but despite their homology, CCK1R shows functional sensitivity to membrane cholesterol while CCK2R does not. Our results offer new insights into how cholesterol modulates GPCR function by showing cholesterol interactions with β2AR that agree with previously published data; additionally, we observe differential and specific cholesterol binding in the CCK receptor subfamily while revealing a previously unreported Cholesterol Recognition Amino-acid Consensus (CRAC) sequence that is also conserved across 38% of class A GPCRs. A thermal denaturation assay (LCP-Tm) shows that mutation of a conserved CRAC sequence on TM7 of the β2AR affects cholesterol stabilization of the receptor in a lipid bilayer. The results of this study provide a better understanding of receptor-cholesterol interactions that can contribute to novel and improved therapeutics for a variety of diseases.

Force Field Parameters for Fe2+4S2-4 Clusters of Dihydropyrimidine Dehydrogenase, the 5-Fluorouracil Cancer Drug Deactivation Protein: A Step towards In Silico Pharmacogenomics Studies

The dimeric dihydropyrimidine dehydrogenase (DPD), metalloenzyme, an adjunct anti-cancer drug target, contains highly specialized 4 × Fe²⁺₄S²⁻₄ clusters per chain. These clusters facilitate the catalysis of the rate-limiting step in the pyrimidine degradation pathway through a harmonized electron transfer cascade that triggers a redox catabolic reaction. In the process, the bulk of the administered 5-fluorouracil (5-FU) cancer drug is inactivated, while a small proportion is activated to nucleic acid antimetabolites. The occurrence of missense mutations in DPD protein within the general population, including those of African descent, has adverse toxicity effects due to altered 5-FU metabolism. Thus, deciphering mutation effects on protein structure and function is vital, especially for precision medicine purposes. We previously proposed combining molecular dynamics (MD) and dynamic residue network (DRN) analysis to decipher the molecular mechanisms of missense mutations in other proteins. However, the presence of Fe²⁺₄S²⁻₄ clusters in DPD poses a challenge for such in silico studies. The existing AMBER force field parameters cannot accurately describe the Fe²⁺ center coordination exhibited by this enzyme. Therefore, this study aimed to derive AMBER force field parameters for DPD enzyme Fe²⁺ centers, using the original Seminario method and the collation features Visual Force Field Derivation Toolkit as a supportive approach. All-atom MD simulations were performed to validate the results. Both approaches generated similar force field parameters, which accurately described the human DPD protein Fe²⁺₄S²⁻₄ cluster architecture. This information is crucial and opens new avenues for in silico cancer pharmacogenomics and drug discovery related research on 5-FU drug efficacy and toxicity issues.

Can molecular dynamics simulations improve the structural accuracy and virtual screening performance of GPCR models?

The determination of G protein-coupled receptor (GPCR) structures at atomic resolution has improved understanding of cellular signaling and will accelerate the development of new drug candidates. However, experimental structures still remain unavailable for a majority of the GPCR family. GPCR structures and their interactions with ligands can also be modelled computationally, but such predictions have limited accuracy. In this work, we explored if molecular dynamics (MD) simulations could be used to refine the accuracy of in silico models of receptor-ligand complexes that were submitted to a community-wide assessment of GPCR structure prediction (GPCR Dock). Two simulation protocols were used to refine 30 models of the D3 dopamine receptor (D3R) in complex with an antagonist. Close to 60 μs of simulation time was generated and the resulting MD refined models were compared to a D3R crystal structure. In the MD simulations, the receptor models generally drifted further away from the crystal structure conformation. However, MD refinement was able to improve the accuracy of the ligand binding mode. The best refinement protocol improved agreement with the experimentally observed ligand binding mode for a majority of the models. Receptor structures with improved virtual screening performance, which was assessed by molecular docking of ligands and decoys, could also be identified among the MD refined models. Application of weak restraints to the transmembrane helixes in the MD simulations further improved predictions of the ligand binding mode and second extracellular loop. These results provide guidelines for application of MD refinement in prediction of GPCR-ligand complexes and directions for further method development.

Protein Structure Modeling with MODELLER

Genome sequencing projects have resulted in a rapid increase in the number of known protein sequences. In contrast, only about one-hundredth of these sequences have been characterized at atomic resolution using experimental structure determination methods. Computational protein structure modeling techniques have the potential to bridge this sequence-structure gap. In the following chapter, we present an example that illustrates the use of MODELLER to construct a comparative model for a protein with unknown structure. Automation of a similar protocol has resulted in models of useful accuracy for domains in more than half of all known protein sequences.

Demystifying the Molecular Basis of Pyrazoloquinolinones Recognition at the Extracellular α1+/β3- Interface of the GABAA Receptor by Molecular Modeling

GABA_A receptors are pentameric ligand-gated ion channels that serve as major inhibitory neurotransmitter receptors in the mammalian brain and the target of numerous clinically relevant drugs interacting with different ligand binding sites. Here, we report an in silico approach to investigate the binding of pyrazoloquinolinones (PQs) that mediate allosteric effects through the extracellular α+/β- interface of GABA_A receptors. First, we docked a potent prototype of PQs into the α1+/β3- site of a homology model of the human α1β3γ2 subtype of the GABA_A receptor. Next, for each docking pose, we computationally derived protein-ligand complexes for 18 PQ analogs with known experimental potency. Subsequently, binding energy was calculated for all complexes using the molecular mechanics-generalized Born surface area method. Finally, docking poses were quantitatively assessed in the light of experimental data to derive a binding hypothesis. Collectively, the results indicate that PQs at the α1+/β3- site likely exhibit a common binding mode that can be characterized by a hydrogen bond interaction with β3Q64 and hydrophobic interactions involving residues α1F99, β3Y62, β3M115, α1Y159, and α1Y209. Importantly, our results are in good agreement with the recently resolved cryo-Electron Microscopy structures of the human α1β3γ2 and α1β2γ2 subtypes of GABA_A receptors.

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.

Structure-based drug repositioning over the human TMPRSS2 protease domain: search for chemical probes able to repress SARS-CoV-2 Spike protein cleavages

In December 2019, a new coronavirus was identified in the Hubei province of central china and named SARS-CoV-2. This new virus induces COVID-19, a severe respiratory disease with high death rate. A putative target to interfere with the virus is the host transmembrane serine protease family member II (TMPRSS2). This enzyme is critical for the entry of coronaviruses into human cells by cleaving and activating the spike protein (S) of SARS-CoV-2. Repositioning approved, investigational and experimental drugs on the serine protease domain of TMPRSS2 could thus be valuable. There is no experimental structure for TMPRSS2 but it is possible to develop quality structural models for the serine protease domain using comparative modeling strategies as such domains are highly structurally conserved. Beside the TMPRSS2 catalytic site, we predicted on our structural models a main exosite that could be important for the binding of protein partners and/or substrates. To block the catalytic site or the exosite of TMPRSS2 we used structure-based virtual screening computations and two different collections of approved, investigational and experimental drugs. We propose a list of 156 molecules that could bind to the catalytic site and 100 compounds that may interact with the exosite. These small molecules should now be tested in vitro to gain novel insights over the roles of TMPRSS2 or as starting point for the development of second generation analogs.

Structural and functional consequences of replacement of His403 with Arg near the catalytic site of Anoxybacillus flavithermus cyclomaltodextrinase

The hydrolytic activity of a thermophilic cyclomaltodextrinase (CMD) from Anoxybacillus flavithermus ZNU-NGA and a representative single mutant were investigated against soluble substrates including α-, β- and γ-cyclomaltodestrines (CDs). Based on the occurrence of arginine (Arg) at position 403 in some homologue proteins, His⁴⁰³ in Wild-type (WT) CMD was replaced with Arg (H403R variant) with site-directed mutagenesis procedures. According to bioinformatics data, Arg⁴⁰³ in mutant protein is located near Glu³⁵⁷ as one of the catalytic residues in a manner that they are able to create a medium-range attractive electrostatistic interaction. Structural studies by Far UV-CD showed that this mutation is accompanied by conversion of a small fraction of α-helix to β-form structure. Fluorescence data reveals that, the hydrophobic regions at the surface of protein, as the binding sites for ANS (8-Anilinonaphthalene-1-sulfonic acid) increase in mutant protein, demonstrating relative inflation of H403R variant compared with WT protein. However, the polarity of microenvironment around chromophores did not change upon mutation. Activity measurement in different ranges of pH and temperatures showed that the optimum values of pH and temperature in mutant enzyme is the same as WT enzyme, however; the activity at optimum points increased in H403R variant. It was also revealed that the H403R variant had slightly improved catalytic efficiency for γ-CD. The same value of activation parameters for both protein variants indicates that mutation does not alter the mechanism of catalysis during enzyme-substrate formation.

Structure-based drug target prioritisation and rational drug design for targeting Chlamydia trachomatis eye infections

Chlamydia trachomatis (C.t) is a major causative of infectious blindness in world. It is a real challenge to combat Chlamydial infection as it is an intracellular pathogen. Hence, it is essential to determine the most potential targets of C.t in order to inhibit or suppress its virulence during its infectious phase. Thus, in this study, the highly expressed-cum-most essential genes reported through our earlier study were reprioritized by structure-based comparative binding site analysis with host proteome. Therefore, computational approaches involving molecular modelling, large-scale binding site prediction and comparison, molecular dynamics simulation studies were performed to narrow down the most potential targets. Furthermore, high-throughput virtual screening and ADMETox were also performed to identify potential hits that shall efficiently inhibit the prioritised targets. Hence, by this study we report Pyruvoyl-dependent arginine decarboxylase (PvlArgDC), DNA-repair protein (RecO) and porin (outer membrane protein) as the most viable targets of C.t which can be potentially targeted by compounds, NSC_13086, MFCD00276409, MFCD05662003, respectively. AbbreviationsC.tChlamydia trachomatisSTDSexually transmitted diseaseHTVSHigh-throughput virtual screeningADMEToxAbsorption, Distribution, Metabolism, Excretion and ToxicityPMPocketMatchMDMolecular Dynamics simulationSPStandard precisionXPExtra precisionMMGBSAMolecular mechanics energies combined with generalised Born and surface area continuum solvationOMPOuter membrane proteinPvlArgDCPyruvoyl-dependent arginine decarboxylaseRecORecombination protein O.Communicated by Ramaswamy H. Sarma.

Structural Insights into Hearing Loss Genetics from Polarizable Protein Repacking

Hearing loss is associated with ∼8100 mutations in 152 genes, and within the coding regions of these genes are over 60,000 missense variants. The majority of these variants are classified as "variants of uncertain significance" to reflect our inability to ascribe a phenotypic effect to the observed amino acid change. A promising source of pathogenicity information is biophysical simulation, although input protein structures often contain defects because of limitations in experimental data and/or only distant homology to a template. Here, we combine the polarizable atomic multipole optimized energetics for biomolecular applications force field, many-body optimization theory, and graphical processing unit acceleration to repack all deafness-associated proteins and thereby improve average structure MolProbity score from 2.2 to 1.0. We then used these optimized wild-type models to create over 60,000 structures for missense variants in the Deafness Variation Database, which are being incorporated into the Deafness Variation Database to inform deafness pathogenicity prediction. Finally, this work demonstrates that advanced polarizable atomic multipole force fields are efficient enough to repack the entire human proteome.

Designing and construction of novel variants of Chondroitinase ABC I to reduce aggregation rate

Chondroitinase ABC I (cABC I) can degrade inhibitory molecules for axon regrowth at the site of damage after spinal cord injury (SCI). One of the main problems in the practical application is the possibility of structural changes that lead to the inactivation of the enzyme. In current work, three variants of cABC I was designed and constructed by manipulation of a short helix conformation (Gln⁶⁷⁸-Leu⁶⁷⁹-Ser⁶⁸⁰-Gln⁶⁸¹); where Gln residues were converted to Glu. According to the enzyme kinetics studies, the catalytic efficiency of the Q681E and double mutant (Q678E/Q681E) increases in comparison with WT enzyme; while that of Q678E decreases. It was also shown that the rate of the inactivation of the enzyme variants over time is greater in WT and Q678E variants than that of the Q681E and double mutant. Negative values of entropy change of thermal inactivation measurements; demonstrate that inactivation of the WT and Q678E variants are mainly originated from aggregation. These observations can be explained by considering the repulsive electrostatic interaction between enzyme molecules that prevents protein aggregation over time. It is concluded that increasing the solubility of the Q681E and double mutant via favorable interactions of surface-exposed charged residues with dipole momentum of water molecules accompanied by the presence of intermolecular repulsive electrostatic interaction leads to decreasing the rate of aggregation in both long-term storage and heat-induced structural changes.

Understanding Ligand Binding to G-Protein Coupled Receptors Using Multiscale Simulations

Human G-protein coupled receptors (GPCRs) convey a wide variety of extracellular signals inside the cell and they are one of the main targets for pharmaceutical intervention. Rational drug design requires structural information on these receptors; however, the number of experimental structures is scarce. This gap can be filled by computational models, based on homology modeling and docking techniques. Nonetheless, the low sequence identity across GPCRs and the chemical diversity of their ligands may limit the quality of these models and hence refinement using molecular dynamics simulations is recommended. This is the case for olfactory and bitter taste receptors, which constitute the first and third largest GPCR groups and show sequence identities with the available GPCR templates below 20%. We have developed a molecular dynamics approach, based on the combination of molecular mechanics and coarse grained (MM/CG), tailored to study ligand binding in GPCRs. This approach has been applied so far to bitter taste receptor complexes, showing significant predictive power. The protein/ligand interactions observed in the simulations were consistent with extensive mutagenesis and functional data. Moreover, the simulations predicted several binding residues not previously tested, which were subsequently verified by carrying out additional experiments. Comparison of the simulations of two bitter taste receptors with different ligand selectivity also provided some insights into the binding determinants of bitter taste receptors. Although the MM/CG approach has been applied so far to a limited number of GPCR/ligand complexes, the excellent agreement of the computational models with the mutagenesis and functional data supports the applicability of this method to other GPCRs for which experimental structures are missing. This is particularly important for the challenging case of GPCRs with low sequence identity with available templates, for which molecular docking shows limited predictive power.

Anti-HIV-1 integrase potency of methylgallate from Alchornea cordifolia using in vitro and in silico approaches

According to the 2018 report of the United Nations Programme on HIV/AIDS (UNAIDS), acquired immune deficiency syndrome (AIDS), a disease caused by the human immunodeficiency virus (HIV), remains a significant public health problem. The non-existence of a cure or effective vaccine for the disease and the associated emergence of resistant viral strains imply an urgent need for the discovery of novel anti-HIV drug candidates. The current study aimed to identify potential anti-retroviral compounds from Alchornea cordifolia. Bioactive compounds were identified using several chromatographic and spectroscopic techniques and subsequently evaluated for cytotoxicity and anti-HIV properties. Molecular modelling studies against HIV-1 integrase (HIV-1 IN) were performed to decipher the mode of action of methylgallate, the most potent compound (IC50 = 3.7 nM) and its analogues from ZINC database. Cytotoxicity assays showed that neither the isolated compounds nor the crude methanolic extract displayed cytotoxicity effects on the HeLa cell line. A strong correlation between the in vitro and in silico results was observed and important HIV-1 IN residues interacting with the different compounds were identified. These current results indicate that methylgallate is the main anti-HIV-1 compound in A. cordifolia stem bark, and could be a potential platform for the development of new HIV-1 IN inhibitors.

An IRAK1-PIN1 signalling axis drives intrinsic tumour resistance to radiation therapy

Drug-based strategies to overcome tumour resistance to radiotherapy (R-RT) remain limited by the single-agent toxicity of traditional radiosensitizers (e.g., platinums) and a lack of targeted alternatives. In a screen for compounds that restore radiosensitivity in p53 mutant zebrafish while tolerated in non-irradiated wild-type animals, we identified the benzimidazole anthelmintic, oxfendazole. Surprisingly, oxfendazole acts via inhibition of IRAK1, a kinase otherwise involved in Interleukin-1 and Toll-like receptor (IL-1R/TLR) immune responses. IRAK1 drives R-RT in a pathway involving IRAK4 and TRAF6 but not the IL-1R/TLR—IRAK adaptor MyD88. Rather than stimulating NF-κB, radiation-activated IRAK1 acts to prevent apoptosis mediated by the PIDDosome complex (PIDD/RAIDD/caspase-2). Countering this pathway with IRAK1 inhibitors suppresses R-RT in tumour models derived from cancers in which TP53 mutations predict R-RT. Lastly, IRAK1 inhibitors synergize with inhibitors of PIN1, a prolyl isomerase essential for IRAK1 activation in response to pathogens and, as shown here, ionizing radiation. These data identify an IRAK1 radiation-response pathway as a rational chemo-RT target.

Discovery of Potent Inhibitors for the Large Neutral Amino Acid Transporter 1 (LAT1) by Structure-Based Methods

The large neutral amino acid transporter 1 (LAT1) is a promising anticancer target that is required for the cellular uptake of essential amino acids that serve as building blocks for cancer growth and proliferation. Here, we report a structure-based approach to identify chemically diverse and potent inhibitors of LAT1. First, a homology model of LAT1 that is based on the atomic structures of the prokaryotic homologs was constructed. Molecular docking of nitrogen mustards (NMs) with a wide range of affinity allowed for deriving a common binding mode that could explain the structure−activity relationship pattern in NMs. Subsequently, validated binding hypotheses were subjected to molecular dynamics simulation, which allowed for extracting a set of dynamic pharmacophores. Finally, a library of ~1.1 million molecules was virtually screened against these pharmacophores, followed by docking. Biological testing of the 30 top-ranked hits revealed 13 actives, with the best compound showing an IC₅₀ value in the sub-μM range.

Enzymatic and non-enzymatic pathways of kynurenines' dimerization: the molecular factors for oxidative stress development

Kynurenines, the products of tryptophan oxidative degradation, are involved in multiple neuropathologies, such as Huntington's chorea, Parkinson's disease, senile dementia, etc. The major cause for hydroxykynurenines's neurotoxicity is the oxidative stress induced by the reactive oxygen species (ROS), the by-products of L-3-hydroxykynurenine (L-3HOK) and 3-hydroxyanthranilic acid (3HAA) oxidative self-dimerization. 2-aminophenol (2AP), a structural precursor of L-3HOK and 3HAA, undergoes the oxidative conjugation to form 2-aminophenoxazinone. There are several modes of 2AP dimerization, including both enzymatic and non-enzymatic stages. In this study, the free energies for 2AP, L-3HOK and 3HAA dimerization stages have been calculated at B3LYP/6-311G(d,p)//6-311+(O)+G(d) level, both in the gas phase and in heptane or water solution. For the intermediates, ionization potentials and electron affinities were calculated, as well as free energy and kinetics of molecular oxygen interaction with several non-enzymatically formed dimers. H-atom donating power of the intermediates increases upon the progress of the oxidation, making possible generation of hydroperoxyl radical or hydrogen peroxide from O2 at the last stages. Among the dimerization intermediates, 2-aminophenoxazinole derivatives have the lowest ionization potential and can reduce O2 to superoxide anion. The rate for O-H homolytic bond dissociation is significantly higher than that for C-H bond in non-enzymatic quinoneimine conjugate. However, the last reaction passes irreversibly, reducing O2 to hydroperoxyl radical. The inorganic ferrous iron and the heme group of Drosophila phenoxazinone synthase significantly reduce the energy cost of 2AP H-atom abstraction by O2. We have also shown experimentally that total antioxidant capacity decreases in Drosophila mutant cardinal with L-3HOK excess relative to the wild type Canton-S, and lipid peroxidation decreases in aged cardinal. Taken together, our data supports the conception of hydroxykynurenines' dual role in neurotoxicity: serving as antioxidants themselves, blocking lipid peroxidation by H-atom donation, they also can easily generate ROS upon dimerization, leading to the oxidative stress development.

Reoptimized UNRES Potential for Protein Model Quality Assessment

Ranking protein structure models is an elusive problem in bioinformatics. These models are evaluated on both the degree of similarity to the native structure and the folding pathway. Here, we simulated the use of the coarse-grained UNited RESidue (UNRES) force field as a tool to choose the best protein structure models for a given protein sequence among a pool of candidate models, using server data from the CASP11 experiment. Because the original UNRES was optimized for Molecular Dynamics simulations, we reoptimized UNRES using a deep feed-forward neural network, and we show that introducing additional descriptive features can produce better results. Overall, we found that the reoptimized UNRES performs better in selecting the best structures and tracking protein unwinding from its native state. We also found a relatively poor correlation between UNRES values and the model’s Template Modeling Score (TMS). This is remedied by reoptimization. We discuss some cases where our reoptimization procedure is useful.

Novel selective thiadiazine DYRK1A inhibitor lead scaffold with human pancreatic β-cell proliferation activity

The Dual-Specificity Tyrosine Phosphorylation-Regulated Kinase 1A (DYRK1A) is an enzyme that has been implicated as an important drug target in various therapeutic areas, including neurological disorders (Down syndrome, Alzheimer's disease), oncology, and diabetes (pancreatic β-cell expansion). Current small molecule DYRK1A inhibitors are ATP-competitive inhibitors that bind to the kinase in an active conformation. As a result, these inhibitors are promiscuous, resulting in pharmacological side effects that limit their therapeutic applications. None are in clinical trials at this time. In order to identify a new DYRK1A inhibitor scaffold, we constructed a homology model of DYRK1A in an inactive, DFG-out conformation. Virtual screening of 2.2 million lead-like compounds from the ZINC database, followed by in vitro testing of selected 68 compounds revealed 8 hits representing 5 different chemical classes. We chose to focus on one of the hits from the computational screen, thiadiazine 1 which was found to inhibit DYRK1A with IC50 of 9.41 μM (Kd = 7.3 μM). Optimization of the hit compound 1, using structure-activity relationship (SAR) analysis and in vitro testing led to the identification of potent thiadiazine analogs with significantly improved binding as compared to the initial hit (Kd = 71-185 nM). Compound 3-5 induced human β-cell proliferation at 5 μM while showing selectivity for DYRK1A over DYRK1B and DYRK2 at 10 μM. This newly developed DYRK1A inhibitor scaffold with unique kinase selectivity profiles has potential to be further optimized as novel therapeutics for diabetes.

Whole grain intake associated molecule 5-aminovaleric acid betaine decreases β-oxidation of fatty acids in mouse cardiomyocytes

Despite epidemiological evidence showing that diets rich in whole grains reduce the risk of chronic life-style related diseases, biological mechanisms for these positive effects are mostly unknown. Increased 5-aminovaleric acid betaine (5-AVAB) levels in plasma and metabolically active tissues such as heart have been associated with consumption of diets rich in whole grains. However, biological effects of 5-AVAB are poorly understood. We evaluated 5-AVAB concentrations in human and mouse heart tissue (3-22 µM and 38-78 µM, respectively) using mass spectrometry. We show that 5-AVAB, at physiological concentration range, dose-dependently inhibits oxygen consumption due to β-oxidation of fatty acids, but does not otherwise compromise mitochondrial respiration, as measured with oxygen consumption rate in cultured mouse primary cardiomyocytes. We also demonstrate that this effect is caused by 5-AVAB induced reduction of cellular L-carnitine. Reduced L-carnitine levels are at least partly mediated by the inhibition of cell membrane carnitine transporter (OCTN2) as evaluated by in silico docking, and by siRNA mediated silencing of OCTN2 in cultured cardiomyocytes. 5-AVAB caused inhibition of β-oxidation of fatty acids is a novel mechanism on how diets rich in whole grains may regulate energy metabolism in the body. Elucidating potentially beneficial effects of 5-AVAB e.g. on cardiac physiology will require further in vivo investigations.

Homology Modeling Informs Ligand Discovery for the Glutamine Transporter ASCT2

The Alanine-Serine-Cysteine transporter (SLC1A5, ASCT2), is a neutral amino acid exchanger involved in the intracellular homeostasis of amino acids in peripheral tissues. Given its role in supplying glutamine to rapidly proliferating cancer cells in several tumor types such as triple-negative breast cancer and melanoma, ASCT2 has been identified as a key drug target. Here we use a range of computational methods, including homology modeling and ligand docking, in combination with cell-based assays, to develop hypotheses for structure-function relationships in ASCT2. We perform a phylogenetic analysis of the SLC1 family and its prokaryotic homologs to develop a useful multiple sequence alignment for this protein family. We then generate homology models of ASCT2 in two different conformations, based on the human EAAT1 structures. Using ligand enrichment calculations, the ASCT2 models are then compared to crystal structures of various homologs for their utility in discovering ASCT2 inhibitors. We use virtual screening, cellular uptake and electrophysiology experiments to identify a non-amino acid ASCT2 inhibitor that is predicted to interact with the ASCT2 substrate binding site. Our results provide insights into the structural basis of substrate specificity in the SLC1 family, as well as a framework for the design of future selective and potent ASCT2 inhibitors as cancer therapeutics.

Allosteric Modulation of Human Hsp90α Conformational Dynamics

Central to Hsp90's biological function is its ability to interconvert between various conformational states. Drug targeting of Hsp90's regulatory mechanisms, including its modulation by cochaperone association, presents as an attractive therapeutic strategy for Hsp90 associated pathologies. In this study, we utilized homology modeling techniques to calculate full-length structures of human Hsp90α in closed and partially open conformations and used these structures as a basis for several molecular dynamics based analyses aimed at elucidating allosteric mechanisms and modulation sites in human Hsp90α. Atomistic simulations demonstrated that bound adenosine triphosphate (ATP) stabilizes the dimer by "tensing" each protomer, while adenosine diphosphate (ADP) and apo configurations "relax" the complex by increasing global flexibility, the former case resulting in a fully open "v-like" conformation. Dynamic residue network analysis revealed regions of the protein involved in intraprotein communication and identified several key communication hubs that correlate with known functional sites. Pairwise comparison of betweenness centrality, shortest path, and residue fluctuations revealed that a proportional relationship exists between the latter two measurables and an inverse relationship between these two and betweenness centrality. This analysis showed how protein flexibility, degree of compactness, and the distance cutoff used for network construction influence the correlations between these metrics. These findings are novel and suggest shortest path and betweenness centrality to be more relevant quantities to follow for detecting functional residues in proteins compared to residue fluctuations. Perturbation response scanning analysis identified several potential residue sites capable of modulating conformational change in favor of interstate conversion. For the ATP-bound open conformation, these sites were found to overlap with known Aha1 and client binding sites, demonstrating how naturally occurring forces associated with cofactor binding could allosterically modulate conformational dynamics.

The phosphorelay signal transduction system in Candida glabrata: an in silico analysis

Signaling systems allow microorganisms to sense and respond to different stimuli through the modification of gene expression. The phosphorelay signal transduction system in eukaryotes involves three proteins: a sensor protein, an intermediate protein and a response regulator, and requires the transfer of a phosphate group between two histidine-aspartic residues. The SLN1-YPD1-SSK1 system enables yeast to adapt to hyperosmotic stress through the activation of the HOG1-MAPK pathway. The genetic sequences available from Saccharomyces cerevisiae were used to identify orthologous sequences in Candida glabrata, and putative genes were identified and characterized by in silico assays. An interactome analysis was carried out with the complete genome of C. glabrata and the putative proteins of the phosphorelay signal transduction system. Next, we modeled the complex formed between the sensor protein CgSln1p and the intermediate CgYpd1p. Finally, phosphate transfer was examined by a molecular dynamic assay. Our in silico analysis showed that the putative proteins of the C. glabrata phosphorelay signal transduction system present the functional domains of histidine kinase, a downstream response regulator protein, and an intermediate histidine phosphotransfer protein. All the sequences are phylogenetically more related to S. cerevisiae than to C. albicans. The interactome suggests that the C. glabrata phosphorelay signal transduction system interacts with different proteins that regulate cell wall biosynthesis and responds to oxidative and osmotic stress the same way as similar systems in S. cerevisiae and C. albicans. Molecular dynamics simulations showed complex formation between the response regulator domain of histidine kinase CgSln1 and intermediate protein CgYpd1 in the presence of a phosphate group and interactions between the aspartic residue and the histidine residue. Overall, our research showed that C. glabrata harbors a functional SLN1-YPD1-SSK1 phosphorelay system.

Structural Characterization of the Trimerization of TRAF6 Protein Through Molecular Dynamics Simulations

The tumour necrosis factor (TNF) receptor-associated factor (TRAF) family of proteins having E3 ligase activity are the key molecules involved in cellular immune response pathways. TRAF6 is a unique member of the TRAF superfamily differing from other members of the family, owing to its specific interactions with molecules outside the TNF receptor superfamily. The C-terminal domain of TRAF proteins contains the catalytic residues and are known to be involved in self-oligomerization forming a mushroom-shaped trimeric structure, which is the functional form of the protein. However, the monomeric crystal structure of TRAF6 C-terminal domain has been already determined, but the trimeric structure of the same is still not available. We here applied computational structural modelling and molecular dynamics simulations studies to get insights into the molecular interactions involved in determining the trimeric structure of the TRAF6 C-terminal domain. The non-availability of the trimeric structure of the TRAF6 C-terminal domain prevented the elucidation of the molecular mechanism of many different biological processes. Our results suggest that the trimer complex is transient in nature. The amino acid residues Lys340 and Glu345 in the coiled coil domain in the C-terminus of TRAF6 play a critical role in trimer structure formation. This structural modelling study may therefore be utilized to obtain the experimentally validated trimeric structure of this important protein.

Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis

Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.

An overview of comparative modelling and resources dedicated to large-scale modelling of genome sequences

Computational modelling of proteins has been a major catalyst in structural biology. Bioinformatics groups have exploited the repositories of known structures to predict high-quality structural models with high efficiency at low cost. This article provides an overview of comparative modelling, reviews recent developments and describes resources dedicated to large-scale comparative modelling of genome sequences. The value of subclustering protein domain superfamilies to guide the template-selection process is investigated. Some recent cases in which structural modelling has aided experimental work to determine very large macromolecular complexes are also cited.

Discovery of Competitive and Noncompetitive Ligands of the Organic Cation Transporter 1 (OCT1; SLC22A1)

Organic cation transporter 1 (OCT1) plays a critical role in the hepatocellular uptake of structurally diverse endogenous compounds and xenobiotics. Here we identified competitive and noncompetitive OCT1-interacting ligands in a library of 1780 prescription drugs by combining in silico and in vitro methods. Ligands were predicted by docking against a comparative model based on a eukaryotic homologue. In parallel, high-throughput screening (HTS) was conducted using the fluorescent probe substrate ASP⁺ in cells overexpressing human OCT1. Thirty competitive OCT1 ligands, defined as ligands predicted in silico as well as found by HTS, were identified. Of the 167 ligands identified by HTS, five were predicted to potentially cause clinical drug interactions. Finally, virtual screening of 29 332 metabolites predicted 146 competitive OCT1 ligands, of which an endogenous neurotoxin, 1-benzyl-1,2,3,4-tetrahydroisoquinoline, was experimentally validated. In conclusion, by combining docking and in vitro HTS, competitive and noncompetitive ligands of OCT1 can be predicted.

Protein Structure Modeling with MODELLER

Genome sequencing projects have resulted in a rapid increase in the number of known protein sequences. In contrast, only about one-hundredth of these sequences have been characterized at atomic resolution using experimental structure determination methods. Computational protein structure modeling techniques have the potential to bridge this sequence-structure gap. In the following chapter, we present an example that illustrates the use of MODELLER to construct a comparative model for a protein with unknown structure. Automation of a similar protocol has resulted in models of useful accuracy for domains in more than half of all known protein sequences.

Predicting protein conformational changes for unbound and homology docking: learning from intrinsic and induced flexibility

Predicting protein conformational changes from unbound structures or even homology models to bound structures remains a critical challenge for protein docking. Here we present a study directly addressing the challenge by reducing the dimensionality and narrowing the range of the corresponding conformational space. The study builds on cNMA-our new framework of partner- and contact-specific normal mode analysis that exploits encounter complexes and considers both intrinsic and induced flexibility. First, we established over a CAPRI (Critical Assessment of PRedicted Interactions) target set that the direction of conformational changes from unbound structures and homology models can be reproduced to a great extent by a small set of cNMA modes. In particular, homology-to-bound interface root-mean-square deviation (iRMSD) can be reduced by 40% on average with the slowest 30 modes. Second, we developed novel and interpretable features from cNMA and used various machine learning approaches to predict the extent of conformational changes. The models learned from a set of unbound-to-bound conformational changes could predict the actual extent of iRMSD with errors around 0.6 Å for unbound proteins in a held-out benchmark subset, around 0.8 Å for unbound proteins in the CAPRI set, and around 1 Å even for homology models in the CAPRI set. Our results shed new insights into origins of conformational differences between homology models and bound structures and provide new support for the low-dimensionality of conformational adjustment during protein associations. The results also provide new tools for ensemble generation and conformational sampling in unbound and homology docking. Proteins 2017; 85:544-556. © 2016 Wiley Periodicals, Inc.

Resurrecting ancestral structural dynamics of an antiviral immune receptor: adaptive binding pocket reorganization repeatedly shifts RNA preference

Background

Although resurrecting ancestral proteins is a powerful tool for understanding the molecular-functional evolution of gene families, nearly all studies have examined proteins functioning in relatively stable biological processes. The extent to which more dynamic systems obey the same ‘rules’ governing stable processes is unclear. Here we present the first detailed investigation of the functional evolution of the RIG-like receptors (RLRs), a family of innate immune receptors that detect viral RNA in the cytoplasm.

Results

Using kinetic binding assays and molecular dynamics simulations of ancestral proteins, we demonstrate how a small number of adaptive protein-coding changes repeatedly shifted the RNA preference of RLRs throughout animal evolution by reorganizing the shape and electrostatic distribution across the RNA binding pocket, altering the hydrogen bond network between the RLR and its RNA target. In contrast to observations of proteins involved in metabolism and development, we find that RLR-RNA preference ‘flip flopped’ between two functional states, and shifts in RNA preference were not always coupled to gene duplications or speciation events. We demonstrate at least one reversion of RLR-RNA preference from a derived to an ancestral function through a novel structural mechanism, indicating multiple structural implementations of similar functions.

Conclusions

Our results suggest a model in which frequent shifts in selection pressures imposed by an evolutionary arms race preclude the long-term functional optimization observed in stable biological systems. As a result, the evolutionary dynamics of immune receptors may be less constrained by structural epistasis and historical contingency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0818-6) contains supplementary material, which is available to authorized users.