Molecular dynamics simulations of glycoproteins using CHARMM

Docking and Molecular Dynamics Simulations Clarify Binding Sites for Interactions of Novel Marine Sulfated Glycans with SARS-CoV-2 Spike Glycoprotein

The entry of SARS-CoV-2 into the host cell is mediated by its S-glycoprotein (SGP). Sulfated glycans bind to the SGP receptor-binding domain (RBD), which forms a ternary complex with its receptor angiotensin converting enzyme 2. Here, we have conducted a thorough and systematic computational study of the binding of four oligosaccharide building blocks from novel marine sulfated glycans (isolated from Pentacta pygmaea and Isostichopus badionotus) to the non-glycosylated and glycosylated RBD. Blind docking studies using three docking programs identified five potential cryptic binding sites. Extensive site-targeted docking and molecular dynamics simulations using two force fields confirmed only two binding sites (Sites 1 and 5) for these novel, highly charged sulfated glycans, which were also confirmed by previously published reports. This work showed the structural features and key interactions driving ligand binding. A previous study predicted Site 2 to be a potential binding site, which was not observed here. The use of several molecular modeling approaches gave a comprehensive assessment. The detailed comparative study utilizing multiple modeling approaches is the first of its kind for novel glycan-SGP interaction characterization. This study provided insights into the key structural features of these novel glycans as they are considered for development as potential therapeutics.

Orientational Preferences of GPI-Anchored Ly6/uPAR Proteins

Ly6/uPAR proteins regulate many essential functions in the nervous and immune systems and epithelium. Most of these proteins contain single β-structural LU domains with three protruding loops and are glycosylphosphatidylinositol (GPI)-anchored to a membrane. The GPI-anchor role is currently poorly studied. Here, we investigated the positional and orientational preferences of six GPI-anchored proteins in the receptor-unbound state by molecular dynamics simulations. Regardless of the linker length between the LU domain and GPI-anchor, the proteins interacted with the membrane by polypeptide parts and N-/O-glycans. Lynx1, Lynx2, Lypd6B, and Ly6H contacted the membrane by the loop regions responsible for interactions with nicotinic acetylcholine receptors, while Lypd6 and CD59 demonstrated unique orientations with accessible receptor-binding sites. Thus, GPI-anchoring does not guarantee an optimal 'pre-orientation' of the LU domain for the receptor interaction.

Molecular simulations of complex carbohydrates and glycoconjugates

Complex carbohydrates (glycans) are the most abundant and versatile biopolymers in nature. The broad diversity of biochemical functions that carbohydrates cover is a direct consequence of the variety of 3D architectures they can adopt, displaying branched or linear arrangements, widely ranging in sizes, and with the highest diversity of building blocks of any other natural biopolymer. Despite this unparalleled complexity, a common denominator can be found in the glycans' inherent flexibility, which hinders experimental characterization, but that can be addressed by high-performance computing (HPC)-based molecular simulations. In this short review, I present and discuss the state-of-the-art of molecular simulations of complex carbohydrates and glycoconjugates, highlighting methodological strengths and weaknesses, important insights through emblematic case studies, and suggesting perspectives for future developments.

An Integrated Mass Spectrometry-Based Glycomics-Driven Glycoproteomics Analytical Platform to Functionally Characterize Glycosylation Inhibitors

