Single Residue Substitution Research Articles

Strategic Single-Residue Substitution in the Antimicrobial Peptide Esc(1-21) Confers Activity against Staphylococcus aureus, Including Drug-Resistant and Biofilm Phenotype.

Staphylococcus aureus, a bacterium resistant to multiple drugs, is a significant cause of illness and death worldwide. Antimicrobial peptides (AMPs) provide an excellent potential strategy to cope with this threat. Recently, we characterized a derivative of the frog-skin AMP esculentin-1a, Esc(1-21) (1) that is endowed with potent activity against Gram-negative bacteria but poor efficacy against Gram-positive strains. In this study, three analogues of peptide 1 were designed by replacing Gly8 with α-aminoisobutyric acid (Aib), Pro, and dPro (2-4, respectively). The single substitution Gly8 → Aib8 in peptide 2 makes it active against the planktonic form of Gram-positive bacterial strains, especially Staphylococcus aureus, including multidrug-resistant clinical isolates, with an improved biostability without resulting in cytotoxicity to mammalian cells. Moreover, peptide 2 showed a higher antibiofilm activity than peptide 1 against both reference and clinical isolates of S. aureus. Peptide 2 was also able to induce rapid bacterial killing, suggesting a membrane-perturbing mechanism of action. Structural analysis of the most active peptide 2 evidenced that the improved biological activity of peptide 2 is the consequence of a combination of higher biostability, higher α helical content, and ability to reduce membrane fluidity and to adopt a distorted helix, bent in correspondence of Aib8. Overall, this study has shown how a strategic single amino acid substitution is sufficient to enlarge the spectrum of activity of the original peptide 1, and improve its biological properties for therapeutic purposes, thus paving the way to optimize AMPs for the development of new broad-spectrum anti-infective agents.

Structural and functional insights into the 2′-O-methyltransferase of SARS-CoV-2

A unique feature of coronaviruses is their utilization of self-encoded nonstructural protein 16 (nsp16), 2′-O-methyltransferase (2′-O-MTase), to cap their RNAs through ribose 2′-O-methylation modification. This process is crucial for maintaining viral genome stability, facilitating efficient translation, and enabling immune escape. Despite considerable advances in the ultrastructure of SARS-CoV-2 nsp16/nsp10, insights into its molecular mechanism have so far been limited. In this study, we systematically characterized the 2′-O-MTase activity of nsp16 in SARS-CoV-2, focusing on its dependence on nsp10 stimulation. We observed cross-reactivity between nsp16 and nsp10 in various coronaviruses due to a conserved interaction interface. However, a single residue substitution (K58T) in SARS-CoV-2 nsp10 restricted the functional activation of MERS-CoV nsp16. Furthermore, the cofactor nsp10 effectively enhanced the binding of nsp16 to the substrate RNA and the methyl donor S-adenosyl-l-methionine (SAM). Mechanistically, His-80, Lys-93, and Gly-94 of nsp10 interacted with Asp-102, Ser-105, and Asp-106 of nsp16, respectively, thereby effectively stabilizing the SAM binding pocket. Lys-43 of nsp10 interacted with Lys-38 and Gly-39 of nsp16 to dynamically regulate the RNA binding pocket and facilitate precise binding of RNA to the nsp16/nsp10 complex. By assessing the conformational epitopes of nsp16/nsp10 complex, we further determined the critical residues involved in 2′-O-MTase activity. Additionally, we utilized an in vitro biochemical platform to screen potential inhibitors targeting 2′-O-MTase activity. Overall, our results significantly enhance the understanding of viral 2′-O methylation process and mechanism, providing valuable targets for antiviral drug development.

CRISPR/Cas9-mediated genome editing of the thymidine kinase gene in a clinical HSV-1 isolate identifies F289S as novel acyclovir-resistant mutation

