sAMP-VGG16: Force-field assisted image-based deep neural network prediction model for short antimicrobial peptides.

During the last three decades, antimicrobial peptides (AMPs) have emerged as a promising therapeutic alternative to antibiotics. The approaches for designing AMPs span from experimental trial-and-error methods to synthetic hybrid peptide libraries. To overcome the exceedingly expensive and time-consuming process of designing effective AMPs, many computational and machine-learning tools for AMP prediction have been recently developed. In general, to encode the peptide sequences, featurization relies on approaches based on (a) amino acid (AA) composition, (b) physicochemical properties, (c) sequence similarity, and (d) structural properties. In this work, we present an image-based deep neural network model to predict AMPs, where we are using feature encoding based on Drude polarizable force-field atom types, which can capture the peptide properties more efficiently compared to conventional feature vectors. The proposed prediction model identifies short AMPs (≤30 AA) with promising accuracy and efficiency and can be used as a next-generation screening method for predicting new AMPs. The source code is publicly available at the Figshare server sAMP-VGG16.

LncRNA Malat1 suppresses pyroptosis and T cell-mediated killing of incipient metastatic cells.

The contribution of antitumor immunity to metastatic dormancy is poorly understood. Here we show that the long noncoding RNA Malat1 is required for tumor initiation and metastatic reactivation in mouse models of breast cancer and other tumor types. Malat1 localizes to nuclear speckles to couple transcription, splicing and mRNA maturation. In metastatic cells, Malat1 induces WNT ligands, autocrine loops to promote self-renewal and the expression of Serpin protease inhibitors. Through inhibition of caspase-1 and cathepsin G, SERPINB6B prevents gasdermin D-mediated induction of pyroptosis. In this way, SERPINB6B suppresses immunogenic cell death and confers evasion of T cell-mediated tumor lysis of incipient metastatic cells. On-target inhibition of Malat1 using therapeutic antisense nucleotides suppresses metastasis in a SERPINB6B-dependent manner. These results suggest that Malat1-induced expression of SERPINB6B can titrate pyroptosis and immune recognition at metastatic sites. Thus, Malat1 is at the nexus of tumor initiation, reactivation and immune evasion and represents a tractable and clinically relevant drug target.

Combining SILCS and Artificial Intelligence for High-Throughput Prediction of the Passive Permeability of Drug Molecules.

Membrane permeability of drug molecules plays a significant role in the development of new therapeutic agents. Accordingly, methods to predict the passive permeability of drug candidates during a medicinal chemistry campaign offer the potential to accelerate the drug design process. In this work, we combine the physics-based site identification by ligand competitive saturation (SILCS) method and data-driven artificial intelligence (AI) to create a high-throughput predictive model for the passive permeability of druglike molecules. In this study, we present a comparative analysis of four regression models to predict membrane permeabilities of small druglike molecules; of the tested models, Random Forest was the most predictive yielding an R2 of 0.81 for the independent data set. The input feature vector used to train the developed prediction model includes absolute free energy profiles of ligands through a POPC-cholesterol bilayer based on ligand grid free energy (LGFE) profiles obtained from the SILCS approach. The use of the membrane free energy profiles from SILCS offers information on the physical forces contributing to ligand permeability, while the use of AI yields a more predictive model trained on experimental PAMPA permeability data for a collection of 229 molecules. This combination allows for rapid estimations of ligand permeability at a level of accuracy beyond currently available predictive models while offering insights into the contributions of the functional groups in the ligands to the permeability barrier, thereby offering quantitative information to facilitate rational ligand design.

Enhanced myogenesis through lncFAM-mediated recruitment of HNRNPL to the MYBPC2 promoter.

The mammalian transcriptome comprises a vast family of long noncoding (lnc)RNAs implicated in physiologic processes such as myogenesis, through which muscle forms during embryonic development and regenerates in the adult. However, the specific molecular mechanisms by which lncRNAs regulate human myogenesis are poorly understood. Here, we identified a novel muscle-specific lncRNA, lncFAM71E1-2:2 (lncFAM), which increased robustly during early human myogenesis. Overexpression of lncFAM promoted differentiation of human myoblasts into myotubes, while silencing lncFAM suppressed this process. As lncFAM resides in the nucleus, chromatin isolation by RNA purification followed by mass spectrometry (ChIRP-MS) analysis was employed to identify the molecular mechanisms whereby it might promote myogenesis. Analysis of lncFAM-interacting proteins revealed that lncFAM recruited the RNA-binding protein HNRNPL to the promoter of MYBPC2, in turn increasing MYBPC2 mRNA transcription and enhancing production of the myogenic protein MYBPC2. These results highlight a mechanism whereby a novel ribonucleoprotein complex, lncFAM-HNRNPL, elevates MYBPC2 expression transcriptionally to promote myogenesis.

