Thermodynamics of Hydrogen Bonding in Hydrophilic and Hydrophobic Media

Effects of an external static EF on the conformational transition of 5-HT1A receptor: A molecular dynamics simulation study

The serotonin receptor subtype 1A (5-HT1AR), one of the G-protein-coupled receptor (GPCR) family, has been implicated in several neurological conditions. Understanding the activation and inactivation mechanism of 5-HT1AR at the molecular level is critical for discovering novel therapeutics in many diseases. Recently there has been a growing appreciation for the role of external electric fields (EFs) in influencing the structure and activity of biomolecules. In this study, we used molecular dynamics (MD) simulations to examine conformational features of active states of 5-HT1AR and investigate the effect of an external static EF with 0.02 V/nm applied on the active state of 5-HT1AR. Our results showed that the active state of 5-HT1AR maintained the native structure, while the EF led to structural modifications in 5-HT1AR, particularly inducing the inward movement of transmembrane helix 6 (TM6). Furthermore, it disturbed the conformational switches associated with activation in the CWxP, DRY, PIF, and NPxxY motifs, consequently predisposing an inclination towards the inactive-like conformation. We also found that the EF led to an overall increase in the dipole moment of 5-HT1AR, encompassing TM6 and pivotal amino acids. The analyses of conformational properties of TM6 showed that the changed secondary structure and decreased solvent exposure occurred upon the EF condition. The interaction of 5-HT1AR with the membrane lipid bilayer was also altered under the EF. Our findings reveal the molecular mechanism underlying the transition of 5-HT1AR conformation induced by external EFs, which offer potential novel insights into the prospect of employing structure-based EF applications for GPCRs.

Influence of the Water Model on the Structure and Interactions of the GPR40 Protein with the Lipid Membrane and the Solvent: Rigid versus Flexible Water Models

G protein-coupled receptors (GPCR) are responsible for modulating various physiological functions and are thus related to the pathophysiology of different diseases. Being potential therapeutic targets, multiple computational methodologies have been developed to analyze their behavior and interactions with other species. The solvent, on the other hand, has received much less attention. In this work, we analyzed the effect of four explicit water models on the structure and interactions of the GPR40 receptor in its apo form. We employed the rigid SPC/E and TIP4P models, and their flexible versions, the FBA/ϵ and TIP4P/ϵflex. We explored the structural changes and their correlation with some bulk dynamic properties of water. Our results showed an adverse effect on the conservation of the secondary structure of the receptor with all the models due to the breaking of the intramolecular hydrogen bond network, being more evident for the TIP4P models. Notably, all four models brought the receptor to states similar to the active one, modifying the intracellular part of the TM5 and TM6 domains in a "hinge" type movement, allowing the opening of the structure. Regarding the dynamic properties, the rigid models showed results comparable to those obtained in other studies on membrane systems. However, flexible models exhibit disparities in the molecular representation of systems. Surprisingly, the FBA/ϵ model improves the molecular picture of several properties, even though their agreement with bulk diffusion is poorer. These findings reinforce our idea that exploring other water models or improving the current ones, to better represent the membrane interface, can lead to a positive impact on the description of the signal transduction mechanisms and the search of new drugs by targeting these receptors.

Exploring How System Dimensions and Periodic Boundary Conditions Influence the Molecular Dynamics Simulation of A6H Peptide Self-Assembly Nanostructures

This work presents a study on the effects of periodic boundary conditions (PBC) on the energetic/structural properties and hydrogen bond dynamics (HB) using molecular dynamics (MD) simulations of peptide membranes composed of alanine and histidine. Our results highlight that simulations using small surface areas for the peptide membrane may result in nonconvergent values for membrane properties, which are only observed in regions simulated at a certain distance from the PBCs. Specifically, regarding hydrogen bonds, a property pervasive in peptide membranes, our findings indicate a significant increase in the lifetime of these interactions, reaching values ∼19% higher when observed in structures free from PBCs. For peptide mobility in these nanomembranes, our results compare regions simulated directly under the influence of PBCs with regions free from these conditions, emphasizing greater mobility of amino acid psi/phi angles in the latter model.

Assessing the Efficacy of Acanthoic Acid Isolated from Acanthopanax koreanum Nakai in Male Infertility: An In Vivo and In Silico Approach

Acanthoic acid, a diterpene isolated from the root bark of Acanthopanax koreanum Nakai, possesses diverse pharmacological activities, including anti-inflammatory, anti-diabetic, gastrointestinal protection, and cardiovascular protection. This study is the first to investigate the egg-hatching rates of Drosophila melanogaster affected by acanthoic acid. Notably, male flies supplemented with 10 μM acanthoic acid exhibited a strong increase in hatching rates compared with controls under adverse temperature conditions, suggesting a potential protective effect against environmental stressors. Molecular docking simulations revealed the binding affinities and specific interactions between acanthoic acid and proteins related to male infertility, including SHBG, ADAM17, and DNase I, with binding affinity values of -10.2, -6.8, and -5.8 kcal/mol, respectively. Following the docking studies, molecular dynamic simulations were conducted for a duration of 100 ns to examine the stability of these interactions. Additionally, a total binding energy analysis and decomposition analysis offered insights into the underlying energetic components and identified key contributing residues.

Probing the energy landscape of the lipid interactions of the serotonin1A receptor

G protein-coupled receptors (GPCRs) are lipid-regulated transmembrane proteins that play a central role in cell signaling and pharmacology. Although the role of membrane lipids in GPCR function is well established, the underlying GPCR-lipid interactions have not been thermodynamically characterized due to the complexity of these interactions. In this work, we estimate the energetics and dynamics of lipid association from coarse-grain simulations of the serotonin1A receptor embedded in a complex membrane. We show that lipids bind to the receptor with varying energetics of 1-4 kT, and timescales of 1-10 μs. The most favorable energetics and longest residence times are observed for cholesterol, glycosphingolipid GM1, phosphatidylethanolamine (PE) and phosphatidylserine (PS) lipids. Multi-exponential fitting of the contact probability suggests distinct dynamic regimes, corresponding to ps, ns and μs timescales, that we correlate with the annular, intermediate and non-annular lipid sites. The timescales of lipid binding correspond to high barrier heights, despite their relatively weaker energetics. Our results highlight that GPCR-lipid interactions are driven by both thermodynamic interactions and the dynamical features of lipid binding.

A Virtual Drug Discovery Screening Illuminates Campesterol as a Potent Estrogen Receptor Alpha Inhibitor in Breast Cancer

Breast cancer remains a global health challenge, and innovative strategies are required to target estrogen receptor α (ERα), a key player in its development. This study investigates the potential of campesterol, a natural phytosterol, as an ERα inhibitor for breast cancer. Our approach integrates in silico, in vitro, and ex vivo experiments to assess the therapeutic potential of campesterol. In silico analyses highlight campesterol as a promising ERα ligand with favorable binding affinities and dynamic properties. Structural analysis reveals conformational changes in ERα upon campesterol binding. In vitro studies confirm the selective growth inhibition of campesterol against ERα-positive breast cancer cells. This study extends to ER+ breast cancer patient-derived organoids (PDOs), showing the effectiveness of campesterol in ERα-positive breast cancer PDOs. Importantly, it emphasizes the receptor-specific nature of campesterol, providing insights into its context-dependent action. In conclusion, campesterol displays potential as an ERα inhibitor, offering new avenues for ER+ breast cancer treatment.

