All-atom Molecular Dynamics Simulations Research Articles

Enriching Anticancer Drug Pipeline with Potential Inhibitors of Cyclin-Dependent Kinase-8 Identified from Natural Products.

Cyclin-dependent kinase 8 (CDK8) is highly expressed in various cancers and common complex human diseases, and an important therapeutic target for drug discovery and development. The CDK8 inhibitors are actively sought after, especially among natural products. We performed a virtual screening using the ZINC library comprising approximately 90,000 natural compounds. We applied Lipinski's rule of five, absorption, distribution, metabolism, excretion, and toxicity properties, and pan-assay interference compounds filter to eliminate promiscuous binders. Subsequently, the filtered compounds underwent molecular docking to predict their binding affinity and interactions with the CDK8 protein. Interaction analysis were carried out to elucidate the interaction mechanism of the screened hits with binding pockets of the CDK8. The ZINC02152165, ZINC04236005, and ZINC02134595 were selected with appreciable specificity and affinity with CDK8. An all-atom molecular dynamic (MD) simulation followed by essential dynamics was performed for 200 ns. Taken together, the results suggest that ZINC02152165, ZINC04236005, and ZINC02134595 can be harnessed as potential leads in therapeutic development. Moreover, the binding of the molecules brings change in protein conformation in a way that blocks the ATP-binding site of the protein, obstructing its kinase activity. These new findings from natural products offer insights into the molecular mechanisms underlying CDK8 inhibition. CDK8 was previously associated with behavioral and neurological diseases such as autism spectrum disorder, and cancers, for example, colorectal, prostate, breast, and acute myeloid leukemia. Hence, we call for further research and experimental validation, and with an eye to inform future clinical drug discovery and development in these therapeutic fields.

RNA binding tunes the conformational plasticity and intradomain stability of TDP-43 tandem RNA recognition motifs

TAR DNA binding protein 43 (TDP-43) is a nuclear RNA/DNA-binding protein with pivotal roles in RNA-related processes such as splicing, transcription, transport, and stability. The high binding affinity and specificity of TDP-43 toward its cognate RNA sequences (GU-rich) is mediated by highly conserved residues in its tandem RNA recognition motif (RRM) domains (aa: 104–263). Importantly, the loss of RNA binding to the tandem RRMs caused by physiological stressors and chemical modifications promotes cytoplasmic mislocalization and pathological aggregation of TDP-43. Despite the substantial implications of RNA binding in TDP-43 function and pathology, its precise effects on the intradomain stability, and conformational dynamics of the tandem RRMs is not properly understood. Here, we employed all-atom molecular dynamics (MD) simulations to assess the effect of RNA binding on the conformational landscape and intradomain stability of TDP-43 tandem RRMs. RNA limits the overall conformational space of the tandem RRMs and promotes intradomain stability through a combination of specific base stacking interactions and transient electrostatic interactions. In contrast, tandem RRMs exhibit a high intrinsic conformational plasticity in the absence of RNA, which, surprisingly, is accompanied by a tendency of RRM1 to adopt partially unfolded conformations. Overall, our simulations reveal how RNA binding dynamically tunes the structural and conformational landscape of TDP-43 tandem RRMs, contributing to physiological function and mitigating pathological aggregation.

An integrative characterization of proline cis and trans conformers in a disordered peptide

Intrinsically disordered proteins (IDPs) often contain proline residues that undergo cis/trans isomerization. While molecular dynamics (MD) simulations have the potential to fully characterize the proline cis and trans subensembles, they are limited by the slow timescales of isomerization and force field inaccuracies. NMR spectroscopy can report on ensemble-averaged observables for both the cis-proline and trans-proline states, but a full atomistic characterization of these conformers is challenging. Given the importance of proline cis/trans isomerization for influencing the conformational sampling of disordered proteins, we employed a combination of all-atom MD simulations with enhanced sampling (metadynamics), NMR, and small-angle x-ray scattering (SAXS) to characterize the two subensembles of the ORF6 C-terminal region (ORF6CTR) from SARS-CoV-2 corresponding to the proline-57 (P57) cis and trans states. We performed MD simulations in three distinct force fields: AMBER03ws, AMBER99SB-disp, and CHARMM36m, which are all optimized for disordered proteins. Each simulation was run for an accumulated time of 180–220 μs until convergence was reached, as assessed by blocking analysis. A good agreement between the cis-P57 populations predicted from metadynamic simulations in AMBER03ws was observed with populations obtained from experimental NMR data. Moreover, we observed good agreement between the radius of gyration predicted from the metadynamic simulations in AMBER03ws and that measured using SAXS. Our findings suggest that both the cis-P57 and trans-P57 conformations of ORF6CTR are extremely dynamic and that interdisciplinary approaches combining both multiscale computations and experiments offer avenues to explore highly dynamic states that cannot be reliably characterized by either approach in isolation.

