A helix propensity scale based on experimental studies of peptides and proteins

Elements of the C-terminal tail of a C-terminal domain homolog of the Orange Carotenoid Protein determining xanthophyll uptake from liposomes

Carotenoids perform multifaceted roles in life ranging from coloration over light harvesting to photoprotection. The Orange Carotenoid Protein (OCP), a light-driven photoswitch involved in cyanobacterial photoprotection, accommodates a ketocarotenoid vital for its function. OCP extracts its ketocarotenoid directly from membranes, or accepts it from homologs of its C-terminal domain (CTDH). The CTDH from Anabaena (AnaCTDH) was shown to be important for carotenoid transfer and delivery from/to membranes. The C-terminal tail of AnaCTDH is a critical structural element likely serving as a gatekeeper and facilitator of carotenoid uptake from membranes. We investigated the impact of amino acid substitutions within the AnaCTDH-CTT on echinenone and canthaxanthin uptake from DOPC and DMPG liposomes. The transfer rate was uniformly reduced for substitutions of Arg-137 and Arg-138 to Gln or Ala, and depended on the lipid type, indicating a weaker interaction particularly with the lipid head group. Our results further suggest that Glu-132 has a membrane-anchoring effect on the PC lipids, specifically at the choline motif as inferred from the strongly different effects of the CTT variants on the extraction from the two liposome types. The substitution of Pro-130 by Gly suggests that the CTT is perpendicular to both the membrane and the main AnaCTDH protein during carotenoid extraction. Finally, the simultaneous mutation of Leu-133, Leu-134 and Leu-136 for alanines showed that the hydrophobicity of the CTT is crucial for carotenoid uptake. Since some substitutions accelerated carotenoid transfer into AnaCTDH while others slowed it down, carotenoprotein properties can be engineered toward the requirements of applications.

A Novel Beta-Defensin Isoform from Malabar Trevally, Carangoides malabaricus (Bloch & Schneider, 1801), an Arsenal Against Fish Bacterial Pathogens: Molecular Characterization, Recombinant Production, and Mechanism of Action

Antimicrobial peptides (AMPs), including beta-defensin from fish, are a crucial class of peptide medicines. The focus of the current study is the molecular and functional attributes of CmDef, a 63-amino acid beta-defensin AMP from Malabar trevally, Carangoides malabaricus. This peptide demonstrated typical characteristics of AMPs, including hydrophobicity, amphipathic nature, and +2.8 net charge. The CmDef was recombinantly expressed and the recombinant peptide, rCmDef displayed a strong antimicrobial activity against bacterial fish pathogens with an MIC of 8 µM for V. proteolyticus and 32 µM for A. hydrophila. The E. tarda and V. harveyi showed an inhibition of 94% and 54%, respectively, at 32 µM concentration. No activity was observed against V. fluvialis and V. alginolyticus. The rCmDef has a multimode of action that exerts an antibacterial effect by membrane depolarization followed by membrane permeabilization and ROS production. rCmDef also exhibited anti-cancer activities in silico without causing hemolysis. The peptide demonstrated stability under various conditions, including different pH levels, temperatures, salts, and metal ions (KCl and CaCl2), and remained stable in the presence of proteases such as trypsin and proteinase K at concentrations up to 0.2 µg/100 µl. The strong antibacterial efficacy and non-cytotoxic nature suggest that rCmDef is a single-edged sword that can contribute significantly to aquaculture disease management.

Conformational Dynamics of the Nucleosomal Histone H2B Tails Revealed by Molecular Dynamics Simulations

Epigenetic modifications of histone N-terminal tails play a critical role in regulating the chromatin structure and biological processes such as transcription and DNA repair. One of the key post-translational modifications (PTMs) is the acetylation of lysine residues on histone tails. Epigenetic modifications are ubiquitous in the development of diseases, such as cancer and neurological disorders. Histone H2B tails are critical regulators of nucleosome dynamics, biological processes, and certain diseases. Here, we report all-atomistic molecular dynamics (MD) simulations of the nucleosome to demonstrate that acetylation of the histone tails changes their conformational space and interaction with DNA. We perform simulations of H2B tails, critical regulators of gene regulation, in both the lysine-acetylated (ACK) and unacetylated wild type (WT) states. To explore the effects of salt concentration, we use two different NaCl concentrations to perform simulations at microsecond time scales. Salt can modulate the effects of electrostatic interactions between the DNA phosphate backbone and histone tails. Upon acetylation, H2B tails shift their secondary structure helical propensity. The number of contacts between the DNA and the H2B tail decreases. We characterize the conformational dynamics of the H2B tails by principal component analysis (PCA). The ACK tails become more compact at increased salt concentrations, but conformations from the WT tails display the most contacts with DNA at both salt concentrations. Mainly, H2B acetylation may increase the DNA accessibility for regulatory proteins to bind, which can aid in gene regulation and NCP stability.

Employment of mastoparan-like peptides to prevent Staphylococcus aureus associated with bovine mastitis

Bovine mastitis is a frequent infection in lactating cattle, causing great economic losses. Staphylococcus aureus represents the main etiological agent, which causes recurrent and persistent intramammary infections because conventional antibiotics are ineffective against it. Mastoparan-like peptides are multifunctional molecules with broad antimicrobial potential, constituting an attractive alternative. Nevertheless, their toxicity to host cells has hindered their therapeutic application. Previously, our group engineered three mastoparan-L analogs, namely mastoparan-MO, mastoparan-R1, and [I5, R8] MP, to improve cell selectivity and potential. Here, we were interested in comparing the antibacterial efficacy of mastoparan-L and its analogs against bovine mastitis isolates of S. aureus strains, making a correlation with the physicochemical properties and structural arrangement changes promoted by the sequence modifications. As a result, the analog's hemolytic and/or antimicrobial activity was balanced. All the peptides displayed α-helical folding in hydrophobic and membrane-mimetic environments, as determined by circular dichroism. The peptide [I5, R8] MP stood out for its enhanced selectivity and antibacterial features related to mastoparan-L and the other derivatives. Biophysical approaches revealed that [I5, R8] MP rapidly depolarizes the bacterial membrane of S. aureus, causing cell death by subsequent membrane disruption. Our results demonstrated that the [I5, R8] MP peptide could be a starting point for the development of peptide-based drugs for the treatment of bovine mastitis, with the advantage of no residue in milk, which would help reduce the use of classical antibiotics.IMPORTANCEStaphylococcus aureus is a leading cause of mastitis, the world's most important dairy cattle disease. The multidrug resistance and zoonotic potential of S. aureus, besides the likelihood of antibiotic residues in milk, are of critical concern to public and animal health. Antimicrobial peptides offer a novel antimicrobial strategy. Here, we demonstrate that [I5, R8] MP is a potent and selective peptide, which acts on S. aureus by targeting the bacterial membrane. Therefore, understanding the physicochemical determinants and the modes of action of this class of antimicrobials opens novel prospects for peptide development with enhanced activities in the bovine mastitis context.

Effective Synthesis of High-Integrity mRNA Using In Vitro Transcription

mRNA vaccines are entering a period of rapid development. However, their synthesis is still plagued by challenges related to mRNA impurities and fragments (incomplete mRNA). Most impurities of mRNA products transcribed in vitro are mRNA fragments. Only full-length mRNA transcripts containing both a 5'-cap and a 3'-poly(A) structure are viable for in vivo expression. Therefore, RNA fragments are the primary product-related impurities that significantly hinder mRNA efficacy and must be effectively controlled; these species are believed to originate from either mRNA hydrolysis or premature transcriptional termination. In the manufacturing of commercial mRNA vaccines, T7 RNA polymerase-catalyzed in vitro transcription (IVT) synthesis is a well-established method for synthesizing long RNA transcripts. This study identified a pivotal domain on the T7 RNA polymerase that is associated with erroneous mRNA release. By leveraging the advantageous properties of a T7 RNA polymerase mutant and precisely optimized IVT process parameters, we successfully achieved an mRNA integrity exceeding 91%, thereby further unlocking the immense potential of mRNA therapeutics.

Analysis of AlphaMissense data in different protein groups and structural context

