Barrierless reactions of C2 Criegee intermediates with H2SO4 and their implication to oligomers and new particle formation.

The formation of oligomeric hydrogen peroxide triggered by Criegee intermediate maybe contributes significantly to the formation and growth of secondary organic aerosol (SOA). However, to date, the reactivity of C2 Criegee intermediates (CH3CHOO) in areas contaminated with acidic gas remains poorly understood. Herein, high-level quantum chemical calculations and Born-Oppenheimer molecular dynamics (BOMD) simulations are used to explore the reaction of CH3CHOO and H2SO4 both in the gas phase and at the air-water interface. In the gas phase, the addition reaction of CH3CHOO with H2SO4 to generate CH3HC(OOH)OSO3H (HPES) is near-barrierless, regardless of the presence of water molecules. BOMD simulations show that the reaction at the air-water interface is even faster than that in the gas phase. Further calculations reveal that the HPES has a tendency to aggregate with sulfuric acids, ammonias, and water molecules to form stable clusters, meanwhile the oligomerization reaction of CH3CHOO with HPES in the gas phase is both thermochemically and kinetically favored. Also, it is noted that the interfacial HPES- ion can attract H2SO4, NH3, (COOH)2 and HNO3 for particle formation from the gas phase to the water surface. Thus, the results of this work not only elucidate the high atmospheric reactivity of C2 Criegee intermediates in polluted regions, but also deepen our understanding of the formation process of atmospheric SOA induced by Criegee intermediates.

Theoretical calculation on degradation mechanism of novel copolyesters under CALB enzyme.

Poly(butylene succinate-co-furandicarboxylate) (PBSF) and poly(butylene adipate-co-furandicarboxylate) (PBAF) are novel furandicarboxylic acid-based biodegradable copolyesters with great potential to replace fossil-derived terephthalic acid-based copolyesters such as poly(butylene succinate-co-terephthalate) (PBST) and poly(butylene adipate-co-terephthalate) (PBAT). In this study, quantum chemistry techniques after molecular dynamics simulations are employed to investigate the degradation mechanism of PBSF and PBAF catalyzed by Candida antarctica lipase B (CALB). Computational analysis indicates that the catalytic reaction follows a four-step mechanism resembling the ping-pong bibi mechanism, with the initial two steps being acylation reactions and the subsequent two being hydrolysis reactions. Notably, the first step of the hydrolysis is identified as the rate-determining step. Moreover, by introducing single-point mutations to expand the substrate entrance tunnel, the catalytic distance of the first acylation step decreases. Additionally, energy barrier of the rate-determining step is decreased in the PBSF system by site-directed mutations on key residues increasing hydrophobicity of the enzyme's active site. This study unprecedently show the substrate binding pocket and hydrophobicity of the enzyme's active site have the potential to be engineered to enhance the degradation of copolyesters catalyzed by CALB.

Structure and dynamics of cholesterol-mediated aquaporin-0 arrays and implications for lipid rafts.

Aquaporin-0 (AQP0) tetramers form square arrays in lens membranes through a yet unknown mechanism, but lens membranes are enriched in sphingomyelin and cholesterol. Here, we determined electron crystallographic structures of AQP0 in sphingomyelin/cholesterol membranes and performed molecular dynamics (MD) simulations to establish that the observed cholesterol positions represent those seen around an isolated AQP0 tetramer and that the AQP0 tetramer largely defines the location and orientation of most of its associated cholesterol molecules. At a high concentration, cholesterol increases the hydrophobic thickness of the annular lipid shell around AQP0 tetramers, which may thus cluster to mitigate the resulting hydrophobic mismatch. Moreover, neighboring AQP0 tetramers sandwich a cholesterol deep in the center of the membrane. MD simulations show that the association of two AQP0 tetramers is necessary to maintain the deep cholesterol in its position and that the deep cholesterol increases the force required to laterally detach two AQP0 tetramers, not only due to protein-protein contacts but also due to increased lipid-protein complementarity. Since each tetramer interacts with four such 'glue' cholesterols, avidity effects may stabilize larger arrays. The principles proposed to drive AQP0 array formation could also underlie protein clustering in lipid rafts.

Relative genotoxicity of polycyclic aromatic hydrocarbons inferred from free energy perturbation approaches.

