Phosphorylation of the compartmentalized PKA substrate TAF15 regulates RNA-protein interactions.

Spatiotemporal-controlled second messengers alter molecular interactions of central signaling nodes for ensuring physiological signal transmission. One prototypical second messenger molecule which modulates kinase signal transmission is the cyclic-adenosine monophosphate (cAMP). The main proteinogenic cellular effectors of cAMP are compartmentalized protein kinase A (PKA) complexes. Their cell-type specific compositions precisely coordinate substrate phosphorylation and proper signal propagation which is indispensable for numerous cell-type specific functions. Here we present evidence that TAF15, which is implicated in the etiology of amyotrophic lateral sclerosis, represents a novel nuclear PKA substrate. In cross-linking and immunoprecipitation experiments (iCLIP) we showed that TAF15 phosphorylation alters the binding to target transcripts related to mRNA maturation, splicing and protein-binding related functions. TAF15 appears to be one of multiple PKA substrates that undergo RNA-binding dynamics upon phosphorylation. We observed that the activation of the cAMP-PKA signaling axis caused a change in the composition of a collection of RNA species that interact with TAF15. This observation appears to be a broader principle in the regulation of molecular interactions, as we identified a significant enrichment of RNA-binding proteins within endogenous PKA complexes. We assume that phosphorylation of RNA-binding domains adds another layer of regulation to binary protein-RNAs interactions with consequences to RNA features including binding specificities, localization, abundance and composition.

Membrane platform protein PulF of the Klebsiella type II secretion system forms a trimeric ion channel essential for endopilus assembly and protein secretion.

IMPORTANCE: Type IV pili and type II secretion systems are members of the widespread type IV filament (T4F) superfamily of nanomachines that assemble dynamic and versatile surface fibers in archaea and bacteria. The assembly and retraction of T4 filaments with diverse surface properties and functions require the plasma membrane platform proteins of the GspF/PilC superfamily. Generally considered dimeric, platform proteins are thought to function as passive transmitters of the mechanical energy generated by the ATPase motor, to somehow promote insertion of pilin subunits into the nascent pilus fibers. Here, we generate and experimentally validate structural predictions that support the trimeric state of a platform protein PulF from a type II secretion system. The PulF trimers form selective proton or sodium channels which might energize pilus assembly using the membrane potential. The conservation of the channel sequence and structural features implies a common mechanism for all T4F assembly systems. We propose a model of the oligomeric PulF-PulE ATPase complex that provides an essential framework to investigate and understand the pilus assembly mechanism.

Dengue and Zika RNA-RNA interactomes reveal pro- and anti-viral RNA in human cells.

BACKGROUND: Identifying host factors is key to understanding RNA virus pathogenicity. Besides proteins, RNAs can interact with virus genomes to impact replication. RESULTS: Here, we use proximity ligation sequencing to identify virus-host RNA interactions for four strains of Zika virus (ZIKV) and one strain of dengue virus (DENV-1) in human cells. We find hundreds of coding and non-coding RNAs that bind to DENV and ZIKV viruses. Host RNAs tend to bind to single-stranded regions along the virus genomes according to hybridization energetics. Compared to SARS-CoV-2 interactors, ZIKV-interacting host RNAs tend to be downregulated upon virus infection. Knockdown of several short non-coding RNAs, including miR19a-3p, and 7SK RNA results in a decrease in viral replication, suggesting that they act as virus-permissive factors. In addition, the 3'UTR of DYNLT1 mRNA acts as a virus-restrictive factor by binding to the conserved dumbbell region on DENV and ZIKV 3'UTR to decrease virus replication. We also identify a conserved set of host RNAs that interacts with DENV, ZIKV, and SARS-CoV-2, suggesting that these RNAs are broadly important for RNA virus infection. CONCLUSIONS: This study demonstrates that host RNAs can impact virus replication in permissive and restrictive ways, expanding our understanding of host factors and RNA-based gene regulation during viral pathogenesis.

Flavivirus genome recoding by codon optimisation confers genetically stable in vivo attenuation in both mice and mosquitoes.

