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.

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.

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

Using synthetic oligosaccharides, we examined how deoxyfluorination (site and pattern) impact the solubility and aggregation of cellulose and chitin oligomers.

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.

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.

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.

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.

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.

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.

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.

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.

Tailoring carbon nitride properties and photoactivity by interfacial engineering of hydrogen-bonded frameworks

The rational synthesis of carbon nitride materials, ranging from polymeric carbon nitride to nitrogen-doped carbon, by supramolecular preorganization of their monomers is a powerful tool for the design of their morphology and photophysical and catalytic activities. Here we show a new facile and scalable approach for the synthesis of ordered CN materials with excellent photoactivity, which consists of supramolecular interfacial preorganization of CN monomers at the interface of two non-miscible solvents. Molecular dynamic simulations supported by experimental results reveal that an appropriate choice of monomers and solvents leads to the formation of a supramolecular assembly solely at the interface of the solvents. As a proof of concept, we show that the properties of the CN materials after thermal condensation can be tuned by adding an additional monomer to one solvent only. The advantages of the new method are demonstrated here through the tunable morphologies and surface area, the formation of new electronic junctions and high activity as a photocatalyst for hydrogen evolution and pollutant degradation of the CN materials.

Divalent Sialylated Precision Glycooligomers Binding to Polyomaviruses and the Effect of Different Linkers

Divalent precision glycooligomers terminating in N-acetylneuraminic acid (Neu5Ac) or 3'-sialyllactose (3'-SL) with varying linkers between scaffold and the glycan portions are synthesized via solid phase synthesis for co-crystallization studies with the sialic acid-binding major capsid protein VP1 of human Trichodysplasia spinulosa-associated Polyomavirus. High-resolution crystal structures of complexes demonstrate that the compounds bind to VP1 depending on the favorable combination of carbohydrate ligand and linker. It is found that artificial linkers can replace portions of natural carbohydrate linkers as long as they meet certain requirements such as size or flexibility to optimize contact area between ligand and receptor binding sites. The obtained results will influence the design of future high affinity ligands based on the structures presented here, and they can serve as a blueprint to develop multivalent glycooligomers as inhibitors of viral adhesion.

Nanodroplets at Membranes Create Tight-Lipped Membrane Necks via Negative Line Tension

The response of biomembranes to aqueous-phase separation and to the resulting water-in-water droplets has been recently studied on the micrometer scale using optical microscopy and elasticity theory. When such a droplet adheres to the membrane, it forms a contact area that is bounded by a contact line. For a micrometer-sized droplet, the line tension associated with this contact line can usually be ignored compared with the surface tensions. However, for a small nanoscopic droplet, this line tension is expected to affect the membrane-droplet morphology. Here, we use molecular simulations to study nanodroplets at membranes and to gain insight into these line tension effects. The latter effects are shown to depend strongly on another key parameter, the mechanical tension experienced by the membrane. For a large membrane tension, a droplet adhering to the membrane is only partially engulfed by the membrane, and the membrane-droplet system exhibits an axisymmetric morphology. A reduction of the membrane tension leads to an increase in the contact area and a decrease in the interfacial area of the droplet, initially retaining its axisymmetric shape, which implies a circular contact line and a circular membrane neck. However, when the tension falls below a certain threshold value, the system undergoes a morphological transition toward a non-axisymmetric morphology with a non-circular membrane neck. This morphology persists until the nanodroplet is completely engulfed by the membrane and the membrane neck has closed into a tight-lipped shape. The latter morphology is caused by a negative line tension, which is shown to be a robust feature of membrane-droplet systems. A closed membrane neck with a tight-lipped shape suppresses both thermally activated and protein-induced scission of the neck, implying a reduction in the cellular uptake of nanodroplets by pinocytosis and fluid-phase endocytosis. Furthermore, based on our results, we can also draw important conclusions about the time-dependent processes corresponding to the surface nucleation and growth of nanodroplets at membranes.

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.

Area Increase and Budding in Giant Vesicles Triggered by Light: Behind the Scene

Biomembranes are constantly remodeled and in cells, these processes are controlled and modulated by an assortment of membrane proteins. Here, it is shown that such remodeling can also be induced by photoresponsive molecules. The morphological control of giant vesicles in the presence of a water-soluble ortho-tetrafluoroazobenzene photoswitch (F-azo) is demonstrated and it is shown that the shape transformations are based on an increase in membrane area and generation of spontaneous curvature. The vesicles exhibit budding and the buds can be retracted by using light of a different wavelength. In the presence of F-azo, the membrane area can increase by more than 5% as assessed from vesicle electrodeformation. To elucidate the underlying molecular mechanism and the partitioning of F-azo in the membrane, molecular dynamics simulations are employed. Comparison with theoretically calculated shapes reveals that the budded shapes are governed by curvature elasticity, that the spontaneous curvature can be decomposed into a local and a nonlocal contribution, and that the local spontaneous curvature is about 1/(2.5 µm). The results show that exo- and endocytotic events can be controlled by light and that these photoinduced processes provide an attractive method to change membrane area and morphology.

