2StrucCompare: a webserver for visualizing small but noteworthy differences between protein tertiary structures through interrogation of the secondary structure content

Temperature-dependent structural changes in xylanase II from Trichoderma longibrachiatum

Endo-β-1,4-xylanases degrade heteroxylans that constitute the lignocellulosic plant cell wall. This enzyme is widely used in the food, paper, textile, and biorefinery industries. Temperature affects the optimum activity of xylanase and is an important factor in its application. Various structural analyses of xylanase have been performed, but its structural influence by temperature is not fully elucidated. To better understand the structural influence of xylanase due to temperature, the crystal structure of xylanase II from Trichoderma longibrachiatum (TloXynII) at room and cryogenic temperatures was determined at 2.1 and 1.9 Å resolution, respectively. The room-temperature structure of TloXynII (TloXynIIRT) showed a B-factor value 2.09 times higher than that of the cryogenic-temperature structure of TloXynII (TloXynIICryo). Subtle movement of the catalytic and substrate binding residues was observed between TloXynIIRT and TloXynIICryo. In TloXynIIRT, the thumb domain exhibited high flexibility, whereas in TloXynIICryo, the finger domain exhibited high flexibility. The substrate binding cleft of TloXynIIRT was narrower than that of TloXynIICryo, indicating a distinct finger domain conformation. Numerous water molecule networks were observed in the substrate binding cleft of TloXynIICryo, whereas only a few water molecules were observed in TloXynIIRT. These structural analyses expand our understanding of the temperature-dependent conformational changes in xylanase.

Activation of CBASS Cap5 endonuclease immune effector by cyclic nucleotides

The bacterial cyclic oligonucleotide-based antiphage signaling system (CBASS) is similar to the cGAS-STING system in humans, containing an enzyme that synthesizes a cyclic nucleotide on viral infection and an effector that senses the second messenger for the antiviral response. Cap5, containing a SAVED domain coupled to an HNH DNA endonuclease domain, is the most abundant CBASS effector, yet the mechanism by which it becomes activated for cell killing remains unknown. We present here high-resolution structures of full-length Cap5 from Pseudomonas syringae (Ps) with second messengers. The key to PsCap5 activation is a dimer-to-tetramer transition, whereby the binding of second messenger to dimer triggers an open-to-closed transformation of the SAVED domains, furnishing a surface for assembly of the tetramer. This movement propagates to the HNH domains, juxtaposing and converting two HNH domains into states for DNA destruction. These results show how Cap5 effects bacterial cell suicide and we provide proof-in-principle data that the CBASS can be extrinsically activated to limit bacterial infections.

DichroPipeline: A suite of online and downloadable tools and resources for protein circular dichroism spectroscopic data analyses, interpretations, and their interoperability with other bioinformatics tools and resources

