Comparative Protein Structure Modeling Using MODELLER

Protein interactions with metallothionein-3 promote vectorial active transport in human proximal tubular cells

Metallothionein 3 (MT-3) is a small, cysteine-rich protein that binds to essential metals required for homeostasis, as well as to heavy metals that have the potential to exert toxic effects on cells. MT-3 is expressed by epithelial cells of the human kidney, including the cells of the proximal tubule. Our laboratory has previously shown that mortal cultures of human proximal tubular (HPT) cells express MT-3 and form domes in the cell monolayer, a morphological feature indicative of vectorial active transport, an essential function of the proximal tubule. However, an immortalized proximal tubular cell line HK-2 lacks the expression of MT-3 and fails to form domes in the monolayer. Transfection of HK-2 cells with the MT-3 gene restores dome formation in these cells suggesting that MT-3 is required for vectorial active transport. In order to determine how MT-3 imparts this essential feature to the proximal tubule, we sought to identify proteins that interact either directly or indirectly with MT-3. Using a combination of pulldowns, co-immunoprecipitations, and mass spectrometry analysis, putative protein interactants were identified and subsequently confirmed by Western analysis and confocal microscopy, following which proteins with direct physical interactions were investigated through molecular docking. Our data shows that MT-3 interacts with myosin-9, aldolase A, enolase 1, β-actin, and tropomyosin 3 and that these interactions are maximized at the periphery of the apical membrane of doming proximal tubule cells. Together these observations reveal that MT-3 interacts with proteins involved in cytoskeletal organization and energy metabolism, and these interactions at the apical membrane support vectorial active transport and cell differentiation in proximal tubule cultures.

Migraine-Associated Mutation in the Na,K-ATPase Leads to Disturbances in Cardiac Metabolism and Reduced Cardiac Function

Background Mutations in ATP1A2 gene encoding the Na,K-ATPase α2 isoform are associated with familial hemiplegic migraine type 2. Migraine with aura is a known risk factor for heart disease. The Na,K-ATPase is important for cardiac function, but its role for heart disease remains unknown. We hypothesized that ATP1A2 is a susceptibility gene for heart disease and aimed to assess the underlying disease mechanism. Methods and Results Mice heterozygous for the familial hemiplegic migraine type 2-associated G301R mutation in the Atp1a2 gene (α2+/G301R mice) and matching wild-type controls were compared. Reduced expression of the Na,K-ATPase α2 isoform and increased expression of the α1 isoform were observed in hearts from α2+/G301R mice (Western blot). Left ventricular dilation and reduced ejection fraction were shown in hearts from 8-month-old α2+/G301R mice (cardiac magnetic resonance imaging), and this was associated with reduced nocturnal blood pressure (radiotelemetry). Cardiac function and blood pressure of 3-month-old α2+/G301R mice were similar to wild-type mice. Amplified Na,K-ATPase-dependent Src kinase/Ras/Erk1/2 (p44/42 mitogen-activated protein kinase) signaling was observed in hearts from 8-month-old α2+/G301R mice, and this was associated with mitochondrial uncoupling (respirometry), increased oxidative stress (malondialdehyde measurements), and a heart failure-associated metabolic shift (hyperpolarized magnetic resonance). Mitochondrial membrane potential (5,5´,6,6´-tetrachloro-1,1´,3,3´-tetraethylbenzimidazolocarbocyanine iodide dye assay) and mitochondrial ultrastructure (transmission electron microscopy) were similar between the groups. Proteomics of heart tissue further suggested amplified Src/Ras/Erk1/2 signaling and increased oxidative stress and provided the molecular basis for systolic dysfunction in 8-month-old α2+/G301R mice. Conclusions Our findings suggest that ATP1A2 mutation leads to disturbed cardiac metabolism and reduced cardiac function mediated via Na,K-ATPase-dependent reactive oxygen species signaling through the Src/Ras/Erk1/2 pathway.

Pregnane-based steroids are novel positive NMDA receptor modulators that may compensate for the effect of loss-of-function disease-associated GRIN mutations

BACKGROUND AND PURPOSE

N-methyl-D-aspartate receptors (NMDARs) play a critical role in synaptic plasticity, and mutations in human genes encoding NMDAR subunits have been described in individuals with various neuropsychiatric disorders. Compounds with a positive allosteric effect are thought to compensate for reduced receptor function.

EXPERIMENTAL APPROACH

We have used whole-cell patch-clamp electrophysiology on recombinant rat NMDARs and human variants found in individuals with neuropsychiatric disorders, in combination with in silico modelling, to explore the site of action of novel epipregnanolone-based NMDAR modulators.

KEY RESULTS

Analysis of the action of 4-(20-oxo-5β-pregnan-3β-yl) butanoic acid (EPA-But) at the NMDAR indicates that the effect of this steroid with a "bent" structure is different from that of cholesterol and oxysterols and shares a disuse-dependent mechanism of NMDAR potentiation with the "planar" steroid 20-oxo-pregn-5-en-3β-yl sulfate (PE-S). The potentiating effects of EPA-But and PE-S are additive. Alanine scan mutagenesis identified residues that reduce the potentiating effect of EPA-But. No correlation was found between the effects of EPA-But and PE-S at mutated receptors that were less sensitive to either steroid. The relative degree of potentiation induced by the two steroids also differed in human NMDARs carrying rare variants of hGluN1 or hGluN2B subunits found in individuals with neuropsychiatric disorders, including intellectual disability, epilepsy, developmental delay, and autism spectrum disorder.

CONCLUSION AND IMPLICATIONS

Our results show novel sites of action for pregnanolones at the NMDAR and provide an opportunity for the development of new therapeutic neurosteroid-based ligands to treat diseases associated with glutamatergic system hypofunction.

Binding Affinity Characterization of Four Antennae-Enriched Odorant-Binding Proteins From Harmonia axyridis (Coleoptera: Coccinellidae)

Harmonia axyridis is an important natural enemy that consumes many agricultural and forestry pests. It relies on a sensitive olfactory system to find prey and mates. Odorant-binding proteins (OBPs) as the first-step of recognizing volatiles, transport odors through sensillum lymph to odorant receptors (ORs). However, little is known about the molecular mechanisms of H. axyridis olfaction. In this study, four H. axyridis antenna specific OBP genes, HaxyOBP3, 5, 12, and 15, were bacterially expressed and the binding features of the four recombinant proteins to 40 substances were investigated using fluorescence competitive binding assays. Three-dimensional structure modeling and molecular docking analysis predicted the binding sites between HaxyOBPs and candidate volatiles. Developmental expression analyses showed that the four HaxyOBP genes displayed a variety of expression patterns at different development stages. The expression levels of HaxyOBP3 and HaxyOBP15 were higher in the adult stage than in the other developmental stages, and HaxyOBP15 was significantly transcriptionally enriched in adult stage. Ligand-binding analysis demonstrated that HaxyOBP3 and HaxyOBP12 only combined with two compounds, β-ionone and p-anisaldehyde. HaxyOBP5 protein displayed binding affinities with methyl salicylate, β-ionone, and p-anisaldehyde (Ki = 18.15, 11.71, and 13.45 μM). HaxyOBP15 protein had a broad binding profile with (E)-β-farnesene, β-ionone, α-ionone, geranyl acetate, nonyl aldehyde, dihydro-β-ionone, and linalyl acetate (Ki = 4.33-31.01 μM), and hydrophobic interactions played a key role in the binding of HaxyOBP15 to these substances according to molecular docking. Taken together, HaxyOBP15 exhibited a broader ligand-binding spectrum and a higher expression in adult stage than HaxyOBP3, 5, and 12, indicating HaxyOBP15 may play a greater role in binding volatiles than other three HaxyOBPs. The results will increase our understanding of the molecular mechanism of H. axyridis olfaction and may also result in new management strategies (attractants/repellents) that increase the biological control efficacy of H. axyridis.

An in-frame deletion mutation in the degron tail of auxin coreceptor IAA2 confers resistance to the herbicide 2,4-D in Sisymbrium orientale

Synthetic auxin herbicides intersect basic plant developmental biology and applied weed management. We investigated resistance to 2,4-D in the Australian weed Sisymbrium orientale (Indian hedge mustard). We identified a mechanism involving an in-frame 27-bp deletion in the degron tail of auxin coreceptor IAA2, one member of the gene family of Aux/IAA auxin co-receptors. We show that this deletion in IAA2 is a gain-of-function mutation that confers synthetic auxin resistance. This field-evolved mechanism of resistance to synthetic auxin herbicides confirms previous biochemical studies showing the role of the Aux/IAA degron tail in regulating Aux/IAA protein degradation upon auxin perception. The deletion mutation could be generated in crops using gene-editing approaches for cross-resistance to multiple synthetic auxin herbicides.