LncRNA OIP5-AS1-directed miR-7 degradation promotes MYMX production during human myogenesis.

Long noncoding RNAs (lncRNAs) and microRNAs (miRNAs) modulate gene expression programs in physiology and disease. Here, we report a noncoding RNA regulatory network that modulates myoblast fusion into multinucleated myotubes, a process that occurs during muscle development and muscle regeneration after injury. In early stages of human myogenesis, the levels of lncRNA OIP5-AS1 increased, while the levels of miR-7 decreased. Moreover, OIP5-AS1 bound and induced miR-7 decay via target RNA-directed miRNA decay; accordingly, loss of OIP5-AS1 attenuated, while antagonizing miR-7 accelerated, myotube formation. We found that the OIP5-AS1-mediated miR-7 degradation promoted myoblast fusion, as it derepressed the miR-7 target MYMX mRNA, which encodes the fusogenic protein myomixer (MYMX). Remarkably, an oligonucleotide site blocker interfered with the OIP5-AS1-directed miR-7 degradation, allowing miR-7 to accumulate, lowering MYMX production and suppressing myotube formation. These results highlight a mechanism whereby lncRNA OIP5-AS1-mediated miR-7 decay promotes myotube formation by stimulating a myogenic fusion program.

LINC00162 regulates cell proliferation and apoptosis by sponging PAQR4-targeting miR-485-5p.

Growing evidence indicates that long intergenic noncoding RNAs play an important role in cancer progression by affecting gene regulation at the transcriptional and posttranscriptional levels. Recent studies have shown that long intergenic noncoding RNA functions as a competitive endogenous RNA, which can interact with and mitigate the function of microRNA. In this study, we investigated the molecular mechanism by which LINC00162 regulates cell proliferation and apoptotic cell death. By analyzing RNA sequencing data, LINC00162 was identified to be a target of heterogeneous nuclear ribonucleoprotein K (hnRNPK). HnRNPK positively regulated LINC00162 expression through p38 mitogen-activated protein kinase. Lowering the level of either hnRNPK or LINC00162 decreased proliferation and colony formation while it increased apoptotic cell death. Small RNA sequencing followed by the antisense oligonucleotide pulldown, revealed that LINC00162 interacts directly with miR-485-5p which exhibited tumor-suppressing effects by suppressing cell proliferation and colony formation, and increasing apoptotic cell death. Through the bioinformatic approaches, progestin and adipoQ receptor 4 (PAQR4) was selected as a common target of LINC00162 and miR-485-5p. miR-485-5p decreased the expression of PAQR4 by directly binding to the 3'-untranslated region of PAQR4 messenger RNA. Knockdown of hnRNPK and LINC00162 increased the level of functional miR-485-5p, indicating that LINC00162 may compete for miR-485-5p, thereby derepressing PAQR4 expression. Overexpression of either hnRNPK or LINC00162, or inhibition of miR-485-5p, protected cells against etoposide-induced apoptotic death. Our findings demonstrate that a regulatory paradigm implicating hnRNPK, LINC00162, miR-485-5p, and PAQR4 plays an important role in cell proliferation and apoptosis, and is a promising target for cancer therapeutics.

Deep Neural Network Model to Predict the Electrostatic Parameters in the Polarizable Classical Drude Oscillator Force Field.