Utilizing molecular dynamics and free energy perturbation, we examine the relative binding affinity of several covalent polycyclic aromatic hydrocarbon - DNA (PAH-DNA) adducts at the central adenine of NRAS codon-61, a mutational hotspot implicated in cancer risk. Several PAHs classified by the International Agency for Research on Cancer as probable, possible, or unclassifiable as to carcinogenicity are found to have greater binding affinity than the known carcinogen, benzo[a]pyrene (B[a]P). van der Waals interactions between the intercalated PAH and neighboring nucleobases, and minimal disruption of the DNA duplex drive increases in binding affinity. PAH-DNA adducts may be repaired by global genomic nucleotide excision repair (GG-NER), hence we also compute relative free energies of complexation of PAH-DNA adducts with RAD4-RAD23 (the yeast ortholog of human XPC-RAD23) which constitutes the recognition step in GG-NER. PAH-DNA adducts exhibiting the greatest DNA binding affinity also exhibit the least RAD4-RAD23 complexation affinity and are thus predicted to resist the GG-NER machinery, contributing to their genotoxic potential. In particular, the fjord region PAHs dibenzo[a,l]pyrene, benzo[g]chrysene, and benzo[c]phenanthrene are found to have greater binding affinity while having weaker RAD4-RAD23 complexation affinity than their respective bay region analogs B[a]P, chrysene, and phenanthrene. We also find that the bay region PAHs dibenzo[a,j]anthracene, dibenzo[a,c]anthracene, and dibenzo[a,h]anthracene exhibit greater binding affinity and weaker RAD4-RAD23 complexation affinity than B[a]P. Thus, the study of PAH genotoxicity likely needs to be substantially broadened, with implications for public policy and the health sciences. This approach can be broadly applied to assess factors contributing to the genotoxicity of other unclassified compounds.

Balancing thermodynamic stability, dynamics, and kinetics in phase separation of intrinsically disordered proteins.

Intrinsically disordered proteins (IDPs) are prevalent participants in liquid-liquid phase separation due to their inherent potential for promoting multivalent binding. Understanding the underlying mechanisms of phase separation is challenging, as phase separation is a complex process, involving numerous molecules and various types of interactions. Here, we used a simplified coarse-grained model of IDPs to investigate the thermodynamic stability of the dense phase, conformational properties of IDPs, chain dynamics, and kinetic rates of forming condensates. We focused on the IDP system, in which the oppositely charged IDPs are maximally segregated, inherently possessing a high propensity for phase separation. By varying interaction strengths, salt concentrations, and temperatures, we observed that IDPs in the dense phase exhibited highly conserved conformational characteristics, which are more extended than those in the dilute phase. Although the chain motions and global conformational dynamics of IDPs in the condensates are slow due to the high viscosity, local chain flexibility at the short timescales is largely preserved with respect to that at the free state. Strikingly, we observed a non-monotonic relationship between interaction strengths and kinetic rates for forming condensates. As strong interactions of IDPs result in high stable condensates, our results suggest that the thermodynamics and kinetics of phase separation are decoupled and optimized by the speed-stability balance through underlying molecular interactions. Our findings contribute to the molecular-level understanding of phase separation and offer valuable insights into the developments of engineering strategies for precise regulation of biomolecular condensates.

Metadynamics simulations reveal mechanisms of Na+ and Ca2+ transport in two open states of the channelrhodopsin chimera, C1C2.

Cation conducting channelrhodopsins (ChRs) are a popular tool used in optogenetics to control the activity of excitable cells and tissues using light. ChRs with altered ion selectivity are in high demand for use in different cell types and for other specialized applications. However, a detailed mechanism of ion permeation in ChRs is not fully resolved. Here, we use complementary experimental and computational methods to uncover the mechanisms of cation transport and valence selectivity through the channelrhodopsin chimera, C1C2, in the high- and low-conducting open states. Electrophysiology measurements identified a single-residue substitution within the central gate, N297D, that increased Ca2+ permeability vs. Na+ by nearly two-fold at peak current, but less so at stationary current. We then developed molecular models of dimeric wild-type C1C2 and N297D mutant channels in both open states and calculated the PMF profiles for Na+ and Ca2+ permeation through each protein using well-tempered/multiple-walker metadynamics. Results of these studies agree well with experimental measurements and demonstrate that the pore entrance on the extracellular side differs from original predictions and is actually located in a gap between helices I and II. Cation transport occurs via a relay mechanism where cations are passed between flexible carboxylate sidechains lining the full length of the pore by sidechain swinging, like a monkey swinging on vines. In the mutant channel, residue D297 enhances Ca2+ permeability by mediating the handoff between the central and cytosolic binding sites via direct coordination and sidechain swinging. We also found that altered cation binding affinities at both the extracellular entrance and central binding sites underly the distinct transport properties of the low-conducting open state. This work significantly advances our understanding of ion selectivity and permeation in cation channelrhodopsins and provides the insights needed for successful development of new ion-selective optogenetic tools.

