Structure and intermolecular interactions in spheroidal high-density lipoprotein subpopulations

•

High-density lipoprotein subpopulations have unique surface profiles and dynamics.

•

Relative hydrophobic surface area decreases with increasing lipoprotein size.

•

Core lipid exposure at the lipoprotein surface decreases with increasing size.

•

Cholesterol molecules localise near apolipoprotein A-I central helices.

•

Lipid and protein interactions stabilise multifoil models of apolipoprotein A-I.

Human serum high-density lipoproteins (HDLs) are a population of small, dense protein-lipid aggregates that are crucial for intravascular lipid trafficking and are protective against cardiovascular disease. The spheroidal HDL subfraction can be separated by size and density into five major subpopulations with distinct molecular compositions and unique biological functionalities: HDL_3c, HDL_3b, HDL_3a, HDL_2a and HDL_2b. Representative molecular models of these five subpopulations were developed and characterised for the first time in the presence of multiple copies of its primary protein component apolipoprotein A-I (apoA-I) using coarse-grained molecular dynamics simulations. Each HDL model exhibited size, morphological and compositional profiles consistent with experimental observables. With increasing particle size the separation of core and surface molecules became progressively more defined, resulting in enhanced core lipid mixing, reduced core lipid exposure at the surface, and the formation of an interstitial region between core and surface molecules in HDL_2b. Cholesterol molecules tended to localise around the central helix-5 of apoA-I, whilst triglyceride molecules predominantly interacted with aromatic, hydrophobic residues located within the terminal helix-10 across all subpopulation models. The three intermediate HDL models exhibited similar surface profiles despite having distinct molecular compositions. ApoA-I in trefoil, quatrefoil and pentafoil arrangements across the surface of HDL particles exhibited significant warping and twisting, but largely retained intermolecular contacts between adjacent apoA-I chains. Representative HDL subpopulations differed in particle size, morphology, intermolecular interaction profiles and lipid and protein dynamics. These findings reveal how different HDL subpopulations might exhibit distinct functional associations depending on particle size, form and composition.

Heparin Inhibits Cellular Invasion by SARS-CoV-2: Structural Dependence of the Interaction of the Spike S1 Receptor-Binding Domain with Heparin

The dependence of development and homeostasis in animals on the interaction of hundreds of extracellular regulatory proteins with the peri- and extracellular glycosaminoglycan heparan sulfate (HS) is exploited by many microbial pathogens as a means of adherence and invasion. Heparin, a widely used anticoagulant drug, is structurally similar to HS and is a common experimental proxy. Exogenous heparin prevents infection by a range of viruses, including S-associated coronavirus isolate HSR1. Here, we show that heparin inhibits severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) invasion of Vero cells by up to 80% at doses achievable through prophylaxis and, particularly relevant, within the range deliverable by nebulisation. Surface plasmon resonance and circular dichroism spectroscopy demonstrate that heparin and enoxaparin, a low-molecular-weight heparin which is a clinical anticoagulant, bind and induce a conformational change in the spike (S1) protein receptor-binding domain (S1 RBD) of SARS-CoV-2. A library of heparin derivatives and size-defined fragments were used to probe the structural basis of this interaction. Binding to the RBD is more strongly dependent on the presence of 2-O or 6-O sulfate groups than on N-sulfation and a hexasaccharide is the minimum size required for secondary structural changes to be induced in the RBD. It is likely that inhibition of viral infection arises from an overlap between the binding sites of heparin/HS on S1 RBD and that of the angiotensin-converting enzyme 2. The results suggest a route for the rapid development of a first-line therapeutic by repurposing heparin and its derivatives as antiviral agents against SARS-CoV-2 and other members of the Coronaviridae.

Heparin Inhibits Cellular Invasion by SARS-CoV-2: Structural Dependence of the Interaction of the Spike S1 Receptor-Binding Domain with Heparin