Efficient Osmotic Pressure Calculations Using Coarse-Grained Molecular Simulations

Osmotic pressure data is increasingly used to parametrize all-atom simulation Force Fields (FFs), leading to large computational cost for larger molecules. Here, we show that the osmotic pressure can be calculated precisely using transferable coarse-grained FFs obtained from short atomistic simulations using an inhomogeneously regularized coarse-graining procedure. This is demonstrated for carbohydrates, where compared to the equivalent atomistic system, an increase of the computational efficiency by a factor of ≈500 is achieved.

Is the Solution Activity Derivative Sufficient to Parametrize Ion-Ion Interactions? Ions for TIP5P Water

Biomolecular processes involve hydrated ions, and thus molecular simulations of such processes require accurate force-field parameters for these ions. In the best force-fields, both ion-water and anion-cation interactions are explicitly parametrized. First, the ion Lennard-Jones parameters are optimized to reproduce, for example, single ion solvation free energies; then ion-pair interactions are adjusted to match experimental activity or activity derivatives. Here, we apply this approach to derive optimized parameters for concentrated NaCl, KCl, MgCl₂, and CaCl₂ salt solutions, to be used with the TIP5P water model. These parameters are of interest because of a number of desirable properties of the TIP5P water model, especially for the simulation of carbohydrates. The results show, that this state of the art approach is insufficient, because the activity derivative often reaches a plateau near the target experimental value, for a wide range of parameter values. The plateau emerges from the interconversion between different types of ion pairs, so parameters leading to equally good agreement with the target solution activity or activity derivative yield very different solution structures. To resolve this indetermination, a second target property, such as the experimentally determined ion-ion coordination number, is required to uniquely determine anion-cation interactions. Simulations show that combining activity derivatives and coordination number as experimental target properties to parametrize ion-ion interactions, is a powerful method for reliable ion-water force field parametrization, and gives insight into the concentration of contact or solvent shared ion pairs in a wide range of salt concentrations. For the alkali and halide ions Li⁺, Rb⁺, Cs⁺, F^-, Br^-, and I^-, we present ion-water parameters appropriate at infinite dilution only.

Patterns of Flexible Nanotubes Formed by Liquid-Ordered and Liquid-Disordered Membranes

Biological membranes form both intra- and intercellular nanotubes that are used for molecular sorting within single cells and for long-distance connections between different cells. Such nanotubes can also develop from synthetic lipid bilayers in their fluid state. Each nanotube has a large area-to-volume ratio and stably encloses a water channel that is thereby shielded from its surroundings. The tubes are rather flexible and can easily change both their length and their conformation. Here, we study nanotubes formed by liquid-ordered (Lo) and liquid-disordered (Ld) membranes with three lipid components exposed to aqueous mixtures of two polymers, polyethylene glycol (PEG) and dextran. Both types of membranes form striking patterns of nanotubes when we reduce the volume of giant vesicles by osmotic deflation, thereby exposing the two bilayer leaflets of the membranes to polymer solutions of different composition. With decreasing volume, three different patterns are observed corresponding to three distinct vesicle morphologies that reflect the interplay of spontaneous curvature and aqueous phase separation. We show that tube nucleation and growth is governed by two kinetic pathways and that the tubes undergo a novel shape transformation from necklace-like to cylindrical tubes at a certain critical tube length. We deduce the spontaneous curvature generated by the membrane-polymer interactions from the observed vesicle morphologies using three different and independent methods of image analysis. The spontaneous curvature of the Ld membranes is found to be 4.7 times larger than that of the Lo membranes. We also show that these curvatures are generated by weak PEG adsorption onto the membranes, with a binding affinity of about 1.6 kBT per chain. In this way, our study provides a direct connection between nanoscopic membrane shapes and molecular interactions. Our approach is rather general and can be applied to many other systems of interest such as polymersomes or membrane-bound proteins and peptides.

Solution Properties of Hemicellulose Polysaccharides with Four Common Carbohydrate Force Fields