The natural auxin indole-3-acetic acid (IAA) is a key regulator of many aspects of plant growth and development. Synthetic auxin herbicides such as 2,4-D mimic the effects of IAA by inducing strong auxinic-signaling responses in plants. To determine the mechanism of 2,4-D resistance in a Sisymbrium orientale (Indian hedge mustard) weed population, we performed a transcriptome analysis of 2,4-D-resistant (R) and -susceptible (S) genotypes that revealed an in-frame 27-nucleotide deletion removing nine amino acids in the degron tail (DT) of the auxin coreceptor Aux/IAA2 (SoIAA2). The deletion allele cosegregated with 2,4-D resistance in recombinant inbred lines. Further, this deletion was also detected in several 2,4-D-resistant field populations of this species. Arabidopsis transgenic lines expressing the SoIAA2 mutant allele were resistant to 2,4-D and dicamba. The IAA2-DT deletion reduced binding to TIR1 in vitro with both natural and synthetic auxins, causing reduced association and increased dissociation rates. This mechanism of synthetic auxin herbicide resistance assigns an in planta function to the DT region of this Aux/IAA coreceptor for its role in synthetic auxin binding kinetics and reveals a potential biotechnological approach to produce synthetic auxin-resistant crops using gene-editing.

Understanding the influence of AMG 510 on the structure of KRASG12C empowered by molecular dynamics simulation

The KRASG12C mutant is often associated with human cancers, and AMG 510 as a promising covalent inhibitor of KRASG12C has achieved surprising efficacy in clinical trials. However, the interaction mechanism between KRASG12C and AMG 510 is not completely understood. Here, we performed all-atom molecular dynamics simulations on the complex of KRASG12C-AMG 510 to explore the influence of this covalent inhibitor on the conformational change of KRASG12C. A PCA (Principal Component Analysis) model was constructed based on known KRAS crystal structures to distinguish different conformations (active, inactive, and other). By mapping simulation trajectories onto the PCA model, we observed that the conformations of KRASG12C bound with AMG 510 were mainly concentrated in the inactive conformation. Further analysis demonstrated that AMG 510 reduced the flexibility of two switch regions to make the complex of KRASG12C-AMG 510 restricted in the inactive conformation. In the meantime, we also identified key interacting residues between KRASG12C and AMG 510 through the calculation of binding energy. Finally, we built a series of KRAS second-site mutation systems (i.e. KRASG12C/mutations) to conduct large-scale screening of potential resistance mutations. By further combining MD simulations and the PCA model, we not only recapitulated the currently known resistance mutations of AMG 510 successfully but also proposed some novel potential resistant mutations. Taken together, these results broaden our insight into the influence of AMG 510 on the conformational change of the KRASG12C mutant at the atomic level, thereby providing crucial hints for the improvement and optimization of drug candidates.

Structure-Based Design of a Chemical Probe Set for the 5-HT5A Serotonin Receptor

The 5-HT_5A receptor (5-HT_5AR), for which no selective agonists and a few antagonists exist, remains the least understood serotonin receptor. A single commercial antagonist, SB-699551, has been widely used to investigate the 5-HT_5AR function in neurological disorders, including pain, but this molecule has substantial liabilities as a chemical probe. Accordingly, we sought to develop an internally controlled probe set. Docking over 6 million molecules against a 5-HT_5AR homology model identified 5 mid-μM ligands, one of which was optimized to UCSF678, a 42 nM arrestin-biased partial agonist at the 5-HT_5AR with a more restricted off-target profile and decreased assay liabilities versus SB-699551. Site-directed mutagenesis supported the docked pose of UCSF678. Surprisingly, analogs of UCSF678 that lost the 5-HT_5AR activity revealed that 5-HT_5AR engagement is nonessential for alleviating pain, contrary to studies with less-selective ligands. UCSF678 and analogs constitute a selective probe set with which to study the function of the 5-HT_5AR.

Sequence Analysis, Structure Prediction of Receptor Proteins and In Silico Study of Potential Inhibitors for Management of Life Threatening COVID-19

Background:

Treatment of the Covid-19 pandemic caused by the highly contagious and pathogenic SARS-CoV-2 is a global menace. Day by day, this pandemic is getting worse. Doctors, scientists and researchers across the world are urgently scrambling for a cure for novel corona virus and continuously working at break neck speed to develop vaccines or drugs. But to date, there are no specific drugs or vaccines available in the market to cope up with the virus.

Objective:

The present study helps us to elucidate 3D structures of SARS-CoV-2 proteins and also to identify natural compounds as potential inhibitors against COVID-19.

Methods:

The 3D structures of the proteins were constructed using Modeller 9.16 modeling tool. Modelled proteins were validated with PROCHECK by Ramachandran plot analysis. In this study, a small library of natural compounds (fifty compounds) was docked to the hACE2 binding site of the modelled surface glycoprotein of SARS-CoV-2 using AutoDock Vina to repurpose these inhibitors against SARS-CoV-2. Conceptual density functional theory calculations of the best eight compounds had been performed by Gaussian-09. Geometry optimizations for these molecules were done at M06-2X/ def2-TZVP level of theory. ADME parameters, pharmacokinetic properties and drug likeness of the compounds were analyzed using swissADME website.

Results:

In this study, we analysed the sequences of surface glycoprotein, nucleocapsid phosphoprotein and envelope protein obtained from different parts of the globe. We modelled all the different sequences of surface glycoprotein and envelop protein in order to derive 3D structure of a molecular target, which is essential for the development of therapeutics. Different electronic properties of the inhibitors have been calculated using DFT through M06-2X functional with def2-TZVP basis set. Docking result at the hACE2 binding site of all modelled surface glycoproteins of SARSCoV- 2 showed that all the eight inhibitors (actinomycin D, avellanin C, ichangin, kanglemycin A, obacunone, ursolic acid, ansamiotocin P-3 and isomitomycin A) studied here were many folds better compared to hydroxychloroquine which has been found to be effective to treat patients suffering from COVID-19. All the inhibitors meet most of the criteria of drug likeness assessment.

Conclusion:

We expect that eight compounds (actinomycin D, avellanin C, ichangin, kanglemycin A, obacunone, ursolic acid, ansamiotocin P-3 and isomitomycin A) can be used as potential inhibitors against SARS-CoV-2.

Potential Drug Targets Identification Against Clostridioides Difficile (CD) and Characterization of Indispensable Proteins by a Subtractive Genomics Approach Followed by Virtual Screening

Background:

Clostridioides difficile (CD) is an enteric multi-drug resistant pathogenic bacterium. CD-associated infections are the leading cause of nosocomial diarrhea that can further lead to pseudomembranous colitis, toxic mega-colon or sepsis with greater mortality and morbidity risks. CD infection possesses higher rates of recurrence due to its greater resistance to antibiotics. Considering its higher rates of recurrence, it has become a major burden on healthcare facilities. Therefore, there is a dire need to identify novel drug targets to combat antibiotic resistance of Clostridioides difficile.

Objective:

To identify and propose new and novel drug targets against the Clostridioides difficile.

Methods:

In the current study, a computational subtractive genomics approach was applied to obtain a set of potential drug targets that exist in the multi-drug resistant strain of Clostridioides difficile. Here, the uncharacterized proteins were studied as potential drug targets. The methodology involved several bioinformatics databases and tools. The druggable proteins sequences were retrieved based on non-homology with host proteome and essentiality for the survival of the pathogen. The uncharacterized proteins were functionally characterized using different computational tools, and sub-cellular localization was also predicted. The metabolic pathways were analyzed using the KEGG database. Eventually, the druggable proteome has been fetched using sequence similarity with the already available drug targets present in the DrugBank database. These druggable proteins were further explored for the structural details to identify drug candidates.

Results :

A priority list of potential drug targets was provided with the help of the applied method on the complete proteome set of the C. difficile. Moreover, the drug-like compounds have been screened against the potential drug targets to prioritize potential drug candidates. To facilitate the need for drug targets and therapies, the study proposed five potential protein drug targets, out of which three proposed drug targets were subjected to homology modeling to explore their structural and functional activities

Conclusion:

In conclusion, we proposed three unique, unexplored drug targets against C. difficile. The structure-based methods were applied and resulted in a list of top-scoring compounds as potential inhibitors to proposed drug targets.

Structural Elements Directing G Proteins and β-Arrestin Interactions with the Human Melatonin Type 2 Receptor Revealed by Natural Variants