Cancer progression is linked to aberrant protein glycosylation due to the overexpression of several glycosylation enzymes. These enzymes are underexploited as potential anticancer drug targets and the development of rapid-screening methods and identification of glycosylation inhibitors are highly sought. An integrated bioinformatics and mass spectrometry-based glycomics-driven glycoproteomics analysis pipeline was performed to identify an N-glycan inhibitor against lung cancer cells. Combined network pharmacology and in silico screening approaches were used to identify a potential inhibitor, pictilisib, against several glycosylation-related proteins, such as Alpha1-6FucT, GlcNAcT-V, and Alpha2,6-ST-I. A glycomics assay of lung cancer cells treated with pictilisib showed a significant reduction in the fucosylation and sialylation of N-glycans, with an increase in high mannose-type glycans. Proteomics analysis and in vitro assays also showed significant upregulation of the proteins involved in apoptosis and cell adhesion, and the downregulation of proteins involved in cell cycle regulation, mRNA processing, and protein translation. Site-specific glycoproteomics analysis further showed that glycoproteins with reduced fucosylation and sialylation were involved in apoptosis, cell adhesion, DNA damage repair, and chemical response processes. To determine how the alterations in N-glycosylation impact glycoprotein dynamics, modeling of changes in glycan interactions of the ITGA5-ITGB1 (Integrin alpha 5-Integrin beta-1) complex revealed specific glycosites at the interface of these proteins that, when highly fucosylated and sialylated, such as in untreated A549 cells, form greater hydrogen bonding interactions compared to the high mannose-types in pictilisib-treated A549 cells. This study highlights the use of mass spectrometry to identify a potential glycosylation inhibitor and assessment of its impact on cell surface glycoprotein abundance and protein-protein interaction.

Host Cell Glycocalyx Remodeling Reveals SARS-CoV-2 Spike Protein Glycomic Binding Sites

Glycans on the host cell membrane and viral proteins play critical roles in pathogenesis. Highly glycosylated epithelial cells represent the primary boundary separating embedded host tissues from pathogens within the respiratory and intestinal tracts. SARS-CoV-2, the causative agent for the COVID-19 pandemic, reaches into the respiratory tract. We found purified human milk oligosaccharides (HMOs) inhibited the viral binding on cells. Spike (S) protein receptor binding domain (RBD) binding to host cells were partly blocked by co-incubation with exogenous HMOs, most by 2-6-sialyl-lactose (6′SL), supporting the notion that HMOs can function as decoys in defense against SARS-Cov2. To investigate the effect of host cell glycocalyx on viral adherence, we metabolically modified and confirmed with glycomic methods the cell surface glycome to enrich specific N-glycan types including those containing sialic acids, fucose, mannose, and terminal galactose. Additionally, Immunofluorescence studies demonstrated that the S protein preferentially binds to terminal sialic acids with α-(2,6)-linkages. Furthermore, site-specific glycosylation of S protein RBD and its human receptor ACE2 were characterized using LC-MS/MS. We then performed molecular dynamics calculations on the interaction complex to further explore the interactive complex between ACE2 and the S protein. The results showed that hydrogen bonds mediated the interactions between ACE2 glycans and S protein with desialylated glycans forming significantly fewer hydrogen bonds. These results supported a mechanism where the virus binds initially to glycans on host cells preferring α-(2,6)-sialic acids and finds ACE2 and with the proper orientation infects the cell.

Challenges and limitations in the studies of glycoproteins: A computational chemist's perspective

Experimenters face challenges and limitations while analyzing glycoproteins due to their high flexibility, stereochemistry, anisotropic effects, and hydration phenomena. Computational studies complement experiments and have been used in characterization of the structural properties of glycoproteins. However, recent investigations revealed that computational studies face significant challenges as well. Here, we introduce and discuss some of these challenges and weaknesses in the investigations of glycoproteins. We also present requirements of future developments in computational biochemistry and computational biology areas that could be necessary for providing more accurate structural property analyses of glycoproteins using computational tools. Further theoretical strategies that need to be and can be developed are discussed herein.

DNA Aptamers Block the Receptor Binding Domain at the Spike Protein of SARS-CoV-2

DNA aptamers are versatile molecular species obtained by the folding of short single-stranded nucleotide sequences, with highly specific recognition capabilities against proteins. Here we test the ability of DNA aptamers to interact with the spike (S-)protein of the SARS-CoV-2 viral capsid. The S-protein, a trimer made up of several subdomains, develops the crucial function of recognizing the ACE2 receptors on the surface of human cells, and subsequent fusioning of the virus membrane with the host cell membrane. In order to achieve this, the S1 domain of one protomer switches between a closed conformation, in which the binding site is inaccessible to the cell receptors, and an open conformation, in which ACE2 can bind, thereby initiating the entry process of the viral genetic material in the host cell. Here we show, by means of state-of-the-art molecular simulations, that small DNA aptamers experimentally identified can recognize the S-protein of SARS-CoV-2, and characterize the details of the binding process. We find that their interaction with different subdomains of the S-protein can effectively block, or at least considerably slow down the opening process of the S1 domain, thereby significantly reducing the probability of virus-cell binding. We provide evidence that, as a consequence, binding of the human ACE2 receptor may be crucially affected under such conditions. Given the facility and low cost of fabrication of specific aptamers, the present findings could open the way to both an innovative viral screening technique with sub-nanomolar sensitivity, and to an effective and low impact curative strategy.

