Identification of possible binding modes of SARS-CoV-2 spike N-terminal domain for ganglioside GM1.

Coarse-grained molecular dynamics simulations of the lipid bilayer mixture of POPC and cholesterol were carried out in the presence and absence of ganglioside monosialo 1 (GM1) with N - terminal domain (NTD) of SARS-CoV-2 spike glycoprotein. The interactions of GM1 with two different NTD orientations were compared. NTD orientation I compactly bind GM1 predominantly through the sialic acid and the external galactose moieties providing more restriction to GM1 mobility whereas orientation II is more distributed on the lipid surface and due to the relaxed mobility of GM1 there, presumably, the NTD receptor penetrates more through the membrane.

Effect of Transmembrane Electric Field on GM1 Containing DMPC-Cholesterol Monolayer: A Computational Study.

Transmembrane electric potentials and membrane curvature have always provided pathways to mediate different cellular processes. We present results of molecular dynamics (MD) simulations of lipid monolayer composed of 1, 2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and cholesterol (CHOL) under a transverse electric field to monitor the effect of electric field on membrane containing ganglioside monosialo 1 (GM1). Four systems were studied with membrane monolayer in the presence and absence of GM1 with and without applying electric field along the normal of the monolayer. The applied transmembrane electric field was 0.4 mV/Å which corresponds to the action potential of animal cell. Our results indicate that the electric field induces a considerable lateral stress on the monolayer in the presence of GM1, which is evident from the lateral pressure profiles. It was found that due to the application of electric field major perturbation was caused to the system containing GM1, manifested by the bending of the monolayer. We believe this study provides correlation between electric field and spontaneous membrane bending, specially based on the membrane composition. The consequences of these MD simulations provide considerable insights to different biological phenomenon and lipid membrane models.

Phase Behavior of GM1-Containing DMPC-Cholesterol Monolayer: Experimental and Theoretical Study.

Organization and distribution of lipids in cellular membranes play an important role in a diverse range of biological processes, such as membrane trafficking and signaling. Here, we present the combined experimental and simulated results to elucidate the phase behavioral features of ganglioside monosialo 1 (GM1)-containing mixed monolayer of the lipids 1,2-dimyristoyl- sn-glycero-3-phosphocholine (DMPC) and cholesterol (CHOL). Two monolayers having compositions DMPC-CHOL and GM1-DMPC-CHOL are investigated at air-water and air-solid interfaces using Langmuir-Blodgett experiments and scanning electron microscopy (SEM), respectively, to ascertain the phase behavior change of the monolayers. Surface pressure isotherms and SEM imaging of domain formation indicate that addition of GM1 to the monolayer at low surface pressure causes a fluidization of the system but once the system attains the surface pressure corresponding to its liquid-condensed phase, the monolayer becomes more ordered than the system devoid of GM1 and interacts among each other more cooperatively. Besides, the condensing effect of cholesterol on the DMPC monolayer was also verified by our experiments. Apart from these, the effects induced by GM1 on the phase behavior of the binary mixture of DMPC-CHOL were studied with and without applying liquid-expanded (LE)-liquid-condensed (LC) equilibrium surface pressure using molecular dynamics (MD) simulation. Our molecular dynamics (MD) simulation results give an atomistic-level explanation of our experimental findings and furnish a similar conclusion.

Interaction Between Luteinizing Hormone-Releasing Hormone and GM1-Doped Cholesterol/Sphingomyelin Vesicles: A Spectroscopic Study.