Single amino acid substitutions can profoundly affect protein folding, dynamics, and function. The ability to discern between benign and pathogenic substitutions is pivotal for therapeutic interventions and research directions. Given the limitations in experimental examination of these variants, AlphaMissense has emerged as a promising predictor of the pathogenicity of missense variants. Since heterogenous performance on different types of proteins can be expected, we assessed the efficacy of AlphaMissense across several protein groups (e.g. soluble, transmembrane, and mitochondrial proteins) and regions (e.g. intramembrane, membrane interacting, and high confidence AlphaFold segments) using ClinVar data for validation. Our comprehensive evaluation showed that AlphaMissense delivers outstanding performance, with MCC scores predominantly between 0.6 and 0.74. We observed low performance on disordered datasets and ClinVar data related to the CFTR ABC protein. However, a superior performance was shown when benchmarked against the high quality CFTR2 database. Our results with CFTR emphasizes AlphaMissense's potential in pinpointing functional hot spots, with its performance likely surpassing benchmarks calculated from ClinVar and ProteinGym datasets.

Multiscale Simulations of Self-Assembling Peptides: Surface and Core Hydrophobicity Determine Fibril Stability and Amyloid Aggregation

Assemblies of peptides and proteins through specific intermolecular interactions set the basis for macroscopic materials found in nature. Peptides provide easily tunable hydrogen-bonding interactions, which can lead to the formation of ordered structures such as highly stable β-sheets that can form amyloid-like supramolecular peptide nanofibrils (PNFs). PNFs are of special interest, as they could be considered as mimics of various fibrillar structures found in nature. In their ability to serve as supramolecular scaffolds, they could mimic certain features of the extracellular matrix to provide stability, interact with pathogens such as virions, and transduce signals between the outside and inside of cells. Many PNFs have been reported that reveal rich bioactivities. PNFs supporting neuronal cell growth or lentiviral gene transduction have been studied systematically, and their material properties were correlated to bioactivities. However, the impact of the structure of PNFs, their dynamics, and stabilities on their unique functions is still elusive. Herein, we provide a microscopic view of the self-assembled PNFs to unravel how the amino acid sequence of self-assembling peptides affects their secondary structure and dynamic properties of the peptides within supramolecular fibrils. Based on sequence truncation, amino acid substitution, and sequence reordering, we demonstrate that peptide-peptide aggregation propensity is critical to form bioactive β-sheet-rich structures. In contrast to previous studies, a very high peptide aggregation propensity reduces bioactivity due to intermolecular misalignment and instabilities that emerge when fibrils are in close proximity to other fibrils in solution. Our multiscale simulation approach correlates changes in biological activity back to single amino acid modifications. Understanding these relationships could lead to future material discoveries where the molecular sequence predictably determines the macroscopic properties and biological activity. In addition, our studies may provide new insights into naturally occurring amyloid fibrils in neurodegenerative diseases.

Generative artificial intelligence performs rudimentary structural biology modeling

Natural language-based generative artificial intelligence (AI) has become increasingly prevalent in scientific research. Intriguingly, capabilities of generative pre-trained transformer (GPT) language models beyond the scope of natural language tasks have recently been identified. Here we explored how GPT-4 might be able to perform rudimentary structural biology modeling. We prompted GPT-4 to model 3D structures for the 20 standard amino acids and an α-helical polypeptide chain, with the latter incorporating Wolfram mathematical computation. We also used GPT-4 to perform structural interaction analysis between nirmatrelvir and its target, the SARS-CoV-2 main protease. Geometric parameters of the generated structures typically approximated close to experimental references. However, modeling was sporadically error-prone and molecular complexity was not well tolerated. Interaction analysis further revealed the ability of GPT-4 to identify specific amino acid residues involved in ligand binding along with corresponding bond distances. Despite current limitations, we show the capacity of natural language generative AI to perform basic structural biology modeling and interaction analysis with atomic-scale accuracy.

Stability and dynamics of extradenticle modulates its function

Extradenticle (EXD) is a partner protein of the HOX transcription factors and plays an important role in the development of Drosophila. It confers increased affinity and specificity of DNA-binding to the HOX proteins. However, the DNA-binding homeodomain of EXD has a significantly weaker affinity to DNA compared to the HOX homeodomains. Here, we show that a glycine residue (G290) in the middle of the EXD DNA-binding helix primarily results in this weaker binding. Glycine destabilizes helices. To probe its role in the stability and function of the protein, G290 was mutated to alanine. The intrinsic stability of the DNA-binding helix increased in the G290A mutant as observed by NMR studies and molecular dynamics (MD) simulation. Also, NMR dynamics and MD simulation show that dynamic motions present in the wild-type protein are quenched in the mutant. This in turn resulted in increased stability of the entire homeodomain (ΔΔGG→A of -2.6 kcal/mol). Increased protein stability resulted in three-fold better DNA-binding affinity of the mutant as compared to the wild-type protein. Molecular mechanics with generalized Born and surface area solvation (MMGBSA) analysis of our MD simulation on DNA-bound models of both wild-type and mutant proteins shows that the contribution to binding is enhanced for most of the interface residues in the mutant compared to the wild-type. Interestingly, the flexible N-terminal arm makes more stable contact with the DNA minor groove in the mutant. We found that the two interaction sites i.e. the DNA-binding helix and the unstructured N-terminal arm influence each other via the bound DNA. These results provide an interesting conundrum: alanine at position 290 enhances both the stability and the DNA-binding affinity of the protein, however, evolution prefers glycine at this position. We have provided several plausible explanations for this apparent conundrum. The function of the EXD as a HOX co-factor requires its ability to discriminate similar DNA sequences, which is most likely comprom.

Acyl Capping Group Identity Effects on α-Helicity: On the Importance of Amide·Water Hydrogen Bonds to α-Helix Stability

Acyl capping groups stabilize α-helices relative to free N-termini by providing one additional C═Oi···Hi+4-N hydrogen bond. The electronic properties of acyl capping groups might also directly modulate α-helix stability: electron-rich N-terminal acyl groups could stabilize the α-helix by strengthening both i/i + 4 hydrogen bonds and i/i + 1 n → π* interactions. This hypothesis was tested in peptides X-AKAAAAKAAAAKAAGY-NH2, where X = different acyl groups. Surprisingly, the most electron-rich acyl groups (pivaloyl and iso-butyryl) strongly destabilized the α-helix. Moreover, the formyl group induced nearly identical α-helicity to that of the acetyl group, despite being a weaker electron donor for hydrogen bonds and for n → π* interactions. Other acyl groups exhibited intermediate α-helicity. These results indicate that the electronic properties of the acyl carbonyl do not directly determine the α-helicity in peptides in water. In order to understand these effects, DFT calculations were conducted on α-helical peptides. Using implicit solvation, α-helix stability correlated with acyl group electronics, with the pivaloyl group exhibiting closer hydrogen bonds and n → π* interactions, in contrast to the experimental results. However, DFT and MD calculations with explicit water solvation revealed that hydrogen bonding to water was impacted by the sterics of the acyl capping group. Formyl capping groups exhibited the closest water-amide hydrogen bonds, while pivaloyl groups exhibited the longest. In α-helices in the PDB, the highest frequency of close amide-water hydrogen bonds is observed when the N-cap residue is Gly. The combination of experimental and computational results indicates that solvation (hydrogen bonding of water) to the N-terminal amide groups is a central determinant of α-helix stability.

Systematic evaluation of lecithin:cholesterol acyltransferase binding sites in apolipoproteins via peptide based nanodiscs: regulatory role of charged residues at positions 4 and 7

Lecithin:cholesterol acyltransferase (LCAT) exhibits α-activity on high-density and β-activity on low-density lipoproteins. However, the molecular determinants governing LCAT activation by different apolipoproteins remain elusive. Uncovering these determinants would offer the opportunity to design and explore advanced therapies against dyslipidemias. Here, we have conducted coarse-grained and all-atom molecular dynamics simulations of LCAT with nanodiscs made with α-helical amphiphilic peptides either derived from apolipoproteins A1 and E (apoA1 and apoE) or apoA1 mimetic peptide 22A that was optimized to activate LCAT. This study aims to explore what drives the binding of peptides to our previously identified interaction site in LCAT. We hypothesized that this approach could be used to screen for binding sites of LCAT in different apolipoproteins and would provide insights to differently localized LCAT activities. Our screening approach was able to discriminate apoA1 helixes 4, 6, and 7 as key contributors to the interaction with LCAT supporting the previous research data. The simulations provided detailed molecular determinants driving the interaction with LCAT: the formation of hydrogen bonds or salt bridges between peptides E4 or D4 and LCAT S236 or K238 residues. Additionally, salt bridging between R7 and D73 was observed, depending on the availability of R7. Expanding our investigation to diverse plasma proteins, we detected novel LCAT binding helixes in apoL1, apoB100, and serum amyloid A. Our findings suggest that the same binding determinants, involving E4 or D4 -S236 and R7-D73 interactions, influence LCAT β-activity on low-density lipoproteins, where apoE and or apoB100 are hypothesized to interact with LCAT.

