Fighting Celiac Disease: Improvement of pH Stability of Cathepsin L In Vitro by Computational Design.

Roughly 1% of the global population is susceptible to celiac disease (CD)-inheritable autoimmune inflammation of the small intestine caused by intolerance to gliadin proteins present in wheat, rye, and barley grains, and called gluten in wheat. Classical treatment is a life-long gluten-free diet, which is constraining and costly. An alternative approach is based upon the development and oral reception of effective peptidases that degrade in the stomach immunogenic proline- and glutamine-rich gliadin peptides, which are the cause of the severe reaction in the intestine. In previous research, we have established that the major digestive peptidase of an insect Tribolium castaneum-cathepsin L-hydrolyzes immunogenic prolamins after Gln residues but is unstable in the extremely acidic environment (pH 2-4) of the human stomach and cannot be used as a digestive aid. In this work, using molecular dynamics simulations, we discover the probable cause of the pH instability of cathepsin L-loss of the catalytically competent rotameric state of one of the active site residues, His 275. To "fix" the correct orientation of this residue, we designed a V277A mutant variant, which extends the range of stability of the peptidase in the acidic environment while retaining most of its activity. We suggest this protein as a lead glutenase for the development of oral medical preparation that fights CD and gluten intolerance in susceptible people.

Specific Binding of the α-Component of the Lantibiotic Lichenicidin to the Peptidoglycan Precursor Lipid II Predetermines Its Antimicrobial Activity.

To date, a number of lantibiotics have been shown to use lipid II-a highly conserved peptidoglycan precursor in the cytoplasmic membrane of bacteria-as their molecular target. The α-component (Lchα) of the two-component lantibiotic lichenicidin, previously isolated from the Bacillus licheniformis VK21 strain, seems to contain two putative lipid II binding sites in its N-terminal and C-terminal domains. Using NMR spectroscopy in DPC micelles, we obtained convincing evidence that the C-terminal mersacidin-like site is involved in the interaction with lipid II. These data were confirmed by the MD simulations. The contact area of lipid II includes pyrophosphate and disaccharide residues along with the first isoprene units of bactoprenol. MD also showed the potential for the formation of a stable N-terminal nisin-like complex; however, the conditions necessary for its implementation in vitro remain unknown. Overall, our results clarify the picture of two component lantibiotics mechanism of antimicrobial action.

Orientational Preferences of GPI-Anchored Ly6/uPAR Proteins.

Ly6/uPAR proteins regulate many essential functions in the nervous and immune systems and epithelium. Most of these proteins contain single ß-structural LU domains with three protruding loops and are glycosylphosphatidylinositol (GPI)-anchored to a membrane. The GPI-anchor role is currently poorly studied. Here, we investigated the positional and orientational preferences of six GPI-anchored proteins in the receptor-unbound state by molecular dynamics simulations. Regardless of the linker length between the LU domain and GPI-anchor, the proteins interacted with the membrane by polypeptide parts and N-/O-glycans. Lynx1, Lynx2, Lypd6B, and Ly6H contacted the membrane by the loop regions responsible for interactions with nicotinic acetylcholine receptors, while Lypd6 and CD59 demonstrated unique orientations with accessible receptor-binding sites. Thus, GPI-anchoring does not guarantee an optimal 'pre-orientation' of the LU domain for the receptor interaction.

The Mechanism of Selective Recognition of Lipid Substrate by hDHHC20 Enzyme.

S-acylation is a post-translational linkage of long chain fatty acids to cysteines, playing a key role in normal physiology and disease. In human cells, the reaction is catalyzed by a family of 23 membrane DHHC-acyltransferases (carrying an Asp-His-His-Cys catalytic motif) in two stages: (1) acyl-CoA-mediated autoacylation of the enzyme; and (2) further transfer of the acyl chain to a protein substrate. Despite the availability of a 3D-structure of human acyltransferase (hDHHC20), the molecular aspects of lipid selectivity of DHHC-acyltransferases remain unclear. In this paper, using molecular dynamics (MD) simulations, we studied membrane-bound hDHHC20 right before the acylation by C12-, C14-, C16-, C18-, and C20-CoA substrates. We found that: (1) regardless of the chain length, its terminal methyl group always reaches the "ceiling" of the enzyme's cavity; (2) only for C16, an optimal "reactivity" (assessed by a simple geometric criterion) permits the autoacylation; (3) in MD, some key interactions between an acyl-CoA and a protein differ from those in the reference crystal structure of the C16-CoA-hDHHS20 mutant complex (probably, because this structure corresponds to a non-native dimer). These features of specific recognition of full-size acyl-CoA substrates support our previous hypothesis of "geometric and physicochemical selectivity" derived for simplified acyl-CoA analogues.