The dependence of development and homeostasis in animals on the interaction of hundreds of extracellular regulatory proteins with the peri- and extracellular glycosaminoglycan heparan sulfate (HS) is exploited by many microbial pathogens as a means of adherence and invasion. Heparin, a widely used anticoagulant drug, is structurally similar to HS and is a common experimental proxy. Exogenous heparin prevents infection by a range of viruses, including S-associated coronavirus isolate HSR1. Here, we show that heparin inhibits severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) invasion of Vero cells by up to 80% at doses achievable through prophylaxis and, particularly relevant, within the range deliverable by nebulisation. Surface plasmon resonance and circular dichroism spectroscopy demonstrate that heparin and enoxaparin, a low-molecular-weight heparin which is a clinical anticoagulant, bind and induce a conformational change in the spike (S1) protein receptor-binding domain (S1 RBD) of SARS-CoV-2. A library of heparin derivatives and size-defined fragments were used to probe the structural basis of this interaction. Binding to the RBD is more strongly dependent on the presence of 2-O or 6-O sulfate groups than on N-sulfation and a hexasaccharide is the minimum size required for secondary structural changes to be induced in the RBD. It is likely that inhibition of viral infection arises from an overlap between the binding sites of heparin/HS on S1 RBD and that of the angiotensin-converting enzyme 2. The results suggest a route for the rapid development of a first-line therapeutic by repurposing heparin and its derivatives as antiviral agents against SARS-CoV-2 and other members of the Coronaviridae.

GlycoTorch Vina: Docking Designed and Tested for Glycosaminoglycans

Glycosaminoglycans (GAGs) are a family of anionic carbohydrates that play an essential role in the physiology and pathology of all eukaryotic life forms. Experimental determination of GAG-protein complexes is challenging due to their difficult isolation from biological sources, natural heterogeneity, and conformational flexibility-including possible ring puckering of sulfated iduronic acid from ¹C₄ to ²S_O conformation. To overcome these challenges, we present GlycoTorch Vina (GTV), a molecular docking tool based on the carbohydrate docking program VinaCarb (VC). Our program is unique in that it contains parameters to model ²S_O sugars while also supporting glycosidic linkages specific to GAGs. We discuss how crystallographic models of carbohydrates can be biased by the choice of refinement software and structural dictionaries. To overcome these variations, we carefully curated 12 of the best available GAG and GAG-like crystal structures (ranging from tetra- to octasaccharides or longer) obtained from the PDB-REDO server and refined using the same protocol. Both GTV and VC produced pose predictions with a mean root-mean-square deviation (RMSD) of 3.1 Å from the native crystal structure-a statistically significant improvement when compared to AutoDock Vina (4.5 Å) and the commercial software Glide (5.9 Å). Examples of how real-space correlation coefficients can be used to better assess the accuracy of docking pose predictions are given. Comparisons between statistical distributions of empirical "salt bridge" interactions, relevant to GAGs, were compared to density functional theory (DFT) studies of model salt bridges, and water-mediated salt bridges; however, there was generally a poor agreement between these data. Water bridges appear to play an important, yet poorly understood, role in the structures of GAG-protein complexes. To aid in the rapid prototyping of future pose scoring functions, we include a module that allows users to include their own torsional and nonbonded parameters.

Effect of hydroxylysine-O-glycosylation on the structure of type I collagen molecule: A computational study

