On the interaction of hyaluronic acid with synovial fluid lipid membranes

Quantifying Induced Dipole Effects in Small Molecule Permeation in a Model Phospholipid Bilayer

The cell membrane functions as a semipermeable barrier that governs the transport of materials into and out of cells. The bilayer features a distinct dielectric gradient due to the amphiphilic nature of its lipid components. This gradient influences various aspects of small molecule permeation and the folding and functioning of membrane proteins. Here, we employ polarizable molecular dynamics simulations to elucidate the impact of the electronic environment on the permeation process. We simulated eight distinct amino-acid side chain analogs within a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayer using the Drude polarizable force field (FF). Our approach includes both unbiased and umbrella sampling simulations. By using a polarizable FF, we sought to investigate explicit dipole responses in relation to local electric fields along the membrane normal. We evaluate molecular dipole moments, which exhibit variation based on their localization within the membrane, and compare the outcomes with analogous simulations using the nonpolarizable CHARMM36 FF. This comparative analysis aims to discern characteristic differences in the free energy surfaces of permeation for the various amino-acid analogs. Our results provide the first systematic quantification of the impact of employing an explicitly polarizable FF in this context compared to the fixed-charge convention inherent to nonpolarizable FFs, which may not fully capture the influence of the membrane dielectric gradient.

Dynamics of an aqueous suspension of short hyaluronic acid chains near a DPPC bilayer

The synergy between hyaluronic acid (HA) and lipid molecules plays a crucial role in synovial fluids, cell coatings, etc. Diseased cells in cancer and arthritis show changes in HA concentration and chain size, impacting the viscoelastic and mechanical properties of the cells. Although the solution behavior of HA is known in experiments, a molecular-level understanding of the role of HA in the dynamics at the interface of HA-water and the cellular boundary is lacking. Here, we perform atomistic molecular dynamics simulation of short HA chains in an explicit water solvent in the presence of a DPPC bilayer, relevant in pathological cases. We identify a stable interface between HA-water and the bilayer where the water molecules are in contact with the bilayer and the HA chains are located away without any direct contact. Both translation and rotation of the interfacial waters in contact with the lipid bilayer and translation of the HA chains exhibit subdiffusive behavior. The diffusive behavior sets in slightly away from the bilayer, where the diffusion coefficients of water and HA decrease monotonically with increase in HA concentration. On the contrary, the dependence on HA chain size is only marginal due to enhanced chain flexibility as their size increases.

Unraveling the Molecular Dance: Insights into TREM2/DAP12 Complex Formation in Alzheimer's Disease through Molecular Dynamics Simulations

Alzheimer's disease (AD) is a widespread neurodegenerative condition affecting millions globally. Recent research has implicated variants of the triggering receptor expressed in myeloid cells 2 (TREM2) as risk factors for AD. TREM2, an immunomodulatory receptor on microglial surfaces, plays a pivotal role in regulating microglial activation by association with DNAX-activation protein 12 (DAP12). Despite its significance, the mechanism underlying the formation of the complex between the transmembrane domains (TMDs) of TREM2 and DAP12 remains unclear. This study employs multiscale molecular dynamics (MD) simulations to investigate three TMD complex models, including two derived from experiments and one generated by AlphaFold2. Conducted within a lipid membrane consisting of an 80:20 mixture of phosphatidylcholine (POPC) and cholesterol, our analysis reveals hydrogen-bonding interactions between K26 of TREM2 and D16 of DAP12 in all three models, consistent with previous experimental findings. Our results elucidate the different spatial conformations observed in the models and offer insights into the structure of the TREM2/DAP12 TMD complex. Furthermore, we elucidate the role of charged residues in the assembly structure of the complex within the lipid membrane. These findings enhance our understanding of the molecular mechanism governing TREM2/DAP12 complex formation, providing a foundation for designing novel therapeutic strategies to address AD and other neurodegenerative diseases.

Martini-Based Coarse-Grained Soil Organic Matter Model Derived from Atomistic Simulations

The significance of soil organic matter (SOM) in environmental contexts, particularly its role in pollutant adsorption, has prompted an increased utilization of molecular simulations to understand microscopic interactions. This study introduces a coarse-grained SOM model, parametrized within the framework of the versatile Martini 3 force field. Utilizing models generated by the Vienna Soil Organic Matter Modeler 2, which constructs humic substance systems from a fragment database, we employed Swarm-CG to parametrize the fragments and subsequently assembled them into macromolecules. Direct Boltzmann inversion (DBI) facilitated the determination of bonded parameters between fragments. The parametrization yielded favorable agreement in the radius of gyration and solvent-accessible surface area. Transfer free energies exhibited a strong correlation with hexadecane-water and chloroform-water values, albeit deviations were noted for octanol-water values. Comparing densities of modeled Leonardite humic acid systems at coarse-grained and atomistic levels revealed promising agreement, particularly at higher water concentrations. The DBI approach effectively reproduced average values of bonded interactions between fragments. Radial distribution functions between carboxylate groups and calcium ions showed partial agreement, however, reproducing certain peaks was challenging due to fixed bead sizes. Detailed analysis of atomistic systems revealed different configurations between the groups, explaining discrepancies. The present contribution provides a comprehensive insight into the properties, strengths, and weaknesses of the coarse-grained SOM model, serving as a foundation for future investigations encompassing pollutant interactions and varied SOM compositions.