Development of Martini 2.2 parameters for N-glycans: a case study of the HIV-1 Env glycoprotein dynamics

N-linked glycans are ubiquitous in nature and play key roles in biology. For example, glycosylation of pathogenic proteins is a common immune evasive mechanism, hampering the development of successful vaccines. Due to their chemical variability and complex dynamics, an accurate molecular understanding of glycans is still limited by the lack of effective resolution of current experimental approaches. Here, we have developed and implemented a reductive model based on the popular Martini 2.2 coarse-grained force field for the computational study of N-glycosylation. We used the HIV-1 Env as a direct applied example of a highly glycosylated protein. Our results indicate that the model not only reproduces many observables in very good agreement with a fully atomistic force field but also can be extended to study large amount of glycosylation variants, a fundamental property that can aid in the development of drugs and vaccines.

Novel Pharmaceutical Strategy for Selective Abrogation of TSP1-Induced Vascular Dysfunction by Decoy Recombinant CD47 Soluble Receptor in Prophylaxis and Treatment Models

Elevated thrombospondin 1 (TSP1) is a prevalent factor, via cognate receptor CD47, in the pathogenesis of cardiovascular conditions, including ischemia-reperfusion injury (IRI) and pulmonary arterial hypertension (PAH). Moreover, TSP1/CD47 interaction has been found to be associated with platelet hyperaggregability and impaired nitric oxide response, exacerbating progression in IRI and PAH. Pathological TSP1 in circulation arises as a target of our novel therapeutic approach. Our “proof-of-concept” pharmacological strategy relies on recombinant human CD47 peptide (rh-CD47p) as a decoy receptor protein (DRP) to specifically bind TSP1 and neutralize TSP1-impaired vasorelaxation, strongly implicated in IRI and PAH. The binding of rh-CD47p and TSP1 was first verified as the primary mechanism via Western blotting and further quantified with modified ELISA, which also revealed a linear molar dose-dependent interaction. Ex vivo, pretreatment protocol with rh-CD47p (rh-CD47p added prior to TSP1 incubation) demonstrated a prophylactic effect against TSP1-impairment of endothelium-dependent vasodilation. Post-treatment set-up (TSP1 incubation prior to rh-CD47p addition), mimicking pre-existing excessive TSP1 in PAH, reversed TSP1-inhibited vasodilation back to control level. Dose titration identified an effective molar dose range (approx. ≥1:3 of tTSP1:rh-CD47p) for prevention of/recovery from TSP1-induced vascular dysfunction. Our results indicate the great potential for proposed novel decoy rh-CD47p-therapy to abrogate TSP1-associated cardiovascular complications, such as PAH.

Exploring the Role of Glycans in the Interaction of SARS-CoV-2 RBD and Human Receptor ACE2

COVID-19 is a highly infectious respiratory disease caused by the novel coronavirus SARS-CoV-2. It has become a global pandemic and its frequent mutations may pose new challenges for vaccine design. During viral infection, the Spike RBD of SARS-CoV-2 binds the human host cell receptor ACE2, enabling the virus to enter the host cell. Both the Spike and ACE2 are densely glycosylated, and it is unclear how distinctive glycan types may modulate the interaction of RBD and ACE2. Detailed understanding of these determinants is key for the development of novel therapeutic strategies. To this end, we perform extensive all-atom simulations of the (i) RBD-ACE2 complex without glycans, (ii) RBD-ACE2 with oligomannose MAN9 glycans in ACE2, and (iii) RBD-ACE2 with complex FA2 glycans in ACE2. These simulations identify the key residues at the RBD-ACE2 interface that form contacts with higher probabilities, thus providing a quantitative evaluation that complements recent structural studies. Notably, we find that this RBD-ACE2 contact signature is not altered by the presence of different glycoforms, suggesting that RBD-ACE2 interaction is robust. Applying our simulated results, we illustrate how the recently prevalent N501Y mutation may alter specific interactions with host ACE2 that facilitate the virus-host binding. Furthermore, our simulations reveal how the glycan on Asn90 of ACE2 can play a distinct role in the binding and unbinding of RBD. Finally, an energetics analysis shows that MAN9 glycans on ACE2 decrease RBD-ACE2 affinity, while FA2 glycans lead to enhanced binding of the complex. Together, our results provide a more comprehensive picture of the detailed interplay between virus and human receptor, which is much needed for the discovery of effective treatments that aim at modulating the physical-chemical properties of this virus.