Molecular Insights into Single Chain Lipid Modulation of Acid-Sensing Ion Channel 3.

Polyunsaturated fatty acids (PUFAs) and their analogs play a significant role in modulating the activity of diverse ion channels, and recent studies show that these lipids potentiate acid-sensing ion channels (ASICs), leading to increased activity. The potentiation of the channel stems from multiple gating changes, but the exact mechanism of these effects remains uncertain. We posit a mechanistic explanation for one of these changes in channel function, the increase in the maximal current, by applying a combination of electrophysiology and all-atom molecular dynamics simulations on the open-state hASIC3. Microsecond-scale simulations were performed on open-state hASIC3 in the absence and presence of a PUFA, docosahexaenoic acid (DHA), and a PUFA analog, N-arachidonyl glycine (AG). Intriguingly, our simulations in the absence of PUFA or PUFA analogs reveal that a tail from the membrane phospholipid POPC inserts itself into the pore of the channel through lateral fenestrations on the sides of the transmembrane segments, obstructing ion permeation through the channel. The binding of either DHA or AG prevented POPC from accessing the pore in our simulations, relieving the block of ionic conduction by phospholipids. Finally, we use the single-channel recording to show that DHA increases the amplitude of the single-channel currents in ASIC3, which is consistent with our hypothesis that PUFAs relieve the pore block of the channel induced by POPCs. Together, these findings offer a potential mechanistic explanation of how PUFAs modulate ASIC maximal current, revealing a novel mechanism of action for PUFA-induced modulation of ion channels.

Role of the structural order of the hydration layer in regulating the heterogeneous ice nucleation efficiency

Understanding the complex microscopic mechanisms of heterogeneous ice nucleation is crucial across diverse fields from cloud science to microbiology. In this work, we examined the ability of solid surfaces to promote ice nucleation as a function of the structural order of interfacial water, including the first and second hydration layers. Comprehensive all-atom molecular dynamics simulations using the TIP4P/ice water model were conducted on a water film atop a modeled surface inspired by the (0001) surface of AgI crystal, with varying surface charges. These simulations revealed non-monotonic nucleation rates. The first hydration layer formed a hexagonal hydrogen-bond network resembling ice structures, exhibiting low mobility essential for inducing ice nucleation. Notably, the second hydration layer, acting as a key-point bridge, displayed varying degrees of orientational preference facilitated by inter-layer hydrogen bonds. The calculation of fluctuation parameters indicated that the greater structural order of the second layer significantly enhances the surface’s ice-forming capabilities. Unlike previous studies that focused primarily on the role of the first layer of interfacial water, our findings demonstrate that the second layer provides a more accurate indicator of ice nucleation capabilities.

Mechanical response of double-stranded DNA: Bend, twist, and overwind.

We employed all-atom molecular dynamics simulations to explore the mechanical response of bending, twisting, and overwinding for double-stranded DNA (dsDNA). We analyzed the bending and twisting deformations, as well as their stiffnesses, using the tilt, roll, and twist modes under stretching force. Findings indicate that the roll and twist angles vary linearly with the stretching force but show opposite trends. The tilt, roll, and twist elastic moduli are considered constants, while the coupling between roll and twist modes slightly decreases under stretching force. The effect of the stretching force on the roll and twist modes, including both their deformations and elasticities, exhibits sequence-dependence, with symmetry around the base pair step. Furthermore, we examined the overwinding path and mechanism of dsDNA from the perspective of the stiffness matrix, based on the tilt, roll, and twist modes. The correlations among tilt, roll, and twist angles imply an alternative overwinding pathway via twist-roll coupling when dsDNA is stretched, wherein entropic contribution prevails.

Rotation and Self-Assembly Driving NLRP3 Activation.