Microphase Separation Produces Interfacial Environment within Diblock Biomolecular Condensates

The phase separation of intrinsically disordered proteins is emerging as an important mechanism for cellular organization. However, efforts to connect protein sequences to the physical properties of condensates, i.e., the molecular grammar, are hampered by a lack of effective approaches for probing high-resolution structural details. Using a combination of multiscale simulations and fluorescence lifetime imaging microscopy experiments, we systematically explored a series of systems consisting of diblock elastin-like polypeptides (ELP). The simulations succeeded in reproducing the variation of condensate stability upon amino acid substitution and revealed different microenvironments within a single condensate, which we verified with environmentally sensitive fluorophores. The interspersion of hydrophilic and hydrophobic residues and a lack of secondary structure formation result in an interfacial environment, which explains both the strong correlation between ELP condensate stability and interfacial hydrophobicity scales, as well as the prevalence of protein-water hydrogen bonds. Our study uncovers new mechanisms for condensate stability and organization that may be broadly applicable.

S100A9 interacts with a dynamic region on CD14 to activate Toll-like receptor 4

S100A9 is a Damage Associated Molecular Pattern (DAMP) that activates inflammatory pathways via Toll-like receptor 4 (TLR4). This activity plays important homeostatic roles in tissue repair, but can also contribute to inflammatory diseases. The mechanism of activation is unknown. Here, we follow up on a previous observation that the protein CD14 is an important co-receptor that enables S100A9 to activate TLR4. Using cell-based functional assays and a combination of mutations and pharmocological perturbations, we found that CD14 must be membrane bound to potentiate TLR4 activation by S100A9. Additionally, S100A9 is sensitive to inhibitors of pathways downstream of TLR4 internalization. Together, this suggests that S100A9 induces activity via CD14-dependent internalization of TLR4. We then used mutagenesis, structural modeling, and in vitro binding experiments to establish that S100A9 binds to CD14's N-terminus in a region that overlaps with, but is not identical to, the region where CD14 binds its canonical ligand, lipopolysaccharide (LPS). In molecular dynamics simulations, this region of the protein is dynamic, allowing it to reorganize to recognize both S100A9 (a soluble protein) and LPS (a small hydrophobic molecule). Our work is the first attempt at a molecular characterization of the S100A9/CD14 interaction, bringing us one step closer to unraveling the full mechanism by which S100A9 activates TLR4/MD-2.

Tertiary Plasticity Drives the Efficiency of Enterocin 7B Interactions with Lipid Membranes

The ability of antimicrobial peptides to efficiently kill their bacterial targets depends on the efficiency of their binding to the microbial membrane. In the case of enterocins, there is a three-part interaction: initial binding, unpacking of helices on the membrane surface, and permeation of the lipid bilayer. Helical unpacking is driven by disruption of the peptide hydrophobic core when in contact with membranes. Enterocin 7B is a leaderless enterocin antimicrobial peptide produced from Enterococcus faecalis that functions alone, or with its cognate partner enterocin 7A, to efficiently kill a wide variety of Gram-stain positive bacteria. To better characterize the role that tertiary structural plasticity plays in the ability of enterocin 7B to interact with the membranes, a series of arginine single-site mutants were constructed that destabilize the hydrophobic core to varying degrees. A series of experimental measures of structure, stability, and function, including CD spectra, far UV CD melting profiles, minimal inhibitory concentrations analysis, and release kinetics of calcein, show that decreased stabilization of the hydrophobic core is correlated with increased efficiency of a peptide to permeate membranes and in killing bacteria. Finally, using the computational technique of adaptive steered molecular dynamics, we found that the atomistic/energetic landscape of peptide mechanical unfolding leads to free energy differences between the wild type and its mutants, whose trends correlate well with our experiment.

Dynamics of activation in the voltage-sensing domain of Ciona intestinalis phosphatase Ci-VSP

The Ciona intestinalis voltage-sensing phosphatase (Ci-VSP) is a membrane protein containing a voltage-sensing domain (VSD) that is homologous to VSDs from voltage-gated ion channels responsible for cellular excitability. Previously published crystal structures of Ci-VSD in putative resting and active conformations suggested a helical-screw voltage sensing mechanism in which the S4 helix translocates and rotates to enable exchange of salt-bridge partners, but the microscopic details of the transition between the resting and active conformations remained unknown. Here, by combining extensive molecular dynamics simulations with a recently developed computational framework based on dynamical operators, we elucidate the microscopic mechanism of the resting-active transition at physiological membrane potential. Sparse regression reveals a small set of coordinates that distinguish intermediates that are hidden from electrophysiological measurements. The intermediates arise from a noncanonical helical-screw mechanism in which translocation, rotation, and side-chain movement of the S4 helix are only loosely coupled. These results provide insights into existing experimental and computational findings on voltage sensing and suggest ways of further probing its mechanism.