NMR of glycoproteins: profiling, structure, conformation and interactions

In glycoproteins, carbohydrates are responsible for the selective interaction and tight regulation of cellular processes, constituting the main information transducer interface in protein-glycoprotein interactions. Increasing experimental and computational evidence suggest that such interactions often induce allosteric changes in the host protein, underlining the importance of studying intact glycoproteins. Technical issues have precluded such studies for years but, nowadays, a promising era is emerging where NMR spectroscopy, among other techniques, allows the characterization of the composition, structure and segmental dynamics of glycoproteins. In this review, we discuss such advances and highlight some selected examples. This novel technology unravels multiple new functional mechanisms, subtly hidden within the sugar code.

Identification of Novel Scaffolds with Dual Role as Antiepileptic and Anti-Breast Cancer

Aromatase inhibitors with an $\mathrm{IC}_{50}$ IC 50 value ranging from 1.4 to 49.7 µM are known to act as antiepileptic drugs besides being potential breast cancer inhibitors. The aim of the present study is to identify novel antiepileptic aromatase inhibitors with higher activity exploiting the ligand-based pharmacophore approach utilizing the experimentally known inhibitors. The resultant Hypo1 consists of four features and was further validated by using three different strategies. Hypo1 was allowed to screen different databases to identify lead molecules and were further subjected to Lipinski's Rule of Five and ADMET to establish their drug-like properties. Consequently, the obtained 68-screened molecules were subjected to molecular docking by GOLD v5.2.2. Furthermore, the compounds with the highest dock scores were assessed for molecular interactions. Later, the MD simulation was applied to evaluate the protein backbone stabilities and binding energies adapting GROMACS v5.0.6 and MM/PBSA which was followed by the density functional theory (DFT), to analyze their orbital energies, and further the energy gap between them. Eventually, the number of Hit molecules was culled to three projecting Hit1, Hit2, and Hit3 as the potential lead compounds based on their highest dock scores, hydrogen bond interaction, lowest energy gap, and the least binding energies and stable MD results.

Characterization of ELISA Antibody-Antigen Interaction using Footprinting-Mass Spectrometry and Negative Staining Transmission Electron Microscopy

We describe epitope mapping data using multiple covalent labeling footprinting-mass spectrometry (MS) techniques coupled with negative stain transmission electron microscopy (TEM) data to analyze the antibody-antigen interactions in a sandwich enzyme-linked immunosorbant assay (ELISA). Our hydroxyl radical footprinting-MS data using fast photochemical oxidation of proteins (FPOP) indicates suppression of labeling across the antigen upon binding either of the monoclonal antibodies (mAbs) utilized in the ELISA. Combining these data with Western blot analysis enabled the identification of the putative epitopes that appeared to span regions containing N-linked glycans. An additional structural mapping technique, carboxyl group footprinting-mass spectrometry using glycine ethyl ester (GEE) labeling, was used to confirm the epitopes. Deglycosylation of the antigen resulted in loss of potency in the ELISA, supporting the FPOP and GEE labeling data by indicating N-linked glycans are necessary for antigen binding. Finally, mapping of the epitopes onto the antigen crystal structure revealed an approximate 90° relative spatial orientation, optimal for a noncompetitive binding ELISA. TEM data shows both linear and diamond antibody-antigen complexes with a similar binding orientation as predicted from the two footprinting-MS techniques. This study is the first of its kind to utilize multiple bottom-up footprinting-MS techniques and TEM visualization to characterize the monoclonal antibody-antigen binding interactions of critical reagents used in a quality control (QC) lot-release ELISA. Graphical Abstract ᅟ.