Understanding the role of neural membrane in translocation and action of neurohormone is of great importance. Luteinizing hormone-releasing hormone (LHRH) is a neuropeptide hormone and it acts as a final signaling molecule by stimulating the synthesis of LH and FSH to maintain reproduction in all vertebrates. The receptors of LHRH are found in breast tumors and pituitary gland in the brain. Moreover, neural plasma membrane is also found to contain specific binding site for LHRH. The mechanism by which LHRH binds to membrane before it binds to the receptors is a very critical step and can have a profound impact upon the translation of peptide across the membrane. A complex form of glycosphingolipids known as Ganglioside is an important component of plasma membrane of nerve cells and breast tumor tissues. They play an important role in various physiological membrane processes. Therefore, the interaction of ganglioside-containing membrane with LHRH might be crucial in aiding the LHRH to translate through the neural membrane and reach its receptor for binding and activation. Using CD, UV-Absorbance, and fluorescence spectroscopy, the effect of Ganglioside Monosialo 1(GM1)-induced conformational changes of LHRH in the presence of Cholesterol (CHOL)/Sphingomyelin (SM) and GM1/CHOL/SM vesicles was studied. The aforesaid spectroscopic studies show that LHRH is able to bind with both the vesicles, but GM1-containing vesicles interact more effectively than vesicles without GM1. CHOL/SM vesicles partially disturb the conformation of the peptide. Moreover, binding of LHRH to GM1/CHOL/SM vesicles induces loss of conformational rigidity and attainment of a random coil.

Localization and dynamics of the anticarcinogenic curcumin with GM1 and other miceller assemblies.

Structural transitions involving shape changes play an important role in cellular physiology and enhance the bioavailability of the natural food like curcumin in surfactant aggregates. In this work, we have studied the localization, dynamics and stability of curcumin in various miceller assemblies using a combination of absorbance and fluorescence spectroscopic approaches. The measurements of absorption and fluorescence spectra of curcumin revealed that the nature of interactions of ionic and nonionic surfactants and the glycosphingolipid, GM1 with curcumin is significantly different with surfactant concentrations. At low concentrations of SDS and the GM1 the head group of SDS and GM1 binds to the central ß-diketone group of curcumin to form SDS-curcumin or GM1-curcumin complexes. At high concentrations, both formed micelles with curcumin completely solubilized inside. Cucurmin is solubilized in the stern layer of SDS micelles. Compared to spherical micelles, rod shaped micelles allow more curcumin to bind through hydrophobic interactions indicated by higher absorption and fluorescence, enhanced partition coefficient and stability. Whereas curcumin associates with GM1 micelles with lower partition coefficient, solubility and remain closer to aqueous phase decreasing its bioavailability and stability. While cucurmin is solubilized in the palisade layer of deoxycholate and octyl glucopyranoside micelles through the alkyl chains providing more hydrophobic microenvironment to curcumin with enhanced stability and bioavailability. Graphical abstract Schematic diagram of the two different types of detergent micelles and larger GM1 micelles.

Ion channel stability of Gramicidin A in lipid bilayers: effect of hydrophobic mismatch.

Hydrophobic mismatch which is defined as the difference between the lipid hydrophobic thickness and the peptide hydrophobic length is known to be responsible in altering the lipid/protein dynamics. Gramicidin A (gA), a 15 residue ß helical peptide which is well recognized to form ion conducting channels in lipid bilayer, may change its structure and function in a hydrophobic mismatched condition. We have performed molecular dynamics simulations of gA dimer in phospholipid bilayers to investigate whether or not the conversion from channel to non-channel form of gA dimer would occur under extreme negative hydrophobic mismatch. By varying the length of lipid bilayers from DLPC (1, 2-Dilauroyl-sn-glycero-3-phosphocholine) to DAPC (1, 2-Diarachidoyl-sn-glycero-3-phosphocholine), a broad range of mismatch was considered from nearly matching to extremely negative. Our simulations revealed that though the ion-channel conformation is retained by gA under a lesser mismatched situation, in extremely negative mismatched situation, in addition to bilayer thinning, the conformation of gA is changed and converted to a non-channel one. Our results demonstrate that although the channel conformation of Gramicidin A is the most stable structure, it is possible for gA to change its conformation from channel to non-channel depending upon the local environment of host bilayers.

Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations.

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man9GlcNAc2). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct ß sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