An extended Tudor domain within Vreteno interconnects Gtsf1L and Ago3 for piRNA biogenesis in Bombyx mori

Piwi-interacting RNAs (piRNAs) direct PIWI proteins to transposons to silence them, thereby preserving genome integrity and fertility. The piRNA population can be expanded in the ping-pong amplification loop. Within this process, piRNA-associated PIWI proteins (piRISC) enter a membraneless organelle called nuage to cleave their target RNA, which is stimulated by Gtsf proteins. The resulting cleavage product gets loaded into an empty PIWI protein to form a new piRISC complex. However, for piRNA amplification to occur, the new RNA substrates, Gtsf-piRISC, and empty PIWI proteins have to be in physical proximity. In this study, we show that in silkworm cells, the Gtsf1 homolog BmGtsf1L binds to piRNA-loaded BmAgo3 and localizes to granules positive for BmAgo3 and BmVreteno. Biochemical assays further revealed that conserved residues within the unstructured tail of BmGtsf1L directly interact with BmVreteno. Using a combination of AlphaFold modeling, atomistic molecular dynamics simulations, and in vitro assays, we identified a novel binding interface on the BmVreteno-eTudor domain, which is required for BmGtsf1L binding. Our study reveals that a single eTudor domain within BmVreteno provides two binding interfaces and thereby interconnects piRNA-loaded BmAgo3 and BmGtsf1L.

Atomistic simulations of calcium aluminosilicate interfaced with liquid water

The dissolution behavior of calcium aluminosilicate based glass fibers, such as stone wool fibers, is an important consideration in mineral wool applications for both the longevity of the mineral wool products in humid environments and limiting the health impacts of released and inhaled fibers from the mineral wool product. Balancing these factors requires a molecular-level understanding of calcium aluminosilicate glass dissolution mechanisms, details that are challenging to resolve with experiment alone. Molecular dynamics simulations are a powerful tool capable of providing complementary atomistic insights regarding dissolution; however, they require force fields capable of describing not-only the calcium aluminosilicate surface structure but also the interactions relevant to dissolution phenomena. Here, a new force field capable of describing amorphous calcium aluminosilicate surfaces interfaced with liquid water is developed by fitting parameters to experimental and first principles simulation data of the relevant oxide-water interfaces, including ab initio molecular dynamics simulations performed for this work for the wüstite and periclase interfaces. Simulations of a calcium aluminosilicate surface interfaced with liquid water were used to test this new force field, suggesting moderate ingress of water into the porous glass interface. This design of the force field opens a new avenue for the further study of calcium and network-modifier dissolution phenomena in calcium aluminosilicate glasses and stone wool fibers at liquid water interfaces.

Correlation between chemical denaturation and the unfolding energetics of Acanthamoeba actophorin

The actin filament network is in part remodeled by the action of a family of filament severing proteins that are responsible for modulating the ratio between monomeric and filamentous actin. Recent work on the protein actophorin from the amoeba Acanthamoeba castellani identified a series of site-directed mutations that increase the thermal stability of the protein by 22°C. Here, we expand this observation by showing that the mutant protein is also significantly stable to both equilibrium and kinetic chemical denaturation, and employ computer simulations to account for the increase in thermal or chemical stability through an accounting of atomic-level interactions. Specifically, the potential of mean force (PMF) can be obtained from steered molecular dynamics (SMD) simulations in which a protein is unfolded. However, SMD can be inefficient for large proteins as they require large solvent boxes, and computationally expensive as they require increasingly many SMD trajectories to converge the PMF. Adaptive steered molecular dynamics (ASMD) overcomes the second of these limitations by steering the particle in stages, which allows for convergence of the PMF using fewer trajectories compared with SMD. Use of the telescoping water scheme within ASMD partially overcomes the first of these limitations by reducing the number of waters at each stage to only those needed to solvate the structure within a given stage. In the PMFs obtained from ASMD, the work of unfolding Acto-2 was found to be higher than the Acto-WT by approximately 120 kCal/mol and reflects the increased stability seen in the chemical denaturation experiments. The evolution of the average number of hydrogen bonds and number of salt bridges during the pulling process provides a mechanistic view of the structural changes of the actophorin protein as it is unfolded, and how it is affected by the mutation in concert with the energetics reported through the PMF.

Polysaccharides-Naturally Occurring Immune Modulators