Herpes simplex virus type 1 (HSV-1) is a neurotropic alphaherpesvirus that establishes a lifelong infection in sensory neurons of infected individuals, accompanied with intermittent reactivation of latent virus causing (a)symptomatic virus shedding. Whereas acyclovir (ACV) is a safe and highly effective antiviral to treat HSV-1 infections, long-term usage can lead to emergence of ACV resistant (ACVR) HSV-1 and subsequently ACV refractory disease. Here, we isolated an HSV-1 strain from a patient with reactivated herpetic eye disease that did not respond to ACV treatment. The isolate carried a novel non-synonymous F289S mutation in the viral UL23 gene encoding the thymidine kinase (TK) protein. Because ACV needs conversion by viral TK and subsequently cellular kinases to inhibit HSV-1 replication, the UL23 gene is commonly mutated in ACVR HSV-1 strains. The potential role of the F289S mutation causing ACVR was investigated using CRISPR/Cas9-mediated HSV-1 genome editing. Reverting the F289S mutation in the original clinical isolate to the wild-type sequence S289F resulted in an ACV-sensitive (ACVS) phenotype, and introduction of the F289S substitution in an ACVS HSV-1 reference strain led to an ACVR phenotype. In summary, we identified a new HSV-1 TK mutation in the eye of a patient with ACV refractory herpetic eye disease, which was identified as the causative ACVR mutation with the aid of CRISPR/Cas9-mediated genome engineering technology. Direct editing of clinical HSV-1 isolates by CRISPR/Cas9 is a powerful strategy to assess whether single residue substitutions are causative to a clinical ACVR phenotype.

Destabilization and Degradation of a Disease-Linked PGM1 Protein Variant.

PGM1-linked congenital disorder of glycosylation (PGM1-CDG) is an autosomal recessive disease characterized by several phenotypes, some of which are life-threatening. Research focusing on the disease-related variants of the α-D-phosphoglucomutase 1 (PGM1) protein has shown that several are insoluble in vitro and expressed at low levels in patient fibroblasts. Due to these observations, we hypothesized that some disease-linked PGM1 protein variants are structurally destabilized and subject to protein quality control (PQC) and rapid intracellular degradation. Employing yeast-based assays, we show that a disease-associated human variant, PGM1 L516P, is insoluble, inactive, and highly susceptible to ubiquitylation and rapid degradation by the proteasome. In addition, we show that PGM1 L516P forms aggregates in S. cerevisiae and that both the aggregation pattern and the abundance of PGM1 L516P are chaperone-dependent. Finally, using computational methods, we perform saturation mutagenesis to assess the impact of all possible single residue substitutions in the PGM1 protein. These analyses identify numerous missense variants with predicted detrimental effects on protein function and stability. We suggest that many disease-linked PGM1 variants are subject to PQC-linked degradation and that our in silico site-saturated data set may assist in the mechanistic interpretation of PGM1 variants.

Influence of Model Structures on Predictors of Protein Stability Changes from Single-Point Mutations.

Missense variation in genomes can affect protein structure stability and, in turn, the cell physiology behavior. Predicting the impact of those variations is relevant, and the best-performing computational tools exploit the protein structure information. However, most of the current protein sequence variants are unresolved, and comparative or ab initio tools can provide a structure. Here, we evaluate the impact of model structures, compared to experimental structures, on the predictors of protein stability changes upon single-point mutations, where no significant changes are expected between the original and the mutated structures. We show that there are substantial differences among the computational tools. Methods that rely on coarse-grained representation are less sensitive to the underlying protein structures. In contrast, tools that exploit more detailed molecular representations are sensible to structures generated from comparative modeling, even on single-residue substitutions.

Identification of three elevenin receptors and roles of elevenin disulfide bond and residues in receptor activation in Aplysia californica

