Deoxyfluorination tunes the aggregation of cellulose and chitin oligosaccharides and highlights the role of specific hydroxyl groups in the crystallization process.

Cellulose and chitin are abundant structural polysaccharides exploited by nature in a large number of applications thanks to their crystallinity. Chemical modifications are commonly employed to tune polysaccharide physical and mechanical properties, but generate heterogeneous mixtures. Thus, the effect of such modifications is not well understood at the molecular level. In this work, we examined how deoxyfluorination (site and pattern) impact the solubility and aggregation of well-defined cellulose and chitin oligomers. While deoxyfluorination increased solubility in water and lowered the crystallinity of cellulose oligomers, chitin was much less affected by the modification. The OH/F substitution also highlighted the role of specific hydroxyl groups in the crystallization process. This work provides guidelines for the design of cellulose- and chitin-based materials. A similar approach can be imagined to prepare cellulose and chitin analogues capable of withstanding enzymatic degradation.

Systematic Structural Characterization of Chitooligosaccharides Enabled by Automated Glycan Assembly.

Chitin, a polymer composed of ß(1-4)-linked N-acetyl-glucosamine monomers, and its partially deacetylated analogue chitosan, are abundant biopolymers with outstanding mechanical as well as elastic properties. Their degradation products, chitooligosaccharides (COS), can trigger the innate immune response in humans and plants. Both material and biological properties are dependent on polymer length, acetylation, as well as the pH. Without well-defined samples, a complete molecular description of these factors is still missing. Automated glycan assembly (AGA) enabled rapid access to synthetic well-defined COS. Chitin-cellulose hybrid oligomers were prepared as important tools for a systematic structural analysis. Intramolecular interactions, identified by molecular dynamics simulations and NMR analysis, underscore the importance of the chitosan amino group for the stabilization of specific geometries.

A multiscale coarse-grained model to predict the molecular architecture and drug transport properties of modified chitosan hydrogels.

Hydrogels constructed with functionalized polysaccharides are of interest in a multitude of applications, chiefly the design of therapeutic and regenerative formulations. Tailoring the chemical modification of polysaccharide-based hydrogels to achieve specific drug release properties involves the optimization of many tunable parameters, including (i) the type, degree (χ), and pattern of the functional groups, (ii) the water-polymer ratio, and (iii) the drug payload. To guide the design of modified polysaccharide hydrogels for drug release, we have developed a computational toolbox that predicts the structure and physicochemical properties of acylated chitosan chains, and their impact on the transport of drug molecules. Herein, we present a multiscale coarse-grained model to investigate the structure of networks of chitosan chains modified with acetyl, butanoyl, or heptanoyl moieties, as well as the diffusion of drugs doxorubicin (Dox) and gemcitabine (Gem) through the resulting networks. The model predicts the formation of different network structures, in particular the hydrophobically-driven transition from a uniform to a cluster/channel morphology and the formation of fibers of chitin chains. The model also describes the impact of structural and physicochemical properties on drug transport, which was confirmed experimentally by measuring Dox and Gem diffusion through an ensemble of modified chitosan hydrogels.

Automated access to well-defined ionic oligosaccharides.

Ionic polysaccharides are part of many biological events, but lack structural characterisation due to challenging purifications and complex synthesis. Four monosaccharides bearing modifications not found in nature are used for the automated synthesis of a collection of ionic oligosaccharides. Structural analysis reveals how the charge pattern affects glycan conformation.

Spherical Nanovesicles Transform into a Multitude of Nonspherical Shapes.

Nanovesicles are closed, bubblelike surfaces with a diameter between 20 and 200 nm, formed by lipid bilayers and biomembranes. Electron microscopy (EM) studies have shown that these vesicles can attain both spherical and nonspherical shapes. One disadvantage of EM methods is that they provide only a single snapshot of each vesicle. Here, we use molecular dynamics simulations to monitor the morphological transformations of individual nanovesicles. We start with the assembly of spherical vesicles that enclose a certain volume of water and contain a certain total number of lipids. When we reduce their volume, the spherical vesicles are observed to transform into a multitude of nonspherical shapes such as oblates and stomatocytes as well as prolates and dumbbells. This surprising polymorphism can be controlled by redistributing a small fraction of lipids between the inner and outer leaflets of the bilayer membranes. As a consequence, the inner and the outer leaflets experience different mechanical tensions. Small changes in the vesicle volume reduce the overall bilayer tension by 2 orders of magnitude, thereby producing tensionless bilayers. In addition, we show how to determine, for a certain total number of lipids, the unique spherical vesicle for which both leaflet tensions vanish individually. We also compute the local spontaneous curvature of the spherical membranes by identifying the first moment of the spherically symmetric stress profiles across the lipid bilayers with the nanoscopic torque as derived from curvature elasticity. Our study can be extended to other types of lipid membranes and sheds new light on cellular nanovesicles such as exosomes, which are increasingly used as biomarkers and drug delivery systems.