As a critical sensor protein, NLRP3 detects cellular perturbation caused by diverse exogenous and endogenous stimuli. NLRP3 activation requires domain rotation within the NEK7-bound NLRP3 monomer and assembly. However, a detailed molecular mechanism for NLRP3 assembly and activation remains elusive, particularly in terms of dynamics and energetics. In this work, all-atom molecular dynamics (MD) simulations are executed to describe large-amplitude closed-to-open conformational transitions along the rotational pathway. From the MD trajectories, the computed potential of mean force (PMF) shows that NLRP3 activation through monomeric domain rotation is an uphill process, during which the active conformation of the NLRP3-NEK7 monomer cannot be stabilized. Further binding free-energy calculations for two neighboring NLRP3-NEK7 subunits in a disc assembly with the C10 symmetry reveal that the protein self-assembly starts approximately at the 86.5° position on the rotary pathway, along which the NLRP3 activation becomes a downhill process to the active state at 90.5°. The active NLRP3-NEK7 monomeric conformation in the disc assembly is stabilized because of the interactions between the neighboring subunits, involving mainly FISNA loop 1 in one subunit and a "crocodile-clip" structure formed by the NBD helix-loop-strand motif (residues 351-373) and the WHD β-hairpin loop (residues 501-521) in the other. Our simulations also demonstrate that NEK7 plays an important role in the NLRP3 cage dissociation in the centrosome, which is consistent with biological experiments. The computational results provide kinetic, energetic, and structural insights into the molecular mechanisms of the activation of NLRP3 and the NEK7-driven dissociation of inactive NLRP3 cages. The activation mechanism of NLRP3 proposed in this work is significantly different from those of previous structural studies.

NP-Collabs: Investigating Counterion-Mediated Bridges in the Multiply Phosphorylated Tau-R2 Repeat.

Tau is an intrinsically disordered (IDP) microtubule-associated protein (MAP) that plays a key part in microtubule assembly and organization. The function of tau can be regulated by multiple phosphorylation sites. These post-translational modifications are known to decrease the binding affinity of tau for microtubules, and abnormal tau phosphorylation patterns are involved in Alzheimer's disease. Using all-atom molecular dynamics simulations, we compared the conformational landscapes explored by the tau R2 repeat domain (which comprises a strong tubulin binding site) in its native state and with multiple phosphorylations on the S285, S289, and S293 residues, with four different standard force field (FF)/water model combinations. We find that the different parameters used for the phosphate groups (which can be more or less flexible) in these FFs and the specific interactions between bulk cations and water lead to the formation of a specific type of counterion bridge, termed nP-collab (for nphosphate collaboration, with n being an integer), where counterions form stable structures binding with two or three phosphate groups simultaneously. The resulting effect of nP-collabs on the tau-R2 conformational space differs when using sodium or potassium cations and is likely to impact the peptide overall dynamics and how this MAP interacts with tubulins. We also investigated the effect of phosphoresidue spacing and ionic concentration by modeling polyalanine peptides containing two phosphoserines located one-six residues apart. Three new metrics specifically tailored for IDPs (proteic Menger curvature, local curvature, and local flexibility) were introduced, which allow us to fully characterize the impact of nP-collabs on the dynamics of disordered peptides at the residue level.

Ribosome Tunnel Environment Drives the Formation of α-Helix during Cotranslational Folding.

Protein conformations in cells are not solely determined by amino acid sequences; they also depend on cellular environments. For instance, the ribosome tunnel induces its specific α-helix formation during cotranslational folding. Owing to the link between these temporally α-helix and biological functions, the mechanism of α-helix formation inside the ribosome tunnel has been previously explored. Consequently, the conformational restrictions of the tunnel were considered one of the driving forces of α-helix formation. Conversely, the ribosomal tunnel environment, including its chemical properties, appears to influence the α-helix formation. However, a comprehensive analysis of the ribosome tunnel environment's impact on the α-helix formation has not been conducted yet due to challenges in experimentally controlling it. Therefore, as a new computational approach, we proposed a ribosome environment-mimicking model (REMM) based on the radius and components of the experimentally determined ribosome tunnel structures. Using REMM, we assessed the impact of the ribosome tunnel environment on α-helix formation. Herein, we employed carbon nanotubes (CNT) as a reference model alongside REMM because CNT reproduce conformational restrictions rather than the ribosome tunnel environment. Quantitatively, the ability to reproduce the α-helix of nascent peptides in the experimental structure was compared between the CNT and REMM using enhanced all-atom molecular dynamics simulations. Consequently, the REMM more accurately reproduced the α-helix of the nascent peptides than the CNT, highlighting the significance of the ribosome tunnel environment in α-helix formation. Additionally, we analyzed the properties of the peptide inside each model to reveal the mechanism of ribosome tunnel-specific α-helix formation. Consequently, we revealed that the chemical diversities of the tunnel are essential for the formation of backbone-to-backbone hydrogen bonds in the peptides. In conclusion, the ribosome tunnel environment, with the diverse chemical properties, drives its specific α-helix formation. By proposing REMM, we newly provide the technical basis for investigating the protein conformations in various cellular environments.