The prevention of disease and infection requires immune systems that operate effectively. This is accomplished by the elimination of infections and abnormal cells. Immune or biological therapy treats disease by either stimulating or inhibiting the immune system, dependent upon the circumstances. In plants, animals, and microbes, polysaccharides are abundant biomacromolecules. Due to the intricacy of their structure, polysaccharides may interact with and impact the immune response; hence, they play a crucial role in the treatment of several human illnesses. There is an urgent need for the identification of natural biomolecules that may prevent infection and treat chronic disease. This article addresses some of the naturally occurring polysaccharides of known therapeutic potential that have already been identified. This article also discusses extraction methods and immunological modulatory capabilities.

The Colloidal Stability of Apolar Nanoparticles in Solvent Mixtures

Solvent engineering is a powerful and versatile method to tune colloidal stability. Here, we link the molecular structure of apolar ligand shells on gold nanoparticles with their colloidal stability in solvent mixtures. The agglomeration temperature of the particles was measured with small-angle X-ray scattering. It depended on solvent composition and changed linearly for hexane-hexadecane mixtures, but nonlinearly for cyclohexane-hexadecane and hexanol-hexadecane mixtures. Molecular dynamics (MD) simulations indicate that agglomeration is dominated by temperature-dependent ligand order in the alkane mixtures and that the temperature at which the ligand shell orders depends on the solvent composition near the ligands, which can differ substantially from the bulk composition. Small-angle neutron scattering confirmed that, at intermediate solvent compositions above the agglomeration temperature, the fraction of cyclohexane near the ligands was larger than in the bulk. The enrichment of cyclohexane near the ligands stabilized their disordered state, which, consequently, led to the experimentally observed nonlinear trend of the agglomeration temperature. In contrast, hexanol was depleted from the ligand shell at all temperatures. This again stabilized the disordered state. Furthermore, we found that agglomeration at high hexanol fractions was driven by a solvophobic effect that exceeded the influence of ligand order. The results show that strong nonlinearities in the colloidal stability of nanoparticle dispersions in solvent mixtures are directly linked to the molecular details of ligand-solvent and solvent-solvent interactions, which can be used to precisely tune stability.

Insight into the assembly of lipid-hyaluronan complexes in osteoarthritic conditions

Interactions between molecules in the synovial fluid and the cartilage surface may play a vital role in the formation of adsorbed films that contribute to the low friction of cartilage boundary lubrication. Osteoarthritis (OA) is the most common degenerative joint disease. Previous studies have shown that in OA-diseased joints, hyaluronan (HA) not only breaks down resulting in a much lower molecular weight (MW), but also its concentration is reduced ten times. Here, we have investigated the structural changes of lipid-HA complexes as a function of HA concentration and MW to simulate the physiologically relevant conditions that exist in healthy and diseased joints. Small angle neutron scattering and dynamic light scattering were used to determine the structure of HA-lipid vesicles in bulk solution, while a combination of atomic force microscopy and quartz crystal microbalance was applied to study their assembly on a gold surface. We infer a significant influence of both MW and HA concentrations on the structure of HA-lipid complexes in bulk and assembled on a gold surface. Our results suggest that low MW HA cannot form an amorphous layer on the gold surface, which is expected to negatively impact the mechanical integrity and longevity of the boundary layer and could contribute to the increased wear of the cartilage that has been reported in joints diseased with OA.

A Polarized Gel Electrolyte for Wide-Temperature Flexible Zinc-Air Batteries

The development of flexible zinc-air batteries (FZABs) has attracted broad attention in the field of wearable electronic devices. Gel electrolyte is one of the most important components in FZABs, which is urgent to be optimized to match with Zn anode and adapt to severe climates. In this work, a polarized gel electrolyte of polyacrylamide-sodium citric (PAM-SC) is designed for FZABs, in which the SC molecules contain large amount of polarized -COO^- functional groups. The polarized -COO^- groups can form an electrical field between gel electrolyte and Zn anode to suppress Zn dendrite growth. Besides, the -COO^- groups in PAM-SC can fix H₂ O molecules, which prevents water from freezing and evaporating. The polarized PAM-SC hydrogel delivers a high ionic conductivity of 324.68 mS cm^-1 and water retention of 96.85 % after being exposed for 96 h. FZABs with the PAM-SC gel electrolyte exhibit long cycling life of 700 cycles at -40 °C, showing the application prospect under extreme conditions.

Structural Model of a Porphyromonas gingivalis type IX Secretion System Shuttle Complex

Porphyromonas gingivalis is a gram-negative oral anaerobic pathogen and is one of the key causative agents of periodontitis. P. gingivalis utilises a range of virulence factors, including the cysteine protease RgpB, to drive pathogenesis and these are exported and attached to the cell surface via the type IX secretion system (T9SS). All cargo proteins possess a conserved C-terminal signal domain (CTD) which is recognised by the T9SS, and the outer membrane β-barrel protein PorV (PG0027/LptO) can interact with cargo proteins as they are exported to the bacterial surface. Using a combination of solution nuclear magnetic resonance (NMR) spectroscopy, biochemical analyses, machine-learning-based modelling and molecular dynamics (MD) simulations, we present a structural model of a PorV:RgpB-CTD complex from P. gingivalis. This is the first structural insight into CTD recognition by the T9SS and shows how the conserved motifs in the CTD are the primary sites that mediate binding. In PorV, interactions with extracellular surface loops are important for binding the CTD, and together these appear to cradle and lock RgpB-CTD in place. This work provides insight into cargo recognition by PorV but may also have important implications for understanding other aspects of type-IX dependent secretion.