Asymmetrically Branched Precision Glycooligomers Targeting Langerin.

Asymmetrically branched precision glycooligomers are synthesized by solid-phase polymer synthesis for studying multivalent carbohydrate-protein interactions. Through the stepwise assembly of Fmoc-protected oligo(amidoamine) building blocks and Fmoc/Dde-protected lysine, straightforward variation of structural parameters such as the number and length of arms, as well as the number and position of carbohydrate ligands, is achieved. Binding of 1-arm and 3-arm glycooligomers toward lectin receptors langerin and concanavalin A (ConA) was evaluated where the smallest 3-arm glycooligomer shows the highest binding toward langerin, and stepwise elongation of one, two, or all three arms leads to decreased binding. When directly comparing binding toward langerin and ConA, we find that structural variation of the scaffold affects glycomimetic ligand binding differently for the different targets, indicating the potential to tune such ligands not only for their avidity but also for their selectivity toward different lectins.

Tailoring the Chemical Modification of Chitosan Hydrogels to Fine-Tune the Release of a Synergistic Combination of Chemotherapeutics.

Combination chemotherapy with a defined ratio and sequence of drug release is a clinically established and effective route to treat advanced solid tumors. In this context, a growing body of literature demonstrates the potential of hydrogels constructed with chemically modified polysaccharides as depots for controlled release of chemotherapeutics. Identifying the appropriate modification in terms of physicochemical properties of the functional group and its degree of substitution (χ) to achieve the desired release profile for multiple drugs is, however, a complex multivariate problem. To address this issue, we have developed a computational toolbox that models the migration of a drug pair through a hydrated network of polysaccharide chains modified with hydrophobic moieties. In this study, we chose doxorubicin (DOX) and Gemcitabine (GEM) as model drugs, as their synergistic effect against breast cancer has been thoroughly investigated, and chitosan as the model polymer. Our model describes how the modification of chitosan chains with acetyl, butanoyl, and heptanoyl moieties at different values χ governs both the structure of the hydrogel network and drug migration through it. Our experimental data confirm the in silico predictions for both single- and dual-drug release and, most notably, the counterintuitive inversion of release vs χ that occurs when switching from a single- to a dual-drug system. Consensus between predicted and experimental data indicates that acetyl modifications (χ = 32-42%) and butanoyl modifications (χ = 19-24%) provide synergistic GEM/DOX release molar ratios (i.e., 5-10). Collectively, these results demonstrate the potential of this model in guiding the design of chemotherapeutic hydrogels to combat cancer.

Systematic Hydrogen-Bond Manipulations To Establish Polysaccharide Structure-Property Correlations.

A dense hydrogen-bond network is responsible for the mechanical and structural properties of polysaccharides. Random derivatization alters the properties of the bulk material by disrupting the hydrogen bonds, but obstructs detailed structure-function correlations. We have prepared well-defined unnatural oligosaccharides including methylated, deoxygenated, deoxyfluorinated, as well as carboxymethylated cellulose and chitin analogues with full control over the degree and pattern of substitution. Molecular dynamics simulations and crystallographic analysis show how distinct hydrogen-bond modifications drastically affect the solubility, aggregation behavior, and crystallinity of carbohydrate materials. This systematic approach to establishing detailed structure-property correlations will guide the synthesis of novel, tailor-made carbohydrate materials.

Multivalent Binding of Precision Glycooligomers on Soft Glycocalyx Mimicking Hydrogels.

We present a synthetic approach toward soft, glycooligomer-functionalized microgel particles mimicking carbohydrate presenting cell surfaces and analyze their specific binding to a model lectin (Concanavalin A, ConA). Focusing on multivalent presentation, a series of sequence-controlled glycooligomers with varying spacing and number of mannose units was synthesized and analyzed for the resulting glycooligomer-ConA affinity. Both direct binding and inhibition studies show a higher affinity with increasing the number of sugar moieties, but they level off for higher valent systems, indicating steric hindrance. Furthermore, the results suggest that increasing the scaffold length tends to decrease binding due to entropic repulsion, which could be compensated by larger scaffolds able to address multiple ConA binding sites. These findings were consistent in all assays (adhesion, fluorescence, and ITC) regardless of binding partner immobilization, demonstrating that flexible ligands exert similar binding modes in solution and when attached to polymer networks, which is relevant for designing glyco-functionalized materials.