Fracture properties of green nano fibrous network with random and aligned fiber distribution: A hierarchical molecular dynamics and peridynamics approach

Polylactic acid (PLA) nanofibrous networks have gained substantial interest across various engineering and scientific disciplines, such as tissue engineering, drug delivery, and filtration, due to their unique and multifunctional attributes, including biodegradability, tuneable mechanical properties, and surface functionality. However, predicting their mechanical behaviour remains challenging due to their structural complexity, multiscale features, and variability in material properties.This study presents a hierarchical approach to investigate the fracture phenomena in both aligned and randomly oriented nanofibrous networks by integrating atomistic modelling and non-local continuum mechanics, peridynamics. At the nanoscale, all-atom molecular dynamics simulations are employed to apply tensile loads to freestanding pristine and silver-doped PLA nanofibres, where key mechanical properties such as Young's modulus, Poisson's ratio, and critical energy release rate are determined using innovative approaches. A new method is introduced to seamlessly transfer data from molecular dynamics to peridynamics by ensuring the convergence of the tensile response of a single fiber in both frameworks. This nano to micro coupling technique is then utilised to examine the Young's modulus, fracture toughness of mode I and II, and crack propagation in PLA nanofibrous networks. The proposed framework can also incorporate the effects of surface coating and fiber arrangements on the measured properties. The current research paves the way for the development of stronger and more durable eco-friendly nanofibrous networks with optimised performance.

Effect of ethanol on the elasticities of double-stranded RNA and DNA revealed by magnetic tweezers and simulations.

The elasticities of double-stranded (ds) DNA and RNA, which are critical to their biological functions and applications in materials science, can be significantly modulated by solution conditions such as ions and temperature. However, there is still a lack of a comprehensive understanding of the role of solvents in the elasticities of dsRNA and dsDNA in a comparative way. In this work, we explored the effect of ethanol solvent on the elasticities of dsRNA and dsDNA by magnetic tweezers and all-atom molecular dynamics simulations. We found that the bending persistence lengths and contour lengths of dsRNA and dsDNA decrease monotonically with the increase in ethanol concentration. Furthermore, the addition of ethanol weakens the positive twist-stretch coupling of dsRNA, while promotes the negative twist-stretch coupling of dsDNA. Counter-intuitively, the lower dielectric environment of ethanol causes a significant re-distribution of counterions and enhanced ion neutralization, which overwhelms the enhanced repulsion along dsRNA/dsDNA, ultimately leading to the softening in bending for dsRNA and dsDNA. Moreover, for dsRNA, ethanol causes slight ion-clamping across the major groove, which weakens the major groove-mediated twist-stretch coupling, while for dsDNA, ethanol promotes the stretch-radius correlation due to enhanced ion binding and consequently enhances the helical radius-mediated twist-stretch coupling.

Structure-based drug-development study against fibroblast growth factor receptor 2: molecular docking and Molecular dynamics simulation approaches