Collagen type II-hyaluronan interactions - the effect of proline hydroxylation: a molecular dynamics study

Hyaluronan-collagen composites have been employed in numerous biomedical applications. Understanding the interactions between hyaluronan and collagen is particularly important in the context of joint cartilage function and the treatment of joint diseases. Many factors affect the affinity of collagen for hyaluronan. One of the important factors is the ratio of 3- or 4-hydroxy proline to proline residues. This article presents the results from molecular dynamics calculations of HA-collagen type II interactions with hyaluronan. The applied protocol employed docking and geometry optimization of complexes built using collagen structures with different numbers of hydroxyl groups attached to proline moieties. It was established that the hydroxyproline/proline ratio affects both structural and energetic features of the collagen-hyaluronan complex. Proline hydroxylation was found to significantly influence the number of all identified types of molecular forces, hydrophobic interactions, water bridges and hydrogen bonds, which can be formed between collagen and hyaluronan. Importantly, an increase in the hydroxyproline/proline ratio in the collagen chain increases the binding affinity for hyaluronan. This is illustrated by the linear correlation between the binding free energy and the hydroxylation degree. A comparison of the results obtained for 3 and 4 hydroxylation of proline indicates that the hydroxyl group attachment position plays a minor role in complex stabilization. However, a slightly stronger affinity was observed for 4 hydroxylation. In order to evaluate the effect of the aqueous environment on the collagen-hyaluronan complex stability, the enthalpic and entropic contributions to the free energy of solvation were analyzed.

Atomic-Resolution Experimental Structural Biology and Molecular Dynamics Simulations of Hyaluronan and Its Complexes

This review summarizes the atomic-resolution structural biology of hyaluronan and its complexes available in the Protein Data Bank, as well as published studies of atomic-resolution explicit-solvent molecular dynamics simulations on these and other hyaluronan and hyaluronan-containing systems. Advances in accurate molecular mechanics force fields, simulation methods and software, and computer hardware have supported a recent flourish in such simulations, such that the simulation publications now outnumber the structural biology publications by an order of magnitude. In addition to supplementing the experimental structural biology with computed dynamic and thermodynamic information, the molecular dynamics studies provide a wealth of atomic-resolution information on hyaluronan-containing systems for which there is no atomic-resolution structural biology either available or possible. Examples of these summarized in this review include hyaluronan pairing with other hyaluronan molecules and glycosaminoglycans, with ions, with proteins and peptides, with lipids, and with drugs and drug-like molecules. Despite limitations imposed by present-day computing resources on system size and simulation timescale, atomic-resolution explicit-solvent molecular dynamics simulations have been able to contribute significant insight into hyaluronan's flexibility and capacity for intra- and intermolecular non-covalent interactions.

A Mechanoelastic Glimpse on Hyaluronan-Coated Extracellular Vesicles

Cancer cells secrete extracellular vesicles (EVs) covered with a carbohydrate polymer, hyaluronan (HA), linked to tumor malignancy. Herein, we have unravelled the contour lengths of HA on a single cancer cell-derived EV surface using single-molecule force spectroscopy (SMFS), which divulges the presence of low molecular weight HA (LMW-HA < 200 kDa). We also discovered that these LMW-HA-EVs are significantly more elastic than the normal cell-derived EVs. This intrinsic elasticity of cancer EVs could be directly allied to the LMW-HA abundance and associated labile water network on EV surface as revealed by correlative SMFS, hydration dynamics with fluorescence spectroscopy, and molecular dynamics simulations. This method emerges as a molecular biosensor of the cancer microenvironment.

The Effects of Cholesterol Oxidation on Erythrocyte Plasma Membranes: A Monolayer Study

Cholesterol plays a key role in the molecular and mesoscopic organisation of lipid membranes and it is expected that changes in its molecular structure (e.g., through environmental factors such as oxidative stress) may affect adversely membrane properties and function. In this study, we present evidence that oxidation of cholesterol has significant effects on the mechanical properties, molecular and mesoscopic organisation and lipid-sterol interactions in condensed monolayers composed of the main species found in the inner leaflet of the erythrocyte membrane. Using a combination of experimental methods (static area compressibility, surface dilatational rheology, fluorescence microscopy, and surface sensitive X-ray techniques) and atomistic molecular dynamics simulations, we show that oxidation of cholesterol to 7-ketocholesterol leads to stiffening of the monolayer (under both static and dynamic conditions), significant changes in the monolayer microdomain organisation, disruption in the van der Waals, electrostatic and hydrophobic interactions between the sterol and the other lipid species, and the lipid membrane hydration. Surface sensitive X-ray techniques reveal that, whilst the molecular packing mode is not significantly affected by cholesterol oxidation in these condensed phases, there are subtle changes in membrane thickness and a significant decrease in the coherence length in monolayers containing 7-ketocholesterol.