Fucose-Functionalized Precision Glycomacromolecules Targeting Human Norovirus Capsid Protein.

Norovirus infection is the major cause of nonbacterial gastroenteritis in humans and has been the subject of numerous studies investigating the virus's biophysical properties and biochemical function with the aim of deriving novel and highly potent entry inhibitors to prevent infection. Recently, it has been shown that the protruding P domain dimer (P-dimer) of a GII.10 Norovirus strain exhibits two new binding sites for l-fucose in addition to the canonical binding sites. Thus, these sites provide a novel target for the design of multivalent fucose ligands as entry inhibitors of norovirus infections. In this current study, a first generation of multivalent fucose-functionalized glycomacromolecules was synthesized and applied as model structures to investigate the potential targeting of fucose binding sites in human norovirus P-dimer. Following previously established solid phase polymer synthesis, eight precision glycomacromolecules varying in number and position of fucose ligands along an oligo(amidoamine) backbone were obtained and then used in a series of binding studies applying native MS, NMR, and X-ray crystallography. We observed only one fucose per glycomacromolecule binding to one P-dimer resulting in similar binding affinities for all fucose-functionalized glycomacromolecules, which based on our current findings we attribute to the overall size of macromolecular ligands and possibly to steric hindrance.

Generalizing the Concept of Specific Compound Formulation Additives towards Non-Fluorescent Drugs: A Solubilization Study on Potential Anti-Alzheimer-Active Small-Molecule Compounds.

Tailor-made compound formulation additives enable the testing of potential drugs with undesirable pharmacological profiles. A combinatorial approach using Raman microscopy as the readout method is presented to select peptide sequences from large one-bead-one-compound libraries. The resulting peptide-PEG conjugates solubilize potential prophylactic and therapeutic anti-Alzheimer compounds and can be used as specific additives not only for fluorescent but also for non-fluorescent compounds.

The free energy of nanopores in tense membranes.

Membrane nanopores are central players for a range of important cellular membrane remodeling processes as well as membrane rupture. Understanding pore formation in tense membranes requires comprehension of the molecular mechanism of pore formation and the associated free energy change as a function of the membrane tension. Here we propose a scheme to calculate the free energy change associated with the formation of a nanometer sized pore in molecular dynamics simulations as a function of membrane tension, which requires the calculation of only one computationally expensive potential of mean force. We show that membrane elastic theory can be used to estimate the pore formation free energy at different tension values from the free energy change in a relaxed membrane and the area expansion curves of the membranes. We have computed the pore formation free energy for a dipalmitoyl-phosphatidylcholine (DPPC) membrane at two different lateral pressure values, 1 bar and -40 bar, by calculating the potential of mean force acting on the head group of a single lipid molecule. Unrestrained simulations of the closing process confirm that the intermediate states along this reaction coordinate are reasonable and show that hydrophilic indentations spanning half the bilayer connected by a hydrophobic pore segment represent the corresponding high energy transition state. A comparison of the stability of simulated membranes to experiment at high loading rates show that, contrary to expectation, pores form too easily in small simulated membrane patches. This discrepancy originates from a combination of the absence of ions in the simulations and the small membrane size.

Evaluation of enzyme activity predictions for variants of unknown significance in Arylsulfatase A.

Continued advances in variant effect prediction are necessary to demonstrate the ability of machine learning methods to accurately determine the clinical impact of variants of unknown significance (VUS). Towards this goal, the ARSA Critical Assessment of Genome Interpretation (CAGI) challenge was designed to characterize progress by utilizing 219 experimentally assayed missense VUS in the Arylsulfatase A (ARSA) gene to assess the performance of community-submitted predictions of variant functional effects. The challenge involved 15 teams, and evaluated additional predictions from established and recently released models. Notably, a model developed by participants of a genetics and coding bootcamp, trained with standard machine-learning tools in Python, demonstrated superior performance among submissions. Furthermore, the study observed that state-of-the-art deep learning methods provided small but statistically significant improvement in predictive performance compared to less elaborate techniques. These findings underscore the utility of variant effect prediction, and the potential for models trained with modest resources to accurately classify VUS in genetic and clinical research.

Effect of tension and curvature on the chemical potential of lipids in lipid aggregates.