Strength, number, and kinetics of hydrogen bonds for water confined inside boron nitride nanotubes

Water has shown a myriad of highly interesting properties and behaviors, such as very low friction, phase transition under unexpected conditions, massive property alterations, etc. inside strong nanoconfinements of few-nanometer to sub-nanometer diameters. Water-water hydrogen bonding is one of the most important factors dictating such water behavior and properties inside such strong nanoconfinements. In this paper, we employ Reactive Force Field (ReaxFF) molecular dynamics (MD) simulations for studying multiple facets of such water-water hydrogen bonds (HBs) inside boron-nitride nanotubes (BNNTs) having diameters ranging from a few nanometers to sub-nanometers. First, the strength of the water-water HB interactions, as a function of the HB configuration, is quantified by studying the corresponding PMF (potential of mean force). For water present in extreme confinements (BNNTs with sub-nanometric diameters), we see completely isolated HB basins. On the other hand, for bulk water the HB basin is connected (via a saddle point) to a nearby second PMF well. Therefore, our analysis successfully distinguishes the HB characteristics between the cases of water in extreme confinement and bulk water. Second, we study the kinetics of such water-water HBs: HBs formed by a given pair of water molecules in extreme confinements show a much larger probability of remaining intact once formed or re-forming after they have been broken. Both these results, which shed new light on water-water hydrogen bonding inside strong nanoconfinements, can be explained by the single-file structure formed by the water molecules in extreme BNNT confinements.

Rational Approach toward COVID-19's Main Protease Inhibitors: A Hierarchical Biochemoinformatics Analysis

This study investigated the potential of selected compounds as inhibitors of SARS-CoV-2 Mpro through pharmacokinetic and toxicological analyses, molecular docking, and molecular dynamics simulations. In silico molecular docking simulations revealed promising ligands with favorable binding affinities for Mpro, ranging from -6.2 to -9.5 kcal/mol. Moreover, molecular dynamics simulations demonstrated the stability of protein-ligand complexes over 200 ns, maintaining protein secondary structures. MM-PBSA analysis revealed favorable interactions between ligands and Mpro, with negative binding energy values. Hydrogen bond formation capacity during molecular dynamics was confirmed, indicating consistent interactions with Mpro catalytic residues. Based on these findings, selected ligands show promise for future studies in developing COVID-19 treatments.

Unraveling the binding mechanisms of transglutaminase and substrate subjected to microwaves: Molecular docking and molecular dynamic simulations

Previous studies showed that transglutaminase (TGase) and microwaves acted synergistically to improve the functional properties of proteins. The mechanism behind this has yet to be elucidated. In this study, the phenomenon of microwaves enhancing TGase activity was experimentally validated. Molecular docking and molecular dynamics simulations revealed that moderate microwaves (105 and 108 V/m) increased the structural flexibility of TGase and promoted the orientation of the side chain carboxylate anion group on Asp255, driving the reaction forward. Also, TGase underwent partial transformation from α-helix to turns or coils at 105 and 108 V/m, exposing more residues in the active site and facilitating the binding of the substrate (CBZ-Gln-Gly) to TGase. However, 109 V/m microwaves completely destroyed the TGase structure, inactivating the enzyme. This study provides insights into the molecular mechanisms underlying the interactions between TGase and substrate subjected to microwaves, promoting the future applications of TGase and microwaves in food processing.

Molecular Insights into the Differential Effects of Acetylation on the Aggregation of Tau Microtubule-Binding Repeats

Aggregation of tau protein into intracellular fibrillary inclusions is characterized as the hallmark of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The microtubule-binding (MTB) domain of tau, containing either three or four repeats with sequence similarities, plays an important role in determining tau's aggregation. Previous studies have reported that abnormal acetylation of lysine residues displays a distinct effect on the formation of pathological tau aggregates. However, the underlying molecular mechanism remains mostly elusive. In this study, we performed extensive replica exchange molecular dynamics (REMD) simulations of 144 μs in total to systematically investigate the dimerization of four tau MTB repeats and explore the impacts of Lys280 (K280) or Lys321 (K321) acetylation on the conformational ensembles of the R2 or R3 dimer. Our results show that R3 is the most prone to aggregation among the four repeats, followed by R2 and R4, while R1 displays the weakest aggregation propensity with a disordered structure. Acetylation of K280 could promote the aggregation of R2 peptides by increasing the formation of β-sheet structures and strengthening the interchain interaction. However, K321 acetylation decreases the β-sheet content of the R3 dimer, reduces the ability of R3 peptides to form long β-strands, and promotes the stable helix structure formation. The salt bridge and Y310-Y310 π-π stacking interactions of the R3 dimer are greatly weakened by K321 acetylation, resulting in the inhibition of dimerization. This study uncovers the structural ensembles of tau MTB repeats and provides mechanistic insights into the influences of acetylation on tau aggregation, which may deepen the understanding of the pathogenesis of tauopathies.

Pharmacokinetics and drug-likeness of anti-cancer traditional Chinese medicine: molecular docking and molecular dynamics simulation study

MCM7 (Minichromosome Maintenance Complex Component 7) is a component of the DNA replication licensing factor, which controls DNA replication. The MCM7 protein is linked to tumor cell proliferation and has a function in the development of several human cancers. Several types of cancer may be treated by inhibiting the protein, as it is strongly produced throughout this process. Significantly, Traditional Chinese Medicine (TCM), which has a long history of clinical adjuvant use against cancer, is rapidly gaining traction as a valuable medical resource for the development of novel cancer therapies, including immunotherapy. Therefore, the goal of the research was to find small molecular therapeutic candidates against the MCM7 protein that may be used to treat human cancers. A computational-based virtual screening of 36,000 natural TCM libraries is carried out for this goal using a molecular docking and dynamic simulation technique. Thereby, ∼8 novel potent compounds i.e., ZINC85542762, ZINC95911541, ZINC85542617, ZINC85542646, ZINC85592446, ZINC85568676, ZINC85531303, and ZINC95914464 were successfully shortlisted, each having the capacity to penetrate the cell as potent inhibitors for MCM7 to curb this disorder. These selected compounds were found to have high binding affinities compared to the reference (AGS compound) i.e. < -11.0 kcal/mol. ADMET and pharmacological properties showed that none of these 8 compounds poses any toxic property (carcinogenicity) and have anti-metastatic, and anticancer activity. Additionally, MD simulations were run to assess the compounds' stability and dynamic behavior with the MCM7 complex for about 100 ns. Finally, ZINC95914464, ZINC95911541, ZINC85568676, ZINC85592446, ZINC85531303, and ZINC85542646 are identified as highly stable within the complex throughout the 100 ns simulations. Moreover, the results of binding free energy suggested that the selected virtual hits significantly bind to the MCM7 which implied these compounds may act as a potential MCM7 inhibitor. However, in vitro testing protocols are required to further support these results. Further, assessment through various lab-based trial methods can assist with deciding the action of the compound that will give options in contrast to human cancer immunotherapy.Communicated by Ramaswamy H. Sarma.

Second inflection point of supercooled water surface tension induced by hydrogen bonds: A molecular-dynamics study