Virus genome recoding is an attenuation method that confers genetically stable attenuation by rewriting a virus genome with numerous silent mutations. Prior flavivirus genome recoding attempts utilised codon deoptimisation approaches. However, these codon deoptimisation approaches act in a species dependent manner and were unable to confer flavivirus attenuation in mosquito cells or in mosquito animal models. To overcome these limitations, we performed flavivirus genome recoding using the contrary approach of codon optimisation. The genomes of flaviviruses such as dengue virus type 2 (DENV2) and Zika virus (ZIKV) contain functional RNA elements that regulate viral replication. We hypothesised that flavivirus genome recoding by codon optimisation would introduce silent mutations that disrupt these RNA elements, leading to decreased replication efficiency and attenuation. We chose DENV2 and ZIKV as representative flaviviruses and recoded them by codon optimising their genomes for human expression. Our study confirms that this recoding approach of codon optimisation does translate into reduced replication efficiency in mammalian, human, and mosquito cells as well as in vivo attenuation in both mice and mosquitoes. In silico modelling and RNA SHAPE analysis confirmed that DENV2 recoding resulted in the extensive disruption of genomic structural elements. Serial passaging of recoded DENV2 resulted in the emergence of rescue or adaptation mutations, but no reversion mutations. These rescue mutations were unable to rescue the delayed replication kinetics and in vivo attenuation of recoded DENV2, demonstrating that recoding confers genetically stable attenuation. Therefore, our recoding approach is a reliable attenuation method with potential applications for developing flavivirus vaccines.

Whole Genomic Characterization of Streptococcus iniae Isolates from Barramundi (Lates calcarifer) and Preliminary Evidence of Cross-Protective Immunization.

Lates calcarifer, also known as Barramundi or Asian seabass, is a highly productive and fast-growing species that is well suited to large-scale aquaculture due to its attractive harvestable yields (premium fish). This fish has been envisioned as having the potential to be the "Salmon of Tropics". Cultivating Lates calcarifer in aquaculture poses challenges, as the dense populations that make such aquaculture commercially viable facilitate the rapid spread of infectious diseases, which in turn significantly impact yield. Hence, the immunization of juveniles is necessary, and the development of new immunization agents enhances the efficiency of aquaculture and improves food security. In our study, we characterize seven novel strains of the bacterial pathogen Streptococcus iniae that were collected from commercial fish farms in Singapore and Australia. We find that the capsular operon in our strains is highly conserved and identify a number of major surface antigens previously described in Streptococcus. A genome analysis indicates that the present strains are closely related but form distinct strains within the S. iniae species. We then proceed to demonstrate that inoculation with the inactivated strain P3SAB cross-protects Lates calcarifer against S. iniae infections in vivo from a variety of strains found in both Singapore and Australia.

Tracking and blocking interdependencies of cellular BRAF-MEK oncokinase activities.

The selective targeting of mutated kinases in cancer therapies has the potential to improve therapeutic success and thereby the survival of patients. In the case of melanoma, the constitutively active MAPK pathway is targeted by a combinatorial inhibition of BRAF and MEK activities. These MAPK pathway players may display patient-specific differences in the onco-kinase mutation spectrum, which needs to be considered for the design of more efficient personalized therapies. Here, we extend a bioluminescence-based kinase conformation biosensor (KinCon) to allow for live-cell tracking of interconnected kinase activity states. First, we show that common MEK1 patient mutations promote a structural rearrangement of the kinase to an opened and active conformation. This effect was reversible by the binding of MEK inhibitors to mutated MEK1, as shown in biosensor assays and molecular dynamics simulations. Second, we implement a novel application of the KinCon technology for tracking the simultaneous, vertical targeting of the two functionally linked kinases BRAF and MEK1. Thus, we demonstrate that, in the presence of constitutively active BRAF-V600E, specific inhibitors of both kinases are efficient in driving MEK1 into a closed, inactive conformation state. We compare current melanoma treatments and show that combinations of BRAFi and MEKi display a more pronounced structural change of the drug sensor than the respective single agents, thereby identifying synergistic effects among these drug combinations. In summary, we depict the extension of the KinCon biosensor technology to systematically validate, anticipate, and personalize tailored drug arrangements using a multiplexed setup.

