Atomistic Molecular Dynamics Simulations Research Articles

Structural, topological, and rheological characteristics of entangled short-chain branched polymer melts under shear flow in comparison with the linear analog

We present a detailed analysis of the general influence of short branches on the structural, topological, and rheological behaviors of entangled short-chain branched (SCB) polyethylene (PE) melt systems under shear flow via direct comparison with the corresponding linear analogs using extensive atomistic nonequilibrium molecular dynamics (NEMD) simulations, for a wide range of flow strengths. In comparison with the linear melt, the SCB systems generally exhibit more compact chain structures and larger dynamic resistance, in response to an imposed flow field at all flow strengths. These features essentially arise from (i) the increased chain stiffness due to the torsional restriction of backbone atoms around the branch points and (ii) the fast random Brownian motion of short branches via their very short characteristic relaxation time. We analyzed various structural and rheological properties, such as anisotropic chain dimension and orientation and their detailed distributions, topological characteristics of the entanglement network, material functions, chain rotation dynamics, and flow birefringence. Distinctive physical characteristics of the entangled SCB systems exposed by these individual properties can be consistently understood based on the fundamental structural and dynamical roles of short branches. These findings are considered informative in our systematic understanding and prediction for the general rheological behaviors of long entangled SCB polymer systems under flow, and in tuning the material properties of SCB polymers in practical applications.

Nanopore Translocation Reveals Electrophoretic Force on Noncanonical RNA:DNA Double Helix.

Electrophoretic transport plays a pivotal role in advancing sensing technologies. So far, systematic studies have focused on the translocation of canonical B-form or A-form nucleic acids, while direct RNA analysis is emerging as the new frontier for nanopore sensing and sequencing. Here, we compare the less-explored dynamics of noncanonical RNA:DNA hybrids in electrophoretic transport to the well-researched transport of B-form DNA. Using DNA/RNA nanotechnology and solid-state nanopores, the translocation of RNA:DNA (RD) and DNA:DNA (DD) duplexes was examined. Notably, RD duplexes were found to translocate through nanopores faster than DD duplexes, despite containing the same number of base pairs. Our experiments reveal that RD duplexes present a noncanonical helix, with distinct transport properties from B-form DD molecules. We find that RD and DD molecules, with the same contour length, move with comparable velocity through nanopores. We examined the physical characteristics of both duplex forms using atomic force microscopy, atomistic molecular dynamics simulations, agarose gel electrophoresis, and dynamic light scattering measurements. With the help of coarse-grained and molecular dynamics simulations, we find the effective force per unit length applied by the electric field to a fragment of RD or DD duplex in nanopores with various geometries or shapes to be approximately the same. Our results shed light on the significance of helical form in nucleic acid translocation, with implications for RNA sensing, sequencing, and the molecular understanding of electrophoretic transport.

Atomistic mechanism of coupling between cytosolic sensor domain and selectivity filter in TREK K2P channels

The two-pore domain potassium (K2P) channels TREK-1 and TREK-2 link neuronal excitability to a variety of stimuli including mechanical force, lipids, temperature and phosphorylation. This regulation involves the C-terminus as a polymodal stimulus sensor and the selectivity filter (SF) as channel gate. Using crystallographic up- and down-state structures of TREK-2 as a template for full atomistic molecular dynamics (MD) simulations, we reveal that the SF in down-state undergoes inactivation via conformational changes, while the up-state structure maintains a stable and conductive SF. This suggests an atomistic mechanism for the low channel activity previously assigned to the down state, but not evident from the crystal structure. Furthermore, experimentally by using (de-)phosphorylation mimics and chemically attaching lipid tethers to the proximal C-terminus (pCt), we confirm the hypothesis that moving the pCt towards the membrane induces the up-state. Based on MD simulations, we propose two gating pathways by which movement of the pCt controls the stability (i.e., conductivity) of the filter gate. Together, these findings provide atomistic insights into the SF gating mechanism and the physiological regulation of TREK channels by phosphorylation.

Nanocellulose-Bovine Serum Albumin Interactions in an Aqueous Medium: Investigations Using In Situ Nanocolloidal Probe Microscopy and Reactive Molecular Dynamics Simulations.