Organization and dynamics of tryptophan residues in brain spectrin: novel insight into conformational flexibility.

Brain spectrin enjoys overall structural and sequence similarity with erythroid spectrin, but less is known about its function. We utilized the fluorescence properties of tryptophan residues to monitor their organization and dynamics in brain spectrin. Keeping in mind the functional relevance of hydrophobic binding sites in brain spectrin, we monitored the organization and dynamics of brain spectrin bound to PRODAN. Results from red edge excitation shift (REES) indicate that the organization of tryptophans in brain spectrin is maintained to a considerable extent even after denaturation. These results are supported by acrylamide quenching experiments. To the best of our knowledge, these results constitute the first report of the presence of residual structure in urea-denatured brain spectrin. We further show from REES and time-resolved emission spectra that PRODAN bound to brain spectrin is characterized by motional restriction. These results provide useful information on the differences between erythroid spectrin and brain spectrin.

In silico phase separation in the presence of GM1 in ternary and quaternary lipid bilayers.

Cell membranes are multi-component mixtures with structural and compositional heterogeneity exhibiting a complex phase behavior. Domains formed in cell membranes often known as "Rafts" are of immense importance. Using coarse grained molecular dynamics simulations, we have studied the spontaneous phase separation of the ternary (POPC [1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine]/cholesterol/GM1) and quaternary (POPC/PSM[palmitoyl sphingomyelin]/cholesterol/GM1) lipid bilayers into liquid ordered (Lo) and liquid disordered (Ld) domains due to self-aggregation of GM1 molecules and co-localization of cholesterol with GM1 in accordance with experiments. It is found that GM1 molecules have the ability to associate strongly with each other which leads to the formation of ordered domains in the lipid mixture and the interactions are through the head group and unsaturated tails present in GM1. Preference of cholesterol for association with GM1 over PSM is observed, the domains consisting of GM1 and cholesterol are formed even in the presence of PSM. PSM also forms small domains with cholesterol that are randomly distributed in the Ld phase. Estimation of dynamic quantities like diffusion coefficient also shows that cholesterol has the highest diffusion rate in the Ld phase which is further attributed to its flip flop ability. It is found that in the presence of PSM, cholesterol can undergo flip flop even in the Lo phase. This is accredited to the interaction of cholesterol with PSM from which it can be concluded that in the presence of PSM, the domains formed by GM1 are less tightly packed and less stable than that in the ternary mixture.

Capability of ganglioside GM1 in modulating interactions, structure, location and dynamics of peptides/proteins: biophysical approaches: interaction of ganglioside GM1 with peptides/proteins.

Gangliosides, are glycosphingolipids, present in all vertebrate plasma membranes with particular abundance in nerve cell membrane. Gangliosides can act as portals for antimicrobial peptides, hormones, viruses, lectins, toxins and pathogens. They are strategically positioned on the outer membrane and hence can participate in a large number of recognition processes. Their abundance in nerve cell membrane makes them "likely" receptor candidates for neuropeptides. In this review we outline our work in the area of GM1-peptide/protein interaction. We have explored the effect of GM1 containing micelles/bicelles on structures of peptides, proteins as well as on denatured proteins. It has been observed that the peptides that are disordered or having random coil structure in aqueous solution, attained an ordered three-dimensional structure when interact with GM1. It is also observed that denatured proteins undergo refolding in presence of ganglioside. Peptides/proteins show stronger interaction with membrane lipid bilayer in presence of ganglioside than that without ganglioside. This review mainly focuses on capability of ganglioside GM1 in modulating interaction, structural, location and dynamics of peptides/proteins using a number of biophysical techniques-solution NMR, DOSY, CD, fluorescence etc.

Insights into binding of cholera toxin to GM1 containing membrane.