The Drude polarizable force field (FF) captures electronic polarization effects via auxiliary Drude particles that are attached to non-hydrogen atoms, distinguishing it from commonly used additive FFs that rely on fixed charges. The Drude FF currently includes parameters for biomolecules such as proteins, nucleic acids, lipids, and carbohydrates and small-molecule representative of those classes of molecules as well as a range of atomic ions. Extension of the Drude FF to novel small druglike molecules is challenging as it requires the assignment of partial charges, atomic polarizabilities, and Thole scaling factors. In the present article, deep neural network (DNN) models are trained on quantum mechanical (QM)-based partial charges and atomic polarizabilities along with Thole scale factors trained to target QM molecular dipole moments and polarizabilities. Training of the DNN model used a collection of 39 421 molecules with molecular weights up to 200 Da and containing H, C, N, O, P, S, F, Cl, Br, or I atoms. The DNN model utilizes bond connectivity, including 1,2, 1,3, 1,4, and 1,5 terms and distances of Drude FF atom types as the feature vector to build the model, allowing it to capture both local and nonlocal effects in the molecules. Novel methods have been developed to determine restrained electrostatic potential (RESP) charges on atoms and external points representing lone pairs and to determine Thole scale factors, which have no QM analogue. A penalty scheme is devised as a performance predictor of the trained model. Validation studies show that these DNN models can precisely predict molecular dipole and polarizabilities of Food and Drug Administration (FDA)-approved drugs compared to reference MP2 calculations. The availability of the DNN model allowing for the rapid estimation of the Drude electrostatic parameters will facilitate its applicability to a wider range of molecular species.

RNA-mediated immunotherapy regulating tumor immune microenvironment: next wave of cancer therapeutics.

Accumulating research suggests that the tumor immune microenvironment (TIME) plays an essential role in regulation of tumor growth and metastasis. The cellular and molecular nature of the TIME influences cancer progression and metastasis by altering the ratio of immune- suppressive versus cytotoxic responses in the vicinity of the tumor. Targeting or activating the TIME components show a promising therapeutic avenue to combat cancer. The success of immunotherapy is both astounding and unsatisfactory in the clinic. Advancements in RNA-based technology have improved understanding of the complexity and diversity of the TIME and its effects on therapy. TIME-related RNA or RNA regulators could be promising targets for anticancer immunotherapy. In this review, we discuss the available RNA-based cancer immunotherapies targeting the TIME. More importantly, we summarize the potential of various RNA-based therapeutics clinically available for cancer treatment. RNA-dependent targeting of the TIME, as monotherapy or combined with other evolving therapeutics, might be beneficial for cancer patients' treatment in the near future.

Identification of circRNA-Interacting Proteins by Affinity Pulldown.

Circular RNAs (circRNAs) comprise a vast class of covalently closed transcripts, generated primarily via backsplicing. Most circRNAs arise from full or partial exons, but they can also arise from introns, and from combinations of introns and exons. While high-throughput RNA-sequencing analysis has identified tens of thousands of circRNAs expressed in different tissues and growth conditions, the function of circRNAs has only been described for a handful of them. As most circRNAs appear not to encode peptides, their function is presumed to be linked to their interaction with a range of molecules, particularly other nucleic acids (notably microRNAs) and proteins. A major impediment to identifying circRNA-associated molecules is a lack of suitable methodologies capable of analyzing specifically circRNAs and not their linear RNA counterparts with which they share most of their sequence. Here, we describe a flexible and robust method for identifying the proteins that associate with a given circRNA. The affinity pulldown assay is based on the use of a biotinylated antisense oligomer that recognizes the circRNA-specific junction sequence. Following pulldown using streptavidin beads, the proteins are eluted from the circRNP (circribonucleoprotein) complex and identified by mass spectroscopy; validation by Western blot analysis and other methods would then confirm the identity of the circRNA-associated proteins. We present a detailed step-by-step protocol, tips to optimize the analysis, troubleshooting suggestions, and assistance in interpreting the results. In sum, this protocol enables the discovery of proteins present in circRNPs, a critical effort toward elucidating circRNA function.

Acid ceramidase promotes senescent cell survival.