Developing new therapeutic strategies to target specific molecular pathways has become a primary focus in modern drug discovery science. Fibroblast growth factor receptor 2 (FGFR2) is a critical signaling protein involved in various cellular processes and implicated in numerous diseases, including cancer. Existing FGFR2 inhibitors face limitations like drug resistance and specificity issues. In this study, we present an integrated structure-based bioinformatics analysis to explore the potential of FGFR2 inhibitors-like compounds from the PubChem database with the Tanimoto threshold of 80%. We conducted a structure-based virtual screening approach on a dataset comprising 2336 compounds sourced from the PubChem database. Primarily, the selection of promising compounds was based on several criteria, such as drug-likeness, binding affinities, docking scores, and selectivity. Further, we conducted all-atom molecular dynamics (MD) simulations for 200 ns, followed by an essential dynamics analysis. Finally, a promising FGFR2 inhibitor with PubChem CID:507883 (1-[7-(1H-benzimidazol-2-yl)-4-fluoro-1H-indol-3-yl]-2-(4-benzoylpiperazin-1-yl)ethane-1,2-dione) was screened out from the study. This compound indicates a higher potential for inhibiting FGFR2 than the control inhibitor, Zoligratinib. The identified compound, CID:507883 shows >80% structural similarity with Zoligratinib. ADMET analysis showed promising pharmacokinetic potential of the screened compound. Overall, the findings indicate that the compound CID:507883 may have promising potential to serve as a lead candidate against FGFR2 and could be further exploited in therapeutic development.

Single-molecule sensing inside stereo- and regio-defined hetero-nanopores.

Heteromeric pore-forming proteins often contain recognition patterns or stereospecific selection filters. However, the construction of heteromeric pore-forming proteins for single-molecule sensing is challenging due to the uncontrollability of producing position isomers and difficulties in purification of regio-defined products. To overcome these preparation obstacles, we present an in situ strategy involving single-molecule chemical modification of a heptameric pore-forming protein to build a stereo- and regio-specific heteromeric nanopore (hetero-nanopore) with a subunit stoichiometric ratio of 3:4. The steric hindrance inherent in the homo-nanopore of K238C aerolysin directs the stereo- and regio-selective modification of maleimide derivatives. Our method utilizes real-time ionic current recording to facilitate controlled voltage manipulation for stoichiometric modification and position-based side-isomer removal. Single-molecule experiments and all-atom molecular dynamics simulations revealed that the hetero-nanopore features an asymmetric stereo- and regio-defined residue structure. The hetero-nanopore produced was characterized by mass spectrometry and single-particle cryogenic electron microscopy. In a proof-of-concept single-molecule sensing experiment, the hetero-nanopore exhibited 95% accuracy for label-free discrimination of four peptide stereoisomers with single-amino-acid structural and chiral differences in the mixtures. The customized hetero-nanopores could advance single-molecule sensing.

Membrane curvature sensing and symmetry breaking of the M2 proton channel from Influenza A.

The M2 proton channel aids in the exit of mature influenza viral particles from the host plasma membrane through its ability to stabilize regions of high negative Gaussian curvature (NGC) that occur at the neck of budding virions. The channels are homo-tetramers that contain a cytoplasm-facing amphipathic helix (AH) that is necessary and sufficient for NGC generation; however, constructs containing the transmembrane spanning helix, which facilitates tetramerization, exhibit enhanced curvature generation. Here, we used all-atom molecular dynamics (MD) simulations to explore the conformational dynamics of M2 channels in lipid bilayers revealing that the AH is dynamic, quickly breaking the fourfold symmetry observed in most structures. Next, we carried out MD simulations with the protein restrained in four- and twofold symmetric conformations to determine the impact on the membrane shape. While each pattern was distinct, all configurations induced pronounced curvature in the outer leaflet, while conversely, the inner leaflets showed minimal curvature and significant lipid tilt around the AHs. The MD-generated profiles at the protein-membrane interface were then extracted and used as boundary conditions in a continuum elastic membrane model to calculate the membrane-bending energy of each conformation embedded in different membrane surfaces characteristic of a budding virus. The calculations show that all three M2 conformations are stabilized in inward-budding, concave spherical caps and destabilized in outward-budding, convex spherical caps, the latter reminiscent of a budding virus. One of the C2-broken symmetry conformations is stabilized by 4 kT in NGC surfaces with the minimum energy conformation occurring at a curvature corresponding to 33 nm radii. In total, our work provides atomistic insight into the curvature sensing capabilities of M2 channels and how enrichment in the nascent viral particle depends on protein shape and membrane geometry.

Molecular simulations of understanding the Zn2+ ion structure, dynamics and thermodynamic properties in water in ionic liquids