Collagen undergoes many types of post-translational modifications (PTMs), including intracellular modifications and extracellular modifications. Among these PTMs, glycosylation of hydroxylysine (Hyl) is the most complicated. Experimental studies demonstrated that this PTM ceases once the collagen triple helix is formed and that Hyl-O-glycosylation modulates collagen fibrillogenesis. However, the underlying atomic-level mechanisms of these phenomena remain unclear. In this study, we first adapted the force field parameters for O-linkages between Hyl and carbohydrates and then investigated the influence of Hyl-O-glycosylation on the structure of type I collagen molecule, by performing comprehensive molecular dynamic simulations in explicit solvent of collagen molecule segment with and without the glycosylation of Hyl. Data analysis demonstrated that (i) collagen triple helices remain in a triple-helical structure upon glycosylation of Hyl; (ii) glycosylation of Hyl modulates the peptide backbone conformation and their solvation environment in the vicinity and (iii) the attached sugars are arranged such that their hydrophilic faces are well exposed to the solvent, while their hydrophobic faces point towards the hydrophobic portions of collagen. The adapted force field parameters for O-linkages between Hyl and carbohydrates will aid future computational studies on proteins with Hyl-O-glycosylation. In addition, this work, for the first time, presents the detailed effect of Hyl-O-glycosylation on the structure of human type I collagen at the atomic level, which may provide insights into the design and manufacture of collagenous biomaterials and the development of biomedical therapies for collagen-related diseases.

The role of SP-B1–25 peptides in lung surfactant monolayers exposed to gold nanoparticles

Lung surfactant monolayer’s (acts as the first line barrier for inhaled nanoparticles) components (lipids and peptides) rearrange themselves by the influence of exposed gold nanoparticles at various stages of the breathing cycle.

Barrier-to-autointegration factor 1 (Banf1) regulates poly [ADP-ribose] polymerase 1 (PARP1) activity following oxidative DNA damage

The DNA repair capacity of human cells declines with age, in a process that is not clearly understood. Mutation of the nuclear envelope protein barrier-to-autointegration factor 1 (Banf1) has previously been shown to cause a human progeroid disorder, Néstor–Guillermo progeria syndrome (NGPS). The underlying links between Banf1, DNA repair and the ageing process are unknown. Here, we report that Banf1 controls the DNA damage response to oxidative stress via regulation of poly [ADP-ribose] polymerase 1 (PARP1). Specifically, oxidative lesions promote direct binding of Banf1 to PARP1, a critical NAD⁺-dependent DNA repair protein, leading to inhibition of PARP1 auto-ADP-ribosylation and defective repair of oxidative lesions, in cells with increased Banf1. Consistent with this, cells from patients with NGPS have defective PARP1 activity and impaired repair of oxidative lesions. These data support a model whereby Banf1 is crucial to reset oxidative-stress-induced PARP1 activity. Together, these data offer insight into Banf1-regulated, PARP1-directed repair of oxidative lesions.

Mutation of the nuclear envelope protein, barrier-to-autointegration factor 1 (Banf1), has previously been associated with the development of ageing associated diseases in a human progeria syndrome. Here, the authors reveal the functional link between Banf1-regulated, PARP1-directed repair of oxidative lesions.

Heparanase: A Challenging Cancer Drug Target

Heparanase has been viewed as a promising anti-cancer drug target for almost two decades, but no anti-heparanase therapy has yet reached the clinic. This endoglycosidase is highly expressed in a variety of malignancies, and its high expression is associated with greater tumor size, more metastases, and a poor prognosis. It was first described as an enzyme cleaving heparan sulfate chains of proteoglycans located in extracellular matrices and on cell surfaces, but this is not its only function. It is a multi-functional protein with activities that are enzymatic and non-enzymatic and which take place both outside of the cell and intracellularly. Knowledge of the crystal structure of heparanase has assisted the interpretation of earlier structure-function studies as well as in the design of potential anti-heparanase agents. This review re-examines the various functions of heparanase in light of the structural data. The functions of the heparanase variant, T5, and structure and functions of heparanase-2 are also examined as these heparanase related, but non-enzymatic, proteins are likely to influence the in vivo efficacy of anti-heparanase drugs. The anti-heparanase drugs currently under development predominately focus on inhibiting the enzymatic activity of heparanase, which, in the absence of inhibitors with high clinical efficacy, prompts a discussion of whether this is the best approach. The diversity of outcomes attributed to heparanase and the difficulties of unequivocally determining which of these are due to its enzymatic activity is also discussed and leads us to the conclusion that heparanase is a valid, but challenging drug target for cancer.