Neuropeptides are ubiquitous intercellular signaling molecules in the CNS and play diverse roles in modulating physiological functions by acting on specific G-protein coupled receptors (GPCRs). Among them, the elevenin signaling system is now believed to be present primarily in protostomes. Although elevenin was first identified from the L11 neuron of the abdominal ganglion in mollusc Aplysia californica, no receptors have been described in Aplysia, nor in any other molluscs. Here, using two elevenin receptors in annelid Platynereis dumerilii, we found three putative elevenin GPCRs in Aplysia. We cloned the three receptors and tentatively named them apElevR1, apElevR2, and apElevR3. Using an inositol monophosphate (IP1) accumulation assay, we demonstrated that Aplysia elevenin with the disulfide bond activated the three putative receptors with low EC50 values (ranging from 1.2 to 25 nM), supporting that they are true receptors for elevenin. In contrast, elevenin without the disulfide bond could not activate the receptors, indicating that the disulfide bond is required for receptor activity. Using alanine substitution of individual conserved residues other than the two cysteines, we showed that these residues appear to be critical to receptor activity, and the three different receptors had different sensitivities to the single residue substitution. Finally, we examined the roles of those residues outside the disulfide bond ring by removing these residues and found that they also appeared to be important to receptor activity. Thus, our study provides an important basis for further study of the functions of elevenin and its receptors in Aplysia and other molluscs.

Inhibition of translation termination by the antimicrobial peptide Drosocin.

The proline-rich antimicrobial peptide (PrAMP) Drosocin (Dro) from fruit flies shows sequence similarity to other PrAMPs that bind to the ribosome and inhibit protein synthesis by varying mechanisms. The target and mechanism of action of Dro, however, remain unknown. Here we show that Dro arrests ribosomes at stop codons, probably sequestering class 1 release factors associated with the ribosome. This mode of action is comparable to that of apidaecin (Api) from honeybees, making Dro the second member of the type II PrAMP class. Nonetheless, analysis of a comprehensive library of endogenously expressed Dro mutants shows that the interactions of Dro and Api with the target are markedly distinct. While only a few C-terminal amino acids of Api are critical for binding, the interaction of Dro with the ribosome relies on multiple amino acid residues distributed throughout the PrAMP. Single-residue substitutions can substantially enhance the on-target activity of Dro.

Amino acid substitutions in human growth hormone affect secondary structure and receptor binding.

The interaction between human Growth Hormone (hGH) and hGH Receptor (hGHR) has basic relevance to cancer and growth disorders, and hGH is the scaffold for Pegvisomant, an anti-acromegaly therapeutic. For the latter reason, hGH has been extensively engineered by early workers to improve binding and other properties. We are particularly interested in E174 which belongs to the hGH zinc-binding triad; the substitution E174A is known to significantly increase binding, but to now no explanation has been offered. We generated this and several computationally-selected single-residue substitutions at the hGHR-binding site of hGH. We find that, while many successfully slow down dissociation of the hGH-hGHR complex once bound, they also slow down the association of hGH to hGHR. The E174A substitution induces a change in the Circular Dichroism spectrum that suggests the appearance of coiled-coiling. Here we show that E174A increases affinity of hGH against hGHR because the off-rate is slowed down more than the on-rate. For E174Y (and certain mutations at other sites) the slowdown in on-rate was greater than that of the off-rate, leading to decreased affinity. The results point to a link between structure, zinc binding, and hGHR-binding affinity in hGH.

Characterization of an Aplysia vasotocin signaling system and actions of posttranslational modifications and individual residues of the ligand on receptor activity.