Circular Dichroism (CD) spectroscopy is a widely-used method for characterizing individual protein structures in solutions, membranes, films and macromolecular complexes, as well as for probing macromolecular interactions, conformational changes associated with binding substrates, and in different functionally-related environments. This paper describes a series of related computational and display tools that have been developed over many years to aid in those characterizations and functional interpretations. The new DichroPipeline described herein links a series of format-compatible data processing, analysis, and display tools to enable users to facilely produce the spectra, which can then be made available in the Protein Circular Dichroism Data Bank (https://pcddb.cryst.bbk.ac.uk/) resource, in which the CD spectral and associated metadata for each entry are linked to other structural and functional data bases including the Protein Data Bank (PDB), and the UniProt sequence data base, amongst others. These tools and resources thus provide the basis for a wide range of traceable structural characterizations of soluble, membrane and intrinsically-disordered proteins.

Slc11 Synapomorphy: A Conserved 3D Framework Articulating Carrier Conformation Switch

Transmembrane carriers of the Slc11 family catalyze proton (H+)-dependent uptake of divalent metal ions (Me2+) such as manganese and iron-vital elements coveted during infection. The Slc11 mechanism of high-affinity Me2+ cell import is selective and conserved between prokaryotic (MntH) and eukaryotic (Nramp) homologs, though processes coupling the use of the proton motive force to Me2+ uptake evolved repeatedly. Adding bacterial piracy of Nramp genes spread in distinct environmental niches suggests selective gain of function that may benefit opportunistic pathogens. To better understand Slc11 evolution, Alphafold (AF2)/Colabfold (CF) 3D predictions for bacterial sequences from sister clades of eukaryotic descent (MCb and MCg) were compared using both native and mutant templates. AF2/CF model an array of native MCb intermediates spanning the transition from outwardly open (OO) to inwardly open (IO) carriers. In silico mutagenesis targeting (i) a set of (evolutionarily coupled) sites that may define Slc11 function (putative synapomorphy) and (ii) residues from networked communities evolving during MCb transition indicates that Slc11 synapomorphy primarily instructs a Me2+-selective conformation switch which unlocks carrier inner gate and contributes to Me2+ binding site occlusion and outer gate locking. Inner gate opening apparently proceeds from interaction between transmembrane helix (h) h5, h8 and h1a. MCg1 xenologs revealed marked differences in carrier shape and plasticity, owing partly to an altered intramolecular H+ network. Yet, targeting Slc11 synapomorphy also converted MCg1 IO models to an OO state, apparently mobilizing the same residues to control gates. But MCg1 response to mutagenesis differed, with extensive divergence within this clade correlating with MCb-like modeling properties. Notably, MCg1 divergent epistasis marks the emergence of the genus Bordetella-Achromobacter. Slc11 synapomorphy localizes to the 3D areas that deviate least among MCb and MCg1 models (either IO or OO) implying that it constitutes a 3D network of residues articulating a Me2+-selective carrier conformation switch which is maintained in fast-evolving clades at the cost of divergent epistatic interactions impacting carrier shape and dynamics.

Structure determination and analysis of titin A-band fibronectin type III domains provides insights for disease-linked variants and protein oligomerisation

Titin is the largest protein found in nature and spans half a sarcomere in vertebrate striated muscle. The protein has multiple functions, including in the organisation of the thick filament and acting as a molecular spring during the muscle contraction cycle. Missense variants in titin have been linked to both cardiac and skeletal myopathies. Titin is primarily composed of tandem repeats of immunoglobulin and fibronectin type III (Fn3) domains in a variety of repeat patterns; however, the vast majority of these domains have not had their high-resolution structure determined experimentally. Here, we present the crystal structures of seven wild type titin Fn3 domains and two harbouring rare missense variants reported in hypertrophic cardiomyopathy (HCM) patients. All domains present the typical Fn3 fold, with the domains harbouring variants reported in HCM patients retaining the wild-type conformation. The effect on domain folding and stability were assessed for five rare missense variants found in HCM patients: four caused thermal destabilization of between 7 and 13 °C and one prevented the folding of its domain. The structures also allowed us to locate the positions of residues whose mutations have been linked to congenital myopathies and rationalise how they convey their deleterious effects. We find no evidence of physiological homodimer formation, excluding one hypothesised mechanism as to how titin variants could exert pathological effects.

Two new Scianna variants causing loss of high prevalence antigens: ERMAP model and 3D analysis of the antigens

BACKGROUND

Scianna (Sc) antigens, seven high and two of low prevalence, are expressed on erythrocyte membrane-associated protein (ERMAP). We investigated SC (ERMAP) in individuals who made antibodies to high prevalence Scianna antigens, and propose a 3D model for ERMAP to precisely localize the residues associated with the known antigens.

METHODS

Serological testing and DNA sequencing was performed by standard methods. A 3D structural model was built using a multi-template homology approach. Protein structures representing missense variants associated with the loss or gain of an antigen were generated. Residue accessibility and intraprotein interactions were compared with the wild-type protein.

RESULTS

Two new SC alleles, one with c.349C > T (p.Arg117Cys) in a woman from South India with anti-Sc3 in her plasma, and a c.217_219delinsTGT (p.Arg73Cys) in an African-American woman with an antibody to a new high prevalence antigen, termed SCAC, were identified. Six structural templates were used to model ERMAP. 3D analysis showed that residues key for Scianna antigen expression were all exposed at the surface of the extracellular domain. The p.Arg117Cys change was predicted to abolish interactions between residues 93 and 117, with no compensating interactions.

CONCLUSION

We confirm the extracellular location of Scianna residues responsible for antigen expression which predicts direct accessibility to antibodies. Loss of intraprotein interactions appear to be responsible for a Sc null and production of anti-Sc3 with p.117Cys, SC*01 N.03, and for loss of a high prevalence antigen with p.73Cys, termed SCAC for Sc Arg to Cys. Comparative modeling aids our understanding of new alleles and Scianna antigen expression.

Atomic Simulation of the Binding of JAK1 and JAK2 with the Selective Inhibitor Ruxolitinib

Rheumatoid arthritis belongs to the group of chronic systemic autoimmune diseases characterized by the development of destructive synovitis and extra-articular manifestations. Cytokines regulate a wide range of inflammatory processes involved in the pathogenesis of rheumatoid arthritis and contribute to the induction of autoimmunity and chronic inflammation. Janus-associated kinase (JAK) and signal transducer and activator of transcription (STAT) proteins mediate cell signaling from cytokine receptors, and are involved in the pathogenesis of autoimmune and inflammatory diseases. Targeted small-molecule drugs that inhibit the functional activity of JAK proteins are used in clinical practice for the treatment of rheumatoid arthritis. In our study, we modeled the interactions of the small-molecule drug ruxolitinib with JAK1 and JAK2 isoforms and determined the binding selectivity using molecular docking. Molecular modeling data show that ruxolitinib selectively binds the JAK1 and JAK2 isoforms with a binding affinity of −8.3 and −8.0 kcal/mol, respectively. The stabilization of ligands in the cavity of kinases occurs primarily through hydrophobic interactions. The amino acid residues of the protein globules of kinases that are responsible for the correct positioning of the drug ruxolitinib and its retention have been determined.

Molecular Dynamics Study of Citrullinated Proteins Associated with the Development of Rheumatoid Arthritis

Biological activity regulation by protein post-translational modification (PTM) is critical for cell function, development, differentiation, and survival. Dysregulation of PTM proteins is present in various pathological conditions, including rheumatoid arthritis (RA). RA is a systemic autoimmune disease that primarily affects joints, and there are three main types of protein PTMs associated with the development of this disease, namely, glycosylation, citrullination, and carbamylation. Glycosylation is important for the processing and presentation of antigen fragments on the cell surface and can modulate immunoglobulin activity. The citrullination of autoantigens is closely associated with RA, as evidenced by the presence of antibodies specific to citrullinated proteins in the serum of patients. Carbamylation and dysregulation have recently been associated with RA development in humans.In this study, we performed an overview analysis of proteins with post-translational modifications associated with the development of RA adverted in peer-reviewed scientific papers for the past 20 years. As a result of the search, a list of target proteins and corresponding amino acid sequences with PTM in RA was formed. Structural characteristics of the listed modified proteins were extracted from the Protein Data Bank. Then, molecular dynamics experiments of intact protein structures and corresponding structures with PTMs were performed regarding structures in the list announced in the ProtDB service. This study aimed to conduct a molecular dynamics study of intact proteins and proteins, including post-translational modification and protein citrullination, likely associated with RA development. We observed another exhibition of the fundamental physics concept, symmetry, at the submolecular level, unveiled as the autonomous repetitions of outside the protein structural motif performance globule corresponding to those in the whole protein molecule.

Validation of determinate (dt) gene-based DNA marker in inter-specific hybrid sesame and in-silico analysis of the predicted dt protein structures

Determinacy is a desirable trait in sesame, an important oilseed crop. We have developed an inter-specific hybrid between basally branched indeterminate cultivated Sesamum indicum genotype and wild S. prostratum with no branching yet synchronous pods on the shoot. The hybrid and a few exotic sesame germplasms were successfully screened with a determinacy (dt) gene-based DNA marker. In-silico translation of the partial coding sequences of the dt gene from the two contrasting parent genotypes revealed an SNP (V159A) in S. prostratum. The predicted cytoplasmic dt protein showed a high resemblance with flowering protein centroradialis.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01135-1.

Effect of setting data collection parameters on the reliability of a circular dichroism spectrum

Circular dichroism (CD) spectroscopy is a well-established biophysical technique used to investigate the structure of molecules. The analysis of a protein CD spectrum depends on the quality of the original CD data, which can be affected by the sample purity, background absorption of the additives/solvent/buffer, the choice of the parameters used for data collection, etc. In this paper, the CD spectrum of myoglobin was used as a model to exploit how variations on each data collection parameter could affect the final protein CD spectrum and, the subsequent effect of them on the quantitative analysis of protein secondary structure. Bioinformatics analysis carried out with SESCA package and PDBMD2CD server predicted a theoretical myoglobin CD spectrum, and a Monte Carlo-like model was implemented to estimate the uncertainty in secondary structure predictions performed with CDSSTR, Selcon 3 and ContinLL algorithms. An inappropriate choice of data collection parameters can lead to a misinterpretation of the CD data in terms of the protein structural content.

Bioinformatic analysis of subfamily-specific regions in 3D-structures of homologs to study functional diversity and conformational plasticity in protein superfamilies

Local 3D-structural differences in homologous proteins contribute to functional diversity observed in a superfamily, but so far received little attention as bioinformatic analysis was usually carried out at the level of amino acid sequences. We have developed Zebra3D - the first-of-its-kind bioinformatic software for systematic analysis of 3D-alignments of protein families using machine learning. The new tool identifies subfamily-specific regions (SSRs) - patterns of local 3D-structure (i.e. single residues, loops, or secondary structure fragments) that are spatially equivalent within families/subfamilies, but are different among them, and thus can be associated with functional diversity and function-related conformational plasticity. Bioinformatic analysis of protein superfamilies by Zebra3D can be used to study 3D-determinants of catalytic activity and specific accommodation of ligands, help to prepare focused libraries for directed evolution or assist development of chimeric enzymes with novel properties by exchange of equivalent regions between homologs, and to characterize plasticity in binding sites. A companion Mustguseal web-server is available to automatically construct a 3D-alignment of functionally diverse proteins, thus reducing the minimal input required to operate Zebra3D to a single PDB code. The Zebra3D + Mustguseal combined approach provides the opportunity to systematically explore the value of SSRs in superfamilies and to use this information for protein design and drug discovery. The software is available open-access at https://biokinet.belozersky.msu.ru/Zebra3D.

LAHMA: structure analysis through local annotation of homology-matched amino acids

Comparison of homologous structure models is a key step in analyzing protein structure. With a wealth of homologous structures, comparison becomes a tedious process, and often only a small (user-biased) selection of data is used. A multitude of structural superposition algorithms are then typically used to visualize the structures together in 3D and to compare them. Here, the Local Annotation of Homology-Matched Amino acids (LAHMA) website (https://lahma.pdb-redo.eu) is presented, which compares any structure model with all of its close homologs from the PDB-REDO databank. LAHMA displays structural features in sequence space, allowing users to uncover differences between homologous structure models that can be analyzed for their relevance to chemistry or biology. LAHMA visualizes numerous structural features, also allowing one-click comparison of structure-quality plots (for example the Ramachandran plot) and `in-browser' structural visualization of 3D models.

Identification of a druggable binding pocket in the spike protein reveals a key site for existing drugs potentially capable of combating Covid-19 infectivity

BACKGROUND

Following the recent outbreak of the new coronavirus pandemic (Covid-19), the rapid determination of the structure of the homo-trimeric spike glycoprotein has prompted the study reported here. The aims were to identify potential "druggable" binding pockets in the protein and, if located, to virtual screen pharmaceutical agents currently in use for predicted affinity to these pockets which might be useful to restrict, reduce, or inhibit the infectivity of the virion.

RESULTS

Our analyses of this structure have revealed a key potentially druggable pocket where it might be viable to bind pharmaceutical agents to inhibit its ability to infect human cells. This pocket is found at the inter-chain interface that exists between two domains prior to the virion binding to human Angiotensin Converting Enzyme 2 (ACE2) protein. One of these domains is the highly mobile receptor binding domain, which must move into position to interact with ACE2, which is an essential feature for viral entry to the host cell. Virtual screening with a library of purchasable drug molecules has identified pharmaceuticals currently in use as prescription and over the counter medications that, in silico, readily bind into this pocket.

CONCLUSIONS

This study highlights possible drugs already in use as pharmaceuticals that may act as agents to interfere with the movements of the domains within this protein essential for the infectivity processes and hence might slow, or even halt, the infection of host cells by this new coronavirus. As these are existing pharmaceuticals already approved for use in humans, this knowledge could accelerate their roll-out, through repurposing, for affected individuals and help guide the efforts of other researchers in finding effective treatments for the disease.

The key role of E418 carboxyl group in the formation of Nt.BspD6I nickase active site: Structural and functional properties of Nt.BspD6I E418A mutant

The mutated nickase Nt.BspD6I E418A has been obtained by site-directed mutagenesis. The purified protein has been crystallized, and its spatial structure has been determined at 2.45 Å resolution. An analysis of the crystal structures of the wild-type and mutated nickase have shown that the elimination of a carboxyl group due to the E418A mutation initiates marked conformational changes in both the N-terminal recognition domain and the C-terminal catalytic domain of nickase and insignificantly affects its linker domain. This is supported by changes in the functional properties of mutated nickase: an increase in the oligomerization capacity in the presence of a substrate, a reduction in the capacity to bind a substrate, and complete loss of catalytic activity.