Cellular senescence is linked to chronic age-related diseases including atherosclerosis, diabetes, and neurodegeneration. Compared to proliferating cells, senescent cells express distinct subsets of proteins. In this study, we used cultured human diploid fibroblasts rendered senescent through replicative exhaustion or ionizing radiation to identify proteins differentially expressed during senescence. We identified acid ceramidase (ASAH1), a lysosomal enzyme that cleaves ceramide into sphingosine and fatty acid, as being highly elevated in senescent cells. This increase in ASAH1 levels in senescent cells was associated with a rise in the levels of ASAH1 mRNA and a robust increase in ASAH1 protein stability. Furthermore, silencing ASAH1 in pre-senescent fibroblasts decreased the levels of senescence proteins p16, p21, and p53, and reduced the activity of the senescence-associated ß-galactosidase. Interestingly, depletion of ASAH1 in pre-senescent cells sensitized these cells to the senolytics Dasatinib and Quercetin (D+Q). Together, our study indicates that ASAH1 promotes senescence, protects senescent cells, and confers resistance against senolytic drugs. Given that inhibiting ASAH1 sensitizes cells towards senolysis, this enzyme represents an attractive therapeutic target in interventions aimed at eliminating senescent cells.

Functional Group Distributions, Partition Coefficients, and Resistance Factors in Lipid Bilayers Using Site Identification by Ligand Competitive Saturation.

Small molecules such as metabolites and drugs must pass through the membrane of the cell, a barrier primarily comprising phospholipid bilayers and embedded proteins. To better understand the process of passive diffusion, knowledge of the ability of various functional groups to partition across bilayers and the associated energetics would be of utility. In the present study, the site identification by ligand competitive saturation (SILCS) methodology has been applied to sample the distributions of a diverse set of chemical solutes representing the functional groups of small molecules across phospholipid bilayers composed of 0.9:0.1 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/cholesterol and a mixture of 0.52:0.18:0.3 1,2-dioleoyl-sn-glycero-3-phospho-l-serine/1,2-dioleoyl-sn-glycero-3-phosphocholine/cholesterol used in parallel artificial membrane permeability assay experiments. A combination of oscillating chemical potential grand canonical Monte Carlo and molecular dynamics in the SILCS simulations was applied to achieve solute sampling through the bilayers and surrounding aqueous environment from which the distribution of solutes and the functional groups they represent were obtained. Results show differential distribution of aliphatic versus aromatic groups with the former having increased sampling in the center of the bilayers versus in the region of the glycerol linker for the latter. Variations in the distribution of different polar groups are evident, with large differences between negative acetate and positive methylammonium with accumulation of the polar-neutral and acetate solutes above the bilayer head groups. Conversion of the distributions to absolute free energies allows for a detailed understanding of energetics of functional groups in different regions of the bilayers and for calculation of absolute free-energy profiles of multifunctional drug-like molecules across the bilayers from which partition coefficients and resistance factors suitable for insertion into the homogenous solubility-diffusion equation for calculation of permeability were obtained. Comparisons of the calculated bilayer/solution partition coefficients with 1-octanol/water experimental data for both drug-like molecules and the solutes show overall good agreement, validating the calculated distributions and associated absolute free-energy profiles.

hnRNPK-regulated LINC00263 promotes malignant phenotypes through miR-147a/CAPN2.

Malignant characteristics of cancers, represented by rapid cell proliferation and high metastatic potential, are a major cause of high cancer-related mortality. As a multifunctional RNA-binding protein, heterogeneous nuclear ribonucleoprotein K (hnRNPK) is closely associated with cancer progression in various types of cancers. In this study, we sought to identify hnRNPK-regulated long intergenic non-coding RNAs (lincRNAs) that play a critical role in the regulation of cancer malignancy. We found that hnRNPK controlled malignant phenotypes including invasiveness, proliferation, and clonogenicity. RNA sequencing and functional studies revealed that LINC00263, a novel target of hnRNPK, is involved in the oncogenic functions of hnRNPK. Knockdown of LINC00263 mitigated the malignant capabilities. Conversely, increased malignant phenotypes were observed in LINC00263-overexpressing cells. Since LINC00263 was mainly localized in the cytosol and highly enriched in Argonaute 2-immunoprecipitation (Ago2-IP), we hypothesized that LINC00263 acts as a competitive endogenous RNA (ceRNA), and thus sought to identify LINC00263-associated microRNAs. Using small RNA sequencing followed by antisense oligonucleotide pull-down, miR-147a was selected for further study. We found that miR-147a negatively regulates LINC00263 via direct interaction, thus suppressing malignant capabilities. Moreover, knockdown of hnRNPK and LINC00263 upregulated miR-147a, indicating that LINC00263 serves as a ceRNA for miR-147a. By analyzing RNA sequencing data and miRNA target prediction, calpain 2 (CAPN2) was identified as a putative target of miR-147a. Ago2-IP and luciferase reporter assay revealed that miR-147a suppressed CAPN2 expression by directly binding to the 3'UTR of CAPN2 mRNA. In addition, we found that the weakened malignant capabilities following knockdown of hnRNPK or LINC00263 were restored by miR-147a inhibition or CAPN2 overexpression. Furthermore, our findings were validated in various other types of cancer cells including lung cancer, colorectal cancer, neuroblastoma, and melanoma. Collectively, we demonstrate that hnRNPK-regulated LINC00263 plays an important role in cancer malignancy by acting as a miR-147a decoy and thus upregulating CAPN2.