Cyclophilin A allows the allosteric regulation of a structural motif in the disordered domain 2 of NS5A and thereby fine-tunes HCV RNA replication

Implicated in numerous human diseases, intrinsically disordered proteins (IDPs) are dynamic ensembles of interconverting conformers that often contain many proline residues. Whether and how proline conformation regulates the functional aspects of IDPs remains an open question, however. Here, we studied the disordered domain 2 of nonstructural protein 5A (NS5A-D2) of hepatitis C virus (HCV). NS5A-D2 comprises a short structural motif (PW-turn) embedded in a proline-rich sequence, whose interaction with the human prolyl isomerase cyclophilin A (CypA) is essential for viral RNA replication. Using NMR, we show here that the PW-turn motif exists in a conformational equilibrium between folded and disordered states. We found that the fraction of conformers in the NS5A-D2 ensemble that adopt the structured motif is allosterically modulated both by the cis/trans isomerization of the surrounding prolines that are CypA substrates and by substitutions conferring resistance to cyclophilin inhibitor. Moreover, we noted that this fraction is directly correlated with HCV RNA replication efficiency. We conclude that CypA can fine-tune the dynamic ensemble of the disordered NS5A-D2, thereby regulating viral RNA replication efficiency.

Analysis of crystallization phenomenon in Indian honey using molecular dynamics simulations and artificial neural network

Molecular dynamics simulation was performed on sugar profile and moisture content-based mixture systems of six Indian honey samples. Comparative studies were performed to understand the interactive effects of fructose, glucose, sucrose, maltose and water on crystallization. All simulations led to formation of stable crystal but with different interaction energies. Post-simulation analysis showed that Fructose:Glucose of 1.18 formed the most stable crystal with highest van der Waals and electrostatic interactions. The stability of crystal was further validated with least gyration radius (209 ± 1.81 nm²), accessible surface area (4.09 ± 0.04 nm) and root mean square displacement (3.51 ± 0.00261 nm). Results indicated that not only Fructose:Glucose ratio but also sucrose, maltose and water had a significant effect on the overall crystallization process. The simulation data was used to train the artificial neural network which predicted the stability of honey crystallization depending on Fructose:Glucose and Glucose:Water ratios.

Novel analogs of sulfasalazine as system xc - antiporter inhibitors: Insights from the molecular modeling studies

System x_c ^- (Sx_c ^- ), a cystine-glutamate antiporter, is established as an interesting target for the treatment of several pathologies including epileptic seizures, glioma, neurodegenerative diseases, and multiple sclerosis. Erastin, sorafenib, and sulfasalazine (SSZ) are a few of the established inhibitors of Sx_c ^- . However, its pharmacological inhibition with novel and potent agents is still very much required due to potential issues, for example, potency, bioavailability, and blood-brain barrier (BBB) permeability, with the current lead molecules such as SSZ. Therefore, in this study, we report the synthesis and structure-activity relationships (SAR) of SSZ derivatives along with molecular docking and dynamics simulations using the developed homology model of xCT chain of Sx_c ^- antiporter. The generated homology model attempted to address the limitations of previously reported comparative protein models, thereby increasing the confidence in the computational modeling studies. The main objective of the present study was to derive a suitable lead structure from SSZ eliminating its potential issues for the treatment of glioblastoma multiforme (GBM), a deadly and malignant grade IV astrocytoma. The designed compounds with favorable Sx_c ^- inhibitory activity following in vitro Sx_c ^- inhibition studies, showed moderately potent cytotoxicity in patient-derived human glioblastoma cells, thereby generating potential interest in these compounds. The xCT-ligand model can be further optimized in search of potent lead molecules for novel drug discovery and development studies.