G protein-coupled receptors (GPCRs) can engage distinct subsets of signaling pathways, but the structural determinants of this functional selectivity remain elusive. The naturally occurring genetic variants of GPCRs, selectively affecting different pathways, offer an opportunity to explore this phenomenon. We previously identified 40 coding variants of the MTNR1B gene encoding the melatonin MT₂ receptor (MT₂). These mutations differently impact the β-arrestin 2 recruitment, ERK activation, cAMP production, and Gα_i1 and Gα_z activation. In this study, we combined functional clustering and structural modeling to delineate the molecular features controlling the MT₂ functional selectivity. Using non-negative matrix factorization, we analyzed the signaling signatures of the 40 MT₂ variants yielding eight clusters defined by unique signaling features and localized in distinct domains of MT₂. Using computational homology modeling, we describe how specific mutations can selectively affect the subsets of signaling pathways and offer a proof of principle that natural variants can be used to explore and understand the GPCR functional selectivity.

Biochemical and structural characterization of quizalofop-resistant wheat acetyl-CoA carboxylase

A novel nucleotide mutation in ACC1 resulting in an alanine to valine amino acid substitution in acetyl-CoA carboxylase (ACCase) at position 2004 of the Alopecurus myosuroides reference sequence (A2004V) imparts quizalofop resistance in wheat. Genotypes endowed with the homozygous mutation in one or two ACC1 homoeologs are seven- and 68-fold more resistant to quizalofop than a wildtype winter wheat in greenhouse experiments, respectively. In vitro ACCase activities in soluble protein extracts from these varieties are 3.8- and 39.4-fold more resistant to quizalofop with the homozygous mutation in either one or two genomes, relative to the wildtype. The A2004V mutation does not alter the specific activity of wheat ACCase, suggesting that this resistance trait does not affect the catalytic functions of ACCase. Modeling of wildtype and quizalofop-resistant wheat ACCase demonstrates that the A2004V amino acid substitution causes a reduction in the volume of the binding pocket that hinders quizalofop’s interaction with ACCase. Docking studies confirm that the mutation reduces the binding affinity of quizalofop. Interestingly, the models suggest that the A2004V mutation does not affect haloxyfop binding. Follow up in vivo and in vitro experiments reveal that the mutation, in fact, imparts negative cross-resistance to haloxyfop, with quizalofop-resistant varieties exhibiting higher sensitivity to haloxyfop than the wildtype winter wheat line.

Predicting Protein-Peptide Complex Structures by Accounting for Peptide Flexibility and the Physicochemical Environment

Predicting protein-peptide complex structures is crucial to the understanding of a vast variety of peptide-mediated cellular processes and to peptide-based drug development. Peptide flexibility and binding mode ranking are the two major challenges for protein-peptide complex structure prediction. Peptides are highly flexible molecules, and therefore, brute-force modeling of peptide conformations of interest in protein-peptide docking is beyond current computing power. Inspired by the fact that the protein-peptide binding process is like protein folding, we developed a novel strategy, named MDockPeP2, which tries to address these challenges using physicochemical information embedded in abundant monomeric proteins with an exhaustive search strategy, in combination with an integrated global search and a local flexible minimization method. Only the peptide sequence and the protein crystal structure are required. The method was systemically assessed using a newly constructed structural database of 89 nonredundant protein-peptide complexes with the peptide sequence length ranging from 5 to 29 in which about half of the peptides are longer than 15 residues. MDockPeP2 yielded a total success rate of 58.4% (70.8, 79.8%) for the bound docking (i.e., with the bound receptor and fully flexible peptides) and 19.0% (44.8, 70.7%) for the challenging unbound docking when top 10 (100, 1000) models were considered for each prediction. MDockPeP2 achieved significantly higher success rates on two other datasets, peptiDB and LEADS-PEP, which contain only short- and medium-size peptides (≤ 15 residues). For peptiDB, our method obtained a success rate of 62.0% for the bound docking and 35.9% for the unbound docking when the top 10 models were considered. For LEADS-PEP, MDockPeP2 achieved a success rate of 69.8% when the top 10 models were considered. The program is available at https://zougrouptoolkit.missouri.edu/mdockpep2/download.html.

Q61 mutant-mediated dynamics changes of the GTP-KRAS complex probed by Gaussian accelerated molecular dynamics and free energy landscapes

Understanding the molecular mechanism of the GTP-KRAS binding is significant for improving the target roles of KRAS in cancer treatment. In this work, multiple replica Gaussian accelerated molecular dynamics (MR-GaMD) simulations were applied to decode the effect of Q61A, Q61H and Q61L on the activity of KRAS. Dynamics analyses based on MR-GaMD trajectory reveal that motion modes and dynamics behavior of the switch domain in KRAS are heavily affected by the three Q61 mutants. Information of free energy landscapes (FELs) shows that Q61A, Q61H and Q61L induce structural disorder of the switch domain and disturb the activity of KRAS. Analysis of the interaction network uncovers that the decrease in the stability of hydrogen bonding interactions (HBIs) of GTP with residues V29 and D30 induced by Q61A, Q61H and Q61L is responsible for the structural disorder of the switch-I and that in the occupancy of the hydrogen bond between GTP and residue G60 leads to the structural disorder of the switch-II. Thus, the high disorder of the switch domain caused by three current Q61 mutants produces a significant effect on binding of KRAS to its effectors. This work is expected to provide useful information for further understanding function and target roles of KRAS in anti-cancer drug development.

Preparation of Immune and Synthetic VNAR Libraries as Sources of High-Affinity Binders

The shark-derived autonomous variable antibody domains known as VNARs are attractive tools for therapeutic and diagnostic applications due to their favorable properties like small size (approximately 12 kDa), high thermal and chemical stability, and good tissue penetration. Currently, different techniques have been reported to generate VNAR domains against targets of therapeutic interest. Here, we describe methods for the preparation of an immune VNAR library based on bacteriophage display, and for the preparation of a synthetic library of VNAR domains using a modified protocol based on Kunkel mutagenesis. Finally, we describe procedures for in silico maturation of a VNAR using a bioinformatic approach to obtain higher affinity binders.

P17 induces chemotaxis and differentiation of monocytes via MRGPRX2-mediated mast cell-line activation

BACKGROUND

P17, a peptide isolated from Tetramorium bicarinatum ant venom, is known to induce an alternative phenotype of human monocyte-derived macrophages via activation of an unknown G protein-coupled receptor (GPCR).

OBJECTIVE

We sought to investigate the mechanism of action and the immunomodulatory effects of P17 mediated through MRGPRX2 (Mas-related G protein-coupled receptor X2).

METHODS

To identify the GPCR for P17, we screened 314 GPCRs. Upon identification of MRGPRX2, a battery of in silico, in vitro, ex vivo, and in vivo assays along with the receptor mutation studies were performed. In particular, to investigate the immunomodulatory actions, we used β-hexosaminidase release assay, cytokine releases, quantification of mRNA expression, cell migration and differentiation assays, immunohistochemical labeling, hematoxylin and eosin, and immunofluorescence staining.

RESULTS

P17 activated MRGPRX2 in a dose-dependent manner in β-arrestin recruitment assay. In LAD2 cells, P17 induced calcium and β-hexosaminidase release. Quercetin- and short hairpin RNA-mediated knockdown of MRGPRX2 reduced P17-evoked β-hexosaminidase release. In silico and in vitro mutagenesis studies showed that residue Lys8 of P17 formed a cation-π interaction with the Phe172 of MRGPRX2 and [Ala8]P17 lost its activity partially. P17 activated LAD2 cells to recruit THP-1 and human monocytes in Transwell migration assay, whereas MRGPRX2-impaired LAD2 cells cannot. In addition, P17-treated LAD2 cells stimulated differentiation of THP-1 and human monocytes, as indicated by the enhanced expression of macrophage markers cluster of differentiation 11b and TNF-α by quantitative RT-PCR. Immunohistochemical and immunofluorescent staining suggested monocyte recruitment in mice ears injected with P17.

CONCLUSIONS

Our data provide novel structural information regarding the interaction of P17 with MRGPRX2 and intracellular pathways for its immunomodulatory action.

Cytokinin Perception in Ancient Plants beyond Angiospermae