Exploration of alpha-glucosidase inhibitors: A comprehensive in silico approach targeting a large set of triazole derivatives.

BACKGROUND: The increasing prevalence of diabetes and the side effects associated with current medications necessitate the development of novel candidate drugs targeting alpha-glucosidase as a potential treatment option. METHODS: This study employed computer-aided drug design techniques to identify potential alpha-glucosidase inhibitors from the PubChem database. Molecular docking was used to evaluate 81,197 compounds, narrowing the set for further analysis and providing insights into ligand-target interactions. An ADMET study assessed the pharmacokinetic properties of these compounds, including absorption, distribution, metabolism, excretion, and toxicity. Molecular dynamics simulations validated the docking results. RESULTS: 9 compounds were identified as potential candidate drugs based on their ability to form stable complexes with alpha-glucosidase and their favorable pharmacokinetic profiles, three of these compounds were subjected to the molecular dynamics, which showed stability throughout the entire 100 ns simulation. CONCLUSION: These findings suggest promising new alpha-glucosidase inhibitors for diabetes treatment. Further validation through in vitro and in vivo studies is recommended to confirm their efficacy and safety.

Syk protein inhibitors treatment for the allergic symptoms associated with hyper immunoglobulin E syndromes: A focused on a computational approach.

Background: Mutations in the Spleen tyrosine kinase (Syk) protein have significant implications for its function and response to treatments. Understanding these mutations and identifying new inhibitors can lead to more effective therapies for conditions like autosomal dominant hyper-IgE syndrome (AD-HIES) and related immunological disorders. Objective: To investigate the impact of mutations in the Syk protein on its function and response to reference treatments, and to explore new inhibitors tailored to the mutational profile of Syk. Methods: We collected and analyzed mutations affecting the Syk protein to assess their functional impact. We screened 94 deleterious mutations in the kinase domain using molecular docking techniques. A library of 997 compounds with potential inhibitory activity against Syk was filtered based on Lipinski and Veber rules and toxicity assessments. We evaluated the binding affinity of reference inhibitors and 14 eligible compounds against wild-type and mutant Syk proteins. Molecular dynamics simulations were conducted to evaluate the interaction of Syk protein complexes with the reference inhibitor and potential candidate inhibitors. Results: Among the analyzed mutations, 60.5% were identified as deleterious, underscoring their potential impact on cellular processes. Virtual screening identified three potential inhibitors (IDs: 118558008, 118558000, and 118558092) with greater therapeutic potential than reference treatments, meeting all criteria and exhibiting lower IC50 values. Ligand 1 (ID: 118558000) demonstrated the most stable binding, favorable compactness, and extensive interaction with solvents. A 3D pharmacophore model was constructed, identifying structural features common to these inhibitors. Conclusion: This study found that 60.5% of reported mutations affecting the Syk protein are deleterious. Virtual screening revealed three top potential inhibitors, with ligand 1 (ID: 118558000) showing the most stable binding and favorable interactions. These inhibitors hold promise for more effective therapies targeting Syk-mediated signaling pathways. The pharmacophore model provides valuable insights for developing targeted therapies for AD-HIES and related disorders, offering hope for patients suffering from Hyper IgE syndrome with allergic symptoms.

Divergent mechanisms of steroid inhibition in the human ρ1 GABAA receptor.

ρ-type γ-aminobutyric acid-A (GABAA) receptors are widely distributed in the retina and brain, and are potential drug targets for the treatment of visual, sleep and cognitive disorders. Endogenous neuroactive steroids including ß-estradiol and pregnenolone sulfate negatively modulate the function of ρ1 GABAA receptors, but their inhibitory mechanisms are not clear. By combining five cryo-EM structures with electrophysiology and molecular dynamics simulations, we characterize binding sites and negative modulation mechanisms of ß-estradiol and pregnenolone sulfate at the human ρ1 GABAA receptor. ß-estradiol binds in a pocket at the interface between extracellular and transmembrane domains, apparently specific to the ρ subfamily, and disturbs allosteric conformational transitions linking GABA binding to pore opening. In contrast, pregnenolone sulfate binds inside the pore to block ion permeation, with a preference for activated structures. These results illuminate contrasting mechanisms of ρ1 inhibition by two different neuroactive steroids, with potential implications for subtype-specific gating and pharmacological design.