Surface tension of supercooled water is a fundamental property in various scientific processes. In this study, we perform molecular dynamics simulations with the TIP4P-2005 model to investigate the surface tension of supercooled water down to 220 K. Our results show a second inflection point (SIP) in the surface tension at temperature TSIP ≈ 267.5 ± 2.3 K. Using an extended IAPWS-E functional fit for the water surface tension, we calculate the surface excess internal-energy and entropy terms of the excess Helmholtz free energy. Similar to prior studies [Wang et al., Phys. Chem. Chem. Phys. 21, 3360 (2019); Gorfer et al., J. Chem. Phys. 158, 054503 (2023)], our results show that the surface tension is governed by two driving forces: a surface excess entropy change above the SIP and a surface excess internal-energy change below it. We study hydrogen-bonding near the SIP because it is the main cause of water's anomalous properties. With decreasing temperature, our results show that the entropy contribution to the surface tension reaches a maximum slightly below the SIP and then decreases. This is because the number of hydrogen bonds increases more slowly below the SIP. Moreover, the strengths and lifetimes of the hydrogen bonds also rise dramatically below the SIP, causing the internal-energy term to dominate the excess surface free energy. Thus, the SIP in the surface tension of supercooled TIP4P-2005 water is associated with an increase in the strengths and lifetimes of hydrogen bonds, along with a decrease in the formation rate (#/K) of new hydrogen bonds.

De novo antioxidant peptide design via machine learning and DFT studies

Antioxidant peptides (AOPs) are highly valued in food and pharmaceutical industries due to their significant role in human function. This study introduces a novel approach to identifying robust AOPs using a deep generative model based on sequence representation. Through filtration with a deep-learning classification model and subsequent clustering via the Butina cluster algorithm, twelve peptides (GP1-GP12) with potential antioxidant capacity were predicted. Density functional theory (DFT) calculations guided the selection of six peptides for synthesis and biological experiments. Molecular orbital representations revealed that the HOMO for these peptides is primarily localized on the indole segment, underscoring its pivotal role in antioxidant activity. All six synthesized peptides exhibited antioxidant activity in the DPPH assay, while the hydroxyl radical test showed suboptimal results. A hemolysis assay confirmed the non-hemolytic nature of the generated peptides. Additionally, an in silico investigation explored the potential inhibitory interaction between the peptides and the Keap1 protein. Analysis revealed that ligands GP3, GP4, and GP12 induced significant structural changes in proteins, affecting their stability and flexibility. These findings highlight the capability of machine learning approaches in generating novel antioxidant peptides.

Kinetic description of water transport during spontaneous emulsification induced by Span 80

Spontaneous emulsification is a phenomenon that forms nanometer-sized droplets (nanodroplets) without the application of any external force, and the mechanism has been actively studied for application to various technologies. In this study, we analyzed the kinetics of spontaneous emulsification induced by Span 80. The measurement of water concentration in Span 80 hexadecane solution indicated that the chemical potential of water in the nanodroplets decreased as the amount of water in the nanodroplets decreased. Based on this result, water transport between the aqueous phase and nanodroplets in which the chemical potential of water was controlled was quantitatively investigated by using a microfluidic device. The results demonstrate that the kinetics of water transport during spontaneous emulsification induced by Span 80 was described by a model of osmotic transport through an organic liquid film between the aqueous phase and nanodroplets.

Atomic insights into the inhibition of R3 domain of tau protein by epigallocatechin gallate, quercetin and gallic acid

Inhibiting tau protein aggregation has become a prospective avenue for the therapeutic development of tauopathies. The third microtubule-binding repeat (R3) domain of tau is confirmed as the most aggregation-favorable fragment of the whole protein. As dimerization is the first step of the aggregation of tau into amyloid fibrils, impeding the dimerization of the R3 domain is critical to prevent the full-length tau aggregation. Natural polyphenol small molecules epigallocatechin gallate (EGCG), quercetin (QE) and gallic acid (GA) are proven to inhibit the aggregation of the full-length recombinant tau (For EGCG and QE) or the R3 domain (For GA) of tau in vitro. However, the underlying molecular mechanisms of the inhibitive effects on the R3 domain of tau remain largely unknown. In this study, we conducted numerous all-atom molecular dynamics simulations on R3 dimers with and without EGCG, QE or GA, respectively. The results reveal that all three molecules can effectively decrease the β structure composition of the R3 dimer, induce the dimer to adopt loosely-packed conformations, and weaken interchain interactions, thus impeding the dimerization of the R3 peptide chains. The specific preferentially binding sites for the three molecules exhibit similarities and differences. Hydrophobic, π-π stacking and hydrogen-bonding interactions collectively drive EGCG, QE and GA respectively binding on the R3 dimer, while QE also binds with the dimer through cation-π interaction. Given the incurable nature of tauopathies hitherto, our research provides helpful knowledge for the development of drugs to treat tauopathies.

N6-methyladenosine in 5' UTR does not promote translation initiation

The most abundant N6-methyladenosine (m6A) modification on mRNAs is installed non-stoichiometrically across transcripts, with 5' untranslated regions (5' UTRs) being the least conductive. 5' UTRs are essential for translation initiation, yet the molecular mechanisms orchestrated by m6A remain poorly understood. Here, we combined structural, biochemical, and single-molecule approaches and show that at the most common position, a single m6A does not affect translation yields, the kinetics of translation initiation complex assembly, or start codon recognition both under permissive growth and following exposure to oxidative stress. Cryoelectron microscopy (cryo-EM) structures of the late preinitiation complex reveal that m6A purine ring established stacking interactions with an arginine side chain of the initiation factor eIF2α, although with only a marginal energy contribution, as estimated computationally. These findings provide molecular insights into m6A interactions with the initiation complex and suggest that the subtle stabilization is unlikely to affect the translation dynamics under homeostatic conditions or stress.

Direct Observation of the Mutual Coupling Effect in the Protein-Water-Glycerol Mixture by Combining Neutron Scattering and Selective Deuteration

Numerous studies have discussed the impact of cosolvents on the structure, dynamics, and stability of proteins in aqueous solutions. However, the dynamics of cosolvents in the protein-water-cosolvent ternary system is largely unexplored in experiments due to technical difficulty. Consequently, a comprehensive understanding of the interplay among proteins, water, and cosolvents is still lacking. Here, we employed selective deuteration and neutron scattering techniques to characterize the individual motions of each component in the protein/water/glycerol (GLY) mixture across various temperatures. The consistent dynamic onset temperatures and the correlation between the MSD of the protein and the viscosity of solvents revealed the mutual coupling effects among the three components. Furthermore, our experimental and simulation results showed that the hydrogen bond relaxation energy barrier in the ternary system is ∼43 kJ/mol, whereas in the protein-water binary system it is merely ∼35 kJ/mol. Therefore, we suggest that GLY can enhance hydrogen bond interactions in the ternary system through the mutual coupling effect, thereby serving as one of the protective mechanisms of protein preservation by GLY.

Theoretical models of staurosporine and analogs uncover detailed structural information in biological solution