Cytokinins (CKs) control many plant developmental processes and responses to environmental cues. Although the CK signaling is well understood, we are only beginning to decipher its evolution. Here, we investigated the CK perception apparatus in early-divergent plant species such as bryophyte Physcomitrium patens, lycophyte Selaginella moellendorffii, and gymnosperm Picea abies. Of the eight CHASE-domain containing histidine kinases (CHKs) examined, two CHKs, PpCHK3 and PpCHK4, did not bind CKs. All other CHK receptors showed high-affinity CK binding (K_D of nM range), with a strong preference for isopentenyladenine over other CK nucleobases in the moss and for trans-zeatin over cis-zeatin in the gymnosperm. The pH dependences of CK binding for these six CHKs showed a wide range, which may indicate different subcellular localization of these receptors at either the plasma- or endoplasmic reticulum membrane. Thus, the properties of the whole CK perception apparatuses in early-divergent lineages were demonstrated. Data show that during land plant evolution there was a diversification of the ligand specificity of various CHKs, in particular, the rise in preference for trans-zeatin over cis-zeatin, which indicates a steadily increasing specialization of receptors to various CKs. Finally, this distinct preference of individual receptors to different CK versions culminated in vascular plants, especially angiosperms.

Conformational exploration of RbgA using molecular dynamics: Possible implications in ribosome maturation and GTPase activity in different nucleotide bound states

Ribosome biogenesis GTPase A (RbgA) is involved in the late steps of the 50S ribosomal subunit maturation by binding into the 45S pre-ribosomal subunit. The association of RbgA to the 45S intermediate subunit depends on its bound nucleotide (GTP/GDP), probably because of the conformational shifts that occur between the GTP and GDP bound states. However, the available crystal structures of Staphylococcus aureus RbgA (SaRbgA) do not show any significant variations between different nucleotide bound states. Therefore, conformational exploration of SaRbgA in different nucleotide bound states was carried out using all-atom molecular dynamics (MD) simulations. Exploration of conformational distribution using cluster analysis and principal component analysis (PCA) revealed that GDP and pppGpp bound systems exhibit a larger distribution. This is majorly due to the fluctuations of the C-terminal tail (C-tail) as a result of the unwinding of α-helical secondary conformations into loop conformations which are observed from RMSF and DSSP analyses. Further investigation of the network of interactions revealed that the GTP and GMPPNP bound systems hold the C-tail in an α-helical form through stronger interactions between the active-site and C-tail. We also find that the presence of Mg²⁺ positions Sw-I loop away from the bound nucleotide and stabilizes the active-site water molecules. This seems to assist SaRbgA GTPase activity. In addition, mutations at the C-terminal and Sw-II conserved residues exhibit a larger conformational distribution majorly due to the C-tail fluctuations suggesting that the C-tail of SaRbgA probably interacts with the rRNA or rprotein in the process of ribosome biogenesis.

Paradigmatic De Novo GRIN1 Variants Recapitulate Pathophysiological Mechanisms Underlying GRIN1-Related Disorder Clinical Spectrum

Background: GRIN-related disorders (GRD), the so-called grinpathies, is a group of rare encephalopathies caused by mutations affecting GRIN genes (mostly GRIN1, GRIN2A and GRIN2B genes), which encode for the GluN subunit of the N-methyl D-aspartate (NMDA) type ionotropic glutamate receptors. A growing number of functional studies indicate that GRIN-encoded GluN1 subunit disturbances can be dichotomically classified into gain- and loss-of-function, although intermediate complex scenarios are often present. Methods: In this study, we aimed to delineate the structural and functional alterations of GRIN1 disease-associated variants, and their correlations with clinical symptoms in a Spanish cohort of 15 paediatric encephalopathy patients harbouring these variants. Results: Patients harbouring GRIN1 disease-associated variants have been clinically deeply-phenotyped. Further, using computational and in vitro approaches, we identified different critical checkpoints affecting GluN1 biogenesis (protein stability, subunit assembly and surface trafficking) and/or NMDAR biophysical properties, and their association with GRD clinical symptoms. Conclusions: Our findings show a strong correlation between GRIN1 variants-associated structural and functional outcomes. This structural-functional stratification provides relevant insights of genotype-phenotype association, contributing to future precision medicine of GRIN1-related encephalopathies.

Phosphorylation activates the yeast small heat shock protein Hsp26 by weakening domain contacts in the oligomer ensemble

Hsp26 is a small heat shock protein (sHsp) from S. cerevisiae. Its chaperone activity is activated by oligomer dissociation at heat shock temperatures. Hsp26 contains 9 phosphorylation sites in different structural elements. Our analysis of phospho-mimetic mutations shows that phosphorylation activates Hsp26 at permissive temperatures. The cryo-EM structure of the Hsp26 40mer revealed contacts between the conserved core domain of Hsp26 and the so-called thermosensor domain in the N-terminal part of the protein, which are targeted by phosphorylation. Furthermore, several phosphorylation sites in the C-terminal extension, which link subunits within the oligomer, are sensitive to the introduction of negative charges. In all cases, the intrinsic inhibition of chaperone activity is relieved and the N-terminal domain becomes accessible for substrate protein binding. The weakening of domain interactions within and between subunits by phosphorylation to activate the chaperone activity in response to proteotoxic stresses independent of heat stress could be a general regulation principle of sHsps.

CCA-addition in the cold: Structural characterization of the psychrophilic CCA-adding enzyme from the permafrost bacterium Planococcus halocryophilus

•

A high-resolution structure of a psychrophilic RNA polymerase contributes to our knowledge of cold adaptation.

•

While catalytic core motifs are conserved, at least one shows cold adaptation.

•

Loss of helix-capping increases structural flexibility in a catalytic core motif.

•

Overall reduction of alpha-helical elements appears as a strategy for cold adaptation.

CCA-adding enzymes are highly specific RNA polymerases that add and maintain the sequence C-C-A at tRNA 3‘-ends. Recently, we could reveal that cold adaptation of such a polymerase is not only achieved at the expense of enzyme stability, but also at the cost of polymerization fidelity. Enzymes from psychrophilic organisms usually show an increased structural flexibility to enable catalysis at low temperatures. Here, polymerases face a dilemma, as there is a discrepancy between the need for a tightly controlled flexibility during polymerization and an increased flexibility as strategy for cold adaptation. Based on structural and biochemical analyses, we contribute to clarify the cold adaptation strategy of the psychrophilic CCA-adding enzyme from Planococcus halocryophilus, a gram-positive bacterium thriving in the arctic permafrost at low temperatures down to −15 °C. A comparison with the closely related enzyme from the thermophilic bacterium Geobacillus stearothermophilus reveals several features of cold adaptation - a significantly reduced amount of alpha-helical elements in the C-terminal tRNA-binding region and a structural adaptation in one of the highly conserved catalytic core motifs located in the N-terminal catalytic core of the enzyme.

TopProperty: Robust Metaprediction of Transmembrane and Globular Protein Features Using Deep Neural Networks

Transmembrane proteins (TMPs) are critical components of cellular life. However, due to experimental challenges, the number of experimentally resolved TMP structures is severely underrepresented in databases compared to their cellular abundance. Prediction of (per-residue) features such as transmembrane topology, membrane exposure, secondary structure, and solvent accessibility can be a useful starting point for experimental design or protein structure prediction but often requires different computational tools for different features or types of proteins. We present TopProperty, a metapredictor that predicts all of these features for TMPs or globular proteins. TopProperty is trained on datasets without bias toward a high number of sequence homologs, and the predictions are significantly better than the evaluated state-of-the-art primary predictors on all quality metrics. TopProperty eliminates the need for protein type- or feature-tailored tools, specifically for TMPs. TopProperty is freely available as a web server and standalone at https://cpclab.uni-duesseldorf.de/topsuite/.

In silico evaluation of molecular virus-virus interactions taking place between Cotton leaf curl Kokhran virus- Burewala strain and Tomato leaf curl New Delhi virus

Background

Cotton leaf curl disease (CLCuD) is a disease of cotton caused by begomoviruses, leading to a drastic loss in the annual yield of the crop. Pakistan has suffered two epidemics of this disease leading to the loss of billions in annual exports. The speculation that a third epidemic of CLCuD may result as consequence of the frequent occurrence of Tomato leaf curl New Delhi virus (ToLCNDV) and Cotton leaf curl Kokhran Virus-Burewala Strain (CLCuKoV-Bu) in CLCuD infected samples, demand that the interactions taking between the two viruses be properly evaluated. This study is designed to assess virus-virus interactions at the molecular level and determine the type of co-infection taking place.

Methods

Based on the amino acid sequences of the gene products of both CLCuKoV-Bu and ToLCNDV, protein structures were generated using different software, i.e., MODELLER, I-TASSER, QUARKS, LOMETS and RAPTORX. A consensus model for each protein was selected after model quality assessment using ERRAT, QMEANDisCo, PROCHECK Z-Score and Ramachandran plot analysis. The active and passive residues in the protein structures were identified using the CPORT server. Protein–Protein Docking was done using the HADDOCK webserver, and 169 Protein–Protein Interaction (PPIs) were performed between the proteins of the two viruses. The docked complexes were submitted to the PRODIGY server to identify the interacting residues between the complexes. The strongest interactions were determined based on the HADDOCK Score, Desolvation energy, Van der Waals Energy, Restraint Violation Energy, Electrostatic Energy, Buried Surface Area and Restraint Violation Energy, Binding Affinity and Dissociation constant (K_d). A total of 50 ns Molecular Dynamic simulations were performed on complexes that exhibited the strongest affinity in order to validate the stability of the complexes, and to remove any steric hindrances that may exist within the structures.