Adsorption of Collagen-like Peptides onto Gold Nanosurfaces

The molecular behavior of proteins in the presence of inorganic surfaces is of fundamental biological significance. Examples include extracellular matrix proteins interacting with gold nanoparticles and metallic implant biomaterials, such as titanium and stainless steels. Uncharged inorganic surfaces that interact strongly with the solution phase (hydrophilic surfaces) have been commonly used in disease treatments. A deep understanding of the molecular behavior of body proteins in the presence of hydrophilic surfaces is important in terms of clinical applications. However, the adsorption mechanism of proteins onto hydrophilic surfaces remains not fully understood. Here, comprehensive molecular dynamics simulations are carried out to study the molecular response of a human collagen molecule segment (CMS) to the presence of a planar gold surface (AuNS) in explicit solvent, aiming to unravel the adsorption mechanism of proteins onto hydrophilic surfaces. The results demonstrate that in the presence of AuNS, the CMS first biasedly diffuses toward AuNS, followed by anchoring to the gold surface, and finally adsorbs stepwise onto AuNS, where the protein adjusts its structure to maximize the interaction with AuNS. We conclude that adsorption of proteins onto hydrophilic surfaces adheres to three steps, namely, biased diffusion, anchoring, and stepwise adsorption accompanied by structural adaptation. The obtained adsorption mechanism provides insights into the development of inorganic surfaces for biomedical and therapeutic applications.

Interaction of gold nanosurfaces/nanoparticles with collagen-like peptides

Nanotechnology has quickly emerged as a promising research field with potential effects in disease treatments. For example, gold nanoparticles (AuNPs) have been extensively used in diagnostics and therapeutics. When administrated into human tissues, AuNPs first encounter extracellular matrix (ECM) molecules. Amongst all the ECM components, collagen is the main tension-resisting constituent, whose biofunctional and mechanical properties are strongly dependent on its hierarchical structure. Therefore, an in-depth understanding of the structural response of collagen to the presence of gold nanosurfaces (AuNS) and AuNPs is crucial in terms of clinical applications of AuNPs. However, detailed understanding of the molecular-level and atomic-level interaction between AuNS/AuNPs and collagen in the ECM is elusive. In this study, comprehensive molecular dynamics (MD) simulations have been performed to investigate the molecular behaviour of a collagen molecule segment (CMS) in the presence of AuNS/AuNPs in explicit water, aiming to explore the interaction of AuNS/AuNPs with collagen triple helices at the molecular and atomic levels. The results show that the CMS forms a rapid association with AuNS/AuNPs and undergoes a severe unfolding upon adsorption on AuNS/AuNPs, indicating an unfolding propensity of gold surfaces. We conclude that collagen triple helices unfold readily on AuNS and bare AuNPs, due to the interaction of gold surfaces with the protein backbone. The revealed clear unfolding nature and the unravelled atomic-level unfolding mechanism of collagen triple helices onto AuNPs contribute to the development of AuNPs for biomedical and therapeutic applications, and the design of gold-binding proteins.

Atomistic molecular dynamics simulations of bioactive engrailed 1 interference peptides (EN1-iPeps)

The neural-specific transcription factor Engrailed 1 - is overexpressed in basal-like breast tumours. Synthetic interference peptides - comprising a cell-penetrating peptide/nuclear localisation sequence and the Engrailed 1-specific sequence from the N-terminus have been engineered to produce a strong apoptotic response in tumour cells overexpressing EN1, with no toxicity to normal or non Engrailed 1-expressing cells. Here scaled molecular dynamics simulations were used to study the conformational dynamics of these interference peptides in aqueous solution to characterise their structure and dynamics. Transitions from disordered to α-helical conformation, stabilised by hydrogen bonds and proline-aromatic interactions, were observed throughout the simulations. The backbone of the wild-type peptide folds to a similar conformation as that found in ternary complexes of anterior Hox proteins with conserved hexapeptide motifs important for recognition of pre-B-cell leukemia Homeobox 1, indicating that the motif may possess an intrinsic preference for helical structure. The predicted NMR chemical shifts of these peptides are consistent with the Hox hexapeptides in solution and Engrailed 2 NMR data. These findings highlight the importance of aromatic residues in determining the structure of Engrailed 1 interference peptides, shedding light on the rational design strategy of molecules that could be adopted to inhibit other transcription factors overexpressed in other cancer types, potentially including other transcription factor families that require highly conserved and cooperative protein-protein partnerships for biological activity.