Interactions of cholera toxin (CT) with membrane are associated with the massive secretory diarrhea seen in Asiatic cholera. Ganglioside GM1 has been shown to be responsible for the binding of the B subunit of cholera toxin (CT-B), which then helps CT to pass through the membrane, but the exact mechanism remains to be explored. In this work, we have carried out atomistic scale molecular dynamics simulation to investigate the structural changes of CT upon membrane binding and alteration in membrane structure and dynamics. Starting from the initial structure where the five units of B subunit bind with five GM1, only three of five units remain bound and the whole CT is tilted such that the three binding units are deeper in the membrane. The lipids that are in contact with those units of the CT-B behave differently from the rest of the lipids. Altogether, our results demonstrate the atomistic interaction of CT with GM1 containing lipid membrane and provide a probable mechanism of the early stage alteration of lipid structure and dynamics, which can make a passage for penetration of CT on membrane surface.

Microsecond molecular dynamics simulation of guanidinium chloride induced unfolding of ubiquitin.

An all atom molecular dynamics simulation was used to explore the atomic detail mechanism of guanidinium induced unfolding of the protein ubiquitin. Ubiquitin unfolds through pre-unfolded (intermediate) states, i.e. guanidinium induced unfolding of ubiquitin appears to be a multi-step process, and loss of hydrophobic contacts of C-terminal residues is crucial for ubiquitin unfolding. Free-energy landscapes show that barrier separation between folded and unfolded basins is ∼5.0 kcal mol(-1), and both the basins are of comparable energy. It was observed that guanidinium ions interact directly with ubiquitin. Favorable electrostatic interaction is the main driving force for such accumulation of guanidinium ions near protein, but van der Waals energy also contributes. RDF plots show that accumulation of guanidinium ions near specific residues is the main cause for destabilization of intra-residue interactions crucial to maintain the three-dimensional fold of the protein. One salt-bridge interaction between Lys11 and Glu34 appears to be important to maintain the crystal structure of ubiquitin and this salt-bridge can map the unfolding process of ubiquitin.

Computational studies suggest compounds restoring function of p53 cancer mutants can bind SARS-CoV-2 spike protein.

It is reasonable to think that cancer patients undergoing chemotherapy or immunotherapy may have a more aggressive course if they are positive for the novel coronavirus disease. Their compulsive condition requires investigation into effective drugs. We applied computational techniques to a series of compounds known for restoring the function of p53 cancer mutant p53R175H and p53G245S. Two potent inhibitors, 1-(3-chlorophenyl)-3-(1, 3 -thiazol-2-yl) urea (CTU, PubChem NSC321792) with the highest binding affinity -6.92 kcal/mol followed by a thiosemicarbazone compound N'-(1-(Pyridin-2-yl)ethylidene) azetidine - 1 -carbothiohydrazide (NPC, PubChem NSC319726) with -6.75 kcal/mol were subjected to Molecular Dynamics simulation with receptor binding domain (RBD) and compared with control ligand dexamethasone. In particular, CTU adheres to pocket 1 with an average free energy of binding -21.65 ± 2.89 kcal/mol at the RBD - angiotensin-converting enzyme 2 binding region with the highest frequency of amino acid residues after reaching a local equilibrium in 100 ns MD simulation trajectory. A significant enthalpy contribution from the independent simulations unfolds the possibility of dual binding sites for NPC as shifted pocket 1 (-15.59 ± 5.98 kcal/mol) and pocket 2 (-18.90 ± 5.02 kcal/mol). The obtained results for these two compounds are in good agreement with dexamethasone (-18.45 ± 2.42 kcal/mol). Taken together our findings could facilitate the discovery of small molecules that restore the function of p53 cancer mutants newly against COVID-19 in cancer patients.Communicated by Ramaswamy H. Sarma.

Molecular dynamics simulations suggest Thiosemicarbazones can bind p53 cancer mutant R175H.