Manganese-induced oxidative stress, ultrastructural changes, and proteomics studies in rice plants.

Manganese (Mn) is an essential element for plant growth but it becomes phytotoxic at higher concentrations. The effect of Mn-excess in hydroponics medium was examined on growth, oxidative stress, and ultrastructural changes in chloroplasts and mitochondria as well proteomic alterations in rice (Oryza sativa L.) seedlings. Seedlings grown with 1 mM and 2 mM Mn in nutrient medium for 8 days showed decline in length and fresh biomass, and decline in net photosynthetic rate, transpiration rate, and stomatal conductance. Shoots of the seedlings had higher Mn content than roots. Mn-treated seedlings showed increased production of O2·-, H2O2, and .OH, increased lipid peroxidation, increased carbonylation of proteins, and increased proteolytic activity compared to untreated seedlings. Mn-treated seedlings showed disorganization and swelling of chloroplasts with appearance of plastoglobuli in TEM images and deformity in shape of mitochondria. Using confocal microscopy depolarization of mitochondrial membrane was observed marked by green fluorescence of JC-1 dye monomers in Mn-treated roots. Proteomics studies from leaves of Mn-treated seedlings involving 2DE and PDQuest analysis showed differential expression of 23 proteins, among which MALDI-TOF/TOF mass spectrometry analysis revealed Mn-led downregulation of photosynthesis-related proteins, namely oxygen-evolving complex protein associated with PSII, PAP-3, enzyme involved in protein folding peptidyl-prolyl cis-trans isomerase (PPIase) and carbohydrate metabolizing enzymes hydrolase, fructose-bisphosphate aldolase, transketolase, and isocitrate dehydrogenase, whereas ATP-dependent Clp protease, peroxidase, and nucleic acid-binding proteins were downregulated due to Mn treatment. Results indicate that Mn-excess inhibits growth of rice plants with induction of oxidative stress, causing structural alterations in chloroplasts, mitochondria, inhibiting photosynthesis, and downregulating many photosynthesis and carbohydrate metabolism-related proteins.

Interaction of OIP5-AS1 with MEF2C mRNA promotes myogenic gene expression.

Long noncoding (lnc)RNAs potently regulate gene expression programs in physiology and disease. Here, we describe a key function for lncRNA OIP5-AS1 in myogenesis, the process whereby myoblasts differentiate into myotubes during muscle development and muscle regeneration after injury. In human myoblasts, OIP5-AS1 levels increased robustly early in myogenesis, and its loss attenuated myogenic differentiation and potently reduced the levels of the myogenic transcription factor MEF2C. This effect relied upon the partial complementarity of OIP5-AS1 with MEF2C mRNA and the presence of HuR, an RNA-binding protein (RBP) with affinity for both transcripts. Remarkably, HuR binding to MEF2C mRNA, which stabilized MEF2C mRNA and increased MEF2C abundance, was lost after OIP5-AS1 silencing, suggesting that OIP5-AS1 might serve as a scaffold to enhance HuR binding to MEF2C mRNA, in turn increasing MEF2C production. These results highlight a mechanism whereby a lncRNA promotes myogenesis by enhancing the interaction of an RBP and a myogenic mRNA.

Further Optimization and Validation of the Classical Drude Polarizable Protein Force Field.