Understanding the factors that influence the free energy of lipids in bilayer membranes is an essential step toward understanding exchange processes of lipids between membranes. In general, both lipid composition and membrane geometry can affect lipid exchange rates between bilayer membranes. Here, the free energy change ΔG(des) for the desorption of dipalmitoyl-phosphatidylcholine (DPPC) lipids from different lipid aggregates has been computed using molecular dynamics simulations and umbrella sampling. The value of ΔG(des) is found to depend strongly on the local properties of the aggregate, in that both tension and curvature lead to an increase in ΔG(des). A detailed analysis shows that the increased desorption free energy for tense bilayers arises from the increased conformational entropy of the lipid tails, which reduces the favorable component -TΔS(L) of the desorption free energy.

Leaflet Tensions Control the Spatio-Temporal Remodeling of Lipid Bilayers and Nanovesicles.

Biological and biomimetic membranes are based on lipid bilayers, which consist of two monolayers or leaflets. To avoid bilayer edges, which form when the hydrophobic core of such a bilayer is exposed to the surrounding aqueous solution, a single bilayer closes up into a unilamellar vesicle, thereby separating an interior from an exterior aqueous compartment. Synthetic nanovesicles with a size below 100 nanometers, traditionally called small unilamellar vesicles, have emerged as potent platforms for the delivery of drugs and vaccines. Cellular nanovesicles of a similar size are released from almost every type of living cell. The nanovesicle morphology has been studied by electron microscopy methods but these methods are limited to a single snapshot of each vesicle. Here, we review recent results of molecular dynamics simulations, by which one can monitor and elucidate the spatio-temporal remodeling of individual bilayers and nanovesicles. We emphasize the new concept of leaflet tensions, which control the bilayers' stability and instability, the transition rates of lipid flip-flops between the two leaflets, the shape transformations of nanovesicles, the engulfment and endocytosis of condensate droplets and rigid nanoparticles, as well as nanovesicle adhesion and fusion. To actually compute the leaflet tensions, one has to determine the bilayer's midsurface, which represents the average position of the interface between the two leaflets. Two particularly useful methods to determine this midsurface are based on the density profile of the hydrophobic lipid chains and on the molecular volumes.

Photomanipulation of Minimal Synthetic Cells: Area Increase, Softening, and Interleaflet Coupling of Membrane Models Doped with Azobenzene-Lipid Photoswitches.

Light can effectively interrogate biological systems in a reversible and physiologically compatible manner with high spatiotemporal precision. Understanding the biophysics of photo-induced processes in bio-systems is crucial for achieving relevant clinical applications. Employing membranes doped with the photolipid azobenzene-phosphatidylcholine (azo-PC), a holistic picture of light-triggered changes in membrane kinetics, morphology, and material properties obtained from correlative studies on cell-sized vesicles, Langmuir monolayers, supported lipid bilayers, and molecular dynamics simulations is provided. Light-induced membrane area increases as high as ≈25% and a ten-fold decrease in the membrane bending rigidity is observed upon trans-to-cis azo-PC isomerization associated with membrane leaflet coupling and molecular curvature changes. Vesicle electrodeformation measurements and atomic force microscopy reveal that trans azo-PC bilayers are thicker than palmitoyl-oleoyl phosphatidylcholine (POPC) bilayers but have higher specific membrane capacitance and dielectric constant suggesting an increased ability to store electric charges across the membrane. Lastly, incubating POPC vesicles with azo-PC solutions results in the insertion of azo-PC in the membrane enabling them to become photoresponsive. All these results demonstrate that light can be used to finely manipulate the shape, mechanical and electric properties of photolipid-doped minimal cell models, and liposomal drug carriers, thus, presenting a promising therapeutic alternative for the repair of cellular disorders.

Budding and Fission of Nanovesicles Induced by Membrane Adsorption of Small Solutes.

Membrane budding and fission are essential cellular processes that produce new membrane compartments during cell and organelle division, for intracellular vesicle trafficking as well as during endo- and exocytosis. Such morphological transformations have also been observed for giant lipid vesicles with a size of many micrometers. Here, we report budding and fission processes of lipid nanovesicles with a size below 50 nm. We use coarse-grained molecular dynamics simulations, by which we can visualize the morphological transformations of individual vesicles. The budding and fission processes are induced by low concentrations of small solutes that absorb onto the outer leaflets of the vesicle membranes. In addition to the solute concentration, we identify the solvent conditions as a second key parameter for these processes. For good solvent conditions, the budding of a nanovesicle can be controlled by reducing the vesicle volume for constant solute concentration or by increasing the solute concentration for constant vesicle volume. After the budding process is completed, the budded vesicle consists of two membrane subcompartments which are connected by a closed membrane neck. The budding process is reversible as we demonstrate explicitly by reopening the closed neck. For poor solvent conditions, on the other hand, we observe two unexpected morphological transformations of nanovesicles. Close to the binodal line, at which the aqueous solution undergoes phase separation, the vesicle exhibits recurrent shape changes with closed and open membrane necks, reminiscent of flickering fusion pores (kiss-and-run) as observed for synaptic vesicles. As we approach the binodal line even closer, the recurrent shape changes are truncated by the fission of the membrane neck which leads to the division of the nanovesicle into two daughter vesicles. In this way, our simulations reveal a nanoscale mechanism for the budding and fission of nanovesicles, a mechanism that arises from the interplay between membrane elasticity and solute-mediated membrane adhesion.