As a versatile nanomaterial derived from renewable sources, nanocellulose has attracted considerable attention for its potential applications in various sectors, especially those focused on water treatment and remediation. Here, we have combined atomic force microscopy (AFM) and reactive molecular dynamics (RMD) simulations to characterize the interactions between cellulose nanofibers modified with carboxylate or phosphate groups and the protein foulant model bovine serum albumin (BSA) at pH 3.92, which is close to the isoelectric point of BSA. Colloidal probes were prepared by modification of the AFM probes with the nanofibers, and the nanofiber coating on the AFM tip was for the first time confirmed through fluorescence labeling and confocal optical sectioning. We have found that the wet-state normalized adhesion force is approximately 17.87 ± 8.58 pN/nm for the carboxylated cellulose nanofibers (TOCNF) and about 11.70 ± 2.97 pN/nm for the phosphorylated ones (PCNF) at the studied pH. Moreover, the adsorbed protein partially unfolded at the cellulose interface due to the secondary structure's loss of intramolecular hydrogen bonds. We demonstrate that nanocellulose colloidal probes can be used as a sensitive tool to reveal interactions with BSA at nano and molecular scales and under in situ conditions. RMD simulations helped to gain a molecular- and atomistic-level understanding of the differences between these findings. In the case of PCNF, partially solvated metal ions, preferentially bound to the phosphates, reduced the direct protein-cellulose connections. This understanding can lead to significant advancements in the development of cellulose-based antifouling surfaces and provide crucial insights for expanding the pH range of use and suggesting appropriate recalibrations.

Promoting Cross-Link Reaction of Polymers by the Matrix-Filler Interface Effect: Role of Coupling Agents and Intermediate Linkers.

The matrix-filler interface effect plays an important role in determining the structural stability and mechanical properties of polymer composite materials, which remain ambiguous and need to be studied. The network-forming dynamics of poly(3,3-bis (azidomethyl) oxetane-tetrahydrofuran) (PBT) at the ammonium perchlorate (AP) surface was studied by using atomistic molecular dynamics simulation, considering the additives of curing agent toluene diisocyanate (TDI), cross-linker trimethylolpropane (TMP), and coupling agent triethanolamine (TEA). The presence of the AP surface promotes chain cross-link reaction, which is attributed to the increased production of intermediate linkers formed by TDI, TMP, and TEA. The intermediate linker has three reactive sites that can react with PBT main chains to form a cross-linked structure. Owing to the strong interaction with the AP surface, the coupling agent TEA plays a dominant role in forming the intermediate linker. At the early stage of network forming (reaction ratio r < 30%), the AP surface adsorbs TEA, which leads to a maximum contact density to PBT. As r increases to 60%, the density of intermediate linkers near the AP surface reaches a maximum value. Consequently, the chain cross-link reactions between the intermediate linker and PBT main chains are enhanced as r > 60%. This work explains the micromechanism of the promotion of chain cross-link reaction by the interface effect and provides important insights on designing polymer materials with high mechanical properties.

Cryo-EM and Molecular Dynamics Simulations Reveal Hidden Conformational Dynamics Controlling Ammonia Transport in Human Asparagine Synthetase.

How motions in enzymes might be linked to catalytic function is of considerable general interest. Recent advances in X-ray crystallography and cryogenic electron microscopy offer the promise of elucidating functionally relevant motions in proteins that are not easily amenable to study by other biophysical methods. Here we use 3D variability analysis (3DVA) on cryo-EM maps for wild type (WT) human asparagine synthetase (ASNS) and the R142I ASNS variant to identify conformational changes in the Arg-142 side chain, which mediates the formation of a catalytically relevant intramolecular tunnel. Our 3DVA results for WT ASNS are consistent with independent molecular dynamics (MD) simulations on a model generated from the X-ray structure of human ASNS. Moreover, MD simulations of computational models for the ASNS/β-aspartyl-AMP/MgPPi and R142I/β-aspartyl-AMP/MgPPi ternary complexes, suggest that the structural integrity of the tunnel is impaired in the R142I variant when β-aspartyl-AMP is present in the synthetase active site. The kinetic properties of the R142I ASNS variant support the proposed function of Arg-142. These studies illustrate the power of cryo-EM to identify localized motions and dissect the conformational landscape of large proteins. When combined with MD simulations, 3DVA is a powerful approach to understanding how conformational dynamics might regulate function in multi-domain enzymes possessing multiple active sites.

Magnesium Ions Mitigate Metastable States in the Regulatory Landscape of mRNA Elements.