PAbFold: Linear Antibody Epitope Prediction using AlphaFold2

Defining the binding epitopes of antibodies is essential for understanding how they bind to their antigens and perform their molecular functions. However, while determining linear epitopes of monoclonal antibodies can be accomplished utilizing well-established empirical procedures, these approaches are generally labor- and time-intensive and costly. To take advantage of the recent advances in protein structure prediction algorithms available to the scientific community, we developed a calculation pipeline based on the localColabFold implementation of AlphaFold2 that can predict linear antibody epitopes by predicting the structure of the complex between antibody heavy and light chains and target peptide sequences derived from antigens. We found that this AlphaFold2 pipeline, which we call PAbFold, was able to accurately flag known epitope sequences for several well-known antibody targets (HA / Myc) when the target sequence was broken into small overlapping linear peptides and antibody complementarity determining regions (CDRs) were grafted onto several different antibody framework regions in the single-chain antibody fragment (scFv) format. To determine if this pipeline was able to identify the epitope of a novel antibody with no structural information publicly available, we determined the epitope of a novel anti-SARS-CoV-2 nucleocapsid targeted antibody using our method and then experimentally validated our computational results using peptide competition ELISA assays. These results indicate that the AlphaFold2-based PAbFold pipeline we developed is capable of accurately identifying linear antibody epitopes in a short time using just antibody and target protein sequences. This emergent capability of the method is sensitive to methodological details such as peptide length, AlphaFold2 neural network versions, and multiple-sequence alignment database. PAbFold is available at https://github.com/jbderoo/PAbFold.

Unraveling the underlying mechanisms of reduced amyloidogenic properties in human calcitonin via double mutations

The therapeutic efficacy of peptide-based drugs is commonly hampered by the intrinsic propensity to aggregation. A notable example is human calcitonin (hCT), a peptide hormone comprising 32 amino acids, which is synthesized and secreted by thyroid gland parafollicular cells (C cells). This hormone plays a vital role in regulating blood calcium levels and upholding bone integrity. Despite its physiological importance, utilizing hCT as a drug is hampered by its inclination to form amyloid. To address this limitation, an alternative is provided by salmon calcitonin (sCT), which possesses a lower aggregation propensity. Although sharing the same disulfide bond at the N terminus as hCT, sCT differs from hCT at a total of 16 amino acid positions. However, due to the dissimilarity in sequences, using sCT as a clinical replacement occasionally results in adverse side effects in patients. Earlier investigations have highlighted the significant roles of Tyr-12 and Asn-17 in inducing the formation of amyloid fibrils. By introducing double mutations at these sites, the ability to hinder aggregation can be significantly augmented. This study delves into the oligomerization and helical structure formation of the hCT double mutant (Y12LN17H hCT, noted as DM hCT), as well as two single mutants (Y12L and N17H), aiming to elucidate the mechanism behind hCT fibrillization. In addition, computational prediction tools were employed again to identify potential substitutes. Although the results yielded were not entirely satisfactory, a comparison between the newly examined and previously found hCT double mutants provides insights into the reduced aggregation propensity of the latter. This research endeavor holds the promise of informing the design of more effective therapeutic peptide drugs in the future.

The N-terminal intrinsically disordered region of Ncb5or docks with the cytochrome b5 core to form a helical motif that is of ancient origin

NADH cytochrome b5 oxidoreductase (Ncb5or) is a cytosolic ferric reductase implicated in diabetes and neurological conditions. Ncb5or comprises cytochrome b5 (b5 ) and cytochrome b5 reductase (b5 R) domains separated by a CHORD-Sgt1 (CS) linker domain. Ncb5or redox activity depends on proper inter-domain interactions to mediate electron transfer from NADH or NADPH via FAD to heme. While full-length human Ncb5or has proven resistant to crystallization, we have succeeded in obtaining high-resolution atomic structures of the b5 domain and a construct containing the CS and b5 R domains (CS/b5 R). Ncb5or also contains an N-terminal intrinsically disordered region of 50 residues that has no homologs in other protein families in animals but features a distinctive, conserved L34 MDWIRL40 motif also present in reduced lateral root formation (RLF) protein in rice and increased recombination center 21 in baker's yeast, all attaching to a b5 domain. After unsuccessful attempts at crystallizing a human Ncb5or construct comprising the N-terminal region naturally fused to the b5 domain, we were able to obtain a high-resolution atomic structure of a recombinant rice RLF construct corresponding to residues 25-129 of human Ncb5or (52% sequence identity; 74% similarity). The structure reveals Trp120 (corresponding to invariant Trp37 in Ncb5or) to be part of an 11-residue α-helix (S116 QMDWLKLTRT126 ) packing against two of the four helices in the b5 domain that surround heme (α2 and α5). The Trp120 side chain forms a network of interactions with the side chains of four highly conserved residues corresponding to Tyr85 and Tyr88 (α2), Cys124 (α5), and Leu47 in Ncb5or. Circular dichroism measurements of human Ncb5or fragments further support a key role of Trp37 in nucleating the formation of the N-terminal helix, whose location in the N/b5 module suggests a role in regulating the function of this multi-domain redox enzyme. This study revealed for the first time an ancient origin of a helical motif in the N/b5 module as reflected by its existence in a class of cytochrome b5 proteins from three kingdoms among eukaryotes.

Rational Design of TDP-43 Derived α-Helical Peptide Inhibitors: An In Silico Strategy to Prevent TDP-43 Aggregation in Neurodegenerative Disorders

TDP-43, an essential RNA/DNA-binding protein, is central to the pathology of neurodegenerative diseases, such as amyotrophic lateral sclerosis and frontotemporal dementia. Pathological mislocalization and aggregation of TDP-43 disrupt RNA splicing, mRNA stability, and mRNA transport, thereby impairing neuronal function and survival. The formation of amyloid-like TDP-43 filaments is largely facilitated by the destabilization of an α-helical segment within the disordered C-terminal region. In this study, we hypothesized that preventing the destabilization of the α-helical domain could potentially halt the growth of these pathological filaments. To explore this, we utilized a range of in silico techniques to design and evaluate peptide-based therapeutics that bind to pathological TDP-43 amyloid-like filament crystal structures and resist β sheet conversion. Our computational approaches, including biophysical and secondary structure property prediction, molecular docking, 3D structure prediction, and molecular dynamics simulations, were used to assess the structure, stability, and binding affinity of these peptides in relation to pathological TDP-43 filaments. The results of our in silico analyses identified a selection of promising peptides which displayed a stable α-helical structure, exhibited an increased number of intramolecular hydrogen bonds within the helical domain, and demonstrated high binding affinities for pathological TDP-43 amyloid-like filaments. Molecular dynamics simulations provided further support for the structural and thermodynamic stability of these peptides, as they exhibited lower root-mean-square deviation and more favorable free energy landscapes over 300 ns. These findings establish α-helical propensity peptides as potential lead molecules for the development of novel therapeutics against TDP-43 aggregation. This structure-based computational approach for the rational design of peptide inhibitors opens a new direction in the search for effective interventions for ALS, FTD, and other related neurodegenerative diseases. The peptides identified as the most promising candidates in this study are currently subject to further testing and validation through both in vitro and in vivo experiments.

Continuous evolution of compact protein degradation tags regulated by selective molecular glues

Conditional protein degradation tags (degrons) are usually >100 amino acids long or are triggered by small molecules with substantial off-target effects, thwarting their use as specific modulators of endogenous protein levels. We developed a phage-assisted continuous evolution platform for molecular glue complexes (MG-PACE) and evolved a 36-amino acid zinc finger (ZF) degron (SD40) that binds the ubiquitin ligase substrate receptor cereblon in complex with PT-179, an orthogonal thalidomide derivative. Endogenous proteins tagged in-frame with SD40 using prime editing are degraded by otherwise inert PT-179. Cryo-electron microscopy structures of SD40 in complex with ligand-bound cereblon revealed mechanistic insights into the molecular basis of SD40's activity and specificity. Our efforts establish a system for continuous evolution of molecular glue complexes and provide ZF tags that overcome shortcomings associated with existing degrons.