Ginger (Zingiber officinale) phytochemicals-gingerenone-A and shogaol inhibit SaHPPK: molecular docking, molecular dynamics simulations and in vitro approaches

BACKGROUND

Antibiotic resistance is a defense mechanism, harbored by pathogens to survive under unfavorable conditions. Among several antibiotic resistant microbial consortium, Staphylococcus aureus is one of the most havoc microorganisms. Staphylococcus aureus encodes a unique enzyme 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (SaHPPK), against which, none of existing antibiotics have been reported.

METHODS

Computational approaches have been instrumental in designing and discovering new drugs for several diseases. The present study highlights the impact of ginger phytochemicals on Staphylococcus aureus SaHPPK. Herein, we have retrieved eight ginger phytochemicals from published literature and investigated their inhibitory interactions with SaHPPK. To authenticate our work, the investigation proceeds considering the known antibiotics alongside the phytochemicals. Molecular docking was performed employing GOLD and CDOCKER. The compounds with the highest dock score from both the docking programmes were tested for their inhibitory capability in vitro. The binding conformations that were seated within the binding pocket showing strong interactions with the active sites residues rendered by highest dock score were forwarded towards the molecular dynamic (MD) simulation analysis.

RESULTS

Based on molecular dock scores, molecular interaction with catalytic active residues and MD simulations studies, two ginger phytochemicals, gingerenone-A and shogaol have been proposed as candidate inhibitors against Staphylococcus aureus. They have demonstrated higher dock scores than the known antibiotics and have represented interactions with the key residues within the active site. Furthermore, these compounds have rendered considerable inhibitory activity when tested in vitro. Additionally, their superiority was corroborated by stable MD results conducted for 100 ns employing GROMACS package.

CONCLUSIONS

Finally, we suggest that gingerenone-A and shogaol may either be potential SaHPPK inhibitors or can be used as fundamental platforms for novel SaHPPK inhibitor development.

Exploration for novel inhibitors showing back-to-front approach against VEGFR-2 kinase domain (4AG8) employing molecular docking mechanism and molecular dynamics simulations

Background

Angiogenesis is a process of formation of new blood vessels and is an important criteria demonstrated by cancer cells. Over a period of time, these cancer cells infect the other parts of the healthy body by a process called progression. The objective of the present article is to identify a drug molecule that inhibits angiogenesis and progression.

Methods

In this pursuit, ligand based pharmacophore virtual screening was employed, generating a pharmacophore model, Hypo1 consisting of four features. Furthermore, this Hypo1 was validated recruiting, Fischer’s randomization, test set method and decoy set method. Later, Hypo1 was allowed to screen databases such as Maybridge, Chembridge, Asinex and NCI and were further filtered by ADMET filters and Lipinski’s Rule of Five. A total of 699 molecules that passed the above criteria, were challenged against 4AG8, an angiogenic drug target employing GOLD v5.2.2.

Results

The results rendered by molecular docking, DFT and the MD simulations showed only one molecule (Hit) obeyed the back-to-front approach. This molecule displayed a dock score of 89.77, involving the amino acids, Glu885 and Cys919, Asp1046, respectively and additionally formed several important hydrophobic interactions. Furthermore, the identified lead molecule showed interactions with key residues when challenged with CDK2 protein, 1URW.

Conclusion

The lead candidate showed several interactions with the crucial residues of both the targets. Furthermore, we speculate that the residues Cys919 and Leu83 are important in the development of dual inhibitor. Therefore, the identified lead molecule can act as a potential inhibitor for angiogenesis and progression.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4050-1) contains supplementary material, which is available to authorized users.

Antifreeze Glycoproteins Bind Reversibly to Ice via Hydrophobic Groups