Results

Our results indicate significant interactions taking place between the proteins of the two viruses. Out of all the interactions, the strongest were observed between the Replication Initiation protein (Rep) of CLCuKoV-Bu with the Movement protein (MP), Nuclear Shuttle Protein (NSP) of ToLCNDV (DNA-B), while the weakest were seen between the Replication Enhancer protein (REn) of CLCuKoV-Bu with the REn protein of ToLCNDV. The residues identified to be taking a part in interaction belonged to domains having a pivotal role in the viral life cycle and pathogenicity. It maybe deduced that the two viruses exhibit antagonistic behavior towards each other, and the type of infection may be categorised as a type of Super Infection Exclusion (SIE) or homologous interference. However, further experimentation, in the form of transient expression analysis, is needed to confirm the nature of these interactions and increase our understanding of the direct interactions taking place between two viruses.

Three-Component Synthesis of 2-Amino-3-cyano-4H-chromenes, In Silico Analysis of Their Pharmacological Profile, and In Vitro Anticancer and Antifungal Testing

Chromenes are compounds that may be useful for inhibiting topoisomerase and cytochrome, enzymes involved in the growth of cancer and fungal cells, respectively. The aim of this study was to synthesize a series of some novel 2-amino-3-cyano-4-aryl-6,7-methylendioxy-4H-chromenes 4a-o and 2-amino-3-cyano-5,7-dimethoxy-4-aryl-4H-chromenes 6a-h by a three-component reaction, and test these derivatives for anticancer and antifungal activity. Compounds 4a and 4b were more active than cisplatin (9) and topotecan (7) in SK-LU-1 cells, and more active than 9 in PC-3 cells. An evaluation was also made of the series of compounds 4 and 6 as potential antifungal agents against six Candida strains, finding their MIC₅₀ to be less than or equal to that of fluconazole (8). Molecular docking studies are herein reported, for the interaction of 4 and 6 with topoisomerase IB and the active site of CYP51 of Candida spp. Compounds 4a-o and 6a-h interacted in a similar way as 7 with key amino acids of the active site of topoisomerase IB and showed better binding energy than 8 at the active site of CYP51. Hence, 4a-o and 6a-h are good candidates for further research, having demonstrated their dual inhibition of enzymes that participate in the growth of cancer and fungal cells.

Pseudokinases repurpose flexibility signatures associated with the protein kinase fold for noncatalytic roles

The bilobal protein kinase-like fold in pseudokinases lack one or more catalytic residues, conserved in canonical protein kinases, and are considered enzymatically deficient. Tertiary structures of pseudokinases reveal that their loops topologically equivalent to activation segments of kinases adopt contracted configurations, which is typically extended in active conformation of kinases. Herein, anisotropic network model based normal mode analysis (NMA) was conducted on 51 active conformation structures of protein kinases and 26 crystal structures of pseudokinases. Our observations indicate that although backbone fluctuation profiles are similar for individual kinase-pseudokinase families, low intensity mean square fluctuations in pseudo-activation segment and other sub-structures impart rigidity to pseudokinases. Analyses of collective motions from functional modes reveal that pseudokinases, compared to active kinases, undergo distinct conformational transitions using the same structural fold. All-atom NMA of protein kinase-pseudokinase pairs from each family, sharing high amino acid sequence identities, yielded distinct community clusters, partitioned by residues exhibiting highly correlated fluctuations. It appears that atomic fluctuations from equivalent activation segments guide community membership and network topologies for respective kinase and pseudokinase. Our findings indicate that such adaptations in backbone and side-chain fluctuations render pseudokinases competent for catalysis-independent roles.

A missense mutation converts the Na+,K+-ATPase into an ion channel and causes therapy-resistant epilepsy

The ion pump Na⁺,K⁺-ATPase is a critical determinant of neuronal excitability; however, its role in the etiology of diseases of the central nervous system (CNS) is largely unknown. We describe here the molecular phenotype of a Trp931Arg mutation of the Na⁺,K⁺-ATPase catalytic α1 subunit in an infant diagnosed with therapy-resistant lethal epilepsy. In addition to the pathological CNS phenotype, we also detected renal wasting of Mg²⁺. We found that membrane expression of the mutant α1 protein was low, and ion pumping activity was lost. Arginine insertion into membrane proteins can generate water-filled pores in the plasma membrane, and our molecular dynamic (MD) simulations of the principle states of Na⁺,K⁺-ATPase transport demonstrated massive water inflow into mutant α1 and destabilization of the ion-binding sites. MD simulations also indicated that a water pathway was created between the mutant arginine residue and the cytoplasm, and analysis of oocytes expressing mutant α1 detected a nonspecific cation current. Finally, neurons expressing mutant α1 were observed to be depolarized compared with neurons expressing wild-type protein, compatible with a lowered threshold for epileptic seizures. The results imply that Na⁺,K⁺-ATPase should be considered a neuronal locus minoris resistentia in diseases associated with epilepsy and with loss of plasma membrane integrity.

PIP2-induced membrane binding of the vinculin tail competes with its other binding partners

Vinculin plays a key role during the first phase of focal adhesion formation and interacts with the plasma membrane through specific binding of its tail domain to the lipid phosphatidylinositol 4,5-bisphosphate (PIP₂). Our understanding of the PIP₂-vinculin interaction has been hampered by contradictory biochemical and structural data. Here, we used a multiscale molecular dynamics simulation approach, in which unbiased coarse-grained molecular dynamics were used to generate starting structures for subsequent microsecond-long all-atom simulations. This allowed us to map the interaction of the vinculin tail with PIP₂-enriched membranes in atomistic detail. In agreement with experimental data, we have shown that membrane binding is sterically incompatible with the intramolecular interaction between vinculin's head and tail domain. Our simulations further confirmed biochemical and structural results, which identified two positively charged surfaces, the basic collar and the basic ladder, as the main PIP₂ interaction sites. By introducing a valency-disaggregated binding network analysis, we were able to map the protein-lipid interactions in unprecedented detail. In contrast to the basic collar, in which PIP₂ is specifically recognized by an up to hexavalent binding pocket, the basic ladder forms a series of low-valency binding sites. Importantly, many of these PIP₂ binding residues are also involved in maintaining vinculin in a closed, autoinhibited conformation. These findings led us to propose a molecular mechanism for the coupling between vinculin activation and membrane binding. Finally, our refined binding site suggests an allosteric relationship between PIP₂ and F-actin binding that disfavors simultaneous interaction with both ligands, despite nonoverlapping binding sites.

Designing BH3-Mimetic Peptide Inhibitors for the Viral Bcl-2 Homologues A179L and BHRF1: Importance of Long-Range Electrostatic Interactions

Viruses have evolved strategies to prevent apoptosis of infected cells at early stages of infection. The viral proteins (vBcl-2s) from specific viral genes adopt a helical fold that is structurally similar to that of mammalian antiapoptotic Bcl-2 proteins and exhibit little sequence similarity. Hence, vBcl-2 homologues are attractive targets to prevent viral infection. However, very few studies have focused on developing inhibitors for vBcl-2 homologues. In this study, we have considered two vBcl-2 homologues, A179L from African swine fever virus and BHRF1 from Epstein-Barr virus. We generated two sets of 8000 randomized BH3-like sequences from eight wild-type proapoptotic BH3 peptides. During this process, the four conserved hydrophobic residues and an Asp residue were retained at their respective positions, and all other positions were substituted randomly without any bias. We constructed 8000 structures each for A179L and BHRF1 in complex with BH3-like sequences. Histograms of interaction energies calculated between the peptide and the protein resulted in negatively skewed distributions. The BH3-like peptides with high helical propensities selected from the negative tail of the respective interaction energy distributions exhibited more favorable interactions with A179L and BHRF1, and they are rich in basic residues. Molecular dynamics studies and electrostatic potential maps further revealed that both acidic and basic residues favorably interact with A179L, while only basic residues have the most favorable interactions with BHRF1. As in mammalian homologues, the role of long-range interactions and nonhotspot residues has to be taken into account while designing specific BH3-mimetic inhibitors for vBcl-2 homologues.

A resource of high-quality and versatile nanobodies for drug delivery