Structural insights into peptidoglycan glycosidase EtgA binding to the inner rod protein EscI of the type III secretion system via a designed EscI-EtgA fusion protein

Bacteria express lytic enzymes such as glycosidases, which have potentially self-destructive peptidoglycan (PG)-degrading activity and, therefore, require careful regulation in bacteria. The PG glycosidase EtgA is regulated by localization to the assembling type III secretion system (T3SS), generating a hole in the PG layer for the T3SS to reach the outer membrane. The EtgA localization was found to be mediated via EtgA interacting with the T3SS inner rod protein EscI. To gain structural insights into the EtgA recognition of EscI, we determined the 2.01 Å resolution structure of an EscI (51-87)-linker-EtgA fusion protein designed based on AlphaFold2 predictions. The structure revealed EscI residues 72-87 forming an α-helix interacting with the backside of EtgA, distant from the active site. EscI residues 56-71 also were found to interact with EtgA, with these residues stretching across the EtgA surface. The ability of the EscI to interact with EtgA was also probed using an EscI peptide. The EscI peptide comprising residues 66-87, slightly larger than the observed EscI α-helix, was shown to bind to EtgA using microscale thermophoresis and thermal shift differential scanning fluorimetry. The EscI peptide also had a two-fold activity-enhancing effect on EtgA, whereas the EscI-EtgA fusion protein enhanced activity over four-fold compared to EtgA. Our studies suggest that EtgA regulation by EscI could be trifold involving protein localization, protein activation, and protein stabilization components. Analysis of the sequence conservation of the EscI EtgA interface residues suggested a possible conservation of such regulation for related proteins from different bacteria.

Dual mechanisms contribute to enhanced voltage dependence of an electric fish potassium channel

The voltage dependence of different voltage-gated potassium channels, described by the voltage at which half of the channels are open (V1/2), varies over a range of 80 mV and is influenced by factors such as the number of positive gating charges and the identity of the hydrophobic amino acids in the channel's voltage sensor (S4). Here we explore by experimental manipulations and molecular dynamics simulation the contributions of two derived features of an electric fish potassium channel (Kv1.7a) that is among the most voltage-sensitive Shaker family potassium channels known. These are a patch of four contiguous negatively charged glutamates in the S3-S4 extracellular loop and a glutamate in the S3b helix. We find that these negative charges affect V1/2 by separate, complementary mechanisms. In the closed state, the S3-S4 linker negative patch reduces the membrane surface charge biasing the channel to enter the open state while, upon opening, the negative amino acid in the S3b helix faces the second (R2) gating charge of the voltage sensor electrostatically biasing the channel to remain in the open state. This work highlights two evolutionary novelties that illustrate the potential influence of negatively charged amino acids in extracellular loops and adjacent helices to voltage dependence.

Genome-wide identification and characterization of flowering genes in Citrus sinensis (L.) Osbeck: a comparison among C. Medica L., C. Reticulata Blanco, C. Grandis (L.) Osbeck and C. Clementina

BACKGROUND

Flowering plays an important role in completing the reproductive cycle of plants and obtaining next generation of plants. In case of citrus, it may take more than a year to achieve progeny. Therefore, in order to fasten the breeding processes, the juvenility period needs to be reduced. The juvenility in plants is regulated by set of various flowering genes. The citrus fruit and leaves possess various medicinal properties and are subjected to intensive breeding programs to produce hybrids with improved quality traits. In order to break juvenility in Citrus, it is important to study the role of flowering genes. The present study involved identification of genes regulating flowering in Citrus sinensis L. Osbeck via homology based approach. The structural and functional characterization of these genes would help in targeting genome editing techniques to induce mutations in these genes for producing desirable results.

RESULTS

A total of 43 genes were identified which were located on all the 9 chromosomes of citrus. The in-silico analysis was performed to determine the genetic structure, conserved motifs, cis-regulatory elements (CREs) and phylogenetic relationship of the genes. A total of 10 CREs responsible for flowering were detected in 33 genes and 8 conserved motifs were identified in all the genes. The protein structure, protein-protein interaction network and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was performed to study the functioning of these genes which revealed the involvement of flowering proteins in circadian rhythm pathways. The gene ontology (GO) and gene function analysis was performed to functionally annotate the genes. The structure of the genes and proteins were also compared among other Citrus species to study the evolutionary relationship among them. The expression study revealed the expression of flowering genes in floral buds and ovaries. The qRT-PCR analysis revealed that the flowering genes were highly expressed in bud stage, fully grown flower and early stage of fruit development.

CONCLUSIONS

The findings suggested that the flowering genes were highly conserved in citrus species. The qRT-PCR analysis revealed the tissue specific expression of flowering genes (CsFT, CsCO, CsSOC, CsAP, CsSEP and CsLFY) which would help in easy detection and targeting of genes through various forward and reverse genetic approaches.

An alpha-helical lid guides the target DNA toward catalysis in CRISPR-Cas12a

CRISPR-Cas12a is a powerful RNA-guided genome-editing system that generates double-strand DNA breaks using its single RuvC nuclease domain by a sequential mechanism in which initial cleavage of the non-target strand is followed by target strand cleavage. How the spatially distant DNA target strand traverses toward the RuvC catalytic core is presently not understood. Here, continuous tens of microsecond-long molecular dynamics and free-energy simulations reveal that an α-helical lid, located within the RuvC domain, plays a pivotal role in the traversal of the DNA target strand by anchoring the crRNA:target strand duplex and guiding the target strand toward the RuvC core, as also corroborated by DNA cleavage experiments. In this mechanism, the REC2 domain pushes the crRNA:target strand duplex toward the core of the enzyme, while the Nuc domain aids the bending and accommodation of the target strand within the RuvC core by bending inward. Understanding of this critical process underlying Cas12a activity will enrich fundamental knowledge and facilitate further engineering strategies for genome editing.

Leucine Motifs Stabilize Residual Helical Structure in Disordered Proteins

Many examples are known of regions of intrinsically disordered proteins that fold into α-helices upon binding to their targets. These helical binding motifs (HBMs) can be partially helical also in the unbound state, and this so-called residual structure can affect binding affinity and kinetics. To investigate the underlying mechanisms governing the formation of residual helical structure, we assembled a dataset of experimental helix contents of 65 peptides containing HBM that fold-upon-binding. The average residual helicity is 17% and increases to 60% upon target binding. The helix contents of residual and target-bound structures do not correlate, however the relative location of helix elements in both states shows a strong overlap. Compared to the general disordered regions, HBMs are enriched in amino acids with high helix preference and these residues are typically involved in target binding, explaining the overlap in helix positions. In particular, we find that leucine residues and leucine motifs in HBMs are the major contributors to helix stabilization and target-binding. For the two model peptides, we show that substitution of leucine motifs to other hydrophobic residues (valine or isoleucine) leads to reduction of residual helicity, supporting the role of leucine as helix stabilizer. From the three hydrophobic residues only leucine can efficiently stabilize residual helical structure. We suggest that the high occurrence of leucine motifs and a general preference for leucine at binding interfaces in HBMs can be explained by its unique ability to stabilize helical elements.

Phosphomimetic substitutions in TDP-43's transiently α-helical region suppress phase separation

Phosphorylated TAR DNA-binding protein of 43 kDa (TDP-43) is present within the aggregates of several age-related neurodegenerative disorders, such as amyotrophic lateral sclerosis, frontotemporal lobar degeneration, and Alzheimer's disease, to the point that the presence of phosphorylated TDP-43 is considered a hallmark of some of these diseases. The majority of known TDP-43 phosphorylation sites detected in amyotrophic lateral sclerosis and frontotemporal lobar degeneration patients is located in the low-complexity domain (LCD), the same domain that has been shown to be critical for TDP-43 liquid-liquid phase separation (LLPS). However, the effect of these LCD phosphorylation sites on TDP-43 LLPS has been largely unexplored, and any work that has been done has mainly focused on sites near the C-terminal end of the LCD. Here, we used a phosphomimetic approach to explore the impact of phosphorylation at residues S332 and S333, sites located within the transiently α-helical region of TDP-43 that have been observed to be phosphorylated in disease, on protein LLPS. Our turbidimetry and fluorescence microscopy data demonstrate that these phosphomimetic substitutions greatly suppress LLPS, and solution NMR data strongly suggest that this effect is at least in part due to the loss of α-helical propensity of the phosphomimetic protein variant. We also show that the S332D and S333D substitutions slow TDP-43 LCD droplet aging and fibrillation of the protein. Overall, these findings provide a biophysical basis for understanding the effect of phosphorylation within the transiently α-helical region of TDP-43 LCD on protein LLPS and fibrillation, suggesting that phosphorylation at residues 332 and 333 is not necessarily directly related to the pathogenic process.