Implementation of Telescoping Boxes in Adaptive Steered Molecular Dynamics

Long-time dynamical processes, such as those involving protein unfolding and ligand interactions, can be accelerated and realized through steered molecular dynamics (SMD). The challenge has been the extraction of information from such simulations that generalize for complex nonequilibrium processes. The use of Jarzynski's equality opened the possibility of determining the free energy along the steered coordinate, but sampling over the nonequilibrium trajectories is slow to converge. Adaptive steered molecular dynamics (ASMD) and other related techniques have been introduced to overcome this challenge through the use of stages. Here, we take advantage of these stages to address the numerical cost that arises from the required use of very large solvent boxes. We introduce telescoping box schemes within adaptive steered molecular dynamics (ASMD) in which we adjust the solvent box between stages and thereby vary (and optimize) the required number of solvent molecules. We have benchmarked the method on a relatively long α-helical peptide, Ala₃₀, with respect to the potential of mean force and hydrogen bonds. We show that the use of telescoping boxes introduces little numerical error while significantly reducing the computational cost.

General Characteristics, Biomedical and Dental Application, and Usage of Chitosan in the Treatment of Temporomandibular Joint Disorders: A Narrative Review

The aim of this narrative review was to present research investigating chitosan, including its general characteristics, properties, and medical and dental applications, and finally to present the current state of knowledge regarding the efficacy of chitosan in the treatment of temporomandibular disorders (TMDs) based on the literature. The PICO approach was used for the literature search strategy. The PubMed database was analyzed with the following keywords: (“chitosan”[MeSH Terms] OR “chitosan”[All Fields] OR “chitosans”[All Fields] OR “chitosan s”[All Fields] OR “chitosane”[All Fields]) AND (“temporomandibular joint”[MeSH Terms] OR (“tem-poromandibular”[All Fields] AND “joint”[All Fields]) OR “temporomandibular joint”[All Fields] OR (“temporomandibular”[All Fields] AND “joints”[All Fields]) OR “temporo-mandibular joints”[All Fields]). After screening 8 results, 5 studies were included in this review. Chitosan presents many biological properties and therefore it can be widely used in several branches of medicine and dentistry. Chitosan promotes wound healing, helps to control bleeding, and is used in wound dressings, such as sutures and artificial skin. Apart from its antibacterial property, chitosan has many other properties, such as antifungal, mucoadhesive, anti-inflammatory, analgesic, antioxidant, antihyperglycemic, and antitumoral properties. Further clinical studies assessing the efficacy of chitosan in the treatment of TMD are required. According to only one clinical study, chitosan was effective in the treatment of TMD; however, better clinical results were obtained with platelet-rich plasma.

Cryo-EM structure of PepT2 reveals structural basis for proton-coupled peptide and prodrug transport in mammals

The SLC15 family of proton-coupled solute carriers PepT1 and PepT2 play a central role in human physiology as the principal route for acquiring and retaining dietary nitrogen. A remarkable feature of the SLC15 family is their extreme substrate promiscuity, which has enabled the targeting of these transporters for the improvement of oral bioavailability for several prodrug molecules. Although recent structural and biochemical studies on bacterial homologs have identified conserved sites of proton and peptide binding, the mechanism of peptide capture and ligand promiscuity remains unclear for mammalian family members. Here, we present the cryo-electron microscopy structure of the outward open conformation of the rat peptide transporter PepT2 in complex with an inhibitory nanobody. Our structure, combined with molecular dynamics simulations and biochemical and cell-based assays, establishes a framework for understanding peptide and prodrug recognition within this pharmaceutically important transporter family.

On the Structure and Flip-flop of Free Docosahexaenoic Acid in a Model Human Brain Membrane

Among the omega-3 fatty acids, docosahexaenoic acid (DHA, sn22:6) is particularly vital in human brain cell membranes. There is considerable interest in DHA because low-level DHA has been associated with declined cognitive function and poor visual acuity. In this work, atomistic molecular dynamics simulations were used to investigate the effects of free protonated DHA (DHAP) in molar fractions of 0, 17, 30, and 38% in a realistic model of a healthy brain cell membrane comprising 26 lipid types. Numerous flip-flop events of DHAP were observed and categorized as successful or aborted. Novel use of the machine learning technique, density-based spatial clustering of applications with noise (DBSCAN), effectively identified flip-flop events by way of clustering. Our data show that increasing amounts of DHAP in the membrane disorder the bilayer packing, fluidize the membrane, and increase the rates of successful flip-flop from k = 0.2 μs^-1 (17% DHAP) to 0.8 μs^-1 (30% DHAP) and to 1.3 μs^-1 (38% DHAP). In addition, we also provided a detailed understanding of the flip-flop mechanism of DHAP across this complex membrane. Interestingly, we noted the role of hydrogen bonds in two distinct coordinated flip-flop phenomena between two DHAP molecules: double flip-flop and assisted flip-flop. Understanding the effects of various concentrations of DHAP on the dynamics within a lipid membrane and the resulting structural properties of the membrane are important when considering the use of DHAP as a dietary supplement or as a potential therapeutic.