We utilize all-atom molecular dynamics simulations to explore the intermolecular structure, dynamics and thermodynamic properties of Zn2+ ions in water-in-ionic liquid. Two ionic liquid systems featuring the same cation 1-ethyl-3-methyl-imidazolium [EMIM]+ and distinct anions tetrafluoroborate [BF4]− and hexafluorophosphate [PF6]− are considered. We consider the water (H2O) mole fractions (x) varying from 0.33 to 0.71. We observe a significant interactions of Zn2+ ions with water in the case of [BF4]− when compared to [PF6]−. On the other hand Zn2+ ions mobility rises more in [EMIM]+[BF4]− as compared to [EMIM]+[PF6]− with x. Higher self-diffusion (D) of Zn2+ ions is seen in the case of [EMIM]+[BF4]−. The ionic conductivity (σ) of [EMIM]+[BF4]− is greater compared to [EMIM]+[PF6]− with the rise in x. Overall, this article furnishes in-depth molecular-level insights into the behaviour of Zn2+ ions in the presence of water mixed ionic liquid electrolytes.

PyCTRAMER: A Python package for charge transfer rate constant of condensed-phase systems from Marcus theory to Fermi's golden rule.

In this work, we introduce PyCTRAMER, a comprehensive Python package designed for calculating charge transfer (CT) rate constants in disordered condensed-phase systems at finite temperatures, such as organic photovoltaic (OPV) materials. PyCTRAMER is a restructured and enriched version of the CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory) package [Tinnin et al. J. Chem. Phys. 154, 214108 (2021)], enabling the computation of the Marcus CT rate constant and the six levels of the linearized semiclassical approximations of Fermi's golden rule (FGR) rate constant. It supports various types of intramolecular and intermolecular CT transitions from the excitonic states to CT state. Integrating quantum chemistry calculations, all-atom molecular dynamics (MD) simulations, spin-boson model construction, and rate constant calculations, PyCTRAMER offers an automatic workflow for handling photoinduced CT processes in explicit solvent environments and interfacial CT in amorphous donor/acceptor blends. The package also provides versatile tools for individual workflow steps, including electronic state analysis, state-specific force field construction, MD simulations, and spin-boson model construction from energy trajectories. We demonstrate the software's capabilities through two examples, highlighting both intramolecular and intermolecular CT processes in prototypical OPV systems.

Proton-triggered rearrangement of the AMPA receptor N-terminal domains impacts receptor kinetics and synaptic localization.

AMPA glutamate receptors (AMPARs) are ion channel tetramers that mediate the majority of fast excitatory synaptic transmission. They are composed of four subunits (GluA1-GluA4); the GluA2 subunit dominates AMPAR function throughout the forebrain. Its extracellular N-terminal domain (NTD) determines receptor localization at the synapse, ensuring reliable synaptic transmission and plasticity. This synaptic anchoring function requires a compact NTD tier, stabilized by a GluA2-specific NTD interface. Here we show that low pH conditions, which accompany synaptic activity, rupture this interface. All-atom molecular dynamics simulations reveal that protonation of an interfacial histidine residue (H208) centrally contributes to NTD rearrangement. Moreover, in stark contrast to their canonical compact arrangement at neutral pH, GluA2 cryo-electron microscopy structures exhibit a wide spectrum of NTD conformations under acidic conditions. We show that the consequences of this pH-dependent conformational control are twofold: rupture of the NTD tier slows recovery from desensitized states and increases receptor mobility at mouse hippocampal synapses. Therefore, a proton-triggered NTD switch will shape both AMPAR location and kinetics, thereby impacting synaptic signal transmission.

Tailoring Phage Nanosomes for Enhanced Theranostic Properties of Near Infrared Dyes.

Near-infrared (NIR) phototherapies offer noninvasive, cost-effective solutions for treating tumors and microbial infections. However, organic NIR dyes commonly used suffer from solubility and stability issues requiring frequent dosing. We address this challenge by exploring the bacteriophage-mediated enhancement of NIR dye properties. Upon encapsulation within phage nanosomes, IR780 and Indocyanine green (ICG), with similar optical properties but distinct water solubility and exhibit enhanced UV-vis absorbance and photothermal transduction efficacy compared to liposomes. Experimental characterization corroborated with all-atom molecular dynamics simulations imprints the nanoscale structure, solubility, dynamics, and binding of these NIR dye molecules to the membrane and protein molecules present in Phage capsid. These NIR dye-loaded phage nanosomes, coencapsulated with mitoxantrone, demonstrate enhanced anticancer activity, and when combined with amphotericin B, these dye molecules exhibit superior photothermal effects against fungal infections. Our findings present a simple and efficient approach for tuning the photothermal performance of existing NIR dyes through a rational design for enhanced therapeutic outcomes.