Residing in the 5' untranslated region of the mRNA, the 2'-deoxyguanosine (2'-dG) riboswitch mRNA element adopts an alternative structure with the binding of the 2'-dG molecule, which terminates transcription. In general, RNA conformations are strongly affected by positively charged metal ions (especially Mg2+). We have quantitatively explored the combined effect of ligand (2'-dG) and Mg2+ binding on the energy landscape of the aptamer domain of the 2'-dG riboswitch with both explicit solvent all atom molecular dynamics simulations and SHAPE biochemical probing experiments. We show that both ligand and Mg2+ are required for the stabilization of the aptamer domain; however, the two factors function in different ways. While the addition of Mg2+ remodels the energy landscape and reduces its frustration by the formation of additional contacts, the binding of 2'-dG eliminates the metastable states by building a compact core for the aptamer domain. In particular, Mg2+ ions and ligand binding are required to stabilize the most unstable helix P1 (which needs to unfold to activate the transcription platform), and the riboswitch core formed by the backbone of the P2 and P3 helices. They also facilitate a more compact structure in the three-way junction region.

Interaction of Statherin-Derived Peptide with the Surface of Hydroxyapatite: Perspectives Based on Molecular Dynamics Simulations

Introduction: Statherin-derived peptide (StatpSpS) has shown promise against erosive tooth wear. To elucidate its interaction with the hydroxyapatite (HAP) surface, the mechanism related to adsorption of this peptide with HAP was investigated through nanosecond-long all-atom molecular dynamics simulations. Methods: StatpSpS was positioned parallel to the HAP surface in 2 orientations: 1 – neutral and negative residues facing the surface and 2 – positive residues facing the surface. A system containing StatpSpS without HAP was also simulated as control. In the case of systems with HAP, both partially restrained surface and unrestrained surface were constructed. Structural analysis, interaction pattern, and binding-free energy were calculated. Results: In the peptide system without the HAP, there were some conformational changes during the simulation. In the presence of the surface, only moderate changes were observed. Many residues exhibited short and stable distances to the surface, indicating strong interaction. Specially, the residues ASP1 and SER2 have an important role to anchor the peptide to the surface, with positively charged residues, mainly arginine, playing a major role in the further stabilization of the peptide in an extended conformation, with close contacts to the HAP surface. Conclusion: The interaction between StatpSpS and HAP is strong, and the negative charged residues are important to the anchoring of the peptide in the surface, but after the initial placement the peptide rearranges itself to maximize the interactions between positive charged residues.

Molecular dynamics simulations of lipid-protein interactions in SLC4 proteins

The SLC4 family of secondary bicarbonate transporters is responsible for the transport of HCO3−, CO32−, Cl−, Na+, K+, NH3, and H+, which are necessary for regulation of pH and ion homeostasis. They are widely expressed in numerous tissues throughout the body and function in different cell types with different membrane properties. Potential lipid roles in SLC4 function have been reported in experimental studies, focusing mostly on two members of the family: AE1 (Cl−/HCO3− exchanger) and NBCe1 (Na+-CO32−cotransporter). Previous computational studies of the outward-facing state of AE1 with model lipid membranes revealed enhanced protein-lipid interactions between cholesterol (CHOL) and phosphatidylinositol bisphosphate (PIP2). However, the protein-lipid interactions in other members of the family and other conformation states are still poorly understood and this precludes the detailed studies of a potential regulatory role for lipids in the SLC4 family. In this work, we performed coarse-grained and atomistic molecular dynamics simulations on three members of the SLC4 family with different transport modes: AE1, NBCe1, and NDCBE (an Na+-CO32−/Cl− exchanger), in model HEK293 membranes consisting of CHOL, PIP2, phosphatidylcholine, phosphatidylethanolamine, phosphatidylserine, and sphingomyelin. The recently resolved inward-facing state of AE1 was also included in the simulations. Lipid-protein contact analysis of the simulated trajectories was performed with the ProLint server, which provides a multitude of visualization tools for illustration of areas of enhanced lipid-protein contact and identification of putative lipid binding sites within the protein matrix. We observed enrichment of CHOL and PIP2 around all proteins with subtle differences in their distribution depending on the protein type and conformation state. Putative binding sites were identified for CHOL, PIP2, phosphatidylcholine, and sphingomyelin in the three studied proteins, and their potential roles in the SLC4 transport function, conformational transition, and protein dimerization are discussed.

Molecular Insights into the Effects of F16L and F19L Substitutions on the Conformation and Aggregation Dynamics of Human Calcitonin.

Human calcitonin (hCT) regulates calcium-phosphorus metabolism, but its amyloid aggregation disrupts physiological activity, increases thyroid carcinoma risk, and hampers its clinical use for bone-related diseases like osteoporosis and Paget's disease. Improving hCT with targeted modifications to mitigate amyloid formation while maintaining its function holds promise as a strategy. Understanding how each residue in hCT's amyloidogenic core affects its structure and aggregation dynamics is crucial for designing effective analogues. Mutants F16L-hCT and F19L-hCT, where Phe residues in the core are replaced with Leu as in nonamyloidogenic salmon calcitonin, showed different aggregation kinetics. However, the molecular effects of these substitutions in hCT are still unclear. Here, we systematically investigated the folding and self-assembly conformational dynamics of hCT, F16L-hCT, and F19L-hCT through multiple long-time scale independent atomistic discrete molecular dynamics (DMD) simulations. Our results indicated that the hCT monomer primarily assumed unstructured conformations with dynamic helices around residues 4-12 and 14-21. During self-assembly, the amyloidogenic core of hCT14-21 converted from dynamic helices to β-sheets. However, substituting F16L did not induce significant conformational changes, as F16L-hCT exhibited characteristics similar to those of wild-type hCT in both monomeric and oligomeric states. In contrast, F19L-hCT exhibited substantially more helices and fewer β-sheets than did hCT, irrespective of their monomers or oligomers. The substitution of F19L significantly enhanced the stability of the helical conformation for hCT14-21, thereby suppressing the helix-to-β-sheet conformational conversion. Overall, our findings elucidate the molecular mechanisms underlying hCT aggregation and the effects of F16L and F19L substitutions on the conformational dynamics of hCT, highlighting the critical role of F19 as an important target in the design of amyloid-resistant hCT analogs for future clinical applications.

When Tautomers Matter: UV-Vis Absorption Spectra of Hypoxanthine in Aqueous Solution from Fully Atomistic Simulations.

The UV-Vis spectrum of the solvated purine derivative Hypoxanthine (HYX) is investigated using the Quantum Mechanics/Fluctuating Charges (QM/FQ) multiscale approach combined with a sampling of configurations through atomistic Molecular Dynamics (MD) simulations. Keto 1H7H and 1H9H tautomeric forms of HYX are the most stable in aqueous solution and form different stable complexes with the surrounding water molecules, ultimately affecting the electronic absorption spectra. The final simulated spectrum resulting from the combination of the individual spectra of tautomers agrees very well with most of the characteristics in the measured spectrum. The importance of considering the effect of the solute tautomers and, in parallel, the contribution of the different solvent arrangements around the solute when modeling spectral properties, is highlighted. In addition, the high quality of the computed spectra leads to suggesting an alternative way for acquiring tautomeric populations from combined computational/experimental spectra.

Enhanced pyroptosis induction with pore-forming gene delivery for osteosarcoma microenvironment reshaping

Osteosarcoma is the most common malignant bone tumor without efficient management for improving 5-year event-free survival. Immunotherapy is also limited due to its highly immunosuppressive tumor microenvironment (TME). Pore-forming gasdermins (GSDMs)-mediated pyroptosis has gained increasing concern in reshaping TME, however, the expressions and relationships of GSDMs with osteosarcoma remain unclear. Herein, gasdermin E (GSDME) expression is found to be positively correlated with the prognosis and immune infiltration of osteosarcoma patients, and low GSDME expression was observed. A vector termed as LPAD contains abundant hydroxyl groups for hydrating layer formation was then prepared to deliver the GSDME gene to upregulate protein expression in osteosarcoma for efficient TME reshaping via enhanced pyroptosis induction. Atomistic molecular dynamics simulations analysis proved that the hydroxyl groups increased LPAD hydration abilities by enhancing coulombic interaction. The upregulated GSDME expression together with cleaved caspase-3 provided impressive pyroptosis induction. The pyroptosis further initiated proinflammatory cytokines release, increased immune cell infiltration, activated adaptive immune responses and create a favorable immunogenic hot TME. The study not only confirms the role of GSDME in the immune infiltration and prognosis of osteosarcoma, but also provides a promising strategy for the inhibition of osteosarcoma by pore-forming GSDME gene delivery induced enhanced pyroptosis to reshape the TME of osteosarcoma.

Atomistic molecular dynamics simulation and COSMO-SAC approach for enhanced 1,3-propanediol extraction with imidazolium-based ionic liquids.

1,3-Propanediol (1,3-PDO) is a key chemical in various industries, including pharmaceuticals and material sciences, and is projected to see significant market growth. However, the current challenges in its downstream processing, particularly in terms of cost and efficiency, highlight the need for innovative solutions. Our study delves into using ionic liquids (ILs) as a potential alternative, aiming to address these critical separation challenges more sustainably and efficiently. In this study, we utilized molecular dynamics (MD) simulations and the COSMO-SAC to examine 1,3-propanediol (1,3-PDO) extraction using four imidazolium-based ionic liquids with 1-butyl-3-methylimidazolium [Bmim] cation and with different anions bis(pentafluoroethanesulfonyl)imide [NPF2]-, bis(trifluoromethylsulfonyl)imide [NTF2]-, thiocyanate [SCN]-, and trifluoromethanesulfonate [TFO]-. Molecular dynamics simulations, incorporating analysis of radial distribution functions (RDF) and spatial distribution functions (SDF), revealed that [Bmim][SCN] and [Bmim][TFO] exhibit enhanced interactions with 1,3-PDO. Notably, [Bmim][SCN] formed the most hydrogen bonds, averaging 1.639 per molecule, due to its coordinating [SCN]- anion. This was in contrast to the fewer hydrogen bonds formed by non-coordinating anions in [Bmim][NPF2] and [Bmim][NTF2]. In ternary systems, [Bmim][SCN] and [Bmim][TFO] demonstrated superior selectivity for 1,3-PDO extraction compared to the other ionic liquids, with selectivity values around 29. These findings, supported by COSMO-SAC predictive modeling, highlight the potential of [Bmim][SCN] as a promising candidate for 1,3-PDO extraction, emphasizing the importance of anion selection in optimizing ionic liquid properties for this application. In our study, we employed MD simulations, incorporating the OPLS-AA force field, and COSMO-SAC to investigate the extraction of 1,3-PDO using imidazolium-based ionic liquids: [Bmim][NTF2], [Bmim][NPF2], [Bmim][SCN], and [Bmim][TFO]. The MD simulations were conducted using LAMMPS software, focusing on elucidating the RDF, SDF, and hydrogen bonding. Analysis of the distribution coefficient (β) and selectivity (S) for the ternary mixture was also conducted. These aspects of the simulation were analyzed using TRAVIS and VMD software. Additionally, the COSMO-SAC model was employed to determine the activity coefficients of 1,3-PDO in the ionic liquids, with molecular optimization conducted using Gaussian16 and sigma profile calculations performed using COSMO-SAC.

Atomistic insights into the superior performance of C60-decorated graphene supercapacitors

Properties of graphene-based supercapacitors decorated with ▪ fullerenes and three different electrolytes based on ▪ cation and two different anions (▪ and ▪) were investigated by atomistic molecular dynamics simulations. Our simulations reveal that the electrolyte first layer near the electrode, which is made up of the region decorated with fullerenes, is a key factor in explaining the remarkable performance of electrodes coated with ▪ fullerenes. The access of the electrolyte between the interstitial spaces formed between the ▪ molecules provides a network of extended and intense electrostatic interactions that do not exist in a pure graphene electrode. The study demonstrates that ▪ molecules are effective in creating small pores in the graphene sheets, which permit ion mobility while retaining the structural and electrostatic characteristics that are consistent with the experimentally measured values. Notably, this engineered porosity has proven to significantly increase the electrical response at the interface, dramatically improving device performance. These simulation results provide a better understanding of the ▪ porosity of the electrode and the performance of the graphene/▪-based supercapacitor.

A study of hydrophobic domain formation of polymeric drug precipitation inhibitors in aqueous solution

Despite the widespread use of polymers as precipitation inhibitors in supersaturating drug formulations, the current understanding of their mechanisms of action is still incomplete. Specifically, the role of hydrophobic drug interactions with polymers by considering possible supramolecular conformations in aqueous dispersion is an interesting topic. Accordingly, this study investigated the tendency of polymers to create hydrophobic domains, where lipophilic compounds may nest to support drug solubilisation and supersaturation. Fluorescence spectroscopy with the environment-sensitive probe pyrene was compared with atomistic molecular dynamics simulations of the model drug fenofibrate (FENO). Subsequently, kinetic drug supersaturation and thermodynamic solubility experiments were conducted. As a result, the different polymers showed hydrophobic domain formation to a varying degree and the molecular simulations supported interpretation of fluorescence spectroscopy data. Molecular insights were gained into the conformational structure of how the polymers interacted with FENO in solution phase, which apart from nucleation and crystal growth effects, determined drug concentrations in solution. Notable was that even at the lowest polymer concentration of 0.01 %, w/v, there were polymer-specific solubilisation effects of FENO observed and the resulting reduction in apparent drug supersaturation provided relevant knowledge both from a mechanistic and practical perspective.

Simulations predict preferred Mg2+ coordination in a nonenzymatic primer-extension reaction center

The mechanism by which genetic information was copied prior to the evolution of ribozymes is of great interest because of its importance to the origin of life. The most effective known process for the nonenzymatic copying of an RNA template is primer extension by a two-step pathway in which 2-aminoimidazole-activated nucleotides first react with each other to form an imidazolium-bridged intermediate that subsequently reacts with the primer. Reaction kinetics, structure-activity relationships, and X-ray crystallography have provided insight into the overall reaction mechanism, but many puzzles remain. In particular, high concentrations of Mg2+ are required for efficient primer extension, but the mechanism by which Mg2+ accelerates primer extension remains unknown. By analogy with the mechanism of DNA and RNA polymerases, a role for Mg2+ in facilitating the deprotonation of the primer 3′-hydroxyl is often assumed, but no catalytic metal ion is seen in crystal structures of the primer-extension complex. To explore the potential effects of Mg2+ binding in the reaction center, we performed atomistic molecular dynamics simulations of a series of modeled complexes in which a Mg2+ ion was placed in the reaction center with inner-sphere coordination with different sets of functional groups. Our simulations suggest that coordination of a Mg2+ ion with both O3′ of the terminal primer nucleotide and the pro-Sp nonbridging oxygen of the reactive phosphate of an imidazolium-bridged dinucleotide would help to pre-organize the structure of the primer/template substrate complex to favor the primer-extension reaction. Our results suggest that the catalytic metal ion may play an important role in overcoming electrostatic repulsion between a deprotonated O3′ and the reactive phosphate of the bridged dinucleotide and lead to testable predictions of the mode of Mg2+ binding that is most relevant to catalysis of primer extension.

Molecular simulation on the interaction between trehalose and asymmetric lipid bilayer mimicking the membrane of human red blood cells

Trehalose is widely acknowledged for its ability to stabilize plasma membranes during dehydration. However, the exact mechanism by which trehalose interacts with lipid bilayers remains presently unclear. In this study, we conducted atomistic molecular dynamic simulations on asymmetric model bilayers that mimic the membrane of human red blood cells at various trehalose and water contents. We considered three different hydration levels mimicking the full hydration to desiccation scenarios. Results indicate that the asymmetric distribution of lipids did not significantly influence the computed structural characteristics at full and low hydration. At dehydration, however, the order parameter obtained from the symmetric bilayer is significantly higher compared to those obtained from asymmetric ones. Analysis of hydrogen bonds revealed that the protective ability of trehalose is well described by the water replacement hypothesis at full and low hydration, while at dehydration other interaction mechanisms associated with trehalose exclusion from the bilayer may involve. In addition, we found that trehalose exclusion is not attributed to sugar saturation but rather to the reduction in hydration levels. It can be concluded that the protective effect of trehalose is not only related to the hydration level of the bilayer, but also closely tied to the asymmetric distribution of lipids within each leaflet.

Sodium dodecyl sulfate micelles: Accurate analysis of small-angle X-ray scattering data through form factor and atomistic molecular dynamics modelling

The structure of micelles of sodium dodecyl sulfate (SDS) is probed via analysis of small-angle x-ray scattering (SAXS) data with the aim to fit the data over an extended wavenumber q range. This provides detailed information on the micelle shape, which can be described as a polydisperse triaxial ellipsoidal core-shell structure according to model form factor fitting. This model was necessary to fit the data over a wide q range, which is not accurately represented by simpler models such as biaxial ellipsoidal core-shell structures. Data for SDS (at fixed concentration) in a NaCl concentration series revealed increasing structure factor effects with decreasing salt concentration. This reflects decreased charge screening on the headgroups. The structure factor could be modelled using a simple hard sphere structure factor. The analysis of form factor was complemented by atomistic molecular dynamics (MD) simulations, starting from an unbiased initial configuration of a defined number of molecules in a box. The MD configurations were used to calculate the form factor using the software CRYSOL (for small-angle scattering analysis of solution scattering, traditionally for proteins, here for micelles) and accounting for the boundary layer hydration effects. Good agreement with experimental data was found for systems with association numbers close to p = 60. This association number is consistent with that obtained from analysis of the form factor (in the case where structure factor effects could be neglected) and from model-free analysis of the forward scattering intensity. It is also in agreement with prior literature and our findings in regard to form factor parameters are also compared to previous reports.

Nanoactuator for Neuronal Optoporation.

Light-driven modulation of neuronal activity at high spatial-temporal resolution is becoming of high interest in neuroscience. In addition to optogenetics, nongenetic membrane-targeted nanomachines that alter the electrical state of the neuronal membranes are in demand. Here, we engineered and characterized a photoswitchable conjugated compound (BV-1) that spontaneously partitions into the neuronal membrane and undergoes a charge transfer upon light stimulation. The activity of primary neurons is not affected in the dark, whereas millisecond light pulses of cyan light induce a progressive decrease in membrane resistance and an increase in inward current matched to a progressive depolarization and action potential firing. We found that illumination of BV-1 induces oxidation of membrane phospholipids, which is necessary for the electrophysiological effects and is associated with decreased membrane tension and increased membrane fluidity. Time-resolved atomic force microscopy and molecular dynamics simulations performed on planar lipid bilayers revealed that the underlying mechanism is a light-driven formation of pore-like structures across the plasma membrane. Such a phenomenon decreases membrane resistance and increases permeability to monovalent cations, namely, Na+, mimicking the effects of antifungal polyenes. The same effect on membrane resistance was also observed in nonexcitable cells. When sustained light stimulations are applied, neuronal swelling and death occur. The light-controlled pore-forming properties of BV-1 allow performing "on-demand" light-induced membrane poration to rapidly shift from cell-attached to perforated whole-cell patch-clamp configuration. Administration of BV-1 to ex vivo retinal explants or in vivo primary visual cortex elicited neuronal firing in response to short trains of light stimuli, followed by activity silencing upon prolonged light stimulations. BV-1 represents a versatile molecular nanomachine whose properties can be exploited to induce either photostimulation or space-specific cell death, depending on the pattern and duration of light stimulation.

Development, Evaluation, and Molecular Dynamics Study of Ampicillin-Loaded Chitosan-Hyaluronic Acid Films as a Drug Delivery System.

Periodontitis is an inflammatory periodontal disease defined by the progressive loss of tissues surrounding the tooth. Ampicillin is an antibiotic for managing and treating specific bacterial infections, including periodontitis. Periodontal pockets occur due to periodontal disease progression and act as a natural reservoir that is easily reachable for the insertion of a delivery system, and the amount of drug to be released has a major role in the efficiency of treatment of the disease. Polyelectrolyte complexes (PECs), particularly those based on chitosan and hyaluronic acid combinations, offer a promising avenue to overcome the challenges associated with drug delivery. These complexes are both biodegradable and biocompatible, making them an optimal choice for enabling targeted drug delivery. This study centers on developing and assessing the structure and dynamic attributes of a drug-PEC system encompassing ampicillin and chitosan-hyaluronic acid components, which represents a targeted drug delivery system to better alleviate the periodontitis. To achieve this goal, we conducted experiments including weight and drug content uniformity, swelling ındex, drug release %, FT-IR and SEM analyses, and atomistic molecular dynamics simulations on the drug PECs loaded with ampicillin with varying amounts of hyaluronic acid. All simulations and the experimental analysis suggested that increased HA amount resulted in an increase in drug release % and swelling index. The simulation outcomes provide insights into the nature of the drug and PEC interactions alongside transport properties such as drug diffusion coefficients. These coefficients offer valuable insights into the molecular behavior of ampicillin-PEC drug delivery systems, particularly in the context of their application in periodontitis treatment.