The vasopressin/oxytocin signaling system is present in both protostomes and deuterostomes and plays various physiological roles. Although there were reports for both vasopressin-like peptides and receptors in mollusc Lymnaea and Octopus, no precursor or receptors have been described in mollusc Aplysia. Here, through bioinformatics, molecular and cellular biology, we identified both the precursor and two receptors for Aplysia vasopressin-like peptide, which we named Aplysia vasotocin (apVT). The precursor provides evidence for the exact sequence of apVT, which is identical to conopressin G from cone snail venom, and contains 9 amino acids, with two cysteines at position 1 and 6, similar to nearly all vasopressin-like peptides. Through inositol monophosphate (IP1) accumulation assay, we demonstrated that two of the three putative receptors we cloned from Aplysia cDNA are true receptors for apVT. We named the two receptors as apVTR1 and apVTR2. We then determined the roles of post-translational modifications (PTMs) of apVT, i.e., the disulfide bond between two cysteines and the C-terminal amidation on receptor activity. Both the disulfide bond and amidation were critical for the activation of the two receptors. Cross-activity with conopressin S, annetocin from an annelid, and vertebrate oxytocin showed that although all three ligands can activate both receptors, the potency of these peptides differed depending on their residue variations from apVT. We, therefore, tested the roles of each residue through alanine substitution and found that each substitution could reduce the potency of the peptide analog, and substitution of the residues within the disulfide bond tended to have a larger impact on receptor activity than the substitution of those outside the bond. Moreover, the two receptors had different sensitivities to the PTMs and single residue substitutions. Thus, we have characterized the Aplysia vasotocin signaling system and showed how the PTMs and individual residues in the ligand contributed to receptor activity.

Variation in the tonoplast cadmium transporter heavy metal ATPase 3 (HMA3) homolog gene in Aegilops tauschii.

The functionality of HMA3 is a key determinant controlling Cd accumulation in the shoots and grains of plants. Wild relatives of modern crop plants can serve as sources of valuable genetic variation for various traits. Here, resequencing of HMA3 homoeologous genes from Aegilops tauschii (the donor of the wheat D genome) was carried out to identify natural variation at both the nucleotide and polypeptide levels. HMA3 homoeologs are highly conserved, and 10 haplotypes were revealed based on 19 single nucleotide polymorphisms (eight induced single amino acid residue substitutions, including 2 altered amino acids in transmembrane domains) in 80 widely distributed Ae. tauschii accessions. The results provide genetic resources for low/no Cd concentration wheat improvement.

High-throughput investigation of temperature- and ligand-activation of the TRPM8 channel.

The TRPM8 channel functions as the main cold sensor in sensory nerves, where decreases in temperature cause opening of the cation-conducting channel pore. Activation of TRPM8 in neurons is modulated by voltage, lipids, G-proteins, protein kinases, cytosolic calcium and natural compounds such as menthol. These mechanisms tune neural excitability and are likely important for how sensory stimuli are detected, integrated, and perceived as cold, warmth, pain, or analgesia. The molecular mechanisms of activation and regulation of the TRPM8 channel by chemical and physical stimuli are largely unknown, including the mechanism of cold-activation. Here we set out to systematically identify which residues in the transmembrane domain of the TRPM8 channel are determinant for its sensitivity to menthol or cold. We generated a library containing all possible single-residue substitutions and deletions in the transmembrane domain of the channel. This library was benchmarked using patch-clamp electrophysiology and we have developed a high-throughput fluorescence-based assay of ion channel activity. Finally, we will use Sort-seq methods to link genotypes to phenotypes for each of the variants in our libraries. Through this method, it will be possible to link individual residue changes to functional outcomes in TRPM8.

Stabilisation of the RirA [4Fe-4S] cluster results in loss of iron-sensing function.

RirA is a global iron regulator in diverse Alphaproteobacteria that belongs to the Rrf2 superfamily of transcriptional regulators, which can contain an iron-sulfur (Fe-S) cluster. Under iron-replete conditions, RirA contains a [4Fe-4S] cluster, enabling high-affinity binding to RirA-regulated operator sequences, thereby causing the repression of cellular iron uptake. Under iron deficiency, one of the cluster irons dissociates, generating an unstable [3Fe-4S] form that subsequently degrades to a [2Fe-2S] form and then to apo RirA, resulting in loss of high-affinity DNA-binding. The cluster is coordinated by three conserved cysteine residues and an unknown fourth ligand. Considering the lability of one of the irons and the resulting cluster fragility, we hypothesized that the fourth ligand may not be an amino acid residue. To investigate this, we considered that the introduction of an amino acid residue that could coordinate the cluster might stabilize it. A structural model of RirA, based on the Rrf2 family nitrosative stress response regulator NsrR, highlighted residue 8, an Asn in the RirA sequence, as being appropriately positioned to coordinate the cluster. Substitution of Asn8 with Asp, the equivalent, cluster-coordinating residue of NsrR, or with Cys, resulted in proteins that contained a [4Fe-4S] cluster, with N8D RirA exhibiting spectroscopic properties very similar to NsrR. The variant proteins retained the ability to bind RirA-regulated DNA, and could still act as repressors of RirA-regulated genes in vivo. However, they were significantly more stable than wild-type RirA when exposed to O2 and/or low iron. Importantly, they exhibited reduced capacity to respond to cellular iron levels, even abolished in the case of the N8D version, and thus were no longer iron sensing. This work demonstrates the importance of cluster fragility for the iron-sensing function of RirA, and more broadly, how a single residue substitution can alter cluster coordination and functional properties in the Rrf2 superfamily of regulators.