Staurosporine and its analogs (STA-analogs) are indolocarbazoles (ICZs) compounds able to inhibit kinase proteins in a non-specific way, while present antimicrobial and cytostatic properties. The knowledge of molecular features associated to the complexation, including the ligand shape in solution and thermodynamics of complexation, is substantial to the development of new bioactive ICZs with improved therapeutic properties. In this context, the empirical approach of GROMOS force field is able to accurately reproduce condensed phase physicochemical properties of molecular systems after parameterization. Hence, through parameterization under GROMOS force field and molecular simulations, we assessed STA-analogs dynamics in aqueous solution, as well as its interaction with water to probe conformational and structural features involved in complexation to therapeutic targets. The coexistence of multiple conformers observed in simulations, and confirmed by metadynamics calculations, expanding the conformational space knowledge of these ligands with potential implications in understanding the ligand conformational selection during complexation. Also, changes in availability to H-bonding concerning the different substituents and water can reflect on effects at complexation free energy due to variation at the desolvation energetic costs. Based on these results, we expect the obtained structural data provide systemic framework for rational chemical modification of STA-analogs.

300-fold higher neuro- and immunotoxicity from low-redox transformation of carbamazepine

Current challenges in (eco)toxicology are in understanding the transformation of (reactive) substances, and how transformation affects toxic modes of action. Empirical assessment of transformation products of, practically an infinite number of substances, via experimentation, is impossible. Predicting transformation products for (benchmarking) compounds from conditions, facilitates risk analyses. This study applied calculus to predict transformation products of an important environmental and medicinal/toxicological marker, carbamazepine. As radicals are ubiquitous in humans and the environment, we looked into radical-mediated transformations of carbamazepine as a benchmark. We calculated proportions of their speciation states as function of redox conditions, which we took as pH and O2 concentration, describing transformation via covalent and ionic interactions. Formation of ring-contracted products with neuro-immunological activity is thermodynamically favored under anaerobic conditions and at low pH. Experimentally observed product distributions and toxicities reflect that pattern. Our predictive method may support toxicity predictions for other substances and conditions 'similar' to the current case study via interpolation. This paves the way for a more coherent, effective and easier risk assessment of transformation products.

Improving the Thermostability of Thermomyces lanuginosus Lipase by Restricting the Flexibility of N-Terminus and C-Terminus Simultaneously via the 25-Loop Substitutions

(1) Lipases are catalysts widely applied in industrial fields. To sustain the harsh treatments in industries, optimizing lipase activities and thermal stability is necessary to reduce production loss. (2) The thermostability of Thermomyces lanuginosus lipase (TLL) was evaluated via B-factor analysis and consensus-sequence substitutions. Five single-point variants (K24S, D27N, D27R, P29S, and A30P) with improved thermostability were constructed via site-directed mutagenesis. (3) The optimal reaction temperatures of all the five variants displayed 5 °C improvement compared with TLL. Four variants, except D27N, showed enhanced residual activities at 80 °C. The melting temperatures of three variants (D27R, P29S, and A30P) were significantly increased. The molecular dynamics simulations indicated that the 25-loop (residues 24-30) in the N-terminus of the five variants generated more hydrogen bonds with surrounding amino acids; hydrogen bond pair D254-I255 preserved in the C-terminus of the variants also contributes to the improved thermostability. Furthermore, the newly formed salt-bridge interaction (R27…E56) in D27R was identified as a crucial determinant for thermostability. (4) Our study discovered that substituting residues from the 25-loop will enhance the stability of the N-terminus and C-terminus simultaneously, restrict the most flexible regions of TLL, and result in improved thermostability.

Molecular dynamics simulations identify the topological weak spots of a protease CN2S8A

Thermophilic enzymes are highly desired in industrial applications due to their efficient catalytic activity at high temperature. However, most enzymes exhibit inferior thermostability and it remains challenging to identify the optimal sites for designing mutations to improve protein stability. To tackle this issue, we integrated topological analysis and all-atom molecular dynamics simulations to efficiently pinpoint the thermally-unstable regions in protein structures. Using a protease CN2S8A as the model, we analyzed the intramolecular hydrogen bonding interactions between adjacent secondary structure elements, and then identified the topological weak spots of CN2S8A where weak hydrogen bonding interactions were formed. To examine the role of these sites in protein structural stability, we designed three virtual mutations at different weak spots and characterized the effects of these mutations on the structural properties of CN2S8A. The results showed that all three mutations increased the protein structural stability. In conclusion, these findings provide a novel method to identify the topological weak spots of proteins, with implications in the rational design of biocatalysts with superior thermostability.

Complex role of chemical nature and tacticity in the adsorption free energy of carboxylic acid polymers at the oil-water interface: molecular dynamics simulations

Scientific understanding of the molecular structure and adsorption of polymers at oil-water liquid interfaces is very limited. In this study the adsorption free energy at the oil (CCl4)-water interface was estimated using umbrella sampling molecular dynamics simulations for six carboxylate type vinyl polymers differing in hydrophobic nature and tacticity: isotactic and syndiotactic poly(acrylic acid) (i-PAA, s-PAA), isotactic and syndiotactic poly(methacrylic acid) (i-PMA, s-PMA), and atactic and syndiotactic poly(ethylacrylic acid) (a-PEA, s-PEA). ΔGads values are in the order i-PMA < a-PEA < s-PEA < s-PAA < i-PAA < s-PMA. The results show the significant and complex influence of the chemical nature as well as tacticity of the polymer on its adsorption free energy as related to hydrogen bonding and orientation of bonds with respect to oil and water phases. The influence of tacticity is found to be the highest for PMA, which is interpreted to occur due to the balance between interactions among side groups and those occurring between side groups and solvent. Interactions between side-groups are crucial for determining the conformation of PAA (most hydrophilic) and the solvation of the side-group in water is crucial for determining the conformation of PEA (most hydrophobic). The adsorption of PMA represents the transition between these two dominating effects. The molecular contributions to the enthalpy of adsorption indicate that adsorption is favored mainly through two interactions: polymer-CCl4 and water-water.

Deciphering the Inhibitory Mechanism of Naphthoquinone-Dopamine on the Aggregation of Tau Core Fragments PHF6* and PHF6

The formation of neurofibrillary tangles by abnormal aggregation of tau protein is considered to be an important pathological characteristic of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. Two hexapeptides 275VQIINK280 and 306VQIVYK311 in the microtubule binding region, named PHF6* and PHF6, are known to be aggregation-prone and responsible for tau fibrillization. Previous experiments reported that naphthoquinone-dopamine (NQDA) could effectively inhibit the aggregation of PHF6* and PHF6 and disrupt the fibrillar aggregates into nontoxic species, displaying a dual effect on the amyloid aggregation. However, the underlying molecular mechanism remains mostly elusive. Herein, we performed all-atom molecular dynamics (MD) simulations for 114 μs in total to systematically investigate the impacts of NQDA on the oligomerization of PHF6* and PHF6. The conformational ensembles of PHF6* and PHF6 peptides generated by replica exchange MD simulations show that NQDA could effectively prevent the hydrogen bond formation, reduce the ability of peptides to self-assemble into long β-strand and large β-sheets, and induce peptides to form a loosely packed and coil-rich oligomer. The interaction analysis shows that the binding of NQDA to PHF6* is mainly through hydrophobic interactions with residue I277 and hydrogen bonding interactions with Q276; for the PHF6 peptides, NQDA displays a strong π-π stacking interaction with residue Y310, thus impeding the Y310-Y310 π-π stacking and I308-Y310 CH-π interactions. The DA group of NQDA displays a stronger cation-π interaction than the NQ group, while the NQ group exhibits a stronger π-π stacking interaction. MD simulations demonstrate that NQDA prevents the conformational conversion to β-sheet-rich aggregates and displays an inhibitory effect on the oligomerization dynamics of PHF6* and PHF6. Our results provide a complete picture of inhibitory mechanisms of NQDA on PHF6* and PHF6 oligomerization, which may pave the way for designing drug candidates for the treatment of tauopathies.