Therapeutic and diagnostic efficacies of small biomolecules and chemical compounds are hampered by suboptimal pharmacokinetics. Here, we developed a repertoire of robust and high-affinity antihuman serum albumin nanobodies (Nb_HSA) that can be readily fused to small biologics for half-life extension. We characterized the thermostability, binding kinetics, and cross-species reactivity of Nb_HSAs, mapped their epitopes, and structurally resolved a tetrameric HSA-Nb complex. We parallelly determined the half-lives of a cohort of selected Nb_HSAs in an HSA mouse model by quantitative proteomics. Compared to short-lived control nanobodies, the half-lives of Nb_HSAs were drastically prolonged by 771-fold. Nb_HSAs have distinct and diverse pharmacokinetics, positively correlating with their albumin binding affinities at the endosomal pH. We then generated stable and highly bioactive Nb_HSA-cytokine fusion constructs “Duraleukin” and demonstrated Duraleukin's high preclinical efficacy for cancer treatment in a melanoma model. This high-quality and versatile Nb toolkit will help tailor drug half-life to specific medical needs.

•

We provide a resource of high-affinity and versatile albumin nanobodies for drug delivery

•

We systematically map albumin nanobody epitopes by hybrid structural approaches

•

We parallelly measure the pharmacokinetics of nanobodies in a humanized mouse model

•

We develop nanobody-cytokine conjugates “Duraleukin” for cancer immunotherapy

Navigating the amino acid sequence space between functional proteins using a deep learning framework

MOTIVATION

Shedding light on the relationships between protein sequences and functions is a challenging task with many implications in protein evolution, diseases understanding, and protein design. The protein sequence space mapping to specific functions is however hard to comprehend due to its complexity. Generative models help to decipher complex systems thanks to their abilities to learn and recreate data specificity. Applied to proteins, they can capture the sequence patterns associated with functions and point out important relationships between sequence positions. By learning these dependencies between sequences and functions, they can ultimately be used to generate new sequences and navigate through uncharted area of molecular evolution.

RESULTS

This study presents an Adversarial Auto-Encoder (AAE) approached, an unsupervised generative model, to generate new protein sequences. AAEs are tested on three protein families known for their multiple functions the sulfatase, the HUP and the TPP families. Clustering results on the encoded sequences from the latent space computed by AAEs display high level of homogeneity regarding the protein sequence functions. The study also reports and analyzes for the first time two sampling strategies based on latent space interpolation and latent space arithmetic to generate intermediate protein sequences sharing sequential properties of original sequences linked to known functional properties issued from different families and functions. Generated sequences by interpolation between latent space data points demonstrate the ability of the AAE to generalize and produce meaningful biological sequences from an evolutionary uncharted area of the biological sequence space. Finally, 3D structure models computed by comparative modelling using generated sequences and templates of different sub-families point out to the ability of the latent space arithmetic to successfully transfer protein sequence properties linked to function between different sub-families. All in all this study confirms the ability of deep learning frameworks to model biological complexity and bring new tools to explore amino acid sequence and functional spaces.

Empirical Potential Energy Function Toward ab Initio Folding G Protein-Coupled Receptors

Approximately 40-50 percent of all drugs targets are G protein-coupled receptors (GPCRs). Three-dimensional structure of GPCRs is important to probe their biophysical and biochemical functions and their pharmaceutical applications. Lacking reliable and high quality free function is one of the ugent problems of computational predicting the three-dimensional structure in this community. We proposed a GPCR-specified energy function composed of four novel empirical potential energy terms: a two-dimensional contact energy force field, knowledge-based helix pair connection distance energy term, knowledge-based helix pair angle restraint energy term and a disulfide bond energy term. To validate the energy function, we employed an ab initio GPCR three-dimensional structure predictor to test if the energy function improved the accuracy of prediction. We evaluated 28 solved GPCRs and found that 21(75 percent) targets were correctly folded (TM-score>0.5). Also, the average TM-score using the energy function was 0.54, which was improved 134 percent than the TM-score 0.23 for MODELLER energy function and 170 percent than the TM-score 0.20 for Rosetta membrane energy function. The results confirmed that our empirical potential energy function toward ab initio folding is competitive to state-of-the-art solutions for structural prediction of GPCRs.

Characterization of the multidrug efflux transporter styMdtM from Salmonella enterica serovar Typhi

Salmonellae are foodborne pathogens and the major cause of gastroenteritis in humans. Salmonellae express multidrug efflux transporters that play a key role in their drug resistance, which is becoming an increasing problem for therapeutic intervention. Despite their biomedical importance, the mechanisms underlying substrate transport by multidrug efflux transporters remain poorly understood. Here, we describe the first characterization of a multidrug transporter belonging to the major facilitator superfamily from the genus Salmonella. We show that several clinical Salmonella Typhi (S. Typhi) isolates constitutively express the styMdtM (STY4874) gene, which encodes a known multidrug-resistance (MDR) transporter. Guided by the structure of the Escherichia coli (E. coli) homolog, we studied two residues critical for substrate transport, Asp25 and Arg111. Mutation of Asp25 to glutamate did not affect the transport function of styMdtM, whereas mutation to alanine reduced its transport activity, suggesting that a negative charge at this position is critical for substrate translocation across the membrane. Substrate-affinity measurements by intrinsic fluorescence spectroscopy showed that the Asp25Ala mutant retained its capacity to bind substrate, albeit at a lower level. Mutation of Arg111 to alanine resulted in a decrease in secondary structure content of the transporter, and mutation to lysine completely destabilized the structure of the transporter. A homology model of styMdtM suggests that Arg111 is important for stabilizing the transmembrane domain by mediating necessary interactions between neighboring helices. Together, our studies provide new structural and mechanistic insights into the Salmonella MDR transporter styMdtM.

Reconstruction of Three-Dimensional Conformations of Bacterial ClpB from High-Speed Atomic-Force-Microscopy Images

ClpB belongs to the cellular disaggretase machinery involved in rescuing misfolded or aggregated proteins during heat or other cellular shocks. The function of this protein relies on the interconversion between different conformations in its native condition. A recent high-speed-atomic-force-microscopy (HS-AFM) experiment on ClpB from Thermus thermophilus shows four predominant conformational classes, namely, open, closed, spiral, and half-spiral. Analyses of AFM images provide only partial structural information regarding the molecular surface, and thus computational modeling of three-dimensional (3D) structures of these conformations should help interpret dynamical events related to ClpB functions. In this study, we reconstruct 3D models of ClpB from HS-AFM images in different conformational classes. We have applied our recently developed computational method based on a low-resolution representation of 3D structure using a Gaussian mixture model, combined with a Monte-Carlo sampling algorithm to optimize the agreement with target AFM images. After conformational sampling, we obtained models that reflect conformational variety embedded within the AFM images. From these reconstructed 3D models, we described, in terms of relative domain arrangement, the different types of ClpB oligomeric conformations observed by HS-AFM experiments. In particular, we highlighted the slippage of the monomeric components around the seam. This study demonstrates that such details of information, necessary for annotating the different conformational states involved in the ClpB function, can be obtained by combining HS-AFM images, even with limited resolution, and computational modeling.

A KLK4 proteinase substrate capture approach to antagonize PAR1

Proteinase-activated receptor-1 (PAR1), triggered by thrombin and other serine proteinases such as tissue kallikrein-4 (KLK4), is a key driver of inflammation, tumor invasiveness and tumor metastasis. The PAR1 transmembrane G-protein-coupled receptor therefore represents an attractive target for therapeutic inhibitors. We thus used a computational design to develop a new PAR1 antagonist, namely, a catalytically inactive human KLK4 that acts as a proteinase substrate-capture reagent, preventing receptor cleavage (and hence activation) by binding to and occluding the extracellular R41-S42 canonical PAR1 proteolytic activation site. On the basis of in silico site-saturation mutagenesis, we then generated KLK4S207A,L185D, a first-of-a-kind 'decoy' PAR1 inhibitor, by mutating the S207A and L185D residues in wild-type KLK4, which strongly binds to PAR1. KLK4S207A,L185D markedly inhibited PAR1 cleavage, and PAR1-mediated MAPK/ERK activation as well as the migration and invasiveness of melanoma cells. This 'substrate-capturing' KLK4 variant, engineered to bind to PAR1, illustrates proof of principle for the utility of a KLK4 'proteinase substrate capture' approach to regulate proteinase-mediated PAR1 signaling.

Using yeast two-hybrid system and molecular dynamics simulation to detect venom protein-protein interactions

•

The venom protein-protein interactions in snake venom remain largely unknown.

•

Y2H coupled with MD simulations was used to detect venom protein interactions.

•

Venom PLA₂s interact with themselves and Lys49 PLA₂ interacts with venom CRISP.