Investigating the Effect of Chain Extender on the Phase Separation and Mechanical Properties of Polybutadiene-Based Polyurethane.

The thermodynamic incompatibility between the soft and hard segments of thermoplastic polyurethane (TPU) results in a microphase-separated behavior and excellent mechanical properties. However, the effect of the chain extender on the degree of microphase separation (DMS) and the resultant mechanical properties of TPU have not been well studied because of the complex interactions between the soft and hard segments. Herein, hydroxyl-terminated polybutadiene-based TPUs(HTPB-TPUs) without hydrogen bonding between the soft and hard segments are synthesized using hydroxyl-terminated polybutadiene, toluene diisocyanate, and four different chain extenders, and the effect of the chain extender structure on DMS is analyzed experimentally using a combination of analytical techniques. Furthermore, the solubility parameters of the soft and hard segments, glass transition temperatures, and hydrogen-bond density of the HTPB-TPUs, are computed using all-atom molecular dynamics simulations. The results clearly reveal that the chain extender significantly affects the DMS and thus the mechanical properties of HTPB-TPUs. This study paves the way for studying the relationship between the structure and properties of TPU.

Molecular Simulations Study of the Ionic Adsorption on Oscillating Carbon Electrodes

Molecular dynamics simulations are a method of choice to study, at microscopic scales, the adsorption of electrolyte ions onto the electrodes of capacitive storage systems such as supercapacitors1.Here, we use all-atom molecular dynamics simulations to model a pure ionic liquid electrolyte ([EMIM][TFSI], 1-ethyl-3-methylimidazolium trifluoromethanesulfonylimide) in contact with graphite electrodes (Figure: Supercapacitor composed of graphite electrodes (in light blue) in contact with a pure ionic liquid electrolyte ([EMIM][TFSI], atoms colored according to the atom type. The arrows illustrate the oscillations of electrodes used to probe the influence of such a motion on the interfacial liquid.)When applying a potential difference between the two carbon electrodes, a migration of counter ions is observed, ultimately leading to a compensation of the charges carried by the electrode atoms. The adsorption dynamics depend on several parameters (electrostatic interactions, coordination number, viscosity of the electrolyte, size of the ions etc.) and influence the power of these capacitive systems2. In order to improve the performance of supercapacitors, it is then important to describe and understand interfacial properties at a molecular scale.Experimentally, an Electrochemical Quartz Crystal Microbalance (EQCM) can be used to study ion flow on top of the electrode. The frequency variations recorded by the EQCM are linked to variations in the mass adsorbed on the quartz surface, thus making it possible to study the dynamics of ion adsorption during the charge of a supercapacitor.3 However, variations in quartz frequency can influence the stability of the electrode-electrolyte interface, thus modifying ionic interactions and adsorption dynamics at the electrode surface4.Here, thanks to molecular simulations, we have access to the adsorption rate of ions, to their re-orientation, to the variation of charge of the electrode atoms during adsorption etc. This allows us to study in detail the mechanisms governing the electrolyte-electrode interface during the charging of supercapacitors as part of EQCM. And study the extent to which perturbations caused by quartz oscillations in the ionic liquid affect electrode-electrolyte interfacial properties. References 1 K. Xu, H. Shao, Z. Lin, C. Merlet, G. Feng, et al.. ENERGY & ENVIRONMENTAL MATERIALS, 3, 235-246 (2020) 2C. Péan, B. Rotenberg, P. Simon, M. Salanne, ELECTROCHIM. ACTA, 206, 504-512 (2016) 3Y-C. Wu, J. Ye, G. Jiang, K. Ni, N. Shu, P-L. Taberna, Y. Zhu, P. Simon, ANGEW. CHEM. INT. ED., 60, 13317-13322 (2021) 4Nikolai V. Priezjev, MICROFLUIDICS AND NANOFLUIDICS, 14, 225-233 (2013) Figure 1