Effects and mechanisms of anion and cation on the gelation of nanosilica sol by all-atom molecular dynamics simulation: promotion or inhibition?

Nanosilica sol (NSS) is prone to gelation due to the condensation of silicon hydroxyl at normal temperature and pressure, which is further exacerbated by the addition of electrolytes during production. Therefore, the effects of ions and the mechanism of gelation of NSS are crucial for its stability. Herein, all-atom molecular dynamics (AAMD) was carried out to explore the effects and mechanisms of cations (K+, Na+, Ca2+) and anions (Cl-, NO3-, SO42-, PO43-) on the sol-gel transition. Results indicated that highly electrophilic cations (e.g., Ca2+) and anions with slightly stronger nucleophilicity than Si(OH)3O- (e.g., NO3-) could inhibit gelation by preventing Si(OH)4 and Si(OH)3O- from approaching the silica surface. Such inhibition is more pronounced in NSS with larger particle sizes. Our findings offer some critical insights into the effects of ions on the gel stability of NSS, which also contributes significantly to screening suitable electrolytes for the production of NSS.

Impact of hierarchical water dipole orderings on the dynamics of aqueous salt solutions

Ions exhibit highly ion-specific complex behaviours when solvated in water, which remains a mystery despite the fundamental importance of ion solvation in nature, science, and technology. Here we explain these ion-specific properties by the ion-induced hierarchical dipolar, translational, and bond-orientational orderings of ion hydration shell under the competition between ion-water electrostatic interactions and inter-water hydrogen bonding. We first characterise this competition by a new length λHB(q), explaining the ion-specific effects on solution dynamics. Then, by continuously tuning ion size and charge, we find that the bond-orientational order of the ion hydration shell highly develops for specific ion size and charge combinations. This ordering drastically stabilises the hydration shell; its degree changes the water residence time around ions by 11 orders of magnitude for main-group ions. These findings are fundamental to ionic processes in aqueous solutions, providing a physical principle for electrolyte design and application.

Molecular Dynamics Simulations on the Adsorbed Monolayers of N-Dodecyl Betaine at the Air-Water Interface

Betaine is a kind of zwitterionic surfactant with both positive and negative charge groups on the polar head, showing good surface activity and aggregation behaviors. The interfacial adsorption, structures and properties of n-dodecyl betaine (NDB) at different surface coverages at the air-water interface are studied through molecular dynamics (MD) simulations. Interactions between the polar heads and water molecules, the distribution of water molecules around polar heads, the tilt angle of the NDB molecule, polar head and tail chain with respect to the surface normal, the conformations and lengths of the tail chain, and the interfacial thickness of the NDB monolayer are analyzed. The change of surface coverage hardly affects the locations and spatial distributions of the water molecules around the polar heads. As more NDB molecules are adsorbed at the air-water interface, the number of hydrogen bonds between polar heads and water molecules slightly decreases, while the lifetimes of hydrogen bonds become larger. With the increase in surface coverage, less gauche defects along the alkyl chain and longer NDB chain are obtained. The thickness of the NDB monolayer also increases. At large surface coverages, tilted angles of the polar head, tail chain and whole NDB molecule show little change with the increase in surface area. Surface coverages can change the tendency of polar heads and the tail chain for the surface normal.

Designing Nonconventional Luminescent Materials with Efficient Emission in Dilute Solutions via Modulation of Dynamic Hydrogen Bonds

Nonconventional luminescent materials (NLMs) which do not contain traditional aromatic chromophores are of great interest due to their unique chemical structures, optical properties, and their potential applications in various areas, such as cellular imaging and chemical sensing. However, most reported NLMs show weak or no emission in dilute solutions, which severely limits their applications. In this work, dynamic hydrogen bonds were utilized to design NLMs with efficient emission in dilute solutions. To further validate the results, polymers P1 and P2 were successfully prepared and investigated. It was found that the luminescence quantum efficiency of P1 and P2 at a concentration of 0.1 mg/mL in water solution was 8.9 and 0.6%, respectively. The high efficiency can be attributed to the fact that polymer P1 has more intra- or intermolecular dynamic hydrogen bonds and other short interactions than P2 in dilute solutions, allowing P1 to achieve the through-space conjugation effect to increase the degree of system conjugation, restrict molecular motion, and decrease nonradiative transitions, which can effectively improve luminescence. In addition, polymer P2 exhibits the characteristics of clustering-triggered emission, excitation wavelength-dependent and concentration-dependent fluorescence properties, excellent photobleaching resistance, low cytotoxicity, and selective recognition of Fe³⁺. The present study investigates the manipulation of luminescence properties of NLMs in dilute solutions through the modulation of dynamic hydrogen bonds. This approach can serve as a semi-empirical technique for designing and building innovative NLMs in the times ahead.

Insights from molecular dynamics and DFT calculations into the interaction of 1,4-benzodiazepines with 2-hydroxypropyl-βCD in a theoretical study

This study delves into the interaction between benzodiazepine (BZD) drugs and 2-hydroxypropyl-β-cyclodextrin (2HPβCD), a cyclodextrin (CD) known to improve drug delivery and enhance therapeutic outcomes. We find that the 2HPβCD's atoms become more rigid in the presence of chlordiazepoxide (CDP), clonazepam (CLZ), and diazepam (DZM), whereas they become more flexible in the presence of nordazepam (NDM) and nitrazepam (NZP). We also investigated the structure of 2HPβCD and found that loading these drugs increases both the area and volume of the 2HPβCD cavity, making it more suitable for drug delivery. Moreover, this research found that all drugs exhibited negative values for the binding free energy, indicating thermodynamic favorability and improved solubility. The binding free energy order of the BZDs was consistent in both molecular dynamics and Monte Carlo methods, with CDP and DZM having the highest affinity for binding. We also analyzed the contribution of different interaction energies in binding between the carrier and the drugs and found that Van der Waals energy is the primary component. Our results indicate that the number of hydrogen bonds between 2HPβCD/water slightly decreases in the presence of BZDs, but the hydrogen bond's quality remains constant.

Naphthoquinone-dopamine hybrids disrupt α-synuclein fibrils by their intramolecular synergistic interactions with fibrils and display a better effect on fibril disruption