Molecular Dynamics of DHHC20 Acyltransferase Suggests Principles of Lipid and Protein Substrate Selectivity.

Lipid modification of viral proteins with fatty acids of different lengths (S-acylation) is crucial for virus pathogenesis. The reaction is catalyzed by members of the DHHC family and proceeds in two steps: the autoacylation is followed by the acyl chain transfer onto protein substrates. The crystal structure of human DHHC20 (hDHHC20), an enzyme involved in the acylation of S-protein of SARS-CoV-2, revealed that the acyl chain may be inserted into a hydrophobic cavity formed by four transmembrane (TM) α-helices. To test this model, we used molecular dynamics of membrane-embedded hDHHC20 and its mutants either in the absence or presence of various acyl-CoAs. We found that among a range of acyl chain lengths probed only C16 adopts a conformation suitable for hDHHC20 autoacylation. This specificity is altered if the small or bulky residues at the cavity's ceiling are exchanged, e.g., the V185G mutant obtains strong preferences for binding C18. Surprisingly, an unusual hydrophilic ridge was found in TM helix 4 of hDHHC20, and the responsive hydrophilic patch supposedly involved in association was found in the 3D model of the S-protein TM-domain trimer. Finally, the exchange of critical Thr and Ser residues in the spike led to a significant decrease in its S-acylation. Our data allow further development of peptide/lipid-based inhibitors of hDHHC20 that might impede replication of Corona- and other enveloped viruses.

Potassium channel blocker crafted by α-hairpinin scaffold engineering.

The α-Hairpinins are a family of plant defense peptides with a common fold presenting two short α-helices stabilized by two invariant S-S-bridges. We have shown previously that substitution of just two amino acid residues in a wheat α-hairpinin Tk-AMP-X2 leads to Tk-hefu-2 that features specific affinity to voltage-gated potassium channels KV1.3. Here, we utilize a combined molecular modeling approach based on molecular dynamics simulations and protein surface topography technique to improve the affinity of Tk-hefu-2 to KV1.3 while preserving its specificity. An important advance of this work compared with our previous studies is transition from the analysis of various physicochemical properties of an isolated toxin molecule to its consideration in complex with its target, a membrane-bound ion channel. As a result, a panel of computationally designed Tk-hefu-2 derivatives was synthesized and tested against KV1.3. The most active mutant Tk-hefu-10 showed a half-maximal inhibitory concentration of ∼150 nM being >10 times more active than Tk-hefu-2 and >200 times more active than the original Tk-hefu. We conclude that α-hairpinins provide an attractive disulfide-stabilized scaffold for the rational design of ion channel inhibitors. Furthermore, the success rate can be considerably increased by the proposed "target-based" iterative strategy of molecular design.

Conserved Structure and Evolution of DPF Domain of PHF10-The Specific Subunit of PBAF Chromatin Remodeling Complex.

Transcription activation factors and multisubunit coactivator complexes get recruited at specific chromatin sites via protein domains that recognize histone modifications. Single PHDs (plant homeodomains) interact with differentially modified H3 histone tails. Double PHD finger (DPF) domains possess a unique structure different from PHD and are found in six proteins: histone acetyltransferases MOZ and MORF; chromatin remodeling complex BAF (DPF1-3); and chromatin remodeling complex PBAF (PHF10). Among them, PHF10 stands out due to the DPF sequence, structure, and functions. PHF10 is ubiquitously expressed in developing and adult organisms as four isoforms differing in structure (the presence or absence of DPF) and transcription regulation functions. Despite the importance of the DPF domain of PHF10 for transcription activation, its structure remains undetermined. We performed homology modeling of the human PHF10 DPF domain and determined common and distinct features in structure and histone modifications recognition capabilities, which can affect PBAF complex chromatin recruitment. We also traced the evolution of DPF1-3 and PHF10 genes from unicellular to vertebrate organisms. The data reviewed suggest that the DPF domain of PHF10 plays an important role in SWI/SNF-dependent chromatin remodeling during transcription activation.