Proteins and peptides are major components of snake venom. Venom protein transcriptomes and proteomes of many snake species have been reported; however, snake venom complexity (i.e., the venom protein-protein interactions, PPIs) remains largely unknown. To detect the venom protein interactions, we used the most common snake venom component, phospholipase A₂s (PLA₂s) as a “bait” to identify the interactions between PLA₂s and 14 of the most common proteins in Western diamondback rattlesnake (Crotalus atrox) venom by using yeast two-hybrid (Y2H) analysis, a technique used to detect PPIs. As a result, we identified PLA₂s interacting with themselves, and lysing-49 PLA₂ (Lys49 PLA₂) interacting with venom cysteine-rich secretory protein (CRISP). To reveal the complex structure of Lys49 PLA₂-CRISP interaction at the structural level, we first built the three-dimensional (3D) structures of Lys49 PLA₂ and CRISP by a widely used computational program-MODELLER. The binding modes of Lys49 PLA₂-CRISP interaction were then predicted through three different docking programs including ClusPro, ZDOCK and HADDOCK. Furthermore, the most likely complex structure of Lys49 PLA₂-CRISP was inferred by molecular dynamic (MD) simulations with GROMACS software. The techniques used and results obtained from this study strengthen the understanding of snake venom protein interactions and pave the way for the study of animal venom complexity.

Exome sequencing identifies a disease variant of the mitochondrial ATP-Mg/Pi carrier SLC25A25 in two families with kidney stones

BACKGROUND

Calcium kidney stones are common and recurrences are often not preventable by available empiric remedies. Their etiology is multifactorial and polygenic, and an increasing number of genes are implicated. Their identification will enable improved management.

METHODS

DNA from three stone-formers in a Southampton family (UK) and two from an Italian family were analyzed independently by whole exome sequencing and selected variants were genotyped across all available members of both pedigrees. A disease variant of SLC25A25 (OMIM 608745), encoding the mitochondrial ATP-Mg/Pi carrier 3 (APC3) was identified, and analyzed structurally and functionally with respect to its calcium-regulated transport activity.

RESULTS

All five patients had a heterozygous dominant SLC25A25 variant (rs140777921; GRCh37.p13: chr 9 130868670 G>C; p.Gln349His; Reference Sequence NM_001006641.3). Non-stone formers also carried the variant indicating incomplete penetrance. Modeling suggests that the variant lacks a conserved polar interaction, which may cause structural instability. Calcium-regulated ATP transport was reduced to ~20% of the wild type, showing a large reduction in function.

CONCLUSION

The transporter is important in regulating mitochondrial ATP production. This rare variant may increase urine lithogenicity through impaired provision of ATP for solute transport processes in the kidney, and/or for purinergic signaling. Variants found in other genes may compound this abnormality.

Capabilities of bioinformatics tools for optimizing physicochemical features of proteins used in Nano biosensors: A short overview of the tools related to bioinformatics

Protein-protein ligand is one of the most detection methods used in Nano biosensors. Based on the advantage of specific docking between two special 3D structures, they have become a potent candidate in bioanalysis and Nanodiagnostic tools. These tools lease users to do a simple, fast, cost-effective, sensitive, and specific detection of molecular biomarkers in real samples. Recent advantages of using protein-protein ligand Nano-biosensors application is remarkable due to its special docking that refers to each protein unique 3D conformation. However, it challenges different problems such as low rate of docking and hard process for fixation on the basic layer. These challenges make developers to optimize the structure and functions of proteins. The process has different Nano scale calculation that could be done with algorithms and solutions are available as bioinformatics tools. This article aimed to have a short overview of the abilities of bioinformatics tools for modeling and optimization of physiochemical features of proteins in Nano scale.

•

Nano biosensors use different strategies which based on docking between two molecules to detect and identify different proteins.

•

Molecular docking between transducer in Nano biosensors and proteins rely on physicochemical features of transducer, protein and docking strategy.

•

Nano bioinformatics use bioinformatics tools and algorithms as a collective solution for developing functional structure in Nano scale.

•

Nano bioinformatics use different bioinformatics tools to optimize physicochemical features of proteins as a new approach in Nano biosensors and drug discovery.

Understanding the Folding Mediated Assembly of the Bacteriophage MS2 Coat Protein Dimers

The capsids of RNA viruses such as MS2 are great models for studying protein self-assembly because they are made almost entirely of multiple copies of a single coat protein (CP). Although CP is the minimal repeating unit of the capsid, previous studies have shown that CP exists as a homodimer (CP2) even in an acid-disassembled system, indicating that CP2 is an obligate dimer. Here, we investigate the molecular basis of this obligate dimerization using coarse-grained structure-based models and molecular dynamics simulations. We find that, unlike monomeric proteins of similar size, CP populates a single partially folded ensemble whose "foldedness" is sensitive to denaturing conditions. In contrast, CP2 folds similarly to single-domain proteins populating only the folded and the unfolded ensembles, separated by a prominent folding free energy barrier. Several intramonomer contacts form early, but the CP2 folding barrier is crossed only when the intermonomer contacts are made. A dissection of the structure of CP2 through mutant folding simulations shows that the folding barrier arises both from the topology of CP and the interface contacts of CP2. Together, our results show that CP2 is an obligate dimer because of kinetic stability, that is, dimerization induces a folding barrier and that makes it difficult for proteins in the dimer minimum to partially unfold and access the monomeric state without completely unfolding. We discuss the advantages of this obligate dimerization in the context of dimer design and virus stability.

α-tubulin tail modifications regulate microtubule stability through selective effector recruitment, not changes in intrinsic polymer dynamics

Microtubules are non-covalent polymers of αβ-tubulin dimers. Posttranslational processing of the intrinsically disordered C-terminal α-tubulin tail produces detyrosinated and Δ2-tubulin. Although these are widely employed as proxies for stable cellular microtubules, their effect (and of the α-tail) on microtubule dynamics remains uncharacterized. Using recombinant, engineered human tubulins, we now find that neither detyrosinated nor Δ2-tubulin affect microtubule dynamics, while the α-tubulin tail is an inhibitor of microtubule growth. Consistent with the latter, molecular dynamics simulations show the α-tubulin tail transiently occluding the longitudinal microtubule polymerization interface. The marked differential in vivo stabilities of the modified microtubule subpopulations, therefore, must result exclusively from selective effector recruitment. We find that tyrosination quantitatively tunes CLIP-170 density at the growing plus end and that CLIP170 and EB1 synergize to selectively upregulate the dynamicity of tyrosinated microtubules. Modification-dependent recruitment of regulators thereby results in microtubule subpopulations with distinct dynamics, a tenet of the tubulin code hypothesis.

Polycyclic aromatic hydrocarbons modulate the activity of Atlantic cod (Gadus morhua) vitamin D receptor paralogs in vitro

Vitamin D receptor (VDR) mediates the biological function of the steroid hormone calcitriol, which is the metabolically active version of vitamin D. Calcitriol is important for a wide array of physiological functions, including calcium and phosphate homeostasis. In contrast to mammals, which harbor one VDR encoding gene, teleosts possess two orthologous vdr genes encoding Vdr alpha (Vdra) and Vdr beta (Vdrb). Genome mining identified the vdra and vdrb paralogs in the Atlantic cod (Gadus morhua) genome, which were further characterized regarding their phylogeny, tissue-specific expression, and transactivational properties induced by calcitriol. In addition, a selected set of polycyclic aromatic hydrocarbons (PAHs), including naphthalene, phenanthrene, fluorene, pyrene, chrysene, benzo[a]pyrene (BaP), and 7-methylbenzo[a]pyrene, were assessed for their ability to modulate the transcriptional activity of gmVdra and gmVdrb in vitro. Both gmVdra and gmVdrb were activated by calcitriol with similar potencies, but gmVdra produced significantly higher maximal fold activation. Notably, none of the tested PAHs showed agonistic properties towards the Atlantic cod Vdrs. However, binary exposures of calcitriol together with phenanthrene, fluorene, or pyrene, antagonized the activation of gmVdra, while chrysene and BaP significantly potentiated the calcitriol-mediated activity of both receptors. Homology modeling, solvent mapping, and docking analyses complemented the experimental data, and revealed a putative secondary binding site in addition to the canonical ligand-binding pocket (LBP). Calcitriol was predicted to interact with both binding sites, whereas PAHs docked primarily to the LBP. Importantly, our in vitro data suggest that PAHs can interact with the paralogous gmVdrs and interfere with their transcriptional activities, and thus potentially modulate the vitamin D signaling pathway and contribute to adverse effects of crude oil and PAH exposures on cardiac development and bone deformities in fish.