α-Synuclein (αSyn) is an intrinsically disordered protein and its abnormal aggregation into amyloid fibrils is the main hallmark of Parkinson's disease (PD). The disruption of preformed αSyn fibrils using small molecules is considered as a potential strategy for PD treatment. Recent experiments have reported that naphthoquinone-dopamine hybrids (NQDA), synthesized by naphthoquinone (NQ) and dopamine (DA) molecules, can significantly disrupt αSyn fibrils and cross the blood-brain barrier. To unravel the fibril-disruptive mechanisms at the atomic level, we performed microsecond molecular dynamics simulations of αSyn fibrils in the absence and presence of NQDA, NQ, DA, or NQ+DA molecules. Our simulations showed that NQDA reduces the β-sheet content, disrupts K45-E57 and E46-K80 salt-bridges, weakens the inter-protofibril interaction, and thus destabilizes the αSyn fibril structure. NQDA has the ability to form cation-π and H-bonding interactions with K45/K80, and form π-π stacking interactions with Y39/F94. Those interactions between NQDA and αSyn fibrils play a crucial role in disaggregating αSyn fibrils. Moreover, we found that NQDA has a better fibril destabilization effect than that of NQ, DA, and NQ+DA molecules. This is attributed to the synergistic fibril-binding effect between NQ and DA groups in NQDA molecules. The DA group can form strong π-π stacking interactions with aromatic residues Y39/F94 of the αSyn fibril, while the DA molecule cannot. In addition, NQDA can form stronger cation-π interactions with residues K45/K80 than those of both NQ and DA molecules. Our results provide the molecular mechanism underlying the disaggregation of the αSyn fibril by NQDA and its better performance in fibril disruption than NQ, DA, and NQ+DA molecules, which offers new clues for the screening and development of promising drug candidates to treat PD.

Relaxation time scales of interfacial water upon fluid to ripple to gel phase transitions of bilayers

The slow relaxation of interface water (IW) across three primary phases of membranes is relevant to understand the influence of IW on membrane functions at supercooled conditions. To this objective, a total of ∼16.26μs all-atom molecular dynamics simulations of 1,2-dimyristoyl-sn-glycerol-3-phosphocholine lipid membranes are carried out. A supercooling-driven drastic slow-down in heterogeneity time scales of the IW is found at the fluid to the ripple to the gel phase transitions of the membranes. At both fluid-to-ripple-to-gel phase transitions, the IW undergoes two dynamic crossovers in Arrhenius behavior with the highest activation energy at the gel phase due to the highest number of hydrogen bonds. Interestingly, the Stokes-Einstein (SE) relation is conserved for the IW near all three phases of the membranes for the time scales derived from the diffusion exponents and the non-Gaussian parameters. However, the SE relation breaks for the time scale obtained from the self-intermediate scattering functions. The behavioral difference in different time scales is universal and found to be an intrinsic property of glass. The first dynamical transition in the α relaxation time of the IW is associated with an increase in the Gibbs energy of activation of hydrogen bond breaking with locally distorted tetrahedral structures, unlike the bulk water. Thus, our analyses unveil the nature of the relaxation time scales of the IW across membrane phase transitions in comparison with the bulk water. The results will be useful to understand the activities and survival of complex biomembranes under supercooled conditions in the future.

Light, Water, and Melatonin: The Synergistic Regulation of Phase Separation in Dementia

The swift rise in acceptance of molecular principles defining phase separation by a broad array of scientific disciplines is shadowed by increasing discoveries linking phase separation to pathological aggregations associated with numerous neurodegenerative disorders, including Alzheimer’s disease, that contribute to dementia. Phase separation is powered by multivalent macromolecular interactions. Importantly, the release of water molecules from protein hydration shells into bulk creates entropic gains that promote phase separation and the subsequent generation of insoluble cytotoxic aggregates that drive healthy brain cells into diseased states. Higher viscosity in interfacial waters and limited hydration in interiors of biomolecular condensates facilitate phase separation. Light, water, and melatonin constitute an ancient synergy that ensures adequate protein hydration to prevent aberrant phase separation. The 670 nm visible red wavelength found in sunlight and employed in photobiomodulation reduces interfacial and mitochondrial matrix viscosity to enhance ATP production via increasing ATP synthase motor efficiency. Melatonin is a potent antioxidant that lowers viscosity to increase ATP by scavenging excess reactive oxygen species and free radicals. Reduced viscosity by light and melatonin elevates the availability of free water molecules that allow melatonin to adopt favorable conformations that enhance intrinsic features, including binding interactions with adenosine that reinforces the adenosine moiety effect of ATP responsible for preventing water removal that causes hydrophobic collapse and aggregation in phase separation. Precise recalibration of interspecies melatonin dosages that account for differences in metabolic rates and bioavailability will ensure the efficacious reinstatement of the once-powerful ancient synergy between light, water, and melatonin in a modern world.

Forced Interactions: Ionic Polymers at Charged Surfactant Interfaces

Characterizing electrostatic interactions at heterogeneous interfaces is critical for developing a fundamental description of the dynamic processes at charged interfaces. Water-in-oil reverse micelles (RMs) offer a high degree of tunability across composition, polarity, and temperature, making them ideal systems for studying interactions at heterogeneous liquid-liquid interfaces. In the present study, we use a combination of ultrafast two-dimensional infrared spectroscopy and molecular dynamics (MD) simulations to determine the picosecond interfacial dynamics in RMs containing binary compositions of sorbitan monostearate and anionic or cationic cosurfactants, which are used to tune the ratio of charged to nonionic surfactants at the interface. The positively charged polyethylenimine (PEI) polymer is encapsulated within the RMs, and the carbonyl stretching mode of sorbitan monostearate reports on the interfacial hydrogen-bond populations and dynamics. The results show that hydrogen-bond populations are altered through the inclusion of both negatively and positively charged cosurfactants. Charged surfactants increase interfacial water penetration into the surfactant layer, and the surface localization of polymers decreases water penetration. Local hydrogen-bond dynamics undergo a slowdown with the inclusion of charged surfactants, and the encapsulation of polymers results in similar effects, irrespective of the charge.

Long-range communication between transmembrane- and nucleotide-binding domains does not depend on drug binding to mutant P-glycoprotein

In this study, the impact of four P-gp mutations (G185V, G830V, F978A and ΔF335) on drug-binding and efflux-related signal-transmission mechanism was comprehensively evaluated in the presence of ligands within the drug-binding pocket (DBP), experimentally related with changes in their drug efflux profiles. The severe repacking of the transmembrane helices (TMH), induced by mutations and exacerbated by the presence of ligands, indicates that P-gp is sensitive to perturbations in the transmembrane region. Alterations on drug-binding were also observed as a consequence of the TMH repacking, but were not always correlated with alterations on ligands binding mode and/or binding affinity. Finally, and although all P-gp variants holo systems showed considerable changes in the intracellular coupling helices/nucleotide-binding domain (ICH-NBD) interactions, they seem to be primarily induced by the mutation itself rather than by the presence of ligands within the DBP. The data further suggest that the changes in drug efflux experimentally reported are mostly related with changes on drug specificity rather than effects on signal-transmission mechanism. We also hypothesize that an increase in the drug-binding affinity may also be related with the decreased drug efflux, while minor changes in binding affinities are possibly related with the increased drug efflux observed in transfected cells.Communicated by Ramaswamy H. Sarma.

Effect of Charge State on the Equilibrium and Kinetic Properties of Mechanically Interlocked [5]Rotaxane: A Molecular Dynamics Study