Integrating bioinformatics and ferroptosis to reveal the protective mechanism of Astragaloside IV on chronic heart failure rats.

Ferroptosis is an important pathological mechanism of chronic heart failure (CHF). This study aimed to investigate the protective mechanism of Astragaloside IV (AS-IV) on CHF rats by integrating bioinformatics and ferroptosis. CHF-related targets and ferroptosis-related targets were collected. After the intersection, the common targets were obtained. The PPI network of the common targets was constructed, and topological analysis of the network was carried out. The target with the highest topological parameter values was selected as the key target. The key target p53 was obtained through bioinformatics analysis, and its molecular docking model with AS-IV was obtained, as well as molecular dynamics simulation analysis. The rat models of CHF after myocardial infarction were established by ligation of left coronary artery and treated with AS-IV for 4 weeks. AS-IV treatment significantly improved cardiac function in CHF rats, improved cardiomyocyte morphology and myocardial fibrosis, reduced mitochondrial damage, decreased myocardial MDA and Fe2+ content, increased GSH content, inhibited the expression of p53 and p-p53, and up-regulated the expression of SLC7A11 and GPX4. In conclusion, AS-IV improved cardiac function in CHF rats, presumably by regulating p53/SLC7A11/GPX4 signaling pathway and inhibiting myocardial ferroptosis.

Exploring the molecular mechanism of ginseng against anthracycline-induced cardiotoxicity based on network pharmacology, molecular docking and molecular dynamics simulation.

BACKGROUND: Previous clinical and basic studies have revealed that ginseng might have cardioprotective properties against anthracycline-induced cardiotoxicity (AIC). However, the underlying mechanism of ginseng action against AIC remains insufficiently understood. The aim of this study was to explore the related targets and pathways of ginseng against AIC using network pharmacology, molecular docking, cellular thermal shift assay (CETSA) and molecular dynamics (MD) simulations. RESULTS: Fourteen drug-disease common targets were identified. Enrichment analysis showed that the AGE-RAGE in diabetic complications, fluid shear stress and atherosclerosis, and TNF signaling pathway were potentially involved in the action of ginseng against AIC. Molecular docking demonstrated that the core components including Kaempferol, beta-Sitosterol, and Fumarine had notable binding activity with the three core targets CCNA2, STAT1, and ICAM1. Furthermore, the stable complex of STAT1 and Kaempferol with favorable affinity was further confirmed by CETSA and MD simulation. CONCLUSIONS: This study suggested that ginseng might exert their protective effects against AIC through the derived effector compounds beta-Sitosterol, Kaempferol and Fumarine by targeting CCNA2, STAT1, and ICAM1, and modulating AGE-RAGE in diabetic complications, fluid shear stress and atherosclerosis, and TNF signaling pathways.

Integrative experimental and computational analysis of the impact of KGM's polymerization degree on wheat starch's pasting and retrogradation characteristics.

This study investigated the influence of Konjac Glucomannan (KGM) with varying degrees of polymerization (DKGMx) on the gelatinization and retrogradation characteristics of wheat starch, providing new insights into starch-polysaccharide interactions. This research uniquely focuses on the effects of DKGMx, utilizing multidisciplinary approaches including Rapid Visco Analysis (RVA), Differential Scanning Calorimetry (DSC), rheological testing, Low-Field Nuclear Magnetic Resonance (LF-NMR), and molecular simulations to assess the effects of DKGMx on gelatinization temperature, viscosity, structural changes post-retrogradation, and molecular interactions. Our findings revealed that higher degrees of polymerization (DP) of DKGMx significantly enhanced starch's pasting viscosity and stability, whereas lower DP reduced viscosity and interfered with retrogradation. High DP DKGMx promoted retrogradation by modifying moisture distribution. Molecular simulations revealed the interplay between low DP DKGMx and starch molecules. These interactions, characterized by increased hydrogen bonds and tighter binding to more starch chains, inhibited starch molecular rearrangement. Specifically, low DP DKGMx established a dense hydrogen bond network with starch, significantly restricting molecular mobility and rearrangement. This study provides new insights into the role of the DP of DKGMx in modulating wheat starch's properties, offering valuable implications for the functional improvement of starch-based foods and advancing starch science.

Arabinan branches in the RG-I region of citrus pectin aid acid-induced gelation.