Computational analyses reveal fundamental properties of the AT structure related to thrombosis.

Blood coagulation is a vital process for humans and other species. Following an injury to a blood vessel, a cascade of molecular signals is transmitted, inhibiting and activating more than a dozen coagulation factors and resulting in the formation of a fibrin clot that ceases the bleeding. In this process, antithrombin (AT), encoded by the SERPINC1 gene is a key player regulating the clotting activity and ensuring that it stops at the right time. In this sense, mutations to this factor often result in thrombosis-the excessive coagulation that leads to the potentially fatal formation of blood clots that obstruct veins. Although this process is well known, it is still unclear why even single residue substitutions to AT lead to drastically different phenotypes. In this study, to understand the effect of mutations throughout the AT structure, we created a detailed network map of this protein, where each node is an amino acid, and two amino acids are connected if they are in close proximity in the three-dimensional structure. With this simple and intuitive representation and a machine-learning framework trained using genetic information from more than 130 patients, we found that different types of thrombosis have emerging patterns that are readily identifiable. Together, these results demonstrate how clinical features, genetic data and in silico analysis are converging to enhance the diagnosis and treatment of coagulation disorders. Supplementary data are available at Bioinformatics Advances online.

Structural consequences of sequence variation in mammalian prion β2α2 loop segments.

Sequence variation in the β2α2 loop, residues 165-175 of the mammalian prion protein (PrP), influences its structure. To better understand the consequences of sequence variation in this region of the protein, we biochemically and biophysically interrogate natural and artificial sequence variants of the β2α2 loop of mammalian PrP. Using microcrystal electron diffraction (MicroED), we determine atomic resolution structures of segments encompassing residues 168-176 from the β2α2 loop of PrP with sequences corresponding to human, mouse/cow, bank vole/hamster, rabbit/pig/guinea pig, and naked mole rat (elk-T174S) β2α2 loops, as well as synthetic β2α2 loop sequences. This collection of structures presents two dominant amyloid packing polymorphisms. In the first polymorph, denoted "clasped", side chains within a sheet form polar clasps by facing each other on the same strand, exemplified by the mouse/cow, human, and bank vole/hamster sequences. Because its stability is derived from within a strand and through polar ladders within a sheet, the sequence requirements for the mating strand are less restrictive. A second polymorph, denoted "interdigitated," has sidechains interdigitate across mating sheets, exemplified by the elk, naked mole rat (elk T174S), and rabbit sequences. The two types of packing present distinct networks of stabilizing hydrogen bonds. The identity of residue 174 appears to strongly influence the packing adopted in these peptides, but consideration of the overall sequence of a given segment is needed to understand the stability of its assemblies. Incorporation of these β2α2 loop sequences into an 85 residue recombinant segment encoding wild-type bank vole PrP94-178 demonstrates that even single residue substitutions could impact fibril morphology as evaluated by negative stain electron microscopy. This is in line with recent findings supporting the accessibility of different structural geometries by varied mammalian prion sequences, and indicates that sequence-specific polymorphisms may be influenced by residues in the β2α2 loop.

Role of the reaction-structure coupling in temperature compensation of the KaiABC circadian rhythm.