Rotaxanes can exhibit stimuli-responsive behavior by allowing positional fluctuations of their rota groups in response to physiochemical conditions such as the changes in solution pH. However, ionic strength of the solution also affects the molecular conformation by altering the charge state of the entire molecule, coupling the stimuli-responsiveness of rotaxanes with their conformation. A molecular-scale investigation on a model system can allow the decoupling and identification of various effects and can greatly benefit applications of such molecular switches. By using atomistic molecular dynamics simulations, we study equilibrium and kinetics properties of various charge states of the [5]rotaxane, which is a supramolecular moiety with four rotaxanes bonded to a porphyrin core. We model various physiochemical charge states, each of which can be realized at various solution pH levels as well as several exotic charge distributions. By analyzing molecular configurations, hydrogen bonding, and energetics of single molecules in salt-free water and its polyrotaxanated network at the interface of water and chloroform, we demonstrate that charge-neutral and negatively charged molecules often tend to collapse in a way that they can expose their porphyrin core. Contrarily, positively charged moieties tend to take more extended molecular configurations blocking the core. Further, sudden changes in the charge states emulating the pH alterations in solution conditions lead to rapid, sub-10 ns level, changes in the molecular conformation of [5]rotaxane via shuttling motion of CB6 rings along axles. Finally, simulations of 2D [5]rotaxane network structures support our previous findings on a few nanometer-thick film formation at oil-water interfaces. Overall, our results suggest that rotaxane-based structures can exhibit a rich spectrum of molecular configurations and kinetics depending on the ionic strength of the solution.

Druggable sites identification in Streptococcus mutans VicRK system evaluated by catechols

Dental caries is a global public health problem, being the most common non-communicable disease. Streptococcus mutans, the causative agent of human cariogenic dental biofilms, produce glycosyltransferases (Gtfs) whose gene expression is modulated by the VicRK system, which makes them a promising target for dental biofilm inhibitor developments. Bioinformatics have playing a significant role in drug discovery programs mainly in novel hit identification. In this study, potential inhibitors against the S. mutans VicK system have been identified through Structure-based Virtual Screening performed between the VicK druggable sites followed byMolecular Dynamic simulations (MD) with binding affinity analysis by MM-PBSA approach. First, VicK protein was downloaded from PDB, and druggability analyses were performed by PockDrug and FTMap servers describing three interaction sites (S1, S2, and S3) that covered the most important domains for stability and activity. Next, a catechol virtual screening (n = 383) was performed on AutoDock4.2, and better-docked catechols showed strong binding affinity interaction through hydrogen bonding, hydrophobic interactions, and π-stacking with VicK auto kinase and phosphatase activity sites. Ligand efficiency indexes were also calculated (LE, LELP, LLE, and BEI) and showed optimal values. Furthermore, a 200 ns MD simulation run showed stability (RMSD and RMSF) and a high number of hydrogen bonds into peltatoside and maritimein, the two best VicK complexes. These results supported that catechols could potentially inhibit exopolysaccharides synthesis and be used in the biofilm management of new anti-cariogenic and antimicrobial agents.

Parameterization of a single H-bond in Orange Carotenoid Protein by atomic mutation reveals principles of evolutionary design of complex chemical photosystems

Introduction: Dissecting the intricate networks of covalent and non-covalent interactions that stabilize complex protein structures is notoriously difficult and requires subtle atomic-level exchanges to precisely affect local chemical functionality. The function of the Orange Carotenoid Protein (OCP), a light-driven photoswitch involved in cyanobacterial photoprotection, depends strongly on two H-bonds between the 4-ketolated xanthophyll cofactor and two highly conserved residues in the C-terminal domain (Trp288 and Tyr201). Method: By orthogonal translation, we replaced Trp288 in Synechocystis OCP with 3-benzothienyl-L-alanine (BTA), thereby exchanging the imino nitrogen for a sulphur atom. Results: Although the high-resolution (1.8 Å) crystal structure of the fully photoactive OCP-W288_BTA protein showed perfect isomorphism to the native structure, the spectroscopic and kinetic properties changed distinctly. We accurately parameterized the effects of the absence of a single H-bond on the spectroscopic and thermodynamic properties of OCP photoconversion and reveal general principles underlying the design of photoreceptors by natural evolution. Discussion: Such "molecular surgery" is superior over trial-and-error methods in hypothesis-driven research of complex chemical systems.

An Imbalance in the Force: The Need for Standardized Benchmarks for Molecular Simulation

Force fields (FFs) for molecular simulation have been under development for more than half a century. As with any predictive model, rigorous testing and comparisons of models critically depends on the availability of standardized data sets and benchmarks. While such benchmarks are rather common in the fields of quantum chemistry, this is not the case for empirical FFs. That is, few benchmarks are reused to evaluate FFs, and development teams rather use their own training and test sets. Here we present an overview of currently available tests and benchmarks for computational chemistry, focusing on organic compounds, including halogens and common ions, as FFs for these are the most common ones. We argue that many of the benchmark data sets from quantum chemistry can in fact be reused for evaluating FFs, but new gas phase data is still needed for compounds containing phosphorus and sulfur in different valence states. In addition, more nonequilibrium interaction energies and forces, as well as molecular properties such as electrostatic potentials around compounds, would be beneficial. For the condensed phases there is a large body of experimental data available, and tools to utilize these data in an automated fashion are under development. If FF developers, as well as researchers in artificial intelligence, would adopt a number of these data sets, it would become easier to compare the relative strengths and weaknesses of different models and to, eventually, restore the balance in the force.

In Silico discovery of aptamers with an enhanced library design strategy

With advances in force fields and algorithms, robust tools have been developed for molecular simulation of three-dimensional structures of nucleic acids and investigation of aptamer-target interactions. The traditional aptamer discovery technique, Systematic Evolution of Ligands by EXponential enrichment (SELEX), continues to suffer from high investment and low return, while in vitro screening by simulated SELEX remains a challenging task, where more accurate structural modeling and enhanced sampling limit the large-scale application of the method. Here, we proposed a modified aptamer enhanced library design strategy to facilitate the screening of target-binding aptamers. In this strategy, a comprehensive analysis of the original complexes and the target secondary structure were used to construct an enhanced initial library for screening. Our enhanced sequence library design strategy based on the target secondary structure explored a certain sequence space while ensuring the accuracy of the structural conformation and the calculation method. In an enhanced library of only a few dozen sequences, four sequences showed a similar or better binding free energy than the original aptamer, with consistently high binding stability over three rounds of multi-timescale simulations, ranging from - 30.27 to - 32.25 kcal/mol. Consequently, the enhanced library strategy based on the target secondary structure is shown to have very significant potential as a new aptamer design and optimization strategy.

Effect of Local Stress on Accurate Modeling of Bacterial Outer Membranes Using All-Atom Molecular Dynamics