Gelation is a critical property of citrus pectin. However, the roles played by neutral sugar side-chains on acid-induced pectin gelation remain poorly understood. Herein, galactan- or/and arabinan-eliminated pectins (P-G, P-A, and P-AG) were used to investigate the effects of side-chains on gelation. The gel hardness values of citrus pectin, P-G, P-A, and P-AG were 42.6, 39.9, 5.3, and 2.1 g, respectively, suggesting that arabinan contributed more to gelation than galactan. We next found that arabinan branches promoted pectin chain entanglement more effectively than arabinan backbones. Destabilizer addition experiments showed that hydrogen bonding, electrostatic interaction, and hydrophobic interaction were the main forces affecting pectin gel networks and strength, which was further validated by molecular dynamic simulations. The total number of hydrogen bonds between the arabinan branches and galactan/HG (65.7) was significantly higher than that between the arabinan backbones and galactan/HG (39.1), indicating that arabinan branches predominated in terms of such interactions. This study thus elucidated the roles played by neutral-sugar side-chains, especially the arabinan branches of acid-induced pectin gels, in term of enhancing high-methoxyl pectin gelation, and offers novel insights into the structure-gelling relationships of citrus pectin.

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.

Robustaflavone as a novel scaffold for inhibitors of native and auto-proteolysed human neutrophil elastase.

Human neutrophil elastase (HNE) plays a key role in initiating inflammation in the cardiopulmonary and systemic contexts. Pathological auto-proteolysed two-chain (tc) HNE exhibits reduced binding affinity with inhibitors. Using AutoDock Vina v1.2.0, 66 flavonoid inhibitors, sivelestat and alvelestat were docked with single-chain (sc) HNE and tcHNE. Schrodinger PHASE v13.4.132 was used to generate a 3D-QSAR model. Molecular dynamics (MD) simulations were conducted with AMBER v18. The 3D-QSAR model for flavonoids with scHNE showed r2 = 0.95 and q2 = 0.91. High-activity compounds had hydrophobic A/A2 and C/C2 rings in the S1 subsite, with hydrogen bond donors at C5 and C7 positions of the A/A2 ring, and the C4' position of the B/B1 ring. All flavonoids except robustaflavone occupied the S1'-S2' subsites of tcHNE with decreased AutoDock binding affinities. During MD simulations, robustaflavone remained highly stable with both HNE forms. Principal Component Analysis suggested that robustaflavone binding induced structural stability in both HNE forms. Cluster analysis and free energy landscape plots showed that robustaflavone remained within the sc and tcHNE binding site throughout the 100 ns MD simulation. The robustaflavone scaffold likely inhibits both tcHNE and scHNE. It is potentially superior to sivelestat and alvelestat and can aid in developing therapeutics targeting both forms of HNE.

Mechanism of PARP1 Elongation Reaction Revealed by Molecular Modeling.

Poly(ADP-ribose) polymerase 1 (PARP1) plays a major role in the DNA damage repair and transcriptional regulation, and is targeted by a number of clinical inhibitors. Despite this, catalytic mechanism of PARP1 remains largely underexplored because of the complex substrate/product structure. Using molecular modeling and metadynamics simulations we have described in detail elongation of poly(ADP-ribose) chain in the PARP1 active site. It was shown that elongation reaction proceeds via the SN1-like mechanism involving formation of the intermediate furanosyl oxocarbenium ion. Intriguingly, nucleophilic 2'A-OH group of the acceptor substrate can be activated by the general base Glu988 not directly but through the proton relay system including the adjacent 3'A-OH group.

Interaction Analysis by Fragment Molecular Orbital Method for Drug Discovery Research.

The use of computational methods in drug discovery research has increased substantially in recent years. Computational chemistry techniques, such as quantum chemical calculations and molecular dynamics simulations, continue to be widely used. In this review, we focused on drug discovery-related studies that employ fragment molecular orbital methods. Furthermore, we focused on inhibitor discovery, protein-protein interaction analysis, including antigen-antibody interaction analysis, and integration with molecular dynamics simulations.

Sedative Effects of Daidzin, Possibly Through the GABAA Receptor Interaction Pathway: In Vivo Approach with Molecular Dynamic Simulations.