Antifreeze molecules allow organisms to survive in subzero environments. Antifreeze glycoproteins (AFGPs), produced by polar fish, are the most potent inhibitors of ice recrystallization. To date, the molecular mechanism by which AFGPs bind to ice has not yet been elucidated. Mutation experiments cannot resolve whether the binding occurs through the peptide, the saccharides, or both. Here, we use molecular simulations to determine the mechanism and driving forces for binding of AFGP8 to ice, its selectivity for the primary prismatic plane, and the molecular origin of its exceptional ice recrystallization activity. Consistent with experiments, AFGP8 in simulations preferentially adopts the PPII helix secondary structure in solution. We show that the segregation of hydrophilic and hydrophobic groups in the PPII helix is vital for ice binding. Binding occurs through adsorption of methyl groups of the peptide and disaccharides to ice, driven by the entropy of dehydration of the hydrophobic groups as they nest in the cavities at the ice surface. The selectivity to the primary prismatic plane originates in the deeper cavities it has compared to the basal plane. We estimate the free energy of binding of AFGP8 and the longer AFGPs4-6, and find them to be consistent with the reversible binding demonstrated in experiments. The simulations reveal that AFGP8 binds to ice through a myriad of conformations that it uses to diffuse through the ice surface and find ice steps, to which it strongly adsorbs. We interpret that the existence of multiple, weak binding sites is the key for the exceptional ice recrystallization inhibition activity of AFGPs.

Structural database resources for biological macromolecules

This Briefing reviews the widely used, currently active, up-to-date databases derived from the worldwide Protein Data Bank (PDB) to facilitate browsing, finding and exploring its entries. These databases contain visualization and analysis tools tailored to specific kinds of molecules and interactions, often including also complex metrics precomputed by experts or external programs, and connections to sequence and functional annotation databases. Importantly, updates of most of these databases involves steps of curation and error checks based on specific expertise about the subject molecules or interactions, and removal of sequence redundancy, both leading to better data sets for mining studies compared with the full list of raw PDB entries. The article presents the databases in groups such as those aimed to facilitate browsing through PDB entries, their molecules and their general information, those built to link protein structure with sequence and dynamics, those specific for transmembrane proteins, nucleic acids, interactions of biomacromolecules with each other and with small molecules or metal ions, and those concerning specific structural features or specific protein families. A few webservers directly connected to active databases, and a few databases that have been discontinued but would be important to have back, are also briefly commented on. Along the Briefing, sample cases where these databases have been used to aid structural studies or advance our knowledge about biological macromolecules are referenced. A few specific examples are also given where using these databases is easier and more informative than using raw PDB data.

GENESIS 1.1: A hybrid-parallel molecular dynamics simulator with enhanced sampling algorithms on multiple computational platforms

GENeralized-Ensemble SImulation System (GENESIS) is a software package for molecular dynamics (MD) simulation of biological systems. It is designed to extend limitations in system size and accessible time scale by adopting highly parallelized schemes and enhanced conformational sampling algorithms. In this new version, GENESIS 1.1, new functions and advanced algorithms have been added. The all-atom and coarse-grained potential energy functions used in AMBER and GROMACS packages now become available in addition to CHARMM energy functions. The performance of MD simulations has been greatly improved by further optimization, multiple time-step integration, and hybrid (CPU + GPU) computing. The string method and replica-exchange umbrella sampling with flexible collective variable choice are used for finding the minimum free-energy pathway and obtaining free-energy profiles for conformational changes of a macromolecule. These new features increase the usefulness and power of GENESIS for modeling and simulation in biological research. © 2017 Wiley Periodicals, Inc.

Reparameterization of Solute-Solute Interactions for Amino Acid-Sugar Systems Using Isopiestic Osmotic Pressure Molecular Dynamics Simulations

AMBER/GLYCAM and CHARMM are popular force fields for simulations of amino acids and sugars. Here we report excessively attractive amino acid-sugar interactions in both force fields, and corrections to nonbonded interactions that match experimental osmotic pressures of mixed aqueous solutions of diglycine and sucrose. The modified parameters also improve the ΔGtrans of diglycine from water to aqueous sucrose and, with AMBERff99SB/GLYCAM06, eliminate a caging effect seen in previous simulations of the protein ubiquitin with glucose.