Dengue Virus Capsid Protein Facilitates Genome Compaction and Packaging.

Dengue virus (DENV) is a single-stranded (+)-sense RNA virus that infects humans and mosquitoes, posing a significant health risk in tropical and subtropical regions. Mature virions are composed of an icosahedral shell of envelope (E) and membrane (M) proteins circumscribing a lipid bilayer, which in turn contains a complex of the approximately 11 kb genomic RNA with capsid (C) proteins. Whereas the structure of the envelope is clearly defined, the structure of the packaged genome in complex with C proteins remains elusive. Here, we investigated the interactions of C proteins with viral RNA, in solution and inside mature virions, via footprinting and cross-linking experiments. We demonstrated that C protein interaction with DENV genomes saturates at an RNA:C protein ratio below 1:250. Moreover, we also showed that the length of the RNA genome interaction sites varies, in a multimodal distribution, consistent with the C protein binding to each RNA site mostly in singlets or pairs (and, in some instances, higher numbers). We showed that interaction sites are preferentially sites with low base pairing, as previously measured by 2'-acetylation analyzed by primer extension (SHAPE) reactivity indicating structuredness. We found a clear association pattern emerged: RNA-C protein binding sites are strongly associated with long-range RNA-RNA interaction sites, particularly inside virions. This, in turn, explains the need for C protein in viral genome packaging: the protein has a chief role in coordinating these key interactions, promoting proper packaging of viral RNA. Such sites are, thus, highly consequential for viral assembly, and, as such, may be targeted in future drug development strategies against these and related viruses.

Identification of potential pathways and biomarkers linked to progression in ALS.

OBJECTIVE: To identify potential diagnostic and prognostic biomarkers for clinical management and clinical trials in amyotrophic lateral sclerosis. METHODS: We analysed proteomics data of ALS patient-induced pluripotent stem cell-derived motor neurons available through the AnswerALS consortium. After stratifying patients using clinical ALSFRS-R and ALS-CBS scales, we identified differentially expressed proteins indicative of ALS disease severity and progression rate as candidate ALS-related and prognostic biomarkers. Pathway analysis for identified proteins was performed using STITCH. Protein sets were correlated with the effects of drugs using the Connectivity Map tool to identify compounds likely to affect similar pathways. RNAi screening was performed in a Drosophila TDP-43 ALS model to validate pathological relevance. A statistical classification machine learning model was constructed using ridge regression that uses proteomics data to differentiate ALS patients from controls. RESULTS: We identified 76, 21, 71 and 1 candidate ALS-related biomarkers and 22, 41, 27 and 64 candidate prognostic biomarkers from patients stratified by ALSFRS-R baseline, ALSFRS-R progression slope, ALS-CBS baseline and ALS-CBS progression slope, respectively. Nineteen proteins enhanced or suppressed pathogenic eye phenotypes in the ALS fly model. Nutraceuticals, dopamine pathway modulators, statins, anti-inflammatories and antimicrobials were predicted starting points for drug repurposing using the connectivity map tool. Ten diagnostic biomarker proteins were predicted by machine learning to identify ALS patients with high accuracy and sensitivity. INTERPRETATION: This study showcases the powerful approach of iPSC-motor neuron proteomics combined with machine learning and biological confirmation in the prediction of novel mechanisms and diagnostic and predictive biomarkers in ALS.

Computational and Experimental Approaches to Study the RNA Secondary Structures of RNA Viruses.

Most pandemics of recent decades can be traced to RNA viruses, including HIV, SARS, influenza, dengue, Zika, and SARS-CoV-2. These RNA viruses impose considerable social and economic burdens on our society, resulting in a high number of deaths and high treatment costs. As these RNA viruses utilize an RNA genome, which is important for different stages of the viral life cycle, including replication, translation, and packaging, studying how the genome folds is important to understand virus function. In this review, we summarize recent advances in computational and high-throughput RNA structure-mapping approaches and their use in understanding structures within RNA virus genomes. In particular, we focus on the genome structures of the dengue, Zika, and SARS-CoV-2 viruses due to recent significant outbreaks of these viruses around the world.