The intrinsically disordered transcriptional activation domain of CIITA is functionally tuneable by single substitutions: An exception or a new paradigm?

During protein evolution, some amino acid substitutions modulate protein function ("tuneability"). In most proteins, the tuneable range is wide and can be sampled by a set of protein variants that each contains multiple amino acid substitutions. In other proteins, the full tuneable range can be accessed by a set of variants that each contains a single substitution. Indeed, in some globular proteins, the full tuneable range can be accessed by the set of site-saturating substitutions at an individual "rheostat" position. However, in proteins with intrinsically disordered regions (IDRs), most functional studies-which would also detect tuneability-used multiple substitutions or small deletions. In disordered transcriptional activation domains (ADs), studies with multiple substitutions led to the "acidic exposure" model, which does not anticipate the existence of rheostat positions. In the few studies that did assess effects of single substitutions on AD function, results were mixed: the ADs of two full-length transcription factors did not show tuneability, whereas a fragment of a third AD was tuneable by single substitutions. In this study, we tested tuneability in the AD of full-length human class II transactivator (CIITA). Sequence analyses and experiments showed that CIITA's AD is an IDR. Functional assays of singly-substituted AD variants showed that CIITA's function was highly tuneable, with outcomes not predicted by the acidic exposure model. Four tested positions showed rheostat behavior for transcriptional activation. Thus, tuneability of different IDRs can vary widely. Future studies are needed to illuminate the biophysical features that govern whether an IDR is tuneable by single substitutions.

Biological importance of arginine: A comprehensive review of the roles in structure, disorder, and functionality of peptides and proteins

Arginine shows Jekyll and Hyde behavior in several respects. It participates in protein folding via ionic and H-bonds and cation-pi interactions; the charge and hydrophobicity of its side chain make it a disorder-promoting amino acid. Its methylation in histones; RNA binding proteins; chaperones regulates several cellular processes. The arginine-centric modifications are important in oncogenesis and as biomarkers in several cardiovascular diseases. The cross-links involving arginine in collagen and cornea are involved in pathogenesis of tissues but have also been useful in tissue engineering and wound-dressing materials. Arginine is a part of active site of several enzymes such as GTPases, peroxidases, and sulfotransferases. Its metabolic importance is obvious as it is involved in production of urea, NO, ornithine and citrulline. It can form unusual functional structures such as molecular tweezers in vitro and sprockets which engage DNA chains as part of histones in vivo. It has been used in design of cell-penetrating peptides as drugs. Arginine has been used as an excipient in both solid and injectable drug formulations; its role in suppressing opalescence due to liquid-liquid phase separation is particularly very promising. It has been known as a suppressor of protein aggregation during protein refolding. It has proved its usefulness in protein bioseparation processes like ion-exchange, hydrophobic and affinity chromatographies. Arginine is an amino acid, whose importance in biological sciences and biotechnology continues to grow in diverse ways.

DNA-Mediated Peptide Assembly into Protein Mimics

The design of new protein structures is challenging due to their vast sequence space and the complexity of protein folding. Here, we report a new modular DNA-templated strategy to construct protein mimics. We achieve the spatial control of multiple peptide units by conjugation with DNA and hybridization to a branched DNA trimer template followed by covalent stapling of the preorganized peptides into a single unit. A library of protein mimics with different lengths, sequences, and heptad registers has been efficiently constructed. DNA-templated protein mimics show an α-helix or coiled-coil motif formation even when they are constructed from weakly interacting peptide units. Their attached DNA handles can be used to exert dynamic control over the protein mimics' secondary and tertiary structures. This modular strategy will facilitate the development of DNA-encoded protein libraries for the rapid discovery of new therapeutics, enzymes, and antibody mimics.

Inhibition of Aurora-A/N-Myc Protein-Protein Interaction Using Peptidomimetics: Understanding the Role of Peptide Cyclization

Using N-Myc61-89 as a starting template we showcase the systematic use of truncation and maleimide constraining to develop peptidomimetic inhibitors of the N-Myc/Aurora-A protein-protein interaction (PPI); a potential anticancer drug discovery target. The most promising of these - N-Myc73-94-N85C/G89C-mal - is shown to favour a more Aurora-A compliant binding ensemble in comparison to the linear wild-type sequence as observed through fluorescence anisotropy competition assays, circular dichroism (CD) and nuclear magnetic resonance (NMR) experiments. Further in silico investigation of this peptide in its Aurora-A bound state, by molecular dynamics (MD) simulations, imply (i) the bound conformation is more stable as a consequence of the constraint, which likely suppresses dissociation and (ii) the constraint may make further stabilizing interactions with the Aurora-A surface. Taken together this work unveils the first orthosteric N-Myc/Aurora-A inhibitor and provides useful insights on the biophysical properties and thus design of constrained peptides, an attractive therapeutic modality.

Assembling membraneless organelles from de novo designed proteins

Recent advances in de novo protein design have delivered a diversity of discrete de novo protein structures and complexes. A new challenge for the field is to use these designs directly in cells to intervene in biological processes and augment natural systems. The bottom-up design of self-assembled objects such as microcompartments and membraneless organelles is one such challenge. Here we describe the design of genetically encoded polypeptides that form membraneless organelles in Escherichia coli. To do this, we combine de novo α-helical sequences, intrinsically disordered linkers and client proteins in single-polypeptide constructs. We tailor the properties of the helical regions to shift protein assembly from arrested assemblies to dynamic condensates. The designs are characterized in cells and in vitro using biophysical methods and soft-matter physics. Finally, we use the designed polypeptide to co-compartmentalize a functional enzyme pair in E. coli, improving product formation close to the theoretical limit.

The Effect of Pulling and Twisting Forces on Chameleon Sequence Peptides

Chameleon sequences are amino acid sequences found in several distinct configurations in experiment. They challenge our understanding of the link between sequence and structure, and provide insight into structural competition in proteins. Here, we study the energy landscapes for three such sequences, and interrogate how pulling and twisting forces impact the available structural ensembles. Chameleon sequences do not necessarily exhibit multiple structural ensembles on a multifunnel energy landscape when we consider them in isolation. The application of even small forces leads to drastic changes in the energy landscapes. For pulling forces, we observe transitions from helical to extended structures in a very small span of forces. For twisting forces, the picture is much more complex, and highly dependent on the magnitude and handedness of the applied force as well as the reference angle for the twist. Depending on these parameters, more complex and more simplistic energy landscapes are observed alongside more and less diverse structural ensembles. The impact of even small forces is significant, confirming their likely role in folding events. In addition, small forces exerted by the remaining scaffold of a protein may be sufficient to lead to the adoption of a specific structural ensemble by a chameleon sequence.

Programmable de novo designed coiled coil-mediated phase separation in mammalian cells

Membraneless liquid compartments based on phase-separating biopolymers have been observed in diverse cell types and attributed to weak multivalent interactions predominantly based on intrinsically disordered domains. The design of liquid-liquid phase separated (LLPS) condensates based on de novo designed tunable modules that interact in a well-understood, controllable manner could improve our understanding of this phenomenon and enable the introduction of new features. Here we report the construction of CC-LLPS in mammalian cells, based on designed coiled-coil (CC) dimer-forming modules, where the stability of CC pairs, their number, linkers, and sequential arrangement govern the transition between diffuse, liquid and immobile condensates and are corroborated by coarse-grained molecular simulations. Through modular design, we achieve multiple coexisting condensates, chemical regulation of LLPS, condensate fusion, formation from either one or two polypeptide components or LLPS regulation by a third polypeptide chain. These findings provide further insights into the principles underlying LLPS formation and a design platform for controlling biological processes.

APP substrate ectodomain defines amyloid-β peptide length by restraining γ-secretase processivity and facilitating product release