Cancer pathologies are associated with the unfolding and aggregation of most recurring mutations in the DNA Binding Domain (DBD) of p53 that coordinate the destabilization of protein. Substitution at the 175th codon with arginine to histidine (R175H, a mutation of large to small side-chain amino acid) destabilizes the DBD by 3 kcal/mol and triggers breasts, lung cancer, etc. Stabilizing the p53 mutant by small molecules offers an attractive drug-targeted anti-cancer therapy. The thiosemicarbazone (TSC) molecules NPC and DPT are known to act as zinc-metallochaperones to reactivate p53R175H. Here, a combination of LESMD simulations for 10 TSC conformations with a p53R175H receptor, single ligand-protein conformation MD, and ensemble docking with multiple p53R175H conformations observed during simulations is suggested to identify the potential binding site of the target protein in light of their importance for the direct TSC - p53R175H binding. NPC binds mutant R175H in the loop region L2-L3, forming pivotal hydrogen bonds with HIS175, pi­sulfur bonds with TYR163, and pi-alkyl linkages with ARG174 and PRO190, all of which are contiguous to the zinc-binding native site on p53DBD. DPT, on the other hand, was primarily targeting alternative binding sites such as the loop-helix L1/H2 region and the S8 strand. The similar structural characteristics of TSC-bound p53R175H complexes with wild-type p53DBD are thought to be attributable to involved interactions that favour binding free energy contributions of TSC ligands. Our findings may be useful in the identification of novel pockets with druggable properties.

Comparison and Possible Binding Orientations of SARS-CoV-2 Spike N-Terminal Domain for Gangliosides GM3 and GM1.

SARS-CoV-2 spike glycoprotein is anchored by gangliosides. The sialic acid in the ganglioside headgroup is responsible for virus attachment and entry into host cells. We used coarse-grained (CG) molecular dynamics simulations to expand on our previous study of GM1 interaction with two different orientations of the SARS-CoV-2 S1 subunit N-terminal domain (NTD) and to confirm the role of sialic acid receptors in driving the viral receptor; GM3 was used as another ganglioside on the membrane. Because of the smaller headgroup, sialic acid is crucial in GM3 interactions, whereas GM1 interacts with NTD via both the sialic acid and external galactose. In line with our previous findings for NTD orientations in GM1 binding, we identified two orientations, "compact" and "distributed", comprising sugar receptor-interacting residues in GM3-embedded lipid bilayers. Gangliosides in closer proximity to the compact NTD orientation might cause relatively greater restrictions to penetrate the bilayer. However, the attachment of a distributed NTD orientation with more negative interaction energies appears to facilitate GM1/GM3 to move quickly across the membrane. Our findings likely shed some light on the orientations that the NTD receptor acquires during the early phases of interaction with GM1 and GM3 in a membrane environment.

NMR evidence of GM1-induced conformational change of Substance P using isotropic bicelles.

Substance P (SP) is one of the target neurotransmitters associated with diseases related to chronic inflammation, pain and depression. The selective receptor for SP, NK(1)R is located in the heterogeneous microdomains or caveolae in membrane. Gangliosides, specifically GM1, are markers of these heterogeneous sites. Also, gangliosides are considered as important regulatory elements in cell-cell recognition and cell signaling. In the present work, we describe the conformations of Substance P in the presence of ternary membrane systems containing GM1 at the physiological concentration. SP is mostly unstructured in water, but appears as extended 3(10) helical or turn III in isotropic bicelles, more pronounced in the presence of GM1. NMR results suggest that, in the GM1 containing bicelles, the peptide is more inserted into the membrane with its C-terminus, while N-terminus lies close to the membrane-water interface. The NMR-derived conformation of SP in GM1 bicelles is docked on homology modeled NK(1)R and resulting interactions satisfy reported mutagenesis, fluorescence, photo-affinity labeling and modeling data. The results highlight efficacy of GM1 in membrane in providing structure in an otherwise flexible neurotransmitter Substance P; thus providing indication that it may be useful also for other neurotransmitter peptides/proteins associated with membrane.