Elucidating biophysical basis of binding of inhibitors to SARS-CoV-2 main protease by using molecular dynamics simulations and free energy calculations

The recent outbreak of novel "coronavirus disease 2019" (COVID-19) has spread rapidly worldwide, causing a global pandemic. In the present work, we have elucidated the mechanism of binding of two inhibitors, namely α-ketoamide and Z31792168, to SARS-CoV-2 main protease (Mpro or 3CLpro) by using all-atom molecular dynamics simulations and free energy calculations. We calculated the total binding free energy (ΔGbind) of both inhibitors and further decomposed ΔGbind into various forces governing the complex formation using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method. Our calculations reveal that α-ketoamide is more potent (ΔGbind= - 9.05 kcal/mol) compared to Z31792168 (ΔGbind= - 3.25 kcal/mol) against COVID-19 3CLpro. The increase in ΔGbind for α-ketoamide relative to Z31792168 arises due to an increase in the favorable electrostatic and van der Waals interactions between the inhibitor and 3CLpro. Further, we have identified important residues controlling the 3CLpro-ligand binding from per-residue based decomposition of the binding free energy. Finally, we have compared ΔGbind of these two inhibitors with the anti-HIV retroviral drugs, such as lopinavir and darunavir. It is observed that α-ketoamide is more potent compared to lopinavir and darunavir. In the case of lopinavir, a decrease in van der Waals interactions is responsible for the lower binding affinity compared to α-ketoamide. On the other hand, in the case of darunavir, a decrease in the favorable intermolecular electrostatic and van der Waals interactions contributes to lower affinity compared to α-ketoamide. Our study might help in designing rational anti-coronaviral drugs targeting the SARS-CoV-2 main protease.Communicated by Ramaswamy H. Sarma.

Structure of the mycobacterial ESX-5 type VII secretion system pore complex

The ESX-5 type VII secretion system is a membrane-spanning protein complex key to the virulence of mycobacterial pathogens. However, the overall architecture of the fully assembled translocation machinery and the composition of the central secretion pore have remained unknown. Here, we present the high-resolution structure of the 2.1-megadalton ESX-5 core complex. Our structure captured a dynamic, secretion-competent conformation of the pore within a well-defined transmembrane section, sandwiched between two flexible protein layers at the cytosolic entrance and the periplasmic exit. We propose that this flexibility endows the ESX-5 machinery with large conformational plasticity required to accommodate targeted protein secretion. Compared to known secretion systems, a highly dynamic state of the pore may represent a fundamental principle of bacterial secretion machineries.

Novel antimicrobial anionic cecropins from the spruce budworm feature a poly-L-aspartic acid C-terminus

Cecropins form a family of amphipathic α-helical cationic peptides with broad-spectrum antibacterial properties and potent anticancer activity. The emergence of bacteria and cancer cells showing resistance to cationic antimicrobial peptides (CAMPs) has fostered a search for new, more selective and more effective alternatives to CAMPs. With this goal in mind, we looked for cecropin homologs in the genome and transcriptome of the spruce budworm, Choristoneura fumiferana. Not only did we find paralogs of the conventional cationic cecropins (Cfcec⁺ ), our screening also led to the identification of previously uncharacterized anionic cecropins (Cfcec^- ), featuring a poly-l-aspartic acid C-terminus. Comparative peptide analysis indicated that the C-terminal helix of Cfcec^- is amphipathic, unlike that of Cfcec⁺ , which is hydrophobic. Interestingly, molecular dynamics simulations pointed to the lower conformational flexibility of Cfcec^- peptides, relative to that of Cfcec⁺ . Phylogenetic analysis suggests that the evolution of distinct Cfcec⁺ and Cfcec^- peptides may have resulted from an ancient duplication event within the Lepidoptera. Finally, we found that both anionic and cationic cecropins contain a BH3-like motif (G-[KQR]-[HKQNR]-[IV]-[KQR]) that could interact with Bcl-2, a protein involved in apoptosis; this observation is congruent with previous reports indicating that cecropins induce apoptosis. Altogether, our observations suggest that cecropins may provide templates for the development of new anticancer drugs. We also estimated the antibacterial activity of Cfcec-2 and a ∆Cfce-2 peptide as AMPs by testing directly their ability in inhibiting bacterial growth in a disk diffusion assay and their potential for development of novel therapeutics.

Characterization of recombinant E. coli expressing a novel fucosidase from Bacillus cereus 2-8 belonging to GH95 family

Fucoidan oligosaccharides possesses diverse physicochemical and biological activities. Specific glycoside hydrolases are valuable tools for degrading polysaccharides to produce oligosaccharides. In this study, BcFucA, a novel fucosidase belonging to GH95 family from Bacillus cereus 2-8, was cloned into pET21a vector, expressed in E. coli BL21 (DE3) and characterized. The protein consists of 1136 amino acid residues encoded by 3411 bases and has a molecular weight of 125.35 kDa. The optimal temperature and pH of this enzyme are 50 °C and pH 4.0. In addition, this study showed that the unknown function domain (encoding Lys261-Thr681) defined as a linker is quite important for its activity. The obtained novel enzyme BcFucA will contribute to the effective degradation of fucoidan and future industrial applications.

Structure elements can be predicted using the contact volume among protein residues

Previously, the structure elements of dihydrofolate reductase (DHFR) were determined using comprehen­sive Ala-insertion mutation analysis, which is assumed to be a kind of protein “building blocks.” It is hypo­thesized that our comprehension of the structure elements could lead to understanding how an amino acid sequence dictates its tertiary structure. However, the comprehensive Ala-insertion mutation analysis is a time- and cost-consuming process and only a set of the DHFR structure elements have been reported so far. Therefore, developing a computational method to predict structure elements is an urgent necessity. We focused on intramolecular residue–residue contacts to predict the structure elements. We introduced a simple and effective parameter: the overlapped contact volume (CV) among the residues and calculated the CV along the DHFR sequence using the crystal structure. Our results indicate that the CV profile can recapitulate its precipitate ratio profile, which was used to define the structure elements in the Ala-insertion mutation analysis. The CV profile allowed us to predict structure elements like the experimentally determined structure elements. The strong correlation between the CV and precipitate ratio profiles indicates the importance of the intramolecular residue–residue contact in maintaining the tertiary structure. Additionally, the CVs between the structure elements are considerably more than those between a structure element and a linker or two linkers, indicating that the structure elements play a funda­mental role in increasing the intramolecular adhesion. Thus, we propose that the structure elements can be considered a type of “building blocks” that maintain and dictate the tertiary structures of proteins.

Structural insights into integrin α5β1 opening by fibronectin ligand

Integrin α₅β₁ is a major fibronectin receptor critical for cell migration. Upon complex formation, fibronectin and α₅β₁ undergo conformational changes. While this is key for cell-tissue connections, its mechanism is unknown. Here, we report cryo-electron microscopy structures of native human α₅β₁ with fibronectin to 3.1-angstrom resolution, and in its resting state to 4.6-angstrom resolution. The α₅β₁-fibronectin complex revealed simultaneous interactions at the arginine-glycine-aspartate loop, the synergy site, and a newly identified binding site proximal to adjacent to metal ion-dependent adhesion site, inducing the translocation of helix α1 to secure integrin opening. Resting α₅β₁ adopts an incompletely bent conformation, challenging the model of integrin sharp bending inhibiting ligand binding. Our biochemical and structural analyses showed that affinity of α₅β₁ for fibronectin is increased with manganese ions (Mn²⁺) while adopting the half-bent conformation, indicating that ligand-binding affinity does not depend on conformation, and α₅β₁ opening is induced by ligand-binding.

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.

PCPD: Plant cytochrome P450 database and web-based tools for structural construction and ligand docking

Plant cytochrome P450s play key roles in the diversification and functional modification of plant natural products. Although over 200,000 plant P450 gene sequences have been recorded, only seven crystalized P450 genes severely hampered the functional characterization, gene mining and engineering of important P450s. Here, we combined Rosetta homologous modeling and MD-based refinement to construct a high-resolution P450 structure prediction process (PCPCM), which was applied to 181 plant P450s with identified functions. Furthermore, we constructed a ligand docking process (PCPLD) that can be applied for plant P450s virtual screening. 10 examples of virtual screening indicated the process can reduce about 80% screening space for next experimental verification. Finally, we constructed a plant P450 database (PCPD: http://p450.biodesign.ac.cn/), which includes the sequences, structures and functions of the 181 plant P450s, and a web service based on PCPCM and PCPLD. Our study not only developed methods for the P450-specific structure analysis, but also introduced a universal approach that can assist the mining and functional analysis of P450 enzymes.