Comparative Protein Structure Modeling Using MODELLER

Molecular dynamics of mismatch detection-How MutS uses indirect readout to find errors in DNA

The mismatch repair protein MutS safeguards genomic integrity by finding and initiating repair of basepairing errors in DNA. Single-molecule studies show MutS diffusing on DNA, presumably scanning for mispaired/unpaired bases, and crystal structures show a characteristic "mismatch-recognition" complex with DNA enclosed within MutS and kinked at the site of error. But how MutS goes from scanning thousands of Watson-Crick basepairs to recognizing rare mismatches remains unanswered, largely because atomic-resolution data on the search process are lacking. Here, 10 μs all-atom molecular dynamics simulations of Thermus aquaticus MutS bound to homoduplex DNA and T-bulge DNA illuminate the structural dynamics underlying the search mechanism. MutS-DNA interactions constitute a multistep mechanism to check DNA over two helical turns for its 1) shape, through contacts with the sugar-phosphate backbone, 2) conformational flexibility, through bending/unbending engineered by large-scale motions of the clamp domain, and 3) local deformability, through basepair destabilizing contacts. Thus, MutS can localize a potential target by indirect readout due to lower energetic costs of bending mismatched DNA and identify a site that distorts easily due to weaker base stacking and pairing as a mismatch. The MutS signature Phe-X-Glu motif can then lock in the mismatch-recognition complex to initiate repair.

M, Katona G. Survivin prevents the polycomb repressor complex 2 from methylating histone 3 lysine 27

This study investigates the role of survivin in epigenetic control of gene transcription through interaction with the polycomb repressive complex 2 (PRC2). PRC2 is responsible for silencing gene expression by trimethylating lysine 27 on histone 3. We observed differential expression of PRC2 subunits in CD4+ T cells with varying levels of survivin expression, and ChIP-seq results indicated that survivin colocalizes with PRC2 along DNA. Inhibition of survivin resulted in a significant increase in H3K27 trimethylation, implying that survivin prevents PRC2 from functioning. Peptide microarray showed that survivin interacts with peptides from PRC2 subunits, and machine learning revealed that amino acid composition contains relevant information for predicting survivin interaction. NMR and BLI experiments supported the interaction of survivin with PRC2 subunit EZH2. Finally, protein-protein docking revealed that the survivin-EZH2 interaction interface overlaps with catalytic residues of EZH2, potentially inhibiting its H3K27 methylation activity. These findings suggest that survivin inhibits PRC2 function.

Histone H3 lysine 27 crotonylation mediates gene transcriptional repression in chromatin

Histone lysine acylation, including acetylation and crotonylation, plays a pivotal role in gene transcription in health and diseases. However, our understanding of histone lysine acylation has been limited to gene transcriptional activation. Here, we report that histone H3 lysine 27 crotonylation (H3K27cr) directs gene transcriptional repression rather than activation. Specifically, H3K27cr in chromatin is selectively recognized by the YEATS domain of GAS41 in complex with SIN3A-HDAC1 co-repressors. Proto-oncogenic transcription factor MYC recruits GAS41/SIN3A-HDAC1 complex to repress genes in chromatin, including cell-cycle inhibitor p21. GAS41 knockout or H3K27cr-binding depletion results in p21 de-repression, cell-cycle arrest, and tumor growth inhibition in mice, explaining a causal relationship between GAS41 and MYC gene amplification and p21 downregulation in colorectal cancer. Our study suggests that H3K27 crotonylation signifies a previously unrecognized, distinct chromatin state for gene transcriptional repression in contrast to H3K27 trimethylation for transcriptional silencing and H3K27 acetylation for transcriptional activation.