Protein surface topography as a tool to enhance the selective activity of a potassium channel blocker.

Tk-hefu is an artificial peptide designed based on the α-hairpinin scaffold, which selectively blocks voltage-gated potassium channels Kv1.3. Here we present its spatial structure resolved by NMR spectroscopy and analyze its interaction with channels using computer modeling. We apply protein surface topography to suggest mutations and increase Tk-hefu affinity to the Kv1.3 channel isoform. We redesign the functional surface of Tk-hefu to better match the respective surface of the channel pore vestibule. The resulting peptide Tk-hefu-2 retains Kv1.3 selectivity and displays ∼15 times greater activity compared with Tk-hefu. We verify the mode of Tk-hefu-2 binding to the channel outer vestibule experimentally by site-directed mutagenesis. We argue that scaffold engineering aided by protein surface topography represents a reliable tool for design and optimization of specific ion channel ligands.

C-Terminal residues in small potassium channel blockers OdK1 and OSK3 from scorpion venom fine-tune the selectivity.

We report isolation, sequencing, and electrophysiological characterization of OSK3 (α-KTx 8.8 in Kalium and Uniprot databases), a potassium channel blocker from the scorpion Orthochirus scrobiculosus venom. Using the voltage clamp technique, OSK3 was tested on a wide panel of 11 voltage-gated potassium channels expressed in Xenopus oocytes, and was found to potently inhibit Kv1.2 and Kv1.3 with IC50 values of ~331nM and ~503nM, respectively. OdK1 produced by the scorpion Odontobuthus doriae differs by just two C-terminal residues from OSK3, but shows marked preference to Kv1.2. Based on the charybdotoxin-potassium channel complex crystal structure, a model was built to explain the role of the variable residues in OdK1 and OSK3 selectivity.

Structural Insight into Specificity of Interactions between Nonconventional Three-finger Weak Toxin from Naja kaouthia (WTX) and Muscarinic Acetylcholine Receptors.

Weak toxin from Naja kaouthia (WTX) belongs to the group of nonconventional "three-finger" snake neurotoxins. It irreversibly inhibits nicotinic acetylcholine receptors and allosterically interacts with muscarinic acetylcholine receptors (mAChRs). Using site-directed mutagenesis, NMR spectroscopy, and computer modeling, we investigated the recombinant mutant WTX analogue (rWTX) which, compared with the native toxin, has an additional N-terminal methionine residue. In comparison with the wild-type toxin, rWTX demonstrated an altered pharmacological profile, decreased binding of orthosteric antagonist N-methylscopolamine to human M1- and M2-mAChRs, and increased antagonist binding to M3-mAChR. Positively charged arginine residues located in the flexible loop II were found to be crucial for rWTX interactions with all types of mAChR. Computer modeling suggested that the rWTX loop II protrudes to the M1-mAChR allosteric ligand-binding site blocking the entrance to the orthosteric site. In contrast, toxin interacts with M3-mAChR by loop II without penetration into the allosteric site. Data obtained provide new structural insight into the target-specific allosteric regulation of mAChRs by "three-finger" snake neurotoxins.

PREDDIMER: a web server for prediction of transmembrane helical dimers.