Weak carbohydrate-carbohydrate interactions in membrane adhesion are fuzzy and generic.

Carbohydrates such as the trisaccharide motif LeX are key constituents of cell surfaces. Despite intense research, the interactions between carbohydrates of apposing cells or membranes are not well understood. In this article, we investigate carbohydrate-carbohydrate interactions in membrane adhesion as well as in solution with extensive atomistic molecular dynamics simulations that exceed the simulation times of previous studies by orders of magnitude. For LeX, we obtain association constants of soluble carbohydrates, adhesion energies of lipid-anchored carbohydrates, and maximally sustained forces of carbohydrate complexes in membrane adhesion that are in good agreement with experimental results in the literature. Our simulations thus appear to provide a realistic, detailed picture of LeX-LeX interactions in solution and during membrane adhesion. In this picture, the LeX-LeX interactions are fuzzy, i.e. LeX pairs interact in a large variety of short-lived, bound conformations. For the synthetic tetrasaccharide Lac 2, which is composed of two lactose units, we observe similarly fuzzy interactions and obtain association constants of both soluble and lipid-anchored variants that are comparable to the corresponding association constants of LeX. The fuzzy, weak carbohydrate-carbohydrate interactions quantified in our simulations thus appear to be a generic feature of small, neutral carbohydrates such as LeX and Lac 2.

Sequential Linkage of Carbohydrate Antigens to Mimic Capsular Polysaccharides: Toward Semisynthetic Glycoconjugate Vaccine Candidates against Streptococcus pneumoniae Serotype 14.

Vaccines based on isolated polysaccharides successfully protect humans from bacterial pathogens such as Streptococcus pneumoniae. Because polysaccharide production and isolation can be technically challenging, glycoconjugates containing synthetic antigens are an attractive alternative. Typically, the shortest possible oligosaccharide antigen is preferable as syntheses of longer structures are more difficult and time-consuming. Combining several protective epitopes or polysaccharide repeating units as blocks by bonds other than glycosidic linkages would greatly reduce the synthetic effort if the immunological response to the polysaccharide could be retained. To explore this concept, we bridged the well-understood and immunologically potent RU of S. pneumoniae serotype 14 (ST14) with an aliphatic spacer and conjugated it to the carrier protein CRM197. Mice immunized with the spacer-bridged glycan conjugates produced high levels of specific antibodies after just one or two vaccine doses, while the tetrasaccharide repeating unit alone required three doses. The antibodies recognized specifically ST14 CPS, while no significant antibody levels were raised against the spacer or unrelated CPS. Synthetic vaccines generated antibodies with opsonic activity. Mimicking polysaccharides by coupling repeating unit antigens via an aliphatic spacer may prove useful also for the development of other glycoconjugate vaccine candidates, thereby reducing the synthetic complexity while enhancing a faster immune response.

The fusion of membranes and vesicles: pathway and energy barriers from dissipative particle dynamics.

The fusion of lipid bilayers is studied with dissipative particle dynamics simulations. First, to achieve control over membrane properties, the effects of individual simulation parameters are studied and optimized. Then, a large number of fusion events for a vesicle and a planar bilayer are simulated using the optimized parameter set. In the observed fusion pathway, configurations of individual lipids play an important role. Fusion starts with individual lipids assuming a splayed tail configuration with one tail inserted in each membrane. To determine the corresponding energy barrier, we measure the average work for interbilayer flips of a lipid tail, i.e., the average work to displace one lipid tail from one bilayer to the other. This energy barrier is found to depend strongly on a certain dissipative particle dynamics parameter, and, thus, can be adjusted in the simulations. Overall, three subprocesses have been identified in the fusion pathway. Their energy barriers are estimated to lie in the range 8-15 k(B)T. The fusion probability is found to possess a maximum at intermediate tension values. As one decreases the tension, the fusion probability seems to vanish before the tensionless membrane state is attained. This would imply that the tension has to exceed a certain threshold value to induce fusion.