Genome-wide RNA structure changes during human neurogenesis modulate gene regulatory networks.

The distribution, dynamics, and function of RNA structures in human development are under-explored. Here, we systematically assayed RNA structural dynamics and their relationship with gene expression, translation, and decay during human neurogenesis. We observed that the human ESC transcriptome is globally more structurally accessible than differentiated cells and undergoes extensive RNA structure changes, particularly in the 3' UTR. Additionally, RNA structure changes during differentiation are associated with translation and decay. We observed that RBP and miRNA binding is associated with RNA structural changes during early neuronal differentiation, and splicing is associated during later neuronal differentiation. Furthermore, our analysis suggests that RBPs are major factors in structure remodeling and co-regulate additional RBPs and miRNAs through structure. We demonstrated an example of this by showing that PUM2-induced structure changes on LIN28A enable miR-30 binding. This study deepens our understanding of the widespread and complex role of RNA-based gene regulation during human development.

Comprehensive mapping of SARS-CoV-2 interactions in vivo reveals functional virus-host interactions.

SARS-CoV-2 is a major threat to global health. Here, we investigate the RNA structure and RNA-RNA interactions of wildtype (WT) and a mutant (Δ382) SARS-CoV-2 in cells using Illumina and Nanopore platforms. We identify twelve potentially functional structural elements within the SARS-CoV-2 genome, observe that subgenomic RNAs can form different structures, and that WT and Δ382 virus genomes fold differently. Proximity ligation sequencing identify hundreds of RNA-RNA interactions within the virus genome and between the virus and host RNAs. SARS-CoV-2 genome binds strongly to mitochondrial and small nucleolar RNAs and is extensively 2'-O-methylated. 2'-O-methylation sites are enriched in viral untranslated regions, associated with increased virus pair-wise interactions, and are decreased in host mRNAs upon virus infection, suggesting that the virus sequesters methylation machinery from host RNAs towards its genome. These studies deepen our understanding of the molecular and cellular basis of SARS-CoV-2 pathogenicity and provide a platform for targeted therapy.

The nanotube express: Delivering a stapled peptide to the cell surface.

HYPOTHESIS: Carbon nanotubes (CNTs) represent a novel platform for cellular delivery of therapeutic peptides. Chemically-functionalized CNTs may enhance peptide uptake by improving their membrane targeting properties. EXPERIMENTS: Using coarse-grained (CG) molecular dynamics (MD) simulations, we investigate membrane interactions of a peptide conjugated to pristine and chemically-modified CNTs. As proof of principle, we focus on their interactions with PM2, an amphipathic stapled peptide that inhibits the E3 ubiquitin ligase HDM2 from negatively regulating the p53 tumor suppressor. CNT interaction with both simple planar lipid bilayers as well as spherical lipid vesicles was studied, the latter as a surrogate for curved cellular membranes. FINDINGS: Membrane permeation was rapid and spontaneous for both pristine and oxidized CNTs when unconjugated. This was slowed upon addition of a noncovalently attached peptide surface "sheath", which may be an effective way to slow CNT entry and avert membrane rupture. The CNT conjugates were observed to "desheath" their peptide layer at the bilayer interface upon insertion, leaving their cargo behind in the outer leaflet. This suggests that a synergy may exist to optimize CNT safety whilst enhancing the delivery efficiency of "hitchhiking" therapeutic molecules.

Allosteric Kinase Inhibitors Reshape MEK1 Kinase Activity Conformations in Cells and In Silico.