Molecular Dynamics Simulation of Tau Peptides for the Investigation of Conformational Changes Induced by Specific Phosphorylation Patterns

The Tau protein plays an important role due to its biomolecular interactions in neurodegenerative diseases. The lack of stable structure and various posttranslational modifications such as phosphorylation at various sites in the Tau protein pose a challenge for many experimental methods that are traditionally used to study protein folding and aggregation. Atomistic molecular dynamics (MD) simulations can help around deciphering relationship between phosphorylation and various intermediate and stable conformations of the Tau protein which occur on longer timescales. This chapter outlines protocols for the preparation, execution, and analysis of all-atom MD simulations of a 21-amino acid-long phosphorylated Tau peptide with the aim of generating biologically relevant structural and dynamic information. The simulations are done in explicit solvent and starting from nearly extended configurations of the peptide. The scaled MD method implemented in AMBER14 was chosen to achieve enhanced conformational sampling in addition to a conventional MD approach, thereby allowing the characterization of folding for such an intrinsically disordered peptide at 293 K. Emphasis is placed on the analysis of the simulation trajectories to establish correlations with NMR data (i.e., chemical shifts and NOEs). Finally, in-depth discussions are provided for commonly encountered problems.

A Phosphorylation-Induced Turn Defines the Alzheimer's Disease AT8 Antibody Epitope on the Tau Protein

Post mortem biochemical staging of Alzheimer's disease is currently based on immunochemical analysis of brain slices with the AT8 antibody. The epitope of AT8 is described around the pSer202/pThr205 region of the hyperphosphorylated form of the neuronal protein tau. In this study, NMR spectroscopy was used to precisely map the AT8 epitope on phosphorylated tau, and derive its defining structural features by a combination of NMR analyses and molecular dynamics. A particular turn conformation is stabilized by a hydrogen bond of the phosphorylated Thr205 residue to the amide proton of Gly207, and is further stabilized by the two Arg residues opposing the pSer202/pThr205.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225-250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations.

Prediction of heparin binding sites in bone morphogenetic proteins (BMPs)

Heparin is a glycosaminoglycan known to bind bone morphogenetic proteins (BMPs) and the growth and differentiation factors (GDFs) and has strong and variable effects on BMP osteogenic activity. In this paper we report our predictions of the likely heparin binding sites for BMP-2 and 14. The N-terminal sequences upstream of TGF-β-type cysteine-knot domains in BMP-2, 7 and 14 contain the basic residues arginine and lysine, which are key components of the heparin/HS-binding sites, with these residues being highly non-conserved. Importantly, evolutionary conserved surfaces on the beta sheets are required for interactions with receptors and antagonists. Furthermore, BMP-2 has electropositive surfaces on two sides compared to BMP-7 and BMP-14. Molecular docking simulations suggest the presence of high and low affinity binding sites in dimeric BMP-2. Histidines were found to play a role in the interactions of BMP-2 with heparin; however, a pK(a) analysis suggests that histidines are likely not protonated. This is indicative that interactions of BMP-2 with heparin do not require acidic pH. Taken together, non-conserved amino acid residues in the N-terminus and residues protruding from the beta sheet (not overlapping with the receptor binding sites and the dimeric interface) and not C-terminal are found to be important for heparin-BMP interactions.