Hemicellulose polysaccharides play an important role in the swelling behavior of the primary plant cell wall, and molecular dynamics simulations provide the means of gaining a concise understanding of the interactions of hemicellulose polysaccharides with water. Here, we compare four of the main polysaccharide force fields (CHARMM36 TIP3P, GROMOS56A6(CARBO) SPC, GLYCAM06h TIP3P, and GLYCAM06h TIP5P) for the most abundant hemicellulose backbone components. In particular, we compare aggregation, diffusion coefficients, system density, and investigate the free energy of hydration of saccharides in water. We find that the saccharides show nonphysical aggregation at low concentrations with the GLYCAM06h TIP3P force field, which can be rectified by the use of the TIP5P water model. As a result of the aggregation, GLYCAM06h TIP3P does not lead to reasonable diffusion coefficients whereas the diffusion coefficients, as well as the system density, agrees best with experimental data for the GLYCAM06h TIP5P force field. Overall, GLYCAM06h TIP5P gives good agreement with experimental free energy of hydration data for small saccharides. In addition, the free energy of hydration for short polysaccharides calculated with the GLYCAM06h TIP5P force field is consistent with the radial distribution functions between the polysaccharides and water, the hydration number of the polysaccharides, and the hydrogen bonds formed in the system.

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.

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.

Nucleotide-dependent lateral and longitudinal interactions in microtubules

Microtubule (MT) stability is related to the hydrolysis of the guanosine triphosphate nucleotide (NT) bound to β-tubulin. However, the molecular mechanism by which the NT state influences the stability of the contacts in the MT lattice remains elusive. Here, we present large-scale atomistic simulations of different tubulin aggregates, including individual dimers, short protofilaments, a small lattice patch, and a piece of the MT lattice with two infinite protofilaments in both NT states. Together with a coarse-grained (CG) analysis of the fluctuations, these simulations highlight several regions of the protein where local changes are induced by the NT state or by the lateral and longitudinal contacts in the aggregates. Additionally, the CG analysis provides an indication of how the structural changes affect the bonds between the proteins. The results suggest a consistent picture of a possible molecular mechanism by which the NT state induces changes in the H1-S2 loop and more stable longitudinal bonds, both of which locate the H1-S2 and M-loop in more favorable positions to form lateral contacts.

Intrinsic bending of microtubule protofilaments

The complex polymerization dynamics of the microtubule (MT) plus end are closely linked to the hydrolysis of the GTP nucleotide bound to the β-tubulin. The destabilization is thought to be associated with the conformational change of the tubulin dimers from the straight conformation in the MT lattice to a curved conformation. It remains under debate whether this transformation is directly related to the nucleotide state, or a consequence of the longitudinal or lateral contacts in the MT lattice. Here, we present large-scale atomistic simulations of short tubulin protofilaments with both nucleotide states, starting from both extreme conformations. Our simulations indicate that both interdimer and intradimer contacts in both GDP and GTP-bound tubulin dimers and protofilaments in solution bend. There are no observable differences between the mesoscopic properties of the contacts in GTP and GDP-bound tubulin or the intradime and interdimer interfaces.

Defining coarse-grained representations of large biomolecules and biomolecular complexes from elastic network models

Coarse-grained (CG) models of large biomolecular complexes enable simulations of these systems over long timescales that are not accessible for atomistic molecular dynamics (MD) simulations. A systematic methodology, called essential dynamics coarse-graining (ED-CG), has been developed for defining coarse-grained sites in a large biomolecule. The method variationally determines the CG sites so that key dynamic domains in the protein are preserved in the CG representation. The original ED-CG method relies on a principal component analysis (PCA) of a MD trajectory. However, for many large proteins and multi-protein complexes such an analysis may not converge or even be possible. This work develops a new ED-CG scheme using an elastic network model (ENM) of the protein structure. In this procedure, the low-frequency normal modes obtained by ENM are used to define dynamic domains and to define the CG representation accordingly. The method is then applied to several proteins, such as the HIV-1 CA protein dimer, ATP-bound G-actin, and the Arp2/3 complex. Numerical results show that ED-CG with ENM (ENM-ED-CG) is much faster than ED-CG with PCA because no MD is necessary. The ENM-ED-CG models also capture functional essential dynamics of the proteins almost as well as those using full MD with PCA. Therefore, the ENM-ED-CG method may be better suited to coarse-grain a very large biomolecule or biomolecular complex that is too computationally expensive to be simulated by conventional MD, or when a high resolution atomic structure is not even available.

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.

Pathway of membrane fusion with two tension-dependent energy barriers

Fusion of bilayer membranes is studied via dissipative particle dynamics (DPD) simulations. A new set of DPD parameters is introduced which leads to an energy barrier for flips of lipid molecules between adhering membranes. A large number of fusion events is monitored for a vesicle in contact with a planar membrane. Several time scales of the fusion process are found to depend exponentially on the membrane tension. This implies an energy barrier of about 10k(B)T for intermembrane flips and a second size-dependent barrier for the nucleation of a small hemifused membrane patch.