Design and self-assembly of a leucine-enkephalin analogue in different nanostructures: application of nanovesicles.

An opioid (leucine-enkephalin) conformational analogue forms diverse nanostructures such as vesicles, tubes, and organogels through self-assembly. The nanovesicles encapsulate the natural hydrophobic drug curcumin and allow the controlled release through cation-generated porogens in membrane mimetic solvent.

Correction to "Aggregation of Lysozyme in the Presence of a Mixed Bilayer of POPC and POPG".

[This corrects the article DOI: 10.1021/acsomega.1c01145.].

Molecular dynamics simulations of the interactions of kinin peptides with an anionic POPG bilayer.

We have performed molecular dynamics simulations of peptide hormone bradykinin (BK) and its fragment des-Arg9-BK in the presence of an anionic lipid bilayer, with an aim toward delineating the mechanism of action related to their bioactivity. Starting from the initial aqueous environment, both of the peptides are quickly adsorbed and stabilized on the cell surface. Whereas BK exhibits a stronger interaction with the membrane and prefers to stay on the interface, des-Arg9-BK, with the loss of C-terminal Arg, penetrates further. The heterogeneous lipid-water interface induces ß-turn-like structure in the otherwise inherently flexible peptides. In the membrane-bound state, we observed C-terminal ß-turn formation in BK, whereas for des-Arg9-BK, with the deletion of Arg9, turn formation occurred in the middle of the peptide. The basic Arg residues anchor the peptide to the bilayer by strong electrostatic interactions with charged lipid headgroups. Simulations with different starting orientations of the peptides with respect to the bilayer surface lead to the same observations, namely, the relative positioning of the peptides on the membrane surface, deeper penetration of the des-Arg9-BK, and the formation of turn structures. The lipid headgroups adjacent to the bound peptides become substantially tilted, causing bilayer thinning near the peptide contact region and increase the degree of disorder in nearby lipids. Again, because of hydrogen bonding with the peptide, the neighboring lipid's polar heads exhibit considerably reduced flexibility. Corroborating findings from earlier experiments, our results provide important information about how the lipid environment promotes peptide orientation/conformation and how the peptide adapts to the environment.

Cholesterol driven alteration of the conformation and dynamics of phospholamban in model membranes.

The effects of cholesterol on various membrane proteins are of long-standing interest in membrane biophysics. Here we present systematic molecular dynamics simulations (totaling 1.4 µs) of integral protein phospholamban incorporated in POPC/cholesterol bilayers (containing 0, 11.11, 22.03, 33.33, and 50 mol% of cholesterol). Phospholamban is a key regulator of cardiac contractility and has recently emerged as a potential drug target. In agreement with experiments, our results show that in a cholesterol-free pure POPC bilayer, phospholamban exhibits broad conformational distribution, ranging from the closed T-state to the extended R-state, crucial for its functionality. Increasing cholesterol concentration progressively stabilizes the bent conformers of phospholamban over open structures, and favors extensive interactions of its amphipathic N-terminal helix with the bilayer surface. The interaction energies between the N-terminal helix of PLB and different POPC/cholesterol bilayers quantitatively confirm its stronger interaction with a higher cholesterol-containing membrane. Simulation with 50 mol% of cholesterol further supports the above conclusions, where phospholamban undergoes rapid conformational transition from extended to closed form, which remains stable for the rest of the simulation time and exhibits the strongest interaction with the membrane. Cholesterol participates in hydrogen-bonding and π-stacking interactions with polar and/or aromatic residues and favors membrane association of phospholamban. We observed cholesterol-enrichment in the neighborhood of phospholamban. Moreover, as a modulator of membrane biophysical properties, cholesterol modifies the hydrophobic matching and trans-membrane tilting of phospholamban and also hinders its 2D-lateral mobility. Altogether, our results highlight atomistic details of protein-lipid interplay and provide new insights into the possible effects of cholesterol on conformational dynamics of phospholamban in membrane bilayers.