Mutation of a conserved glutamine residue does not abolish desensitization of acid-sensing ion channel 1

Desensitization is a common feature of ligand-gated ion channels, although the molecular cause varies widely between channel types. Mutations that greatly reduce or nearly abolish desensitization have been described for many ligand-gated ion channels, including glutamate, GABA, glycine, and nicotinic receptors, but not for acid-sensing ion channels (ASICs) until recently. Mutating Gln276 to a glycine (Q276G) in human ASIC1a was reported to mostly abolish desensitization at both the macroscopic and the single channel levels, potentially providing a valuable tool for subsequent studies. However, we find that in both human and chicken ASIC1, the effect of Q276G is modest. In chicken ASIC1, the equivalent Q277G slightly reduces desensitization when using pH 6.5 as a stimulus but desensitizes, essentially like wild-type, when using more acidic pH values. In addition, steady-state desensitization is intact, albeit right-shifted, and recovery from desensitization is accelerated. Molecular dynamics simulations indicate that the Gln277 side chain participates in a hydrogen bond network that might stabilize the desensitized conformation. Consistent with this, destabilizing this network with the Q277N or Q277L mutations largely mimics the Q277G phenotype. In human ASIC1a, the Q276G mutation also reduces desensitization, but not to the extent reported previously. Interestingly, the kinetic consequences of Q276G depend on the human variant used. In the common G212 variant, Q276G slows desensitization, while in the rare D212 variant desensitization accelerates. Our data reveal that while the Q/G mutation does not abolish or substantially impair desensitization as previously reported, it does point to unexpected differences between chicken and human ASICs and the need for careful scrutiny before using this mutation in future studies.

KEAP1 Cancer Mutants: A Large-Scale Molecular Dynamics Study of Protein Stability

We have performed 280 μs of unbiased molecular dynamics (MD) simulations to investigate the effects of 12 different cancer mutations on Kelch-like ECH-associated protein 1 (KEAP1) (G333C, G350S, G364C, G379D, R413L, R415G, A427V, G430C, R470C, R470H, R470S and G476R), one of the frequently mutated proteins in lung cancer. The aim was to provide structural insight into the effects of these mutants, including a new class of ANCHOR (additionally NRF2-complexed hypomorph) mutant variants. Our work provides additional insight into the structural dynamics of mutants that could not be analyzed experimentally, painting a more complete picture of their mutagenic effects. Notably, blade-wise analysis of the Kelch domain points to stability as a possible target of cancer in KEAP1. Interestingly, structural analysis of the R470C ANCHOR mutant, the most prevalent missense mutation in KEAP1, revealed no significant change in structural stability or NRF2 binding site dynamics, possibly indicating an covalent modification as this mutant's mode of action.

A Bittersweet Computational Journey among Glycosaminoglycans

Glycosaminoglycans (GAGs) are linear polysaccharides. In proteoglycans (PGs), they are attached to a core protein. GAGs and PGs can be found as free molecules, associated with the extracellular matrix or expressed on the cell membrane. They play a role in the regulation of a wide array of physiological and pathological processes by binding to different proteins, thus modulating their structure and function, and their concentration and availability in the microenvironment. Unfortunately, the enormous structural diversity of GAGs/PGs has hampered the development of dedicated analytical technologies and experimental models. Similarly, computational approaches (in particular, molecular modeling, docking and dynamics simulations) have not been fully exploited in glycobiology, despite their potential to demystify the complexity of GAGs/PGs at a structural and functional level. Here, we review the state-of-the art of computational approaches to studying GAGs/PGs with the aim of pointing out the “bitter” and “sweet” aspects of this field of research. Furthermore, we attempt to bridge the gap between bioinformatics and glycobiology, which have so far been kept apart by conceptual and technical differences. For this purpose, we provide computational scientists and glycobiologists with the fundamentals of these two fields of research, with the aim of creating opportunities for their combined exploitation, and thereby contributing to a substantial improvement in scientific knowledge.

Membrane Surface Modulates Slow Diffusion in Small Crowded Droplets