The soy isoflavone daidzin (DZN) has been considered a hopeful bioactive compound having diverse biological activities, including anxiolytic, memory-enhancing, and antiepileptic effects, in experimental animals. However, its sedative and hypnotic effects are yet to be discovered. This study aimed to evaluate its sedative/hypnotic effect on Swiss mice. Additionally, in silico studies were also performed to see the possible molecular mechanisms behind the tested neurological effect. For this, male Swiss albino mice were treated with DZN (5, 10, or 20 mg/kg) intraperitoneally (i.p.) with or without the standard GABAergic medication diazepam (DZP) and/or flumazenil (FLU) and checked for the onset and duration of sleeping time using thiopental sodium-induced as well as DZP-induced sleeping tests. A molecular docking study was also performed to check its interaction capacity with the α1 and ß2 subunits of the GABAA receptor. Findings suggest that DZN dose-dependently and significantly reduced the latency while increasing the duration of sleep in animals. In combination therapy, DZN shows synergistic effects with the DZP-2 and DZP-2 + FLU-0.01 groups, resulting in significantly (p < 0.05) reduced latency and increased sleep duration. Further, molecular docking studies demonstrate that DZN has a strong binding affinity of - 7.2 kcal/mol, which is closer to the standard ligand DZP (- 8.3 kcal/mol) against the GABAA (6X3X) receptor. Molecular dynamic simulations indicated stability and similar binding locations for DZP and DZN with 6X3X. In conclusion, DZN shows sedative effects on Swiss mice, possibly through the GABAA receptor interaction pathway.

Biphenylsulfonamides as effective MMP-2 inhibitors with promising antileukemic efficacy: Synthesis, in vitro biological evaluation, molecular docking, and MD simulation analysis.

Overexpression of matrix metalloproteinase-2 (MMP-2) possesses a correlation with leukemia especially chronic myeloid leukemia (CML). However, no such MMP-2 inhibitor has come out in the market to date for treating leukemia. In this study, synthesis, biological evaluation, and molecular modeling studies of a set of biphenylsulfonamide derivatives as promising MMP-2 inhibitors were performed, focusing on their potential applications as antileukemic therapeutics. Compounds DH-18 and DH-19 exerted the most effective MMP-2 inhibition (IC50 of 139.45 nM and 115.16 nM, respectively) with potent antileukemic efficacy against the CML cell line K562 (IC50 of 0.338 µM and 0.398 µM, respectively). The lead molecules DH-18 and DH-19 reduced the MMP-2 expression by 21.3% and 17.8%, respectively with effective apoptotic induction (45.4% and 39.8%, respectively) in the K562 cell line. Moreover, both these compounds significantly arrested different phases of the cell cycle. Again, both these molecules depicted promising antiangiogenic efficacy in the ACHN cell line. Nevertheless, the molecular docking and molecular dynamics (MD) simulation studies revealed that DH-18 formed strong bidentate chelation with the catalytic Zn2+ ion through the hydroxamate zinc binding group (ZBG). Apart from that, the MD simulation study also disclosed stable binding interactions of DH-18 and MMP-2 along with crucial interactions with active site amino acid residues namely His120, Glu121, His124, His130, Pro140, and Tyr142. In a nutshell, this study highlighted the importance of biphenylsulfonamide-based novel and promising MMP-2 inhibitors to open up a new avenue for potential therapy against CML.

Sensing the structural and conformational properties of single-stranded nucleic acids using electrometry and molecular simulations.

Inferring the 3D structure and conformation of disordered biomolecules, e.g., single stranded nucleic acids (ssNAs), remains challenging due to their conformational heterogeneity in solution. Here, we use escape-time electrometry (ETe) to measure with sub elementary-charge precision the effective electrical charge in solution of short to medium chain length ssNAs in the range of 5-60 bases. We compare measurements of molecular effective charge with theoretically calculated values for simulated molecular conformations obtained from Molecular Dynamics simulations using a variety of forcefield descriptions. We demonstrate that the measured effective charge captures subtle differences in molecular structure in various nucleic acid homopolymers of identical length, and also that the experimental measurements can find agreement with computed values derived from coarse-grained molecular structure descriptions such as oxDNA, as well next generation ssNA force fields. We further show that comparing the measured effective charge with calculations for a rigid, charged rod-the simplest model of a nucleic acid-yields estimates of molecular structural dimensions such as linear charge spacings that capture molecular structural trends observed using high resolution structural analysis methods such as X-ray scattering. By sensitively probing the effective charge of a molecule, electrometry provides a powerful dimension supporting inferences of molecular structural and conformational properties, as well as the validation of biomolecular structural models. The overall approach holds promise for a high throughput, microscopy-based biomolecular analytical approach offering rapid screening and inference of molecular 3D conformation, and operating at the single molecule level in solution.