When the mixture solution of cyanobacterial proteins, KaiA, KaiB, and KaiC, is incubated with ATP in vitro, the phosphorylation level of KaiC shows stable oscillations with the temperature-compensated circadian period. Elucidating this temperature compensation is essential for understanding the KaiABC circadian clock, but its mechanism has remained a mystery. We analyzed the KaiABC temperature compensation by developing a theoretical model describing the feedback relations among reactions and structural transitions in the KaiC molecule. The model showed that the reduced structural cooperativity should weaken the negative feedback coupling among reactions and structural transitions, which enlarges the oscillation amplitude and period, explaining the observed significant period extension upon single amino-acid residue substitution. We propose that an increase in thermal fluctuations similarly attenuates the reaction-structure feedback, explaining the temperature compensation in the KaiABC clock. The model explained the experimentally observed responses of the oscillation phase to the temperature shift or the ADP-concentration change and suggested that the ATPase reactions in the CI domain of KaiC affect the period depending on how the reaction rates are modulated. The KaiABC clock provides a unique opportunity to analyze how the reaction-structure coupling regulates the system-level synchronized oscillations of molecules.

Active site remodelling of a cyclodipeptide synthase redefines substrate scope

Cyclodipeptide synthases (CDPSs) generate a wide range of cyclic dipeptides using aminoacylated tRNAs as substrates. Histidine-containing cyclic dipeptides have important biological activities as anticancer and neuroprotective molecules. Out of the 120 experimentally validated CDPS members, only two are known to accept histidine as a substrate yielding cyclo(His-Phe) and cyclo(His-Pro) as products. It is not fully understood how CDPSs select their substrates, and we must rely on bioprospecting to find new enzymes and novel bioactive cyclic dipeptides. Here, we developed an in vitro system to generate an extensive library of molecules using canonical and non-canonical amino acids as substrates, expanding the chemical space of histidine-containing cyclic dipeptide analogues. To investigate substrate selection we determined the structure of a cyclo(His-Pro)-producing CDPS. Three consecutive generations harbouring single, double and triple residue substitutions elucidated the histidine selection mechanism. Moreover, substrate selection was redefined, yielding enzyme variants that became capable of utilising phenylalanine and leucine. Our work successfully engineered a CDPS to yield different products, paving the way to direct the promiscuity of these enzymes to produce molecules of our choosing.

Solubility-Aware Protein Binding Peptide Design Using AlphaFold.

New protein–protein interactions (PPIs) are identified, but PPIs have different physicochemical properties compared with conventional targets, making it difficult to use small molecules. Peptides offer a new modality to target PPIs, but designing appropriate peptide sequences by computation is challenging. Recently, AlphaFold and RoseTTAFold have made it possible to predict protein structures from amino acid sequences with ultra-high accuracy, enabling de novo protein design. We designed peptides likely to have PPI as the target protein using the “binder hallucination” protocol of AfDesign, a de novo protein design method using AlphaFold. However, the solubility of the peptides tended to be low. Therefore, we designed a solubility loss function using solubility indices for amino acids and developed a solubility-aware AfDesign binder hallucination protocol. The peptide solubility in sequences designed using the new protocol increased with the weight of the solubility loss function; moreover, they captured the characteristics of the solubility indices. Moreover, the new protocol sequences tended to have higher affinity than random or single residue substitution sequences when evaluated by docking binding affinity. Our approach shows that it is possible to design peptide sequences that can bind to the interface of PPI while controlling solubility.

A So-Far Overlooked Secondary Conformation State in the Binding Mode of SARS-CoV-2 Spike Protein to Human ACE2 and Its Conversion Rate Are Crucial for Estimating Infectivity Efficacy of the Underlying Virus Variant.