Biological membranes are fundamental components of living organisms that play an undeniable role in their survival. Molecular dynamics (MD) serves as an essential computational tool for studying biomembranes on molecular and atomistic scales. The status quo of MD simulations of biomembranes studies a nanometer-sized membrane patch periodically extended under periodic boundary conditions (PBCs). In nature, membranes are usually composed of different lipids in their two layers (referred to as leaflets). This compositional asymmetry imposes a fixed ratio of lipid numbers between the two leaflets in a periodically constrained membrane, which needs to be set appropriately. The widely adopted methods of defining a leaflet lipid ratio suffer from the lack of control over the mechanical tension of each leaflet, which could significantly influence research findings. In this study, we investigate the role of membrane-building protocol and the resulting initial stress state on the interaction between small molecules and asymmetric membranes. We model the outer membrane of Pseudomonas aeruginosa bacteria using two different building protocols and probe their interactions with the Pseudomonas quinolone signal (PQS). Our results show that differential stress could shift the position of free energy minimum for the PQS molecule between the two leaflets of the asymmetric membrane. This work provides critical insights into the relationship between the initial per-leaflet tension and the spontaneous intercalation of PQS.

Investigating soil properties and vegetation parameters in different biochar-amended vegetated soil at large suction for application in bioengineered structures

Bioengineered structures, such as landfill cover, vegetated slopes or embankments, green roof and turf are comprised of soil and vegetation where vegetation imparts stability and protection through root reinforcement and hydrologic action. Soil in bioengineered structures often compacted and subjected to prolong drying due to irregular irrigation which necessitates the investigation of soil properties and vegetation growth in biochar-amended soil (BAS) under large suction range and it is scarce in the literature. In the present study, the effect of different biochar types on soil properties, and the vegetation growth in compacted soil and under large suction range was investigated for application in bioengineered structures. The results revealed that the biochar amendment decreased the dry density (5-32%) and increased the water retention capacity (θ_s by 15-104%, θ₁₅₀₀ by 82-445% and plant available water content (PAWC) by 22-55%), pH (28-77%) and cation exchange capacity (CEC, 16-723%) of the soil. Further, the vegetation growth i.e., vegetation density, dry root mass and shoot mass increased by 8-13%, 33-108% and 4-157% respectively after biochar amendment. The vegetation wilting was started at a higher suction (~ 900 kPa) relative to bare soil (800 kPa), the permanent wilting point (PWP) increased (by 3-35%) and the complete photosynthetic activity remained unchanged at a higher suction (1600 kPa) relative to bare soil (1050 kPa) after biochar amendment. Among the biochar types i.e., Sugarcane Bagasse biochar (SBB), Mesquite biochar (MB) and Water Hyacinth biochar (WHB) tested, the MB showed the best performance i.e., the suitable vegetation growth and health status. The improved water retention due to increased porosity, specific surface area (SSA) and presence of hydrophilic functional groups, and the higher pH, CEC and lower dry density in BAS attributed to the higher vegetation growth. The findings of the present study suggest the application of BAS in bioengineered structures.

JM-20, a Benzodiazepine-Dihydropyridine Hybrid Molecule, Inhibits the Formation of Alpha-Synuclein-Aggregated Species

Studies showed that JM-20, a benzodiazepine-dihydropyridine hybrid molecule, protects against rotenone and 6-hydroxydopamine neurotoxicity. However, its protective effects against cytotoxicity induced by endogenous neurotoxins involved in Parkinson's disease (PD) pathogenesis have never been investigated. In this study, we evaluated the ability of JM-20 to inhibit alpha-synuclein (aSyn) aggregation. We also evaluated the interactions of JM-20 with aSyn by molecular docking and molecular dynamics and assessed the protective effect of JM-20 against aminochrome cytotoxicity. We demonstrated that JM-20 induced the formation of heterogeneous amyloid fibrils, which were innocuous to primary cultures of mesencephalic cells. Moreover, JM-20 reduced the average size of aSyn positive inclusions in H4 cells transfected with SynT wild-type and synphilin-1-V5, but not in HEK cells transfected with synphilin-1-GFP. In silico studies showed the interaction between JM-20 and the aSyn-binding site. Additionally, we showed that JM-20 protects SH-SY5Y cells against aminochrome cytotoxicity. These results reinforce the potential of JM-20 as a neuroprotective compound for PD and suggest aSyn as a molecular target for JM-20.

Molecular Simulation-Guided Spectroscopy of Imidazolium-Based Ionic Liquids and Effects of Methylation on Ion-Cage and -Pair Dynamics

Classical molecular dynamics simulations were performed to assess an atomistic interpretation of the ion-probe structural interactions in two typical ionic liquids (ILs), 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf₂] and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BDimIm][NTf₂] through computational ultrafast spectroscopy. The nitrile stretching vibrations of the thiocyanate anion, [SCN]^-, serve as the local mode of the ultrafast system dynamics within the imidazolium-based ionic liquid environment. The wavelet transform of classical trajectories determines the time-varying fluctuating frequencies and the stretch spectral signatures of SCN^- in the normalized distribution. However, computational modeling of the two-dimensional (2D) spectra from the wavelet-derived vibrational frequencies yields time evolution of the local molecular structure along with the varied time-dependent dynamics of the spectral diffusion process. We calculated the frequency-frequency correlation functions (FFCFs), time correlations associated with the ion-pair and -cage dynamics, and mean square displacements as a function of time, depicting diffusive dynamics. The calculated results based on the pair correlation functions and the distribution of atomic density suggest that the hydrogen and methylated carbon at the two-position of the imidazolium ring of [BMIm] and [BDimIm] cations, respectively, strongly interact with the probe through the N of the thiocyanate anion rather than the S atom. The center-of-mass center-of-mass (COM-COM) cation-probe radial distribution functions (RDFs) in conjunction with the site-specific structural analysis further reveal well-structured interactions of the thiocyanate ion and [BMIm]⁺ cation rather than the [BDimIm] cation. In contrast, the anion-probe COM-COM RDFs depict weak interactive associations within the vibrational probe [SCN]^- and [NTf₂]^- ions. Methylation at the two-position of the imidazolium ring predicts slower structural reorganization and breaking and reformation dynamics of the ion pairs and cages within the ionic liquid framework.

Hydrogen Bond Thermodynamics in Aqueous Acid Solutions: A Combined DFT and Classical Force-Field Approach

The thermodynamics of hydrogen bonds in aqueous and acidic solutions significantly impacts the kinetics and thermodynamics of acid reaction chemistry. We utilize in this work a multiscale approach, combining density functional theory (DFT) with classical molecular dynamics (MD) to model hydrogen bond thermodynamics in an acidic solution. Using thermodynamic cycles, we split the solution phase free energy into its gas phase counterpart plus solvation free energies. We validate this DFT/MD approach by calculating the aqueous phase hydrogen bond free energy between two water molecules (H₂O-···-H₂O), the free energy to transform an H₃O⁺ cation into an H₅O₂⁺ cation, and the hydrogen bond free energy of protonated water clusters (H₃O⁺-···-H₂O and H₅O₂⁺-···-H₂O). The computed equilibrium hydrogen bond free energy of H₂O-···-H₂O is remarkably accurate, especially considering the large individual contributions to the thermodynamic cycle. Turning to cations, we find the ion to be more stable than H₃O⁺ by roughly 1-2 k_BT. This small free energy difference allows for thermal fluctuation between the two idealized motifs, consistent with spectroscopic and simulation studies. Lastly, hydrogen bonding free energies between either H⁺ cation and H₂O in solution were found to be stronger than between two H₂O, though much less so than in vacuum because of dielectric screening in solution. Altogether, our results suggest the DFT/MD approach is promising for application in modeling hydrogen bonding and proton transfer thermodynamics in condensed phases.