SUMMARY: Here we present PREDDIMER, a web tool for prediction of dimer structure of transmembrane (TM) helices. PREDDIMER allows (i) reconstruction of a number of dimer structures for given sequence(s) of TM protein fragments, (ii) ranking and filtering of predicted structures according to respective values of a scoring function, (iii) visualization of predicted 3D dimer structures and (iv) visualization of surface hydrophobicity of TM helices and their contacting (interface) regions represented as 2D maps. RESULTS: We implemented online the original PREDDIMER algorithm and benchmarked the server on 11 TM sequences, whose 3D dimer conformations were obtained previously by nuclear magnetic resonance spectroscopy. In the most of tested cases backbone root-mean-square deviations of closest predicted conformations from the experimental reference are below 3 Å. A randomization test displays good anticorrelation (-0.82) between values of the scoring function and statistical significance of the prediction 'by chance'. Going beyond a single dimer conformation, our web tool predicts an ensemble of possible conformations, which may be useful for explanation of a functioning of bitopic membrane proteins, e.g. receptor tyrosine kinases. AVAILABILITY AND IMPLEMENTATION: PREDDIMER can be accessed for free on the web at http://model.nmr.ru/preddimer/ CONTACT: newant@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Modular organization of α-toxins from scorpion venom mirrors domain structure of their targets, sodium channels.

To gain success in the evolutionary "arms race," venomous animals such as scorpions produce diverse neurotoxins selected to hit targets in the nervous system of prey. Scorpion α-toxins affect insect and/or mammalian voltage-gated sodium channels (Na(v)s) and thereby modify the excitability of muscle and nerve cells. Although more than 100 α-toxins are known and a number of them have been studied into detail, the molecular mechanism of their interaction with Na(v)s is still poorly understood. Here, we employ extensive molecular dynamics simulations and spatial mapping of hydrophobic/hydrophilic properties distributed over the molecular surface of α-toxins. It is revealed that despite the small size and relatively rigid structure, these toxins possess modular organization from structural, functional, and evolutionary perspectives. The more conserved and rigid "core module" is supplemented with the "specificity module" (SM) that is comparatively flexible and variable and determines the taxon (mammal versus insect) specificity of α-toxin activity. We further show that SMs in mammal toxins are more flexible and hydrophilic than in insect toxins. Concomitant sequence-based analysis of the extracellular loops of Na(v)s suggests that α-toxins recognize the channels using both modules. We propose that the core module binds to the voltage-sensing domain IV, whereas the more versatile SM interacts with the pore domain in repeat I of Na(v)s. These findings corroborate and expand the hypothesis on different functional epitopes of toxins that has been reported previously. In effect, we propose that the modular structure in toxins evolved to match the domain architecture of Na(v)s.

Deciphering fine molecular details of proteins' structure and function with a Protein Surface Topography (PST) method.

Molecular surfaces are the key players in biomolecular recognition and interactions. Nowadays, it is trivial to visualize a molecular surface and surface-distributed properties in three-dimensional space. However, such a representation trends to be biased and ambiguous in case of thorough analysis. We present a new method to create 2D spherical projection maps of entire protein surfaces and manipulate with them--protein surface topography (PST). It permits visualization and thoughtful analysis of surface properties. PST helps to easily portray conformational transitions, analyze proteins' properties and their dynamic behavior, improve docking performance, and reveal common patterns and dissimilarities in molecular surfaces of related bioactive peptides. This paper describes basic usage of PST with an example of small G-proteins conformational transitions, mapping of caspase-1 intersubunit interface, and intrinsic "complementarity" in the conotoxin-acetylcholine binding protein complex. We suggest that PST is a beneficial approach for structure-function studies of bioactive peptides and small proteins.

Structure-based prediction of T cell receptor recognition of unseen epitopes using TCRen.

T cell receptor (TCR) recognition of foreign peptides presented by major histocompatibility complex protein is a major event in triggering the adaptive immune response to pathogens or cancer. The prediction of TCR-peptide interactions has great importance for therapy of cancer as well as infectious and autoimmune diseases but remains a major challenge, particularly for novel (unseen) peptide epitopes. Here we present TCRen, a structure-based method for ranking candidate unseen epitopes for a given TCR. The first stage of the TCRen pipeline is modeling of the TCR-peptide-major histocompatibility complex structure. Then a TCR-peptide residue contact map is extracted from this structure and used to rank all candidate epitopes on the basis of an interaction score with the target TCR. Scoring is performed using an energy potential derived from the statistics of TCR-peptide contact preferences in existing crystal structures. We show that TCRen has high performance in discriminating cognate versus unrelated peptides and can facilitate the identification of cancer neoepitopes recognized by tumor-infiltrating lymphocytes.