Membranes are ubiquitous structures in cells. The effects of membranes on various functional molecules have been reported, but their behaviors under macromolecular crowding and cell-sized confinement have not fully been understood. In this study, we model an intracellular environment by crowding micrometer-sized droplets and investigate the effects of membrane properties on molecular diffusion. The molecular diffusion inside small droplets covered with a lipid layer of phosphatidylcholine (PC) becomes slower compared with that of the corresponding bulk solutions under a crowding condition of polysaccharide dextran but not of its monomer unit, glucose. The addition of a poly(ethylene glycol) conjugated lipid (PEGylated lipid) to the PC membrane significantly alters the degree of slow diffusion observed inside small droplets of concentrated dextran. Interestingly, the change is not monotonic against dextran concentration; that is, the PEGylated membrane increases and decreases the degree of slow diffusion with increasing dextran concentration. We explain the nonmonotonic alternation from the increase in effective dextran concentration and the hindered temporal adsorption of dextran to the membrane. Because diffusion alteration by adding PEGylated lipid is observed for condensed small droplets of linear polymer PEG and hydrophilic protein bovine serum albumin, the phenomenon is general for other polymer systems as well. Furthermore, our findings may facilitate the understanding of intracellular molecular behaviors based on membrane effects as well as the development of numerous applications using polymer droplets.

When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins

Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction.

Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.

Two Coexisting Membrane Structures Are Defined by Lateral and Transbilayer Interactions between Sphingomyelin and Cholesterol

The structure of fully hydrated bilayers composed of equimolar proportions of palmitoylsphingomyelin (PSM) and cholesterol has been examined by synchrotron X-ray powder diffraction and atomistic molecular dynamics (MD) simulations. Two coexisting bilayer structures, which are distinguished by the transbilayer phosphate-phosphate distance of coupled PSM molecules, are observed by diffraction at 37 °C. The MD simulations reveal that PSM molecules in the thicker membrane are characterized by more ordered, more extended, and less interdigitated hydrocarbon tails compared to those in the thinner membrane. Intermolecular hydrogen bonds further distinguish the two bilayer structures, and we observe the disruption of a sphingomyelin intermolecular hydrogen bond network induced by the proximity of cholesterol. Through an unsupervised clustering of interatomic distances, we show for the first time that the asymmetry of phospholipids is important in driving their interactions with cholesterol. We identify four distinct modes of interaction, two of which lead to the dehydration of cholesterol. These two modes of interaction provide the first description of precise physical mechanisms underlying the umbrella model, which itself explains how phospholipids may shield cholesterol from water. The most dehydrating mode of interaction is particular to the N-acylated fatty acid moiety of PSM and thus may explain the long-held observation that cholesterol preferentially mixes with sphingomyelins over glycerophospholipids.

The Role of Hyaluronic Acid in Cartilage Boundary Lubrication

Hydration lubrication has emerged as a new paradigm for lubrication in aqueous and biological media, accounting especially for the extremely low friction (friction coefficients down to 0.001) of articular cartilage lubrication in joints. Among the ensemble of molecules acting in the joint, phosphatidylcholine (PC) lipids have been proposed as the key molecules forming, in a complex with other molecules including hyaluronic acid (HA), a robust layer on the outer surface of the cartilage. HA, ubiquitous in synovial joints, is not in itself a good boundary lubricant, but binds the PC lipids at the cartilage surface; these, in turn, massively reduce the friction via hydration lubrication at their exposed, highly hydrated phosphocholine headgroups. An important unresolved issue in this scenario is why the free HA molecules in the synovial fluid do not suppress the lubricity by adsorbing simultaneously to the opposing lipid layers, i.e., forming an adhesive, dissipative bridge between them, as they slide past each other during joint articulation. To address this question, we directly examined the friction between two hydrogenated soy PC (HSPC) lipid layers (in the form of liposomes) immersed in HA solution or two palmitoyl-oleoyl PC (POPC) lipid layers across HA-POPC solution using a surface force balance (SFB). The results show, clearly and surprisingly, that HA addition does not affect the outstanding lubrication provided by the PC lipid layers. A possible mechanism indicated by our data that may account for this is that multiple lipid layers form on each cartilage surface, so that the slip plane may move from the midplane between the opposing surfaces, which is bridged by the HA, to an HA-free interface within a multilayer, where hydration lubrication is freely active. Another possibility suggested by our model experiments is that lipids in synovial fluid may complex with HA, thereby inhibiting the HA molecules from adhering to the lipids on the cartilage surfaces.

Interactions of a short hyaluronan chain with a phospholipid membrane

Hyaluronic acid and phospholipids are two components that are present in the synovial fluid, and both are implicated as important facilitators of joint lubrication. In this work we aim to clarify how hyaluronic acid interacts with a phospholipid bilayer through their molecular interactions at the bilayer surface. To this end we performed molecular dynamics simulations of one hyaluronic acid molecule at a phospholipid bilayer in aqueous solution. The simulations were carried out for two aqueous solutions of equal concentrations, containing either NaCl or CaCl₂. We analyzed hydrogen bonds, hydrophobic contacts and cation mediated bridges to clarify how hyaluoronic acid binds to a phospholipid bilayer. The analysis shows that calcium ions promote longer lasting bonds between the species as they create calcium ion bridges between the carboxylate group of hyaluronic acid and the phosphate group of the phospholipid. This type of additional bonding does not significantly influence the total number of contact created, but rather stabilizes the contact. The presented results can facilitate understanding of the role of hyaluronic acid and phospholipid interactions in terms of lubrication of articular cartilage.