Sequential proteolysis of the amyloid precursor protein (APP) by γ-secretases generates amyloid-β (Aβ) peptides and defines the proportion of short-to-long Aβ peptides, which is tightly connected to Alzheimer's disease (AD) pathogenesis. Here, we study the mechanism that controls substrate processing by γ-secretases and Aβ peptide length. We found that polar interactions established by the APPC99 ectodomain (ECD), involving but not limited to its juxtamembrane region, restrain both the extent and degree of γ-secretases processive cleavage by destabilizing enzyme-substrate interactions. We show that increasing hydrophobicity, via mutation or ligand binding, at APPC99 -ECD attenuates substrate-driven product release and rescues the effects of Alzheimer's disease-associated pathogenic γ-secretase and APP variants on Aβ length. In addition, our study reveals that APPC99 -ECD facilitates the paradoxical production of longer Aβs caused by some γ-secretase inhibitors, which act as high-affinity competitors of the substrate. These findings assign a pivotal role to the substrate ECD in the sequential proteolysis by γ-secretases and suggest it as a sweet spot for the potential design of APP-targeting compounds selectively promoting its processing by these enzymes.

Effects of Replacing a Central Glycine Residue in GLP-1 on Receptor Affinity and Signaling Profile

Agonists of the glucagon-like peptide-1 receptor (GLP-1R) are used to treat diabetes and obesity. Cryo-EM structures indicate that GLP-1 is completely α-helical when bound to the GLP-1R. The mature form of this hormone, GLP-1(7-36), contains a glycine residue near the center (Gly22). Since glycine has the second-lowest α-helix propensity among the proteinogenic α-amino acid residues, and Gly22 does not appear to make direct contact with the receptor, we were motivated to explore the impact on agonist activity of altering the α-helix propensity at this position. We examined GLP-1 analogues in which Gly22 was replaced with L-Ala, D-Ala, or β-amino acid residues with varying helix propensities. The results suggest that the receptor is reasonably tolerant of variations in helix propensity, and that the functional receptor-agonist complex may comprise a conformational spectrum rather than a single fixed structure.

In silico protein dynamics in the human cytoplasm: Partial folding, misfolding, fold switching, and non-native interactions

We examine the influence of cellular interactions in all-atom models of a section of the Homo sapiens cytoplasm on the early folding events of the three-helix bundle protein B (PB). While genetically engineered PB is known to fold in dilute water box simulations in three microseconds, the three initially unfolded PB copies in our two cytoplasm models using a similar force field did not reach the native state during 30-microsecond simulations. We did however capture the formation of all three helices in a compact native-like topology. Folding in vivo is delayed because intramolecular contact formation within PB is in direct competition with intermolecular contacts between PB and surrounding macromolecules. In extreme cases, intermolecular beta-sheets are formed. Interactions with other macromolecules are also observed to promote structure formation, for example when a PB helix in our simulations is shielded from solvent by macromolecular crowding. Sticking and crowding in our models initiate sampling of helix/sheet structural plasticity of PB. Relatedly, in past in vitro experiments, similar GA domains were shown to switch between two different folds. Finally, we also observed that stickiness between PB and the cellular environment can be modulated in our simulations through the reduction in protein hydrophobicity when we reversed PB back to the wild-type sequence. This study demonstrates that even fast-folding proteins can get stuck in non-native states in the cell, making them useful models for protein-chaperone interactions and early stages of aggregate formation relevant to cellular disease.

Identifying Selectivity Filters in Protein Biosensor for Ligand Screening

Specialized sensing mechanisms in bacteria enable the identification of cognate ligands with remarkable selectivity in highly xenobiotic-polluted environments where these ligands are utilized as energy sources. Here, via integrating all-atom computer simulation, biochemical assay, and isothermal titration calorimetry measurements, we determine the molecular basis of MopR, a phenol biosensor's complex selection process of ligand entry. Our results reveal a set of strategically placed selectivity filters along the ligand entry pathway of MopR. These filters act as checkpoints, screening diverse aromatic ligands at the protein surface based on their chemical features and sizes. Ligands meeting specific criteria are allowed to enter the sensing site in an orientation-dependent manner. Sequence and structural analyses demonstrate the conservation of this ligand entry mechanism across the sensor class, with individual amino acids along the selectivity filter path playing a critical role in ligand selection. Together, this investigation highlights the importance of interactions with the ligand entry pathway, in addition to interactions within the binding pocket, in achieving ligand selectivity in biological sensing. The findings enhance our understanding of ligand selectivity in bacterial phenol biosensors and provide insights for rational expansion of the biosensor repertoire, particularly for the biotechnologically relevant class of aromatic pollutants.

Pervasive, conserved secondary structure in highly charged protein regions

Understanding how protein sequences confer function remains a defining challenge in molecular biology. Two approaches have yielded enormous insight yet are often pursued separately: structure-based, where sequence-encoded structures mediate function, and disorder-based, where sequences dictate physicochemical and dynamical properties which determine function in the absence of stable structure. Here we study highly charged protein regions (>40% charged residues), which are routinely presumed to be disordered. Using recent advances in structure prediction and experimental structures, we show that roughly 40% of these regions form well-structured helices. Features often used to predict disorder-high charge density, low hydrophobicity, low sequence complexity, and evolutionarily varying length-are also compatible with solvated, variable-length helices. We show that a simple composition classifier predicts the existence of structure far better than well-established heuristics based on charge and hydropathy. We show that helical structure is more prevalent than previously appreciated in highly charged regions of diverse proteomes and characterize the conservation of highly charged regions. Our results underscore the importance of integrating, rather than choosing between, structure- and disorder-based approaches.

Ligation of Gold Nanoparticles with Self-Assembling, Coiled-Coil Peptides

The surface of gold nanoparticles (AuNPs) can be conjugated with a wide range of highly functional biomolecules. A common pitfall when utilizing AuNPs is their tendency to aggregate, especially when their surface is functionalized with ligands of low molecular weight (no steric repulsion) or ligands of neutral charge (no electrostatic repulsion). For biomedical applications, AuNPs that are colloidally stable are desirable because they have a high surface area and thus reactivity, resist sedimentation, and exhibit uniform optical properties. Here, we engineer the surface of AuNPs so that they remain stable when decorated with coiled-coil (CC) peptides while preserving the native polypeptide properties. We achieve this by using a neutral, mixed ligand layer composed of lipoic acid poly(ethylene glycol) and lipoic acid poly(ethylene glycol) maleimide to attach the CCs. Tuning the surface fraction of each component within the mixed ligand layer also allowed us to control the degree of AuNP labeling with CCs. We demonstrate the dynamic surface properties of these CC-AuNPs by performing a place-exchange reaction and their utility by designing an energy-transfer-based caspase-3 sensor. Overall, this study optimizes the surface chemistry of AuNPs to quantitatively present functional biomolecules while maintaining colloid stability.

Cholesterol-dependent amyloid β production: space for multifarious interactions between amyloid precursor protein, secretases, and cholesterol

Amyloid β is considered a key player in the development and progression of Alzheimer's disease (AD). Many studies investigating the effect of statins on lowering cholesterol suggest that there may be a link between cholesterol levels and AD pathology. Since cholesterol is one of the most abundant lipid molecules, especially in brain tissue, it affects most membrane-related processes, including the formation of the most dangerous form of amyloid β, Aβ42. The entire Aβ production system, which includes the amyloid precursor protein (APP), β-secretase, and the complex of γ-secretase, is highly dependent on membrane cholesterol content. Moreover, cholesterol can affect amyloidogenesis in many ways. Cholesterol influences the stability and activity of secretases, but also dictates their partitioning into specific cellular compartments and cholesterol-enriched lipid rafts, where the amyloidogenic machinery is predominantly localized. The most complicated relationships have been found in the interaction between cholesterol and APP, where cholesterol affects not only APP localization but also the precise character of APP dimerization and APP processing by γ-secretase, which is important for the production of Aβ of different lengths. In this review, we describe the intricate web of interdependence between cellular cholesterol levels, cholesterol membrane distribution, and cholesterol-dependent production of Aβ, the major player in AD.

The Research Progress of Bioactive Peptides Derived from Traditional Natural Products in China

Traditional natural products in China have a long history and a vast pharmacological repertoire that has garnered significant attention due to their safety and efficacy in disease prevention and treatment. Among the bioactive components of traditional natural products in China, bioactive peptides (BPs) are specific protein fragments that have beneficial effects on human health. Despite many of the traditional natural products in China ingredients being rich in protein, BPs have not received sufficient attention as a critical factor influencing overall therapeutic efficacy. Therefore, the purpose of this review is to provide a comprehensive summary of the current methodologies for the preparation, isolation, and identification of BPs from traditional natural products in China and to classify the functions of discovered BPs. Insights from this review are expected to facilitate the development of targeted drugs and functional foods derived from traditional natural products in China in the future.