The scorpion toxin BeKm-1 blocks hERG cardiac potassium channels using an indispensable arginine residue.

BeKm-1 is a peptide toxin from scorpion venom that blocks the pore of the potassium channel hERG (Kv11.1) in the human heart. Although individual protein structures have been resolved, the structure of the complex between hERG and BeKm-1 is unknown. Here, we used molecular dynamics and ensemble docking, guided by previous double-mutant cycle analysis data, to obtain an in silico model of the hERG-BeKm-1 complex. Adding to the previous mutagenesis study of BeKm-1, our model uncovers the key role of residue Arg20, which forms three interactions (a salt bridge and hydrogen bonds) with the channel vestibule simultaneously. Replacement of this residue even by lysine weakens the interactions significantly. In accordance, the recombinantly produced BeKm-1R20K mutant exhibited dramatically decreased activity on hERG. Our model may be useful for future drug design attempts.

A single missense mutation in a coiled-coil domain of Escherichia coli ribosomal protein S2 confers a thermosensitive phenotype that can be suppressed by ribosomal protein S1.

Ribosomal protein S2 is an essential component of translation machinery, and its viable mutated variants conferring distinct phenotypes serve as a valuable tool in studying the role of S2 in translation regulation. One of a few available rpsB mutants, rpsB1, shows thermosensitivity and ensures enhanced expression of leaderless mRNAs. In this study, we identified the nature of the rpsB1 mutation. Sequencing of the rpsB1 allele revealed a G-to-A transition in the part of the rpsB gene which encodes a coiled-coil domain of S2. The resulting E132K substitution resides in a highly conserved site, TKKE, a so-called N-terminal capping box, at the beginning of the second alpha helix. The protruding coiled-coil domain of S2 is known to provide binding with 16S rRNA in the head of the 30S subunit and, in addition, to interact with a key mRNA binding protein, S1. Molecular dynamics simulations revealed a detrimental impact of the E132K mutation on the coiled-coil structure and thereby on the interactions between S2 and 16S rRNA, providing a clue for the thermosensitivity of the rpsB1 mutant. Using a strain producing a leaderless lacZ transcript from the chromosomal lac promoter, we demonstrated that not only the rpsB1 mutation generating S2/S1-deficient ribosomes but also the rpsA::IS10 mutation leading to partial deficiency in S1 alone increased translation efficiency of the leaderless mRNA by about 10-fold. Moderate overexpression of S1 relieved all these effects and, moreover, suppressed the thermosensitive phenotype of rpsB1, indicating the role of S1 as an extragenic suppressor of the E132K mutation.

Molecular Basis for Mambalgin-2 Interaction with Heterotrimeric α-ENaC/ASIC1a/γ-ENaC Channels in Cancer Cells.

Cancer progression is characterized by microenvironmental acidification. Tumor cells adapt to low environmental pH by activating acid-sensing trimeric ion channels of the DEG/ENaC family. The α-ENaC/ASIC1a/γ-ENaC heterotrimeric channel is a tumor-specific acid-sensing channel, and its targeting can be considered a new strategy for cancer therapy. Mambalgin-2 from the Dendroaspis polylepis venom inhibits the α-ENaC/ASIC1a/γ-ENaC heterotrimer more effectively than the homotrimeric ASIC1a channel, initially proposed as the target of mambalgin-2. Although the molecular basis of such mambalgin selectivity remained unclear. Here, we built the models of the complexes of mambalgin-2 with the α-ENaC/ASIC1a/γ-ENaC and ASIC1a channels, performed MD and predicted the difference in the binding modes. The importance of the 'head' loop region of mambalgin-2 for the interaction with the hetero-, but not with the homotrimeric channel was confirmed by site-directed mutagenesis and electrophysiology. A new mode of allosteric regulation of the ENaC channels by linking the thumb domain of the ASIC1a subunit with the palm domain of the γ-ENaC subunit was proposed. The data obtained provide new insights into the regulation of various types of acid-sensing ion channels and the development of new strategies for cancer treatment.