ABSTRACTSince its outbreak in 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread with high transmission efficiency across the world, putting health care as well as economic systems under pressure. During the course of the pandemic, the originally identified SARS-CoV-2 variant has been multiple times replaced by various mutant versions, which showed enhanced fitness due to increased infection and transmission rates. In order to find an explanation for why SARS-CoV-2 and its emerging mutated versions showed enhanced transmission efficiency compared with SARS-CoV (2002), an enhanced binding affinity of the spike protein to human angiotensin converting enzyme 2 (hACE2) has been proposed by crystal structure analysis and was identified in cell culture models. Kinetic analysis of the interaction of various spike protein constructs with hACE2 was considered to be best described by a Langmuir-based 1:1 stoichiometric interaction. However, we demonstrate in this report that the SARS-CoV-2 spike protein interaction with hACE2 is best described by a two-step interaction, which is defined by an initial binding event followed by a slower secondary rate transition that enhances the stability of the complex by a factor of ~190 (primary versus secondary state) with an overall equilibrium dissociation constant (KD) of 0.20 nM. In addition, we show that the secondary rate transition is not only present in SARS-CoV-2 wild type (“wt”; Wuhan strain) but also found in the B.1.1.7 variant, where its transition rate is 5-fold increased.IMPORTANCE The current SARS-CoV-2 pandemic is characterized by the high infectivity of SARS-CoV-2 and its derived variants of concern (VOCs). It has been widely assumed that the reason for its increased cell entry compared with SARS-CoV (2002) is due to alterations in the viral spike protein, where single amino acid residue substitutions can increase affinity for hACE2. So far, the interaction of a single unit of the CoV-2 spike protein has been described using the 1:1 Langmuir interaction kinetic. However, we demonstrate here that there is a secondary state binding step that may be essential for novel VOCs in order to further increase their infectivity. These findings are important for quantitatively understanding the infection process of SARS-CoV-2 and characterization of emerging SARS-CoV-2 variants of spike proteins. Thus, they provide a tool for predicting the potential infectivity of the respective viral variants based on secondary rate transition and secondary complex stability.

Conformational spread drives the evolution of the calcium–calmodulin protein kinase II

The calcium calmodulin (Ca2+/CaM) dependent protein kinase II (CaMKII) decodes Ca2+ frequency oscillations. The CaMKIIα isoform is predominantly expressed in the brain and has a central role in learning. I matched residue and organismal evolution with collective motions deduced from the atomic structure of the human CaMKIIα holoenzyme to learn how its ring architecture abets function. Protein dynamic simulations showed its peripheral kinase domains (KDs) are conformationally coupled via lateral spread along the central hub. The underlying β-sheet motions in the hub or association domain (AD) were deconvolved into dynamic couplings based on mutual information. They mapped onto a coevolved residue network to partition the AD into two distinct sectors. A second, energetically stressed sector was added to ancient bacterial enzyme dimers for assembly of the ringed hub. The continued evolution of the holoenzyme after AD–KD fusion targeted the sector’s ring contacts coupled to the KD. Among isoforms, the α isoform emerged last and, it alone, mutated rapidly after the poikilotherm–homeotherm jump to match the evolution of memory. The correlation between dynamics and evolution of the CaMKII AD argues single residue substitutions fine-tune hub conformational spread. The fine-tuning could increase CaMKIIα Ca2+ frequency response range for complex learning functions.

ProteinsPlus: a comprehensive collection of web-based molecular modeling tools.

Upon the ever-increasing number of publicly available experimentally determined and predicted protein and nucleic acid structures, the demand for easy-to-use tools to investigate these structural models is higher than ever before. The ProteinsPlus web server (https://proteins.plus) comprises a growing collection of molecular modeling tools focusing on protein–ligand interactions. It enables quick access to structural investigations ranging from structure analytics and search methods to molecular docking. It is by now well-established in the community and constantly extended. The server gives easy access not only to experts but also to students and occasional users from the field of life sciences. Here, we describe its recently added new features and tools, beyond them a novel method for on-the-fly molecular docking and a search method for single-residue substitutions in local regions of a protein structure throughout the whole Protein Data Bank. Finally, we provide a glimpse into new avenues for the annotation of AlphaFold structures which are directly accessible via a RESTful service on the ProteinsPlus web server.