The CHARMM Drude-2013 polarizable force field (FF) was developed to include the explicit treatment of induced electronic polarizability, resulting in a more accurate description of the electrostatic interactions in molecular dynamics (MD) simulations. While the Drude-2013 protein FF has shown success in improving the folding properties of α-helical peptides and to reproduce experimental observables in simulations up to 1 µs, some limitations were noted regarding the stability of ß-sheet structures in simulations longer than 100 ns as well as larger deviations from crystal structures in simulations of a number of proteins compared to the additive CHARMM36 protein FF. The origin of the instability has been identified and appears to be primarily due to overestimated atomic polarizabilities and induced dipole-dipole interactions on the Cß, Cγ, and Cδ side chain atoms. To resolve this and other issues, a number of aspects of the model were revisited, resulting in Drude-2019 protein FF. Backbone parameters were optimized targeting the conformational properties of the (Ala)5 peptide in solution along with gas phase properties of the alanine dipeptide. Dipeptides that contain N-acetylated and N'-methylamidated termini, excluding Gly, Pro, and Ala, were used as models to optimize the atomic polarizabilities and Thole screening factors on selected Cß, Cγ, and Cδ carbons by targeting quantum mechanical (QM) dipole moments and molecular polarizabilities. In addition, to obtain better conformational properties, side chain χ1 and χ2 dihedral parameters were optimized targeting QM data for the respective side chain dipeptide conformations as well as Protein Data Bank survey data based on the χ1, χ2 sampling from Hamiltonian replica-exchange MD simulations of (Ala)4-X-(Ala)4 in solution, where X is the amino acid of interest. Further improvements include optimizing nonbonded interactions between charged residues to reproduce QM interaction energies of the charged-protein model compounds and experimental osmotic pressures. Validation of the optimized Drude protein FF includes MD simulations of a collection of peptides and proteins including ß-sheet structures, as well as transmembrane ion channels. Results showed that the updated Drude-2019 protein FF yields smaller overall root-mean-square differences of proteins as compared to the additive CHARMM36m and Drude-2013 FFs as well as similar or improved agreement with experimental NMR properties, allowing for long time scale simulation studies of proteins and more complex biomolecular systems in conjunction with the remainder of the Drude polarizable FF.

Translation of insulin granule proteins are regulated by PDI and PABP.

Glucose mediated insulin biosynthesis is tightly regulated and shared between insulin granule proteins such as its processing enzymes, prohormone convertases, PC1/3 and PC2. However, the molecular players involved in the co-ordinated translation remain elusive. The trans-acting factors like PABP (Poly A Binding Protein) and PDI (Protein Disulphide Isomerize) binds to a conserved sequence in the 5'UTR of insulin mRNA and regulates its translation. Here, we demonstrate that 5'UTR of PC1/3 and PC2 also associate with PDI and PABP. We show that a' and RRM 3-4 domains of PDI and PABP respectively, are necessary for RNA binding activity to the 5'UTRs of insulin and its processing enzymes.

circSamd4 represses myogenic transcriptional activity of PUR proteins.

By interacting with proteins and nucleic acids, the vast family of mammalian circRNAs is proposed to influence many biological processes. Here, RNA sequencing analysis of circRNAs differentially expressed during myogenesis revealed that circSamd4 expression increased robustly in mouse C2C12 myoblasts differentiating into myotubes. Moreover, silencing circSamd4, which is conserved between human and mouse, delayed myogenesis and lowered the expression of myogenic markers in cultured myoblasts from both species. Affinity pulldown followed by mass spectrometry revealed that circSamd4 associated with PURA and PURB, two repressors of myogenesis that inhibit transcription of the myosin heavy chain (MHC) protein family. Supporting the hypothesis that circSamd4 might complex with PUR proteins and thereby prevent their interaction with DNA, silencing circSamd4 enhanced the association of PUR proteins with the Mhc promoter, while overexpressing circSamd4 interfered with the binding of PUR proteins to the Mhc promoter. These effects were abrogated when using a mutant circSamd4 lacking the PUR binding site. Our results indicate that the association of PUR proteins with circSamd4 enhances myogenesis by contributing to the derepression of MHC transcription.

Methods for analysis of circular RNAs.