Membrane-mediated interaction of non-conventional snake three-finger toxins with nicotinic acetylcholine receptors.

Nicotinic acetylcholine receptor of α7 type (α7-nAChR) presented in the nervous and immune systems and epithelium is a promising therapeutic target for cognitive disfunctions and cancer treatment. Weak toxin from Naja kaouthia venom (WTX) is a non-conventional three-finger neurotoxin, targeting α7-nAChR with weak affinity. There are no data on interaction mode of non-conventional neurotoxins with nAChRs. Using α-bungarotoxin (classical three-finger neurotoxin with high affinity to α7-nAChR), we showed applicability of cryo-EM to study complexes of α7-nAChR extracellular ligand-binding domain (α7-ECD) with toxins. Using cryo-EM structure of the α7-ECD/WTX complex, together with NMR data on membrane active site in the WTX molecule and mutagenesis data, we reconstruct the structure of α7-nAChR/WTX complex in the membrane environment. WTX interacts at the entrance to the orthosteric site located at the receptor intersubunit interface and simultaneously forms the contacts with the membrane surface. WTX interaction mode with α7-nAChR significantly differs from α-bungarotoxin's one, which does not contact the membrane. Our study reveals the important role of the membrane for interaction of non-conventional neurotoxins with the nicotinic receptors.

Molecular Dynamics Insight into the Lipid II Recognition by Type A Lantibiotics: Nisin, Epidermin, and Gallidermin.

Lanthionine-containing peptides (lantibiotics) have been considered as pharmaceutical candidates for decades, although their clinical application has been restricted. Most lantibiotics kill bacteria via targeting and segregating of the cell wall precursor-membrane-inserted lipid II molecule-in some cases accompanied by pores formation. Nisin-like lantibiotics specifically bind to pyrophosphate (PPi) moiety of lipid II with their structurally similar N-terminal thioether rings A and B. Although possessing higher pore-forming capability, nisin, in some cases, is 10-fold less efficient in vivo as compared to related epidermin and gallidermin peptides, differing just in a few amino acid residues within their target-binding regions. Here, using molecular dynamics simulations, we investigated atomistic details of intermolecular interactions between the truncated analogues of these peptides (residues 1-12) and lipid II mimic (dimethyl pyrophosphate, DMPPi). The peptides adopt similar conformation upon DMPPi binding with backbone amide protons orienting into a single center capturing PPi moiety via simultaneous formation of up to seven hydrogen bonds. Epidermin and gallidermin adopt the complex-forming conformation twice as frequent as nisin does, enhancing the binding by the lysine 4 side chain. Introduction of the similar residue to nisin in silico improves the binding, providing ideas for further design of prototypic antibiotics.

Biochemical Basis of Skin Disease Mal de Meleda: SLURP-1 Mutants Differently Affect Keratinocyte Proliferation and Apoptosis.

Mal de Meleda is an autosomal recessive palmoplantar keratoderma associated with mutations in a gene encoding SLURP-1. SLURP-1 controls growth, differentiation, and apoptosis of keratinocytes by interaction with α7-type nicotinic acetylcholine receptors. SLURP-1 has a three-finger structure with a ß-structural core (head) and three prolonged loops (fingers). To determine the role of SLURP-1 mutations, we produced 22 mutant variants of the protein, including those involved in Mal de Meleda pathogenesis. All mutants except R71H, R71P, T52A, R96P, and L98P were produced in the folded form. SLURP-1 reduces the growth of Het-1A keratinocytes; thus, we studied the influence of the mutations on its antiproliferative activity. Mutations in loops I and III led to the protein inactivation, whereas most mutations in loop II increased SLURP-1 antiproliferative activity. Alanine substitutions of R96 and L98 residues located in the protein head resulted in the appearance of additional pro-apoptotic activity. Our results agree with the diversity of Mal de Meleda phenotypes. Using obtained functional data, the SLURP-1/α7 type nicotinic acetylcholine receptor complex was modeled in silico. Our study provides functional and structural information about the role of the SLURP-1 mutations in Mal de Meleda pathogenesis and predicts SLURP-1 variants, which could drive the disease.