Conserved class B GPCR activation by a biased intracellular agonist

Class B G-protein-coupled receptors (GPCRs), including glucagon-like peptide 1 receptor (GLP1R) and parathyroid hormone 1 receptor (PTH1R), are important drug targets1-5. Injectable peptide drugs targeting these receptors have been developed, but orally available small-molecule drugs remain under development6,7. Here we report the high-resolution structure of human PTH1R in complex with the stimulatory G protein (Gs) and a small-molecule agonist, PCO371, which reveals an unexpected binding mode of PCO371 at the cytoplasmic interface of PTH1R with Gs. The PCO371-binding site is totally different from all binding sites previously reported for small molecules or peptide ligands in GPCRs. The residues that make up the PCO371-binding pocket are conserved in class B GPCRs, and a single alteration in PTH2R and two residue alterations in GLP1R convert these receptors to respond to PCO371. Functional assays reveal that PCO371 is a G-protein-biased agonist that is defective in promoting PTH1R-mediated arrestin signalling. Together, these results uncover a distinct binding site for designing small-molecule agonists for PTH1R and possibly other members of the class B GPCRs and define a receptor conformation that is specific only for G-protein activation but not arrestin signalling. These insights should facilitate the design of distinct types of class B GPCR small-molecule agonist for various therapeutic indications.

IFF: Identifying key residues in intrinsically disordered regions of proteins using machine learning