Eukaryotic cells express a myriad of circular RNAs (circRNAs), many of them displaying tissue-specific expression patterns. They arise from linear precursor RNAs in which 5' and 3' ends become covalently ligated. Given these features, biochemical and computational approaches traditionally used to study linear RNA must be adapted for analysis of circular RNAs. Such circRNA-specific methodologies are allowing the systematic identification of circRNAs and the analysis of their biological functions. Here, we review the resources and molecular methods currently utilized to quantify circRNAs, visualize their distribution, identify interacting partners, and elucidate their function. We discuss the challenges of analyzing circRNAs and propose alternative approaches for studying this unique class of transcripts. This article is characterized under: RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Methods > RNA Analyses in vitro and In Silico RNA Methods > RNA Analyses in Cells.

Drude Polarizable Force Field Parametrization of Carboxylate and N-Acetyl Amine Carbohydrate Derivatives.

In this work, we report the development of Drude polarizable force field parameters for the carboxylate and N-acetyl amine derivatives, extending the functionality of the existing Drude polarizable carbohydrate force field. The force field parameters have been developed in a hierarchical manner, reproducing the quantum mechanical gas-phase properties of small model compounds representing the key functional group in the carbohydrate derivatives, including optimization of the electrostatic and bonded parameters. The optimized parameters were then used to generate the models for carboxylate and N-acetyl amine carbohydrate derivatives. The transferred parameters were further tested and optimized to reproduce crystal geometries and J-coupling data from nuclear magnetic resonance experiments. The parameter development resulted in the incorporation of d-glucuronate, l-iduronate, N-acetyl-d-glucosamine (GlcNAc), and N-acetyl-d-galactosamine (GalNAc) sugars into the Drude polarizable force field. The parameters developed in this study were then applied to study the conformational properties of glycosaminoglycan polymer hyaluronan, composed of d-glucuronate and N-acetyl-d-glucosamine, in aqueous solution. Upon comparing the results from the additive and polarizable simulations, it was found that the inclusion of polarization improved the description of the electrostatic interactions observed in hyaluronan, resulting in enhanced conformational flexibility. The developed Drude polarizable force field parameters in conjunction with the remainder of the Drude polarizable force field parameters can be used for future studies involving carbohydrates and their conjugates in complex, heterogeneous systems.

Elucidating the role of key structural motifs in antifreeze glycoproteins.

Antifreeze glycoproteins (AFGPs) are distinctively riveting class of bio-macromolecules, which endows the survival of organisms inhabiting polar and subpolar regions. These proteins are believed to hinder microscopic freezing by interacting with embryonic ice crystals and precluding their further growth. The underlying molecular mechanism by which AFGPs bind to ice has remained elusive due to insufficient structural characterization, with conflicting hypotheses on the possible binding mode of AFGPs - either via the hydrophobic peptide backbone or via the hydrophilic carbohydrate side chains - when interacting with ice. Chemical synthesis has allowed researchers to access synthetic variants of natural AFGPs. These studies revealed that AFGPs exhibit huge variations in their thermal hysteresis and ice shaping behavior with only slight structural variations, especially to the carbohydrate side chains. Four key structural motifs were identified as crucial to AFGP activity: the presence of a threonine γ-methyl group, an α-glycosidic carbohydrate-protein linkage, an acetylamide group (-NHCOCH3) at the C2 position of the carbohydrate linked to the protein, and the presence of carbohydrate hydroxyl groups. In this study, we use molecular dynamics (MD) simulations to probe the microscopic properties of water accompanying these structural variations of AFGPs. We find that these variations primarily influence the conformation space of AFGPs and also crucially control their hydration dynamics. Owing to the disordered nature of AFGPs we use Markov-state modeling to identify the conformational preferences of AFGPs. The simulations reveal the importance of steric bulk, intra-molecular carbohydrate-protein H-bonds and conformational preferences (α- vs. ß-linkages) in controlling the spatial segregation of the hydrophilic and hydrophobic regions of AFGPs. We hypothesize that the hydrophobic component of AFGPs is crucial to their binding to ice, which determines the ice shaping ability of AFGPs. However, the hydrophilic carbohydrate hydroxyl groups and their ability to form water bridges control the subsequent hydration dynamics, which is key to the antifreeze properties. Investigating the tetrahedral order parameter of water molecules around the carbohydrates revealed competition between solute- and bulk-influenced solvent structures, with maximum restructuring being observed in the interfacial region 2.5-4.5 Å away from the AFGPs.