Coarse-grained molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO 3 -, CO 3 2- , Cl - , Na + , K + , NH 3 and H + necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl - /HCO 3 - exchanger) and NBCe1 (Na + -CO 3 2- cotransporter). Previous computational studies of the outward facing (OF) state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed multiple 50 µs coarse-grained molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1 and NDCBE (a Na + -CO 3 2- /Cl - exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine (POPC), phosphatidylethanolamine (POPE), phosphatidylserine (POPS), and sphingomyelin (POSM). The recently resolved inward-facing (IF) state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, POPC, and POSM in the three studied proteins and their potential roles in the SLC4 transport function, conformational transition and protein dimerization were discussed.

Statement of significance

The SLC4 protein family is involved in critical physiological processes like pH and blood pressure regulation and maintenance of ion homeostasis. Its members can be found in various tissues. A number of studies suggest possible lipid regulation of the SLC4 function. However, the protein-lipid interactions in the SLC4 family are still poorly understood. Here we make use of long coarse-grained molecular dynamics simulations to assess the protein-lipid interactions in three SLC4 proteins with different transport modes, AE1, NBCe1, and NDCBE. We identify putative lipid binding sites for several lipid types of potential mechanistic importance, discuss them in the framework of the known experimental data and provide a necessary basis for further studies on lipid regulation of SLC4 function.

Substrate binding-induced conformational transitions in the omega-3 fatty acid transporter MFSD2A

Major Facilitator Superfamily Domain containing 2 A (MFSD2A) is a transporter that is highly enriched at the blood-brain and blood-retinal barriers, where it mediates Na+-dependent uptake of ω-3 fatty acids in the form of lysolipids into the brain and eyes, respectively. Despite recent structural insights, it remains unclear how this process is initiated, and driven by Na+. Here, we perform Molecular Dynamics simulations which demonstrate that substrates enter outward facing MFSD2A from the outer leaflet of the membrane via lateral openings between transmembrane helices 5/8 and 2/11. The substrate headgroup enters first and engages in Na+ -bridged interactions with a conserved glutamic acid, while the tail is surrounded by hydrophobic residues. This binding mode is consistent with a "trap-and-flip" mechanism and triggers transition to an occluded conformation. Furthermore, using machine learning analysis, we identify key elements that enable these transitions. These results advance our molecular understanding of the MFSD2A transport cycle.

Whole exome sequencing identified five novel variants in CNTN2, CARS2, ARSA, and CLCN4 leading to epilepsy in consanguineous families

Introduction: Epilepsy is a group of neurological disorders characterized by recurring seizures and fits. The Epilepsy genes can be classified into four distinct groups, based on involvement of these genes in different pathways leading to Epilepsy as a phenotype. Genetically the disease has been associated with various pathways, leading to pure epilepsy-related disorders caused by CNTN2 variations, or involving physical or systemic issues along with epilepsy caused by CARS2 and ARSA, or developed by genes that are putatively involved in epilepsy lead by CLCN4 variations. Methods: In this study, five families of Pakistani origin (EP-01, EP-02, EP-04, EP-09, and EP-11) were included for molecular diagnosis. Results: Clinical presentations of these patients included neurological symptoms such as delayed development, seizures, regression, myoclonic epilepsy, progressive spastic tetraparesis, vision and hearing impairment, speech problems, muscle fibrillation, tremors, and cognitive decline. Whole exome sequencing in index patients and Sanger sequencing in all available individuals in each family identified four novel homozygous variants in genes CARS2: c.655G>A p.Ala219Thr (EP-01), ARSA: c.338T>C: p.Leu113Pro (EP-02), c.938G>T p.Arg313Leu (EP-11), CNTN2: c.1699G>T p.Glu567Ter (EP-04), and one novel hemizygous variant in gene CLCN4: c.2167C>T p.Arg723Trp (EP-09). Conclusion: To the best of our knowledge these variants were novel and had not been reported in familial epilepsy. These variants were absent in 200 ethnically matched healthy control chromosomes. Three dimensional protein analyses revealed drastic changes in the normal functions of the variant proteins. Furthermore, these variants were designated as "pathogenic" as per guidelines of American College of Medical Genetics 2015. Due to overlapping phenotypes, among the patients, clinical subtyping was not possible. However, whole exome sequencing successfully pinpointed the molecular diagnosis which could be helpful for better management of these patients. Therefore, we recommend that exome sequencing be performed as a first-line molecular diagnostic test in familial cases.

Neoantigen-specific CD8 T cells with high structural avidity preferentially reside in and eliminate tumors

The success of cancer immunotherapy depends in part on the strength of antigen recognition by T cells. Here, we characterize the T cell receptor (TCR) functional (antigen sensitivity) and structural (monomeric pMHC-TCR off-rates) avidities of 371 CD8 T cell clones specific for neoantigens, tumor-associated antigens (TAAs) or viral antigens isolated from tumors or blood of patients and healthy donors. T cells from tumors exhibit stronger functional and structural avidity than their blood counterparts. Relative to TAA, neoantigen-specific T cells are of higher structural avidity and, consistently, are preferentially detected in tumors. Effective tumor infiltration in mice models is associated with high structural avidity and CXCR3 expression. Based on TCR biophysicochemical properties, we derive and apply an in silico model predicting TCR structural avidity and validate the enrichment in high avidity T cells in patients' tumors. These observations indicate a direct relationship between neoantigen recognition, T cell functionality and tumor infiltration. These results delineate a rational approach to identify potent T cells for personalized cancer immunotherapy.

Gαs slow conformational transition upon GTP binding and a novel Gαs regulator

G proteins are major signaling partners for G protein-coupled receptors (GPCRs). Although stepwise structural changes during GPCR-G protein complex formation and guanosine diphosphate (GDP) release have been reported, no information is available with regard to guanosine triphosphate (GTP) binding. Here, we used a novel Bayesian integrative modeling framework that combines data from hydrogen-deuterium exchange mass spectrometry, tryptophan-induced fluorescence quenching, and metadynamics simulations to derive a kinetic model and atomic-level characterization of stepwise conformational changes incurred by the β2-adrenergic receptor (β2AR)-Gs complex after GDP release and GTP binding. Our data suggest rapid GTP binding and GTP-induced dissociation of Gαs from β2AR and Gβγ, as opposed to a slow closing of the Gαs α-helical domain (AHD). Yeast-two-hybrid screening using Gαs AHD as bait identified melanoma-associated antigen D2 (MAGE D2) as a novel AHD-binding protein, which was also shown to accelerate the GTP-induced closing of the Gαs AHD.

Phylogenetic and Evolutionary Analysis of Plant Small RNA 2'-O-Methyltransferase (HEN1) Protein Family

HUA ENHANCER 1 (HEN1) is a pivotal mediator in protecting sRNAs from 3'-end uridylation and 3' to 5' exonuclease-mediated degradation in plants. Here, we investigated the pattern of the HEN1 protein family evolutionary history and possible relationships in the plant lineages using protein sequence analyses and conserved motifs composition, functional domain identification, architecture, and phylogenetic tree reconstruction and evolutionary history inference. According to our results, HEN1 protein sequences bear several highly conserved motifs in plant species retained during the evolution from their ancestor. However, several motifs are present only in Gymnosperms and Angiosperms. A similar trend showed for their domain architecture. At the same time, phylogenetic analysis revealed the grouping of the HEN1 proteins in the three main super clads. In addition, the Neighbor-net network analysis result provides some nodes have multiple parents indicating a few conflicting signals in the data, which is not the consequence of sampling error, the effect of the selected model, or the estimation method. By reconciling the protein and species tree, we considered the gene duplications in several given species and found 170 duplication events in the evolution of HEN1 in the plant lineages. According to our analysis, the main HEN1 superclass mostly showed orthologous sequences that illustrate the vertically transmitting of HEN1 to the main lines. However, in both orthologous and paralogs, we predicted insignificant structural deviations. Our analysis implies that small local structural changes that occur continuously during the folds can moderate the changes created in the sequence. According to our results, we proposed a hypothetical model and evolutionary trajectory for the HEN1 protein family in the plant kingdom.

Rapid protein stability prediction using deep learning representations

Predicting the thermodynamic stability of proteins is a common and widely used step in protein engineering, and when elucidating the molecular mechanisms behind evolution and disease. Here, we present RaSP, a method for making rapid and accurate predictions of changes in protein stability by leveraging deep learning representations. RaSP performs on-par with biophysics-based methods and enables saturation mutagenesis stability predictions in less than a second per residue. We use RaSP to calculate ∼ 230 million stability changes for nearly all single amino acid changes in the human proteome, and examine variants observed in the human population. We find that variants that are common in the population are substantially depleted for severe destabilization, and that there are substantial differences between benign and pathogenic variants, highlighting the role of protein stability in genetic diseases. RaSP is freely available-including via a Web interface-and enables large-scale analyses of stability in experimental and predicted protein structures.

Protocol for predicting drug-resistant protein mutations to an ERK2 inhibitor using RESISTOR

Prospective predictions of drug-resistant protein mutants could improve the design of therapeutics less prone to resistance. Here, we describe RESISTOR, an algorithm that uses structure- and sequence-based criteria to predict resistance mutations. We demonstrate the process of using RESISTOR to predict ERK2 mutants likely to arise in melanoma ablating the efficacy of the ERK1/2 inhibitor SCH779284. RESISTOR is included in the free and open-source computational protein design software OSPREY. For complete details on the use and execution of this protocol, please refer to Guerin et al..1.

Deimmunization of flagellin adjuvant for clinical application

Flagellin is the cognate ligand for host pattern recognition receptors, toll-like receptor 5 (TLR5) in the cell surface, and NAIP5/NLRC4 inflammasome in the cytosol. TLR5-binding domain is located in D1 domain, where crucial amino acid sequences are conserved among diverse bacteria. The highly conserved C-terminal 35 amino acids of flagellin were proved to be responsible for the inflammasome activation by binding to NAIP5. D2/D3 domains, located in the central region and exposed to the outside surface of flagellar filament, are heterogeneous across bacterial species and highly immunogenic. Taking advantage of TLR5- and NLRC4-stimulating activities, flagellin has been actively developed as a vaccine adjuvant and immunotherapeutic. Because of its immunogenicity, there exist worries concerning diminished efficacy and possible reactogenicity after repeated administration. Deimmunization of flagellin derivatives while preserving the TLR5/NLRC4-mediated immunomodulatory activity should be the most reasonable option for clinical application. This review describes strategies and current achievements in flagellin deimmunization.

Intrinsically disordered region of talin's FERM domain functions as an initial PIP2 recognition site

Focal adhesions (FAs) mediate the interaction of the cytoskeleton with the extracellular matrix in a highly dynamic fashion. Talin is a central regulator, adaptor protein, and mechano-sensor of FA complexes. For recruitment and firm attachment at FAs, talin's N-terminal FERM domain binds to phosphatidylinositol 4,5-bisphosphate (PIP2)-enriched membranes. A newly published autoinhibitory structure of talin-1, where the known PIP2 interaction sites are covered up, lead us to hypothesize that a hitherto less examined loop insertion of the FERM domain acts as an additional and initial site of contact. We evaluated direct interactions of talin-1 with a PIP2 membrane by means of atomistic molecular dynamics simulations. We show that this unstructured, 33-residue-long loop strongly interacts with PIP2 and can facilitate further membrane contacts, including the canonical PIP2 interactions, by serving as a flexible membrane anchor. Under force as present at FAs, the extensible FERM loop ensures talin maintains membrane contacts when pulled away from the membrane by up to 7 nm. We identify key basic residues of the anchor mediating the highly dynamic talin-membrane interaction. Our results put forward an intrinsically disordered loop as a key and highly adaptable PIP2 recognition site of talin and potentially other PIP2-binding mechano-proteins.

In silico characterization of cysteine-stabilized αβ defensins from neglected unicellular microeukaryotes

BACKGROUND

The emergence of multi-resistant pathogens have increased dramatically in recent years, becoming a major public-health concern. Among other promising antimicrobial molecules with potential to assist in this worldwide struggle, cysteine-stabilized αβ (CS-αβ) defensins are attracting attention due their efficacy, stability, and broad spectrum against viruses, bacteria, fungi, and protists, including many known human pathogens.

RESULTS

Here, 23 genomes of ciliated protists were screened and two CS-αβ defensins with a likely antifungal activity were identified and characterized, using bioinformatics, from a culturable freshwater species, Laurentiella sp. (LsAMP-1 and LsAMP-2). Although any potential cellular ligand could be predicted for LsAMP-2; evidences from structural, molecular dynamics, and docking analyses suggest that LsAMP-1 may form stably associations with phosphatidylinositol 4,5-bisphosphates (PIP2), a phospholipid found on many eukaryotic cells, which could, in turn, represent an anchorage mechanism within plasma membrane of targeted cells.

CONCLUSION

These data stress that more biotechnology-oriented studies should be conducted on neglected protists, such ciliates, which could become valuable sources of novel bioactive molecules for therapeutic uses.

Structural modeling and analyses of genetic variations in the human XPC nucleotide excision repair protein

Xeroderma pigmentosum C (XPC) is a key initiator in the global genome nucleotide excision repair pathway in mammalian cells. Inherited mutations in the XPC gene can cause xeroderma pigmentosum (XP) cancer predisposition syndrome that dramatically increases the susceptibility to sunlight-induced cancers. Various genetic variants and mutations of the protein have been reported in cancer databases and literature. The current lack of a high-resolution 3-D structure of human XPC makes it difficult to assess the structural impact of the mutations/genetic variations. Using the available high-resolution crystal structure of its yeast ortholog, Rad4, we built a homology model of human XPC protein and compared it with a model generated by AlphaFold. The two models are largely consistent with each other in the structured domains. We have also assessed the degree of conservation for each residue using 966 sequences of XPC orthologs. Our structure- and sequence conservation-based assessments largely agree with the variant's impact on the protein's structural stability, computed by FoldX and SDM. Known XP missense mutations such as Y585C, W690S, and C771Y are consistently predicted to destabilize the protein's structure. Our analyses also reveal several highly conserved hydrophobic regions that are surface-exposed, which may indicate novel intermolecular interfaces that are yet to be characterized.Communicated by Ramaswamy H. Sarma.

Computational and artificial intelligence-based methods for antibody development

Due to their high target specificity and binding affinity, therapeutic antibodies are currently the largest class of biotherapeutics. The traditional largely empirical antibody development process is, while mature and robust, cumbersome and has significant limitations. Substantial recent advances in computational and artificial intelligence (AI) technologies are now starting to overcome many of these limitations and are increasingly integrated into development pipelines. Here, we provide an overview of AI methods relevant for antibody development, including databases, computational predictors of antibody properties and structure, and computational antibody design methods with an emphasis on machine learning (ML) models, and the design of complementarity-determining region (CDR) loops, antibody structural components critical for binding.

Structural basis of lysophosphatidylserine receptor GPR174 ligand recognition and activation

Lysophosphatidylserine (LysoPS) is a lipid mediator that induces multiple cellular responses through binding to GPR174. Here, we present the cryo-electron microscopy (cryo-EM) structure of LysoPS-bound human GPR174 in complex with G_s protein. The structure reveals a ligand recognition mode, including the negatively charged head group of LysoPS forms extensive polar interactions with surrounding key residues of the ligand binding pocket, and the L-serine moiety buries deeply into a positive charged cavity in the pocket. In addition, the structure unveils a partially open pocket on transmembrane domain helix (TM) 4 and 5 for a lateral entry of ligand. Finally, the structure reveals a G_s engaging mode featured by a deep insertion of a helix 5 (αH5) and extensive polar interactions between receptor and αH5. Taken together, the information revealed by our structural study provides a framework for understanding LysoPS signaling and a rational basis for designing LysoPS receptor-targeting drugs.

Characterization of the coupling mechanism of scorpion β-neurotoxins on the voltage-gated sodium channel hNav1.6

Scorpion β-neurotoxins represent a pharmacological group that affects voltage-gated sodium channels (Nav). Despite knowing the electrophysiological effect of these toxins on Nav channels, the molecular mechanism by which the union is carried out is still undetermined. In this study, computational techniques such as modeling, docking and molecular dynamics were used to elucidate the mechanism of interaction between scorpion β-neurotoxins using the neurotoxin nCssII and its recombinant variant CssII-RCR, which bind to the site-4, an extracellular receptor, of the human sodium channel hNav1.6. Different modes of interaction were observed for both toxins, where the main distinguishing feature was the interaction generated by the residue E15 on such site-4; that is, E15 in nCssII exhibits an interaction with the voltage-sensing domain II, and the same residue E15 of CssII-RCR exhibits an interaction with domain III. Despite this difference in interaction by E15, it is seen that both neurotoxins interact with similar regions of the voltage sensing domain such as the S3-S4 connecting loop (L834-E838) of the hNav1.6. Our simulations present a first approach to the mode of interaction of scorpion beta-neurotoxins in toxin-receptor complexes, being able to explain at the molecular level the phenomenon of voltage sensor entrapment generated by these toxins.Communicated by Ramaswamy H. Sarma.

Bacterial origins of thymidylate metabolism in Asgard archaea and Eukarya

Asgard archaea include the closest known archaeal relatives of eukaryotes. Here, we investigate the evolution and function of Asgard thymidylate synthases and other folate-dependent enzymes required for the biosynthesis of DNA, RNA, amino acids and vitamins, as well as syntrophic amino acid utilization. Phylogenies of Asgard folate-dependent enzymes are consistent with their horizontal transmission from various bacterial groups. We experimentally validate the functionality of thymidylate synthase ThyX of the cultured 'Candidatus Prometheoarchaeum syntrophicum'. The enzyme efficiently uses bacterial-like folates and is inhibited by mycobacterial ThyX inhibitors, even though the majority of experimentally tested archaea are known to use carbon carriers distinct from bacterial folates. Our phylogenetic analyses suggest that the eukaryotic thymidylate synthase, required for de novo DNA synthesis, is not closely related to archaeal enzymes and might have been transferred from bacteria to protoeukaryotes during eukaryogenesis. Altogether, our study suggests that the capacity of eukaryotic cells to duplicate their genetic material is a sum of archaeal (replisome) and bacterial (thymidylate synthase) characteristics. We also propose that recent prevalent lateral gene transfer from bacteria has markedly shaped the metabolism of Asgard archaea.

Discovery, Characterization, and Engineering of LvIC, an α4/4-Conotoxin That Selectively Blocks Rat α6/α3β4 Nicotinic Acetylcholine Receptors

α6β4 nicotinic acetylcholine receptors (nAChRs) are expressed in the central and peripheral nervous systems, but their functions are not fully understood, largely because of a lack of specific ligands. Here, we characterized a novel α-conotoxin, LvIC, and designed a series of analogues to probe structure-activity relationships at the α6β4 nAChR. The potency and selectivity of these conotoxins were tested using two-electrode voltage-clamp recording on nAChR subtypes expressed in Xenopus laevis oocytes. One of the analogues, [D1G,ΔQ14]LvIC, potently blocked α6/α3β4 nAChRs (α6/α3 is a chimera) with an IC₅₀ of 19 nM, with minimal activity at other nAChR subtypes, including the structurally similar α6/α3β2β3 and α3β4 subtypes. Using NMR, molecular docking, and receptor mutation, structure-activity relationships of [D1G,ΔQ14]LvIC at the α6/α3β4 nAChR were defined. It is a potent and specific antagonist of α6β4 nAChRs that could potentially serve as a novel molecular probe to explore α6β4 nAChR-related neurophysiological and pharmacological functions.

Early-Stage Oligomerization of Prion-like Polypeptides Reveals the Molecular Mechanism of Amyloid-Disrupting Capacity by Proline Residues

Proline cis/trans isomerization governs protein local conformational changes via its local mechanical rigidity. The amyloid-disrupting capacity of proline is widely acknowledged; however, the molecular mechanism is still not clear. To understand how proline residues in polypeptide chains influence amyloid propensity, we study several truncated sequences of the TDP-43 C-terminal region (287-322) and their triple proline variants (308PPP310). We use coarse-grained molecular simulation to study the time evolution of the process of aggregation in the early stages in an effective high-concentration condition (∼25 mM). This ensures the long time scales for protein association at laboratory concentrations. We use several experimentally determined structure templates as initial structures of monomer conformations. We carry out oligomer size analysis and cluster analysis, along with several structural measures, to characterize the size distributions of oligomers and their morphological/structural properties. We show that average oligomer size is not a good indicator of amyloid propensity. Structural order and/or morphological properties are better alternatives. We show that proline variants can efficiently maintain the formation of large "ordered" oligomers of shorter truncated sequences, i.e., 307-322. This "order" maintenance is weakened when using longer truncated sequences (i.e., 287-322), leading to the formation of "disordered" oligomers. From an energy trade-off perspective, if the entropic effect is weak (short sequence length), the shape-complementarity of proline variants effectively guides the oligomerization process to form "ordered" oligomer intermediates. This leads to a distinct aggregation pathway that promotes amyloid formation (on-pathway). Strong entropic effects (long sequence length), however, would cause the formation of "disordered" oligomers. This in turn will suppress amyloid formation (off-pathway). The proline shape-complementary effects provide a guided morphological restraint to facilitate the pathways of amyloid formation. Our study supports the importance of structure-based kinetic heterogeneity of prion-like sequence fragments in driving different aggregation pathways. This work sheds light on the role of morphological and structural order of early-stage oligomeric species in regulating amyloid-disrupting capacity by prolines.

Toxicological impacts and likely protein targets of bisphenol a in Paramecium caudatum

Bisphenol A (BPA) is a widely used plasticizer agent and a well-known ubiquitous endocrine disruptor, which is frequently associated with a series of reproductive, developmental, and transgenerational effects over wildlife, livestocks, and humans. Although extensive toxicological data is available for metazoans, the impact of BPA over unicellular eukaryotes, which represents a considerable proportion of eukaryotic diversity, remains largely overlooked. Here, we used acute end-point toxicological assay and an inverted virtual-screening (IVS) approach to evaluate cellular impairments infringed by BPA over the cosmopolitan ciliated protist, Paramecium caudatum. Our data indicate a clear time-dependent effect over P. caudatum survival, which seems to be a consequence of disruptions to multiple core cellular functions, such as DNA and cell replication, transcription, translation and signaling pathways. Finally, the use of this ciliate as a biosensor to monitor BPA within environments and the relevance of bioinformatic methods to leverage our current knowledge on the impacts of emerging contaminants to biological systems are discussed.

Adrenomedullin 2/intermedin is a slow off-rate, long-acting endogenous agonist of the adrenomedullin 2 G protein-coupled receptor

The signaling peptides adrenomedullin 2/intermedin (AM2/IMD), adrenomedullin (AM), and CGRP have overlapping and distinct functions in the cardiovascular, lymphatic, and nervous systems by activating three shared receptors comprised of the class B GPCR CLR in complex with a RAMP1, -2, or -3 modulatory subunit. Here, we report that AM2/IMD, which is thought to be a non-selective agonist, is kinetically selective for CLR-RAMP3, known as the AM 2 R. AM2/IMD-AM 2 R elicited substantially longer duration cAMP signaling than the eight other peptide-receptor combinations due to AM2/IMD slow off-rate binding kinetics. The regions responsible for the slow off-rate were mapped to the AM2/IMD mid-region and the RAMP3 extracellular domain. MD simulations revealed how these bestow enhanced stability to the complex. Our results uncover AM2/IMD-AM 2 R as a cognate pair with unique temporal features, define the mechanism of kinetic selectivity, and explain how AM2/IMD and RAMP3 collaborate to shape the signaling output of a clinically important GPCR.

Insights from incorporating quantum computing into drug design workflows

MOTIVATION

While many quantum computing (QC) methods promise theoretical advantages over classical counterparts, quantum hardware remains limited. Exploiting near-term QC in computer-aided drug design (CADD) thus requires judicious partitioning between classical and quantum calculations.

RESULTS

We present HypaCADD, a hybrid classical-quantum workflow for finding ligands binding to proteins, while accounting for genetic mutations. We explicitly identify modules of our drug-design workflow currently amenable to replacement by QC: non-intuitively, we identify the mutation-impact predictor as the best candidate. HypaCADD thus combines classical docking and molecular dynamics with quantum machine learning (QML) to infer the impact of mutations. We present a case study with the coronavirus (SARS-CoV-2) protease and associated mutants. We map a classical machine-learning module onto QC, using a neural network constructed from qubit-rotation gates. We have implemented this in simulation and on two commercial quantum computers. We find that the QML models can perform on par with, if not better than, classical baselines. In summary, HypaCADD offers a successful strategy for leveraging QC for CADD.

AVAILABILITY AND IMPLEMENTATION

Jupyter Notebooks with Python code are freely available for academic use on GitHub: https://www.github.com/hypahub/hypacadd_notebook.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Conformational Change of Transcription Factors from Search to Specific Binding: A lac Repressor Case Study

In a process known as facilitated diffusion, DNA-binding proteins find their target sites by combining three-dimensional diffusion and one-dimensional scanning of the DNA. Following the trade-off between speed and stability, agile exploration of DNA requires loose binding, whereas, at the DNA target site, the searching protein needs to establish tight interactions with the DNA. To enable both efficient search and stable binding, DNA-binding proteins and DNA often switch conformations upon recognition. Here, we study the one-dimensional diffusion and DNA binding of the dimeric lac repressor (LacI), which was reported to adopt two different conformations when binding different conformations of DNA. Using coarse-grained molecular dynamic simulations, we studied the diffusion and the sequence-specific binding of these conformations of LacI, as well as their truncated or monomeric variants, with two DNA conformations: straight and bent. The simulations were compared to experimental observables. This study supports that linear diffusion along DNA combines tight rotation-coupled groove tracking and rotation-decoupled hopping, where the protein briefly dissociates and reassociates just a few base pairs away. Tight groove tracking is crucial for target-site recognition, while hopping speeds up the overall search process. We investigated the diffusion of different LacI conformations on DNA and show how the flexibility of LacI's hinge regions ensures agility on DNA as well as faithful groove tracking. If the hinge regions instead form α-helices at the protein-DNA interface, tight groove tracking is not possible. On the contrary, the helical hinge region is essential for tight binding to bent, specific DNA, for the formation of the specific complex. Based on our study of different encounter complexes, we argue that the conformational change in LacI and DNA bending are somewhat coupled. Our findings underline the importance of two distinct protein conformations for facilitated diffusion and specific binding, respectively.

Collagen a1(XI) structure prediction by Alphafold 2

Collagen α1(XI) is a minor fibrillar collagen involved in the critical regulation of collagen fibrils such as nucleation, assembly, and regulation of fibril diameter. The amino propeptide domain of the collagen α1(XI) is retained on the surface of the collagen fibril for an extended period of time and may play a crucial role in the interaction with extracellular matrix glycosaminoglycans and other proteins during the process of fibrillogenesis. Understanding the mechanism of action of this protein will ultimately help us understand the organization and assembly of the extracellular matrix that underlies the structural integrity of connective tissues.

CaV1.2 channelopathic mutations evoke diverse pathophysiological mechanisms

The first pathogenic mutation in CaV1.2 was identified in 2004 and was shown to cause a severe multisystem disorder known as Timothy syndrome (TS). The mutation was localized to the distal S6 region of the channel, a region known to play a major role in channel activation. TS patients suffer from life-threatening cardiac symptoms as well as significant neurodevelopmental deficits, including autism spectrum disorder (ASD). Since this discovery, the number and variety of mutations identified in CaV1.2 have grown tremendously, and the distal S6 regions remain a frequent locus for many of these mutations. While the majority of patients harboring these mutations exhibit cardiac symptoms that can be well explained by known pathogenic mechanisms, the same cannot be said for the ASD or neurodevelopmental phenotypes seen in some patients, indicating a gap in our understanding of the pathogenesis of CaV1.2 channelopathies. Here, we use whole-cell patch clamp, quantitative Ca2+ imaging, and single channel recordings to expand the known mechanisms underlying the pathogenesis of CaV1.2 channelopathies. Specifically, we find that mutations within the S6 region can exert independent and separable effects on activation, voltage-dependent inactivation (VDI), and Ca2+-dependent inactivation (CDI). Moreover, the mechanisms underlying the CDI effects of these mutations are varied and include altered channel opening and possible disruption of CDI transduction. Overall, these results provide a structure-function framework to conceptualize the role of S6 mutations in pathophysiology and offer insight into the biophysical defects associated with distinct clinical manifestations.

The human RNA polymerase I structure reveals an HMG-like docking domain specific to metazoans

We characterize the human RNA polymerase I by evolutionary biochemistry and cryo-EM revealing a built-in structural domain that apparently serves as transcription factor–binding platform in metazoans.

Transcription of the ribosomal RNA precursor by RNA polymerase (Pol) I is a major determinant of cellular growth, and dysregulation is observed in many cancer types. Here, we present the purification of human Pol I from cells carrying a genomic GFP fusion on the largest subunit allowing the structural and functional analysis of the enzyme across species. In contrast to yeast, human Pol I carries a single-subunit stalk, and in vitro transcription indicates a reduced proofreading activity. Determination of the human Pol I cryo-EM reconstruction in a close-to-native state rationalizes the effects of disease-associated mutations and uncovers an additional domain that is built into the sequence of Pol I subunit RPA1. This “dock II” domain resembles a truncated HMG box incapable of DNA binding which may serve as a downstream transcription factor–binding platform in metazoans. Biochemical analysis, in situ modelling, and ChIP data indicate that Topoisomerase 2a can be recruited to Pol I via the domain and cooperates with the HMG box domain–containing factor UBF. These adaptations of the metazoan Pol I transcription system may allow efficient release of positive DNA supercoils accumulating downstream of the transcription bubble.

Descriptive and functional analyses of four cyclin proteins in Trichomonas vaginalis

Trichomonas vaginalis is an early divergent protozoan parasite that causes trichomoniasis, the most common non-viral sexually transmitted infection. In metazoans, there is abundant and detailed research on the cell cycle and the components involved in the regulation mechanisms. Regulators such as the cyclin-dependent kinases (CDKs) and cyclins activate the highly regulated processes of cell division. While CDKs have important roles in the phosphorylation of specific substrates, cyclins are important activating-components of CDKs that allow orderly passage through the different stages of the cell cycle. Cell cycle cyclins are characterized by showing drastic changes in their concentration during the cell cycle progression. However, in protists such as T. vaginalis, some biological processes such as cell cycle regulation remain less well studied. In an attempt to gain insight into cell cycle regulation in T. vaginalis, as an initial approach we characterized four proteins with features of cyclins. The genes encoding these putative cyclins were cloned to produce the recombinant proteins TvCYC1, TvCYC2, TvCYC3, and TvCYC4. The functional activity of TvCYC2, TvCYC3, and TvCYC4 was assessed through their complementation of a yeast cln1,2,3Δ mutant strain; TvCYC1 was not able to complement this mutant. Furthermore, our results suggest that TvCYC1, TvCYC2, and TvCYC3, are able to interact with and activate the kinase activity of TvCRK1, a kinase previously characterized by our group. The present study represents the first characterization of cyclins potentially involved in cell cycle regulation in T. vaginalis.

Structural and Functional Impact of Damaging Nonsynonymous Single Nucleotide Polymorphisms (nsSNPs) on Human VPS35 Protein Using Computational Approaches

Parkinson's disease is the second most common progressive neurodegenerative movement disorder. Mutations in retromer complex subunit and VPS35 represent the second most common cause of late-onset familial Parkinson's disease. The mutation in VPS35 can disrupt the normal protein functions resulting in Parkinson's disease. The aim of this study was the identification of deleterious missense Single Nucleotide Polymorphisms (nsSNPs) and their structural and functional impact on the VPS35 protein. In this study, several insilico tools were used to identify deleterious and disease-associated nsSNPs. 3D structure of VPS35 protein was constructed through MODELLER 9.2, normalized using FOLDX, and evaluated through RAMPAGE and ERRAT whereas, FOLDX was used for mutagenesis. 25 ligands were obtained from literature and docked using PyRx 0.8 software. Based on the binding affinity, five ligands i.e., PG4, MSE, GOL, EDO, and CAF were further analyzed. Molecular Dynamic simulation analysis was performed using GROMACS 5.1.4, where temperature, pressure, density, RMSD, RMSF, Rg, and SASA graphs were analyzed. The results showed that the mutations Y67H, R524W, and D620N had a structural and functional impact on the VPS35 protein. The current findings will help in appropriate drug design against the disease caused by these mutations in a large population using in-vitro study.

Molecular docking and biochemical validation of (-)-syringaresinol-4-O-β-D-apiofuranosyl-(1→2)-β-D-glucopyranoside binding to an allosteric site in monoamine transporters

Albizia julibrissin Durazz is one of the most common herbs used for depression and anxiety treatment, but its mechanism of action as an antidepressant or anxiolytic drug have not been fully understood. We previously isolated and identified one lignan glycoside compound from Albizia Julibrissin Durazz, (-)-syringaresinol-4-O-β-D-apiofuranosyl-(1→2)-β-D-glucopyranoside (SAG), that inhibited all three monoamine transporters with a mechanism of action different from that of the conventional antidepressants. In this study, we generated homology models for human dopamine transporter and human norepinephrine transporter, based on the X-ray structure of Drosophila dopamine transporter, and conducted the molecular docking of SAG to all three human monoamine transporters. Our computational results indicated that SAG binds to an allosteric site (S2) that has been demonstrated to be formed by an aromatic pocket positioned in the scaffold domain in the extracellular vestibule connected to the central site (S1) in these monoamine transporters. In addition, we demonstrated that SAG stabilizes a conformation of serotonin transporter with both the extracellular and cytoplasmic pathways closed. Furthermore, we performed mutagenesis of the residues in both the allosteric and orthosteric sites to biochemically validate SAG binding in all three monoamine transporters. Our results are consistent with the molecular docking calculation and support the association of SAG with the allosteric site. We expect that this herbal molecule could become a lead compound for the development of new therapeutic agents with a novel mechanism of action.

The ABC transporter MsbA adopts the wide inward-open conformation in E. coli cells

Membrane proteins are currently investigated after detergent extraction from native cellular membranes and reconstitution into artificial liposomes or nanodiscs, thereby removing them from their physiological environment. However, to truly understand the biophysical properties of membrane proteins in a physiological environment, they must be investigated within living cells. Here, we used a spin-labeled nanobody to interrogate the conformational cycle of the ABC transporter MsbA by double electron-electron resonance. Unexpectedly, the wide inward-open conformation of MsbA, commonly considered a nonphysiological state, was found to be prominently populated in Escherichia coli cells. Molecular dynamics simulations revealed that extensive lateral portal opening is essential to provide access of its large natural substrate core lipid A to the binding cavity. Our work paves the way to investigate the conformational landscape of membrane proteins in cells.

A New Structural Model of Apolipoprotein B100 Based on Computational Modeling and Cross Linking

ApoB-100 is a member of a large lipid transfer protein superfamily and is one of the main apolipoproteins found on low-density lipoprotein (LDL) and very low-density lipoprotein (VLDL) particles. Despite its clinical significance for the development of cardiovascular disease, there is limited information on apoB-100 structure. We have developed a novel method based on the “divide and conquer” algorithm, using PSIPRED software, by dividing apoB-100 into five subunits and 11 domains. Models of each domain were prepared using I-TASSER, DEMO, RoseTTAFold, Phyre2, and MODELLER. Subsequently, we used disuccinimidyl sulfoxide (DSSO), a new mass spectrometry cleavable cross-linker, and the known position of disulfide bonds to experimentally validate each model. We obtained 65 unique DSSO cross-links, of which 87.5% were within a 26 Å threshold in the final model. We also evaluated the positions of cysteine residues involved in the eight known disulfide bonds in apoB-100, and each pair was measured within the expected 5.6 Å constraint. Finally, multiple domains were combined by applying constraints based on detected long-range DSSO cross-links to generate five subunits, which were subsequently merged to achieve an uninterrupted architecture for apoB-100 around a lipoprotein particle. Moreover, the dynamics of apoB-100 during particle size transitions was examined by comparing VLDL and LDL computational models and using experimental cross-linking data. In addition, the proposed model of receptor ligand binding of apoB-100 provides new insights into some of its functions.

Structure of the malaria vaccine candidate Pfs48/45 and its recognition by transmission blocking antibodies

An effective malaria vaccine remains a global health priority and vaccine immunogens which prevent transmission of the parasite will have important roles in multi-component vaccines. One of the most promising candidates for inclusion in a transmission-blocking malaria vaccine is the gamete surface protein Pfs48/45, which is essential for development of the parasite in the mosquito midgut. Indeed, antibodies which bind Pfs48/45 can prevent transmission if ingested with the parasite as part of the mosquito bloodmeal. Here we present the structure of full-length Pfs48/45, showing its three domains to form a dynamic, planar, triangular arrangement. We reveal where transmission-blocking and non-blocking antibodies bind on Pfs48/45. Finally, we demonstrate that antibodies which bind across this molecule can be transmission-blocking. These studies will guide the development of future Pfs48/45-based vaccine immunogens.

Discovery and characterization of SARS-CoV-2 reactive and neutralizing antibodies from humanized CAMouseHG mice through rapid hybridoma screening and high-throughput single-cell V(D)J sequencing

The coronavirus disease 2019 pandemic has caused more than 532 million infections and 6.3 million deaths to date. The reactive and neutralizing fully human antibodies of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are effective detection tools and therapeutic measures. During SARS-CoV-2 infection, a large number of SARS-CoV-2 reactive and neutralizing antibodies will be produced. Most SARS-CoV-2 reactive and neutralizing fully human antibodies are isolated from human and frequently encoded by convergent heavy-chain variable genes. However, SARS-CoV-2 viruses can mutate rapidly during replication and the resistant variants of neutralizing antibodies easily survive and evade the immune response, especially in the face of such focused antibody responses in humans. Therefore, additional tools are needed to develop different kinds of fully human antibodies to compensate for current deficiency. In this study, we utilized antibody humanized CAMouse^HG mice to develop a rapid antibody discovery method and examine the antibody repertoire of SARS-CoV-2 RBD-reactive hybridoma cells derived from CAMouse^HG mice by using high-throughput single-cell V(D)J sequencing analysis. CAMouse^HG mice were immunized by 28-day rapid immunization method. After electrofusion and semi-solid medium screening on day 12 post-electrofusion, 171 hybridoma clones were generated based on the results of SARS-CoV-2 RBD binding activity assay. A rather obvious preferential usage of IGHV6-1 family was found in these hybridoma clones derived from CAMouse^HG mice, which was significantly different from the antibodies found in patients with COVID-19. After further virus neutralization screening and antibody competition assays, we generated a noncompeting two-antibody cocktail, which showed a potent prophylactic protective efficacy against SARS-CoV-2 in cynomolgus macaques. These results indicate that humanized CAMouse^HG mice not only provide a valuable platform to obtain fully human reactive and neutralizing antibodies but also have a different antibody repertoire from humans. Thus, humanized CAMouse^HG mice can be used as a good complementary tool in discovery of fully human therapeutic and diagnostic antibodies.

Engineering SARS-CoV-2 specific cocktail antibodies into a bispecific format improves neutralizing potency and breadth

One major limitation of neutralizing antibody-based COVID-19 therapy is the requirement of costly cocktails to reduce emergence of antibody resistance. Here we engineer two bispecific antibodies (bsAbs) using distinct designs and compared them with parental antibodies and their cocktail. Single molecules of both bsAbs block the two epitopes targeted by parental antibodies on the receptor-binding domain (RBD). However, bsAb with the IgG-(scFv)2 design (14-H-06) but not the CrossMAb design (14-crs-06) shows increased antigen-binding and virus-neutralizing activities against multiple SARS-CoV-2 variants as well as increased breadth of neutralizing activity compared to the cocktail. X-ray crystallography and cryo-EM reveal distinct binding models for individual cocktail antibodies, and computational simulations suggest higher inter-spike crosslinking potentials by 14-H-06 than 14-crs-06. In mouse models of infections by SARS-CoV-2 and multiple variants, 14-H-06 exhibits higher or equivalent therapeutic efficacy than the cocktail. Rationally engineered bsAbs represent a cost-effective alternative to antibody cocktails and a promising strategy to improve potency and breadth.

The allosteric mechanism leading to an open-groove lipid conductive state of the TMEM16F scramblase

TMEM16F is a Ca2+-activated phospholipid scramblase in the TMEM16 family of membrane proteins. Unlike other TMEM16s exhibiting a membrane-exposed hydrophilic groove that serves as a translocation pathway for lipids, the experimentally determined structures of TMEM16F shows the groove in a closed conformation even under conditions of maximal scramblase activity. It is currently unknown if/how TMEM16F groove can open for lipid scrambling. Here we describe the analysis of ~400 µs all-atom molecular dynamics (MD) simulations of the TMEM16F revealing an allosteric mechanism leading to an open-groove, lipid scrambling competent state of the protein. The groove opens into a continuous hydrophilic conduit that is highly similar in structure to that seen in other activated scramblases. The allosteric pathway connects this opening to an observed destabilization of the Ca2+ ion bound at the distal site near the dimer interface, to the dynamics of specific protein regions that produces the open-groove state to scramble phospholipids.

Exploring the RNase A scaffold to combine catalytic and antimicrobial activities. Structural characterization of RNase 3/1 chimeras

Design of novel antibiotics to fight antimicrobial resistance is one of the first global health priorities. Novel protein-based strategies come out as alternative therapies. Based on the structure-function knowledge of the RNase A superfamily we have engineered a chimera that combines RNase 1 highest catalytic activity with RNase 3 unique antipathogen properties. A first construct (RNase 3/1-v1) was successfully designed with a catalytic activity 40-fold higher than RNase 3, but alas in detriment of its anti-pathogenic activity. Next, two new versions of the original chimeric protein were created showing improvement in the antimicrobial activity. Both second generation versions (RNases 3/1-v2 and -v3) incorporated a loop characteristic of RNase 3 (L7), associated to antimicrobial activity. Last, removal of an RNase 1 flexible loop (L1) in the third version enhanced its antimicrobial properties and catalytic efficiency. Here we solved the 3D structures of the three chimeras at atomic resolution by X-ray crystallography. Structural analysis outlined the key functional regions. Prediction by molecular docking of the protein chimera in complex with dinucleotides highlighted the contribution of the C-terminal region to shape the substrate binding cavity and determine the base selectivity and catalytic efficiency. Nonetheless, the structures that incorporated the key features related to RNase 3 antimicrobial activity retained the overall RNase 1 active site conformation together with the essential structural elements for binding to the human ribonuclease inhibitor (RNHI), ensuring non-cytotoxicity. Results will guide us in the design of the best RNase pharmacophore for anti-infective therapies.

A combined in silico and in vivo approach to the structure-function annotation of SPD-2 provides mechanistic insight into its functional diversity

Centrosomes are organelles that function as hubs of microtubule nucleation and organization, with key roles in organelle positioning, asymmetric cell division, ciliogenesis, and signaling. Aberrant centrosome number, structure or function is linked to neurodegenerative diseases, developmental abnormalities, ciliopathies, and tumor development. A major regulator of centrosome biogenesis and function in C. elegans is the conserved Spindle-defective protein 2 (SPD-2), a homolog of the human CEP-192 protein. CeSPD-2 is required for centrosome maturation, centriole duplication, spindle assembly and possibly cell polarity establishment. Despite its importance, the specific molecular mechanism of CeSPD-2 regulation and function is poorly understood. Here, we combined computational analysis with cell biology approaches to uncover possible structure-function relationships of CeSPD-2 that may shed mechanistic light on its function. Domain prediction analysis corroborated and refined previously identified coiled-coils and ASH (Aspm-SPD-2 Hydin) domains and identified new domains: a GEF domain, an Ig-like domain, and a PDZ-like domain. In addition to these predicted structural features, CeSPD-2 is also predicted to be intrinsically disordered. Surface electrostatic maps identified a large basic region unique to the ASH domain of CeSPD-2. This basic region overlaps with most of the residues predicted to be involved in protein-protein interactions. In vivo, ASH::GFP localized to centrosomes and centrosome-associated microtubules. Our analysis groups ASH domains, PapD, Usher chaperone domains, and Major Sperm Protein (MSP) domains into a single superfold within the larger Immunoglobulin superfamily. This study lays the groundwork for designing rational hypothesis-based experiments to uncover the mechanisms of CeSPD-2 function in vivo.Abbreviations: AIR, Aurora kinase; ASH, Aspm-SPD-2 Hydin; ASP, Abnormal Spindle Protein; ASPM, Abnormal Spindle-like Microcephaly-associated Protein; CC, coiled-coil; CDK, Cyclin-dependent Kinase; Ce, Caenorhabditis elegans; CEP, Centrosomal Protein; CPAP, centrosomal P4.1-associated protein; D, Drosophila; GAP, GTPase activating protein; GEF, GTPase guanine nucleotide exchange factor; Hs, Homo sapiens/Human; Ig, Immunoglobulin; MAP, Microtubule associated Protein; MSP, Major Sperm Protein; MDP, Major Sperm Domain-Containing Protein; OCRL-1, Golgi endocytic trafficking protein Inositol polyphosphate 5-phosphatase; PAR, abnormal embryonic PARtitioning of the cytosol; PCM, Pericentriolar material; PCMD, pericentriolar matrix deficient; PDZ, PSD95/Dlg-1/zo-1; PLK, Polo like kinase; RMSD, Root Mean Square Deviation; SAS, Spindle assembly abnormal proteins; SPD, Spindle-defective protein; TRAPP, TRAnsport Protein Particle; Xe, Xenopus; ZYG, zygote defective protein.

Highly efficient production of rebaudioside D enabled by structure-guided engineering of bacterial glycosyltransferase YojK

Owing to zero-calorie, high-intensity sweetness and good taste profile, the plant-derived sweetener rebaudioside D (Reb D) has attracted great interest to replace sugars. However, low content of Reb D in stevia rebaudiana Bertoni as well as low soluble expression and enzymatic activity of plant-derived glycosyltransferase in Reb D preparation restrict its commercial usage. To address these problems, a novel glycosyltransferase YojK from Bacillus subtilis 168 with the ability to glycosylate Reb A to produce Reb D was identified. Then, structure-guided engineering was performed after solving its crystal structure. A variant YojK-I241T/G327N with 7.35-fold increase of the catalytic activity was obtained, which allowed to produce Reb D on a scale preparation with a great yield of 91.29%. Moreover, based on the results from molecular docking and molecular dynamics simulations, the improvement of enzymatic activity of YojK-I241T/G327N was ascribed to the formation of new hydrogen bonds between the enzyme and substrate or uridine diphosphate glucose. Therefore, this study provides an engineered bacterial glycosyltransferase YojK-I241T/G327N with high solubility and catalytic efficiency for potential industrial scale-production of Reb D.

Rotavirus VP4 Epitope of a Broadly Neutralizing Human Antibody Defined by Its Structure Bound with an Attenuated-Strain Virion

Rotavirus live-attenuated vaccines, both mono- and pentavalent, generate broadly heterotypic protection. B-cells isolated from adults encode neutralizing antibodies, some with affinity for VP5*, that afford broad protection in mice. We have mapped the epitope of one such antibody by determining the high-resolution cryo-EM structure of its antigen-binding fragment (Fab) bound to the virion of a candidate vaccine strain, CDC-9. The Fab contacts both the distal end of a VP5* β-barrel domain and the two VP8* lectin-like domains at the tip of a projecting spike. Its interactions with VP8* do not impinge on the likely receptor-binding site, suggesting that the mechanism of neutralization is at a step subsequent to initial attachment. We also examined structures of CDC-9 virions from two different stages of serial passaging. Nearly all the VP4 (cleaved to VP8*/VP5*) spikes on particles from the earlier passage (wild-type isolate) had transitioned from the "upright" conformation present on fully infectious virions to the "reversed" conformation that is probably the end state of membrane insertion, unable to mediate penetration, consistent with the very low in vitro infectivity of the wild-type isolate. About half the VP4 spikes were upright on particles from the later passage, which had recovered substantial in vitro infectivity but had acquired an attenuated phenotype in neonatal rats. A mutation in VP4 that occurred during passaging appears to stabilize the interface at the apex of the spike and could account for the greater stability of the upright spikes on the late-passage, attenuated isolate. IMPORTANCE Rotavirus live-attenuated vaccines generate broadly heterotypic protection, and B-cells isolated from adults encode antibodies that are broadly protective in mice. Determining the structural and mechanistic basis of broad protection can contribute to understanding the current limitations of vaccine efficacy in developing countries. The structure of an attenuated human rotavirus isolate (CDC-9) bound with the Fab fragment of a broadly heterotypic protective antibody shows that protection is probably due to inhibition of the conformational transition in the viral spike protein (VP4) critical for viral penetration, rather than to inhibition of receptor binding. A comparison of structures of CDC-9 virus particles at two stages of serial passaging supports a proposed mechanism for initial steps in rotavirus membrane penetration.

Metabolic and Structural Insights into Hydrogen Sulfide Mis-Regulation in Enterococcus faecalis

Hydrogen sulfide (H2S) is implicated as a cytoprotective agent that bacteria employ in response to host-induced stressors, such as oxidative stress and antibiotics. The physiological benefits often attributed to H2S, however, are likely a result of downstream, more oxidized forms of sulfur, collectively termed reactive sulfur species (RSS) and including the organic persulfide (RSSH). Here, we investigated the metabolic response of the commensal gut microorganism Enterococcus faecalis to exogenous Na2S as a proxy for H2S/RSS toxicity. We found that exogenous sulfide increases protein abundance for enzymes responsible for the biosynthesis of coenzyme A (CoA). Proteome S-sulfuration (persulfidation), a posttranslational modification implicated in H2S signal transduction, is also widespread in this organism and is significantly elevated by exogenous sulfide in CstR, the RSS sensor, coenzyme A persulfide (CoASSH) reductase (CoAPR) and enzymes associated with de novo fatty acid biosynthesis and acetyl-CoA synthesis. Exogenous sulfide significantly impacts the speciation of fatty acids as well as cellular concentrations of acetyl-CoA, suggesting that protein persulfidation may impact flux through these pathways. Indeed, CoASSH is an inhibitor of E. faecalis phosphotransacetylase (Pta), suggesting that an important metabolic consequence of increased levels of H2S/RSS may be over-persulfidation of this key metabolite, which, in turn, inhibits CoA and acyl-CoA-utilizing enzymes. Our 2.05 Å crystallographic structure of CoA-bound CoAPR provides new structural insights into CoASSH clearance in E. faecalis.

Functional and structural characteristics of HLA-B*13:01-mediated specific T cells reaction in dapsone-induced drug hypersensitivity

BACKGROUND

Severe cutaneous adverse drug reactions (SCARs) are a group of serious clinical conditions caused by immune reaction to certain drugs. The allelic variance of human leukocyte antigens of HLA-B*13:01 has been strongly associated with hypersensitivities induced by dapsone (DDS). T-cell receptor mediated activation of cytotoxic T lymphocytes (CTLs) has also been suggested to play an essential role in pathogenesis of SCARs. However, HLA-B*13:01-DDS-TCR immune synapse that plays role in drug-induced hypersensitivity syndrome (DIHS) associated T cells activation remains uncharacterized.

METHODS

To investigate the molecular mechanisms for HLA-B*13:01 in the pathogenesis of Dapsone-induced drug hypersensitivity (DDS-DIHS), we performed crystallization and expanded drug-specific CTLs to analyze the pathological role of DDS-DIHS.

RESULTS

Results showed the crystal structure of HLA-B*13:01-beta-2-microglobulin (β2M) complex at 1.5 Å resolution and performed mutation assays demonstrating that I118 or I119, and R121 of HLA-B*13:01 were the key residues that mediate the binding of DDS. Subsequent single-cell TCR and RNA sequencing indicated that TCRs composed of paired TRAV12-3/TRBV28 clonotype with shared CDR3 region specifically recognize HLA-B*13:01-DDS complex to trigger inflammatory cytokines associated with DDS-DIHS.

CONCLUSION

Our study identified the novel p-i-HLA/TCR as the model of interaction between HLA-B*13:01, DDS and the clonotype-specific TCR in DDS-DIHS.

New KCNN4 Variants Associated With Anemia: Stomatocytosis Without Erythrocyte Dehydration

The K+ channel activated by the Ca2+, KCNN4, has been shown to contribute to red blood cell dehydration in the rare hereditary hemolytic anemia, the dehydrated hereditary stomatocytosis. We report two de novo mutations on KCNN4, We reported two de novo mutations on KCNN4, V222L and H340N, characterized at the molecular, cellular and clinical levels. Whereas both mutations were shown to increase the calcium sensitivity of the K+ channel, leading to channel opening for lower calcium concentrations compared to WT KCNN4 channel, there was no obvious red blood cell dehydration in patients carrying one or the other mutation. The clinical phenotype was greatly different between carriers of the mutated gene ranging from severe anemia for one patient to a single episode of anemia for the other patient or no documented sign of anemia for the parents who also carried the mutation. These data compared to already published KCNN4 mutations question the role of KCNN4 gain-of-function mutations in hydration status and viability of red blood cells in bloodstream.

Chaperoning of the histone octamer by the acidic domain of DNA repair factor APLF

Nucleosome assembly requires the coordinated deposition of histone complexes H3-H4 and H2A-H2B to form a histone octamer on DNA. In the current paradigm, specific histone chaperones guide the deposition of first H3-H4 and then H2A-H2B. Here, we show that the acidic domain of DNA repair factor APLF (APLF^AD) can assemble the histone octamer in a single step and deposit it on DNA to form nucleosomes. The crystal structure of the APLF^AD-histone octamer complex shows that APLF^AD tethers the histones in their nucleosomal conformation. Mutations of key aromatic anchor residues in APLF^AD affect chaperone activity in vitro and in cells. Together, we propose that chaperoning of the histone octamer is a mechanism for histone chaperone function at sites where chromatin is temporarily disrupted.

Targeting the YXXΦ Motifs of the SARS Coronaviruses 1 and 2 ORF3a Peptides by In Silico Analysis to Predict Novel Virus-Host Interactions

The emerging SARS-CoV and SARS-CoV-2 belong to the family of "common cold" RNA coronaviruses, and they are responsible for the 2003 epidemic and the current pandemic with over 6.3 M deaths worldwide. The ORF3a gene is conserved in both viruses and codes for the accessory protein ORF3a, with unclear functions, possibly related to viral virulence and pathogenesis. The tyrosine-based YXXΦ motif (Φ: bulky hydrophobic residue-L/I/M/V/F) was originally discovered to mediate clathrin-dependent endocytosis of membrane-spanning proteins. Many viruses employ the YXXΦ motif to achieve efficient receptor-guided internalisation in host cells, maintain the structural integrity of their capsids and enhance viral replication. Importantly, this motif has been recently identified on the ORF3a proteins of SARS-CoV and SARS-CoV-2. Given that the ORF3a aa sequence is not fully conserved between the two SARS viruses, we aimed to map in silico structural differences and putative sequence-driven alterations of regulatory elements within and adjacently to the YXXΦ motifs that could predict variations in ORF3a functions. Using robust bioinformatics tools, we investigated the presence of relevant post-translational modifications and the YXXΦ motif involvement in protein-protein interactions. Our study suggests that the predicted YXXΦ-related features may confer specific-yet to be discovered-functions to ORF3a proteins, significant to the new virus and related to enhanced propagation, host immune regulation and virulence.

Expanded ACE2 dependencies of diverse SARS-like coronavirus receptor binding domains

Viral spillover from animal reservoirs can trigger public health crises and cripple the world economy. Knowing which viruses are primed for zoonotic transmission can focus surveillance efforts and mitigation strategies for future pandemics. Successful engagement of receptor protein orthologs is necessary during cross-species transmission. The clade 1 sarbecoviruses including Severe Acute Respiratory Syndrome-related Coronavirus (SARS-CoV) and SARS-CoV-2 enter cells via engagement of angiotensin converting enzyme-2 (ACE2), while the receptor for clade 2 and clade 3 remains largely uncharacterized. We developed a mixed cell pseudotyped virus infection assay to determine whether various clades 2 and 3 sarbecovirus spike proteins can enter HEK 293T cells expressing human or Rhinolophus horseshoe bat ACE2 proteins. The receptor binding domains from BtKY72 and Khosta-2 used human ACE2 for entry, while BtKY72 and Khosta-1 exhibited widespread use of diverse rhinolophid ACE2s. A lysine at ACE2 position 31 appeared to be a major determinant of the inability of these RBDs to use a certain ACE2 sequence. The ACE2 protein from Rhinolophus alcyone engaged all known clade 3 and clade 1 receptor binding domains. We observed little use of Rhinolophus ACE2 orthologs by the clade 2 viruses, supporting the likely use of a separate, unknown receptor. Our results suggest that clade 3 sarbecoviruses from Africa and Europe use Rhinolophus ACE2 for entry, and their spike proteins appear primed to contribute to zoonosis under the right conditions.

Knowing which viruses are primed for zoonotic transmission can focus surveillance efforts and mitigation strategies for future pandemics. This study shows that SARS-like coronaviruses identified in bats from Europe and Africa can use a range of horseshoe bat ACE2s for entry. In addition, viruses found in Russia and Kenya also have the ability to at least weakly use human ACE2.

Comparative phylogeny and evolutionary analysis of Dicer-like protein family in two plant monophyletic lineages

BACKGROUND

Small RNAs (sRNAs) that do not get untranslated into proteins exhibit a pivotal role in the expression regulation of their cognate gene(s) in almost all eukaryotic lineages, including plants. Hitherto, numerous protein families such as Dicer, a unique class of Ribonuclease III, have been reported to be involved in sRNAs processing pathways and silencing. In this study, we aimed to investigate the phylogenetic relationship and evolutionary history of the DCL protein family.

RESULTS

Our results illustrated the DCL family of proteins grouped into four main subfamilies (DCLs 1-4) presented in either Eudicotyledons or Liliopsids. The accurate observation of the phylogenetic trees supports the independent expansion of DCL proteins among the Eudicotyledons and Liliopsids species. They share the common origin, and the main duplication events for the formation of the DCL subfamilies occurred before the Eudicotyledons/Liliopsids split from their ancestral DCL. In addition, shreds of evidence revealed that the divergence happened when multicellularization started and since the need for complex gene regulation considered being a necessity by organisms. At that time, they have evolved independently among the monophyletic lineages. The other finding was that the combination of DCL protein subfamilies bears several highly conserved functional domains in plant species that originated from their ancestor architecture. The conservation of these domains happens to be both lineage-specific and inter lineage-specific.

CONCLUSIONS

DCL subfamilies (i.e., DCL1-DCL4) distribute in their single clades after diverging from their common ancestor and before emerging into higher plants. Therefore, it seems that the main duplication events for the formation of the DCL subfamilies occurred before the Eudicotyledons/Liliopsida split and before the appearance of moss, and after the single-cell green algae. We also observed the same trends among the main DCL subfamilies from functional unit composition and architecture. Despite the long evolutionary course from the divergence of Liliopsida lineage from the Eudicotyledons, a significant diversifying force to domain composition and orientation was absent. The results of this study provide a deeper insight into DCL protein evolutionary history and possible sequence and structural relationships between DCL protein subfamilies in the main higher plant monophyletic lineages; i.e., Eudicotyledons and Liliopsida.