Conserved residues in protein homolog sequence alignments are structurally or functionally important. For intrinsically disordered proteins or proteins with intrinsically disordered regions (IDRs), however, alignment often fails because they lack a steric structure to constrain evolution. Although sequences vary, the physicochemical features of IDRs may be preserved in maintaining function. Therefore, a method to retrieve common IDR features may help identify functionally important residues. We applied unsupervised contrastive learning to train a model with self-attention neuronal networks on human IDR orthologs. Parameters in the model were trained to match sequences in ortholog pairs but not in other IDRs. The trained model successfully identifies previously reported critical residues from experimental studies, especially those with an overall pattern (e.g., multiple aromatic residues or charged blocks) rather than short motifs. This predictive model can be used to identify potentially important residues in other proteins, improving our understanding of their functions. The trained model can be run directly from the Jupyter Notebook in the GitHub repository using Binder (mybinder.org). The only required input is the primary sequence. The training scripts are available on GitHub (https://github.com/allmwh/IFF). The training datasets have been deposited in an Open Science Framework repository (https://osf.io/jk29b).

SARS-CoV-2 inhibitory activity of a short peptide derived from internal fusion peptide of S2 subunit of spike glycoprotein

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has posed a great concern in human population. To fight coronavirus emergence, we have dissected the conserved amino acid region of the internal fusion peptide in the S2 subunit of Spike glycoprotein of SARS-CoV-2 to design new inhibitory peptides. Among the 11 overlapping peptides (9-23-mer), PN19, a 19-mer peptide, exhibited a powerful inhibitory activity against different SARS-CoV-2 clinical isolate variants in absence of cytotoxicity. The PN19 inhibitory activity was found to be dependent on conservation of the central Phe and C-terminal Tyr residues in the peptide sequence. Circular dichroism spectra of the active peptide exhibited an alpha-helix propensity, confirmed by secondary structure prediction analysis. The PN19 inhibitory activity, exerted in the first step of virus infection, was reduced after peptide adsorption treatment with virus-cell substrate during fusion interaction. Additionally, PN19 inhibitory activity was reduced by adding S2 membrane-proximal region derived peptides. PN19 showed binding ability to the S2 membrane proximal region derived peptides, confirmed by molecular modelling, playing a role in the mechanism of action. Collectively, these results confirm that the internal fusion peptide region is a good candidate on which develop peptidomimetic anti SARS-CoV-2 antivirals.

Disruptive mutations in the serotonin transporter associate serotonin dysfunction with treatment-resistant affective disorder

Affective or mood disorders are a leading cause of disability worldwide. The serotonergic system has been heavily implicated in the complex etiology and serves as a therapeutic target. The serotonin transporter (SERT) is a major regulator of serotonin neurotransmission, yet the disease-relevance of impaired SERT function remains unknown. Here, we present the first identification and functional characterization of disruptive coding SERT variants found in patients with psychiatric diseases. In a unique cohort of 144 patients characterized by treatment-resistant chronic affective disorders with a lifetime history of electroconvulsive therapy, we identified two previously uncharacterized coding SERT variants: SERT-N217S and SERT-A500T. Both variants were significantly enriched in the patient cohort compared to GnomAD (SERT-N217S: OR = 151, P = 0.0001 and SERT-A500T: OR = 1348, P = 0.0022) and ethnicity-matched healthy controls (SERT-N217S: OR ≥ 17.7, P ≤ 0.013 and SERT-A500T: OR = ∞, P = 0.029). Functional investigations revealed that the mutations exert distinct perturbations to SERT function, but their overall effects converge on a partial loss-of-function molecular phenotype. Thus, the SERT-A500T variant compromises the catalytic activity, while SERT-N217S disrupts proper glycosylation of SERT with a resulting dominant-negative trafficking deficiency. Moreover, we demonstrate that the trafficking deficiency of SERT-N217S is amenable to pharmacochaperoning by noribogaine. Collectively, our findings describe the first disease-associated loss-of-function SERT variants and implicate serotonergic disturbances arising from SERT dysfunction as a risk factor for chronic affective disorders.

Instantaneous Inactivation of Herpes Simplex Virus by Silicon Nitride Bioceramics

Hydrolytic reactions taking place at the surface of a silicon nitride (Si3N4) bioceramic were found to induce instantaneous inactivation of Human herpesvirus 1 (HHV-1, also known as Herpes simplex virus 1 or HSV-1). Si3N4 is a non-oxide ceramic compound with strong antibacterial and antiviral properties that has been proven safe for human cells. HSV-1 is a double-stranded DNA virus that infects a variety of host tissues through a lytic and latent cycle. Real-time reverse transcription (RT)-polymerase chain reaction (PCR) tests of HSV-1 DNA after instantaneous contact with Si3N4 showed that ammonia and its nitrogen radical byproducts, produced upon Si3N4 hydrolysis, directly reacted with viral proteins and fragmented the virus DNA, irreversibly damaging its structure. A comparison carried out upon testing HSV-1 against ZrO2 particles under identical experimental conditions showed a significantly weaker (but not null) antiviral effect, which was attributed to oxygen radical influence. The results of this study extend the effectiveness of Si3N4's antiviral properties beyond their previously proven efficacy against a large variety of single-stranded enveloped and non-enveloped RNA viruses. Possible applications include the development of antiviral creams or gels and oral rinses to exploit an extremely efficient, localized, and instantaneous viral reduction by means of a safe and more effective alternative to conventional antiviral creams. Upon incorporating a minor fraction of micrometric Si3N4 particles into polymeric matrices, antiherpetic devices could be fabricated, which would effectively impede viral reactivation and enable high local effectiveness for extended periods of time.

Structure-function-guided design of synthetic peptides with anti-infective activity derived from wasp venom

Antimicrobial peptides (AMPs) derived from natural toxins and venoms offer a promising alternative source of antibiotics. Here, through structure-function-guided design, we convert two natural AMPs derived from the venom of the solitary eumenine wasp Eumenes micado into α-helical AMPs with reduced toxicity that kill Gram-negative bacteria in vitro and in a preclinical mouse model. To identify the sequence determinants conferring antimicrobial activity, an alanine scan screen and strategic single lysine substitutions are made to the amino acid sequence of these natural peptides. These efforts yield a total of 34 synthetic derivatives, including alanine substituted and lysine-substituted sequences with stabilized α-helical structures and increased net positive charge. The resulting lead synthetic peptides kill the Gram-negative pathogens Escherichia coli and Pseudomonas aeruginosa (PAO1 and PA14) by rapidly permeabilizing both their outer and cytoplasmic membranes, exhibit anti-infective efficacy in a mouse model by reducing bacterial loads by up to three orders of magnitude, and do not readily select for bacterial resistance.

A novel computationally engineered collagenase reduces the force required for tooth extraction in an ex-situ porcine jaw model

The currently employed tooth extraction methods in dentistry involve mechanical disruption of the periodontal ligament fibers, leading to inevitable trauma to the bundle bone comprising the socket walls. In our previous work, we have shown that a recombinantly expressed truncated version of clostridial collagenase G (ColG) purified from Escherichia coli efficiently reduced the force needed for tooth extraction in an ex-situ porcine jaw model, when injected into the periodontal ligament. Considering that enhanced thermostability often leads to higher enzymatic activity and to set the basis for additional rounds of optimization, we used a computational protein design approach to generate an enzyme to be more thermostable while conserving the key catalytic residues. This process generated a novel collagenase (ColG-variant) harboring sixteen mutations compared to ColG, with a nearly 4℃ increase in melting temperature. Herein, we explored the potential of ColG-variant to further decrease the physical effort required for tooth delivery using our established ex-situ porcine jaw model. An average reduction of 11% was recorded in the force applied to extract roots of mandibular split first and second premolar teeth treated with ColG-variant, relative to those treated with ColG. Our results show for the first time the potential of engineering enzyme properties for dental medicine and further contribute to minimally invasive tooth extraction.

Antiviral Activity against SARS-CoV-2 of Conformationally Constrained Helical Peptides Derived from Angiotensin-Converting Enzyme 2

Despite the availability of vaccines, COVID-19 continues to be aggressive, especially in immunocompromised individuals. Therefore, the development of a specific therapeutic agent with antiviral activity against SARS-CoV-2 is necessary. The infection pathway starts when the receptor binding domain of the viral spike protein interacts with the angiotensin converting enzyme 2 (ACE2), which acts as a host receptor for the RBD expressed on the host cell surface. In this scenario, ACE2 analogs binding to the RBD and preventing the cell entry can be promising antiviral agents. Most of the ACE2 residues involved in the interaction belong to the α1 helix, more specifically to the minimal fragment ACE2(24-42). In order to increase the stability of the secondary structure and thus antiviral activity, we designed different triazole-stapled analogs, changing the position and the number of bridges. The peptide called P3, which has the triazole-containing bridge in the positions 36-40, showed promising antiviral activity at micromolar concentrations assessed by plaque reduction assay. On the other hand, the double-stapled peptide P4 lost the activity, showing that excessive rigidity disfavors the interaction with the RBD.

Probing the Conformational Preference to β-Strand during Peptide Self-Assembly

Alanine-rich tetrapeptides like A₃K dominantly exist as polyproline II helices in dilute aqueous solutions. However, during self-assembly, based on the free energy calculation in implicit solvent for various peptide conformations, only the peptides in the β-strand conformation can be packed closely. This necessitates the conformational transition to the β-strand commonly observed during peptide self-assembly such as in amyloid fibril formation. In fact, the closest interpeptide distance of 4.8 Å is consistent with the interstrand distance determined from the X-ray diffraction pattern of many amyloid fibrils. The position of free energy minimum obtained from implicit solvent calculation matches exactly with the explicit solvent simulation through umbrella sampling when the peptide conformations are restrained, demonstrating the applicability of the former for rapid screening of peptide configurations favorable for self-assembly. The barrier in the free energy profile in the presence of water arises out of the entropic restriction on the interstitial water molecules while satisfying the hydrogen bonding of both the peptides by forming water mediated hydrogen bond bridge. Further, the high energy barrier observed for the β-strand suggests that peptides initially tend to self-assemble in the polyproline II structure to mitigate the desolvation energy cost; the transition to the β-strand would happen only in the later stage after crossing the barrier. The umbrella sampling simulations with peptides allowed to change conformations, relative to each other, confirm the dynamic conformational transition during the course of the self-assembly supporting the "dock and lock" mechanism suggested for amyloid fibrillar growth.

Differential Responses of the GLP-1 and GLP-2 Receptors to N-Terminal Modification of a Dual Agonist

Class B1 G protein-coupled receptors (GPCRs), collectively, respond to a diverse repertoire of extracellular polypeptide agonists and transmit the encoded messages to cytosolic partners. To fulfill these tasks, these highly mobile receptors must interconvert among conformational states in response to agonists. We recently showed that conformational mobility in polypeptide agonists themselves plays a role in activation of one class B1 GPCR, the receptor for glucagon-like peptide-1 (GLP-1). Exchange between helical and nonhelical conformations near the N-termini of agonists bound to the GLP-1R was revealed to be critical for receptor activation. Here, we ask whether agonist conformational mobility plays a role in the activation of a related receptor, the GLP-2R. Using variants of the hormone GLP-2 and the designed clinical agonist glepaglutide (GLE), we find that the GLP-2R is quite tolerant of variations in α-helical propensity near the agonist N-terminus, which contrasts with signaling at the GLP-1R. A fully α-helical conformation of the bound agonist may be sufficient for GLP-2R signal transduction. GLE is a GLP-2R/GLP-1R dual agonist, and the GLE system therefore enables direct comparison of the responses of these two GPCRs to a single set of agonist variants. This comparison supports the conclusion that the GLP-1R and GLP-2R differ in their response to variations in helical propensity near the agonist N-terminus. The data offer a basis for development of new hormone analogues with distinctive and potentially useful activity profiles; for example, one of the GLE analogues is a potent agonist of the GLP-2R but also a potent antagonist of the GLP-1R, a novel form of polypharmacology.

Yeast Svf1 binds ceramides and contributes to sphingolipid metabolism at the ER cis-Golgi interface

Ceramides are essential precursors of complex sphingolipids and act as potent signaling molecules. Ceramides are synthesized in the endoplasmic reticulum (ER) and receive their head-groups in the Golgi apparatus, yielding complex sphingolipids (SPs). Transport of ceramides between the ER and the Golgi is executed by the essential ceramide transport protein (CERT) in mammalian cells. However, yeast cells lack a CERT homolog, and the mechanism of ER to Golgi ceramide transport remains largely elusive. Here, we identified a role for yeast Svf1 in ceramide transport between the ER and the Golgi. Svf1 is dynamically targeted to membranes via an N-terminal amphipathic helix (AH). Svf1 binds ceramide via a hydrophobic binding pocket that is located in between two lipocalin domains. We showed that Svf1 membrane-targeting is important to maintain flux of ceramides into complex SPs. Together, our results show that Svf1 is a ceramide binding protein that contributes to sphingolipid metabolism at Golgi compartments.

Uncovering the Contributions of Charge Regulation to the Stability of Single Alpha Helices

The single alpha helix (SAH) is a recurring motif in biology. The consensus sequence has a di-block architecture that includes repeats of four consecutive glutamate residues followed by four consecutive lysine residues. Measurements show that the overall helicity of sequences with consensus E4 K4 repeats is insensitive to a wide range of pH values. Here, we use the recently introduced q-canonical ensemble, which allows us to decouple measurements of charge state and conformation, to explain the observed insensitivity of SAH helicity to pH. We couple the outputs from separate measurements of charge and conformation with atomistic simulations to derive residue-specific quantifications of preferences for being in an alpha helix and for the ionizable residues to be charged vs. uncharged. We find a clear preference for accommodating uncharged Glu residues within internal positions of SAH-forming sequences. The stabilities of alpha helical conformations increase with the number of E4 K4 repeats and so do the numbers of accessible charge states that are compatible with forming conformations of high helical content. There is conformational buffering whereby charge state heterogeneity buffers against large-scale conformational changes thus making the overall helicity insensitive to large changes in pH. Further, the results clearly argue against a single, rod-like alpha helical conformation being the only or even dominant conformation in the ensembles of so-called SAH sequences.

An engineered T7 RNA polymerase that produces mRNA free of immunostimulatory byproducts

In vitro transcription (IVT) is a DNA-templated process for synthesizing long RNA transcripts, including messenger RNA (mRNA). For many research and commercial applications, IVT of mRNA is typically performed using bacteriophage T7 RNA polymerase (T7 RNAP) owing to its ability to produce full-length RNA transcripts with high fidelity; however, T7 RNAP can also produce immunostimulatory byproducts such as double-stranded RNA that can affect protein expression. Such byproducts require complex purification processes, using methods such as reversed-phase high-performance liquid chromatography, to yield safe and effective mRNA-based medicines. To minimize the need for downstream purification processes, we rationally and computationally engineered a double mutant of T7 RNAP that produces substantially less immunostimulatory RNA during IVT compared with wild-type T7 RNAP. The resulting mutant allows for a simplified production process with similar mRNA potency, lower immunostimulatory content and quicker manufacturing time compared with wild-type T7 RNAP. Herein, we describe the computational design and development of this improved T7 RNAP variant.