Mutations at different stages of the mitogen-activated protein kinase (MAPK) signaling pathway lead to aberrant activation of the involved protein kinase entities. These oncogenic modifications alter signal propagation which converge on the gatekeeper kinases MEK1/2, transmitting the input signal to ERK1/2. Thus, targeted MEK inhibition causes qualitative alterations of carcinogenic MAPK signals. Phosphorylation of the MEK1 activation loop at the positions S218 and S222 by RAF kinases triggers the conformational alignment of MEK's catalytic pocket to enable ATP-binding and substrate phosphorylation. We have extended a kinase conformation (KinCon) biosensor platform to record MEK1 activity dynamics. In addition to MEK phosphorylation by BRAF, the integration of the phosphorylation-mimetic mutations S218D/S222D triggered opening of the kinase. Structural rearrangement may involve the flexibility of the N terminal MEK1 A-helix. Application of the allosterically acting MEK inhibitors (MEKi) trametinib, cobimentinib, refametinib, and selumetinib converted activated MEK1 KinCon reporters back into a more closed inactive conformation. We confirmed MEK1 KinCon activity dynamics upon drug engagement using the patient-derived melanoma cell line A2058, which harbors the V600E hotspot BRAF mutation. In order to confirm biosensor dynamics, we simulated structure dynamics of MEK1 kinase in the presence and absence of mutations and/or MEKi binding. We observed increased dynamics for the S218D/S222D double mutant particularly in the region of the distal A-helix and alpha-C helix. These data underline that MEK1 KinCon biosensors have the potential to be subjected to MEKi efficacy validations in an intact cell setting.

The Pseudo-Circular Genomes of Flaviviruses: Structures, Mechanisms, and Functions of Circularization.

The circularization of viral genomes fulfills various functions, from evading host defense mechanisms to promoting specific replication and translation patterns supporting viral proliferation. Here, we describe the genomic structures and associated host factors important for flaviviruses genome circularization and summarize their functional roles. Flaviviruses are relatively small, single-stranded, positive-sense RNA viruses with genomes of approximately 11 kb in length. These genomes contain motifs at their 5' and 3' ends, as well as in other regions, that are involved in circularization. These motifs are highly conserved throughout the Flavivirus genus and occur both in mature virions and within infected cells. We provide an overview of these sequence motifs and RNA structures involved in circularization, describe their linear and circularized structures, and discuss the proteins that interact with these circular structures and that promote and regulate their formation, aiming to clarify the key features of genome circularization and understand how these affect the flaviviruses life cycle.

Antibody affinity versus dengue morphology influences neutralization.

Different strains within a dengue serotype (DENV1-4) can have smooth, or "bumpy" surface morphologies with different antigenic characteristics at average body temperature (37°C). We determined the neutralizing properties of a serotype cross-reactive human monoclonal antibody (HMAb) 1C19 for strains with differing morphologies within the DENV1 and DENV2 serotypes. We mapped the 1C19 epitope to E protein domain II by hydrogen deuterium exchange mass spectrometry, cryoEM and molecular dynamics simulations, revealing that this epitope is likely partially hidden on the virus surface. We showed the antibody has high affinity for binding to recombinant DENV1 E proteins compared to those of DENV2, consistent with its strong neutralizing activities for all DENV1 strains tested regardless of their morphologies. This finding suggests that the antibody could out-compete E-to-E interaction for binding to its epitope. In contrast, for DENV2, HMAb 1C19 can only neutralize when the epitope becomes exposed on the bumpy-surfaced particle. Although HMAb 1C19 is not a suitable therapeutic candidate, this study with HMAb 1C19 shows the importance of choosing a high-affinity antibody that could neutralize diverse dengue virus morphologies for therapeutic purposes.

Computational modelling of flavivirus dynamics: The ins and outs.

Enveloped viruses such as the flaviviruses represent a significant burden to human health around the world, with hundreds of millions of people each year affected by dengue alone. In an effort to improve our understanding of the molecular basis for the infective mechanisms of these viruses, extensive computational modelling approaches have been applied to elucidate their conformational dynamics. Multiscale protocols have been developed to simulate flavivirus envelopes in close accordance with biophysical data, in particular derived from cryo-electron microscopy, enabling high-resolution refinement of their structures and elucidation of the conformational changes associated with adaptation both to host environments and to immunological factors such as antibodies. Likewise, integrative modelling efforts combining data from biophysical experiments and from genome sequencing with chemical modification are providing unparalleled insights into the architecture of the previously unresolved nucleocapsid complex. Collectively, this work provides the basis for the future rational design of new antiviral therapeutics and vaccine development strategies targeting enveloped viruses.

Mutation-oriented profiling of autoinhibitory kinase conformations predicts RAF inhibitor efficacies.

Kinase-targeted therapies have the potential to improve the survival of patients with cancer. However, the cancer-specific spectrum of kinase alterations exhibits distinct functional properties and requires mutation-oriented drug treatments. Besides post-translational modifications and diverse intermolecular interactions of kinases, it is the distinct disease mutation which reshapes full-length kinase conformations, affecting their activity. Oncokinase mutation profiles differ between cancer types, as it was shown for BRAF in melanoma and non-small-cell lung cancers. Here, we present the target-oriented application of a kinase conformation (KinCon) reporter platform for live-cell measurements of autoinhibitory kinase activity states. The bioluminescence-based KinCon biosensor allows the tracking of conformation dynamics of full-length kinases in intact cells and real time. We show that the most frequent BRAF cancer mutations affect kinase conformations and thus the engagement and efficacy of V600E-specific BRAF inhibitors (BRAFi). We illustrate that the patient mutation harboring KinCon reporters display differences in the effectiveness of the three clinically approved BRAFi vemurafenib, encorafenib, and dabrafenib and the preclinical paradox breaker PLX8394. We confirmed KinCon-based drug efficacy predictions for BRAF mutations other than V600E in proliferation assays using patient-derived lung cancer cell lines and by analyzing downstream kinase signaling. The systematic implementation of such conformation reporters will allow to accelerate the decision process for the mutation-oriented RAF-kinase cancer therapy. Moreover, we illustrate that the presented kinase reporter concept can be extended to other kinases which harbor patient mutations. Overall, KinCon profiling provides additional mechanistic insights into full-length kinase functions by reporting protein-protein interaction (PPI)-dependent, mutation-specific, and drug-driven changes of kinase activity conformations.

Molecular simulations unravel the molecular principles that mediate selective permeability of carboxysome shell protein.

Bacterial microcompartments (BMCs) are nanoscale proteinaceous organelles that encapsulate enzymes from the cytoplasm using an icosahedral protein shell that resembles viral capsids. Of particular interest are the carboxysomes (CBs), which sequester the CO2-fixing enzymes ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) to enhance carbon assimilation. The carboxysome shell serves as a semi-permeable barrier for passage of metabolites in and out of the carboxysome to enhance CO2 fixation. How the protein shell directs influx and efflux of molecules in an effective manner has remained elusive. Here we use molecular dynamics and umbrella sampling calculations to determine the free-energy profiles of the metabolic substrates, bicarbonate, CO2 and ribulose bisphosphate and the product 3-phosphoglycerate associated with their transition through the major carboxysome shell protein CcmK2. We elucidate the electrostatic charge-based permeability and key amino acid residues of CcmK2 functioning in mediating molecular transit through the central pore. Conformational changes of the loops forming the central pore may also be required for transit of specific metabolites. The importance of these in-silico findings is validated experimentally by site-directed mutagenesis of the key CcmK2 residue Serine 39. This study provides insight into the mechanism that mediates molecular transport through the shells of carboxysomes, applicable to other BMCs. It also offers a predictive approach to investigate and manipulate the shell permeability, with the intent of engineering BMC-based metabolic modules for new functions in synthetic biology.

The Molecular Basis for Purine Binding Selectivity in the Bacterial ATP Synthase ϵ Subunit.

The ϵ subunit of ATP synthases has been proposed to regulate ATP hydrolysis in bacteria. Prevailing evidence supports the notion that when the ATP concentration falls below a certain threshold, the ϵ subunit changes its conformation from a non-inhibitory down-state to an extended up-state that then inhibits enzymatic ATP hydrolysis by binding to the catalytic domain. It has been demonstrated that the ϵ subunit from Bacillus PS3 is selective for ATP over other nucleotides, including GTP. In this study, the purine triphosphate selectivity is rationalized by using results from MD simulations and free energy calculations for the R103A/R115A mutant of the ϵ subunit from Bacillus PS3, which binds ATP more strongly than the wild-type protein. Our results are in good agreement with experimental data, and the elucidated molecular basis for selectivity could help to guide the design of novel GTP sensors.

Extending the Martini Coarse-Grained Force Field to N-Glycans.

Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.