CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields

Effect of Cations on ATP Binding to the N-domain of Na+, K+-ATPase

The nucleotide-binding domain (N-domain) of the Na+, K+-ATPase (NKA) is physicochemically characterized by a high content of Glu and Asp residues, resulting in a low isoelectric point (pI = 5.0). Acidic proteins are known to interact with cations. The analysis in silico revealed potential cation interaction sites in the NKA N-domain structure. The interaction with cations was tested in vitro by using a recombinant NKA N-domain. The N-domain contains two Trp residues at the protein surface, as determined by acrylamide-mediated fluorescence quenching, that are useful for structural studies through fluorescence changes. Intrinsic fluorescence of the N-domain was decreased by the presence of cations (Na+, K+, Ca2+) indicating an effect on the protein structure. ATP binding also decreased the N-domain intrinsic fluorescence, which allowed nucleotide affinity determination. In the presence of cations, the N-domain affinity for ATP was increased. Molecular docking of fluorescein isothiocyanate (FITC) with the N-domain showed two binding modes with the isothiocyanate group located 5-6 Å close to Lys485 and Lys506 in the nucleotide-binding site. The presence of ATP prevented the FITC covalent labeling of the N-domain demonstrating the competitive behavior for the binding site. It is proposed that cations interact with the N-domain structure and thereby modulate nucleotide (ATP) affinity and possibly affecting NKA catalysis.

Forecasting the Pharmacological Mechanisms of Plumbago zeylanica and Solanum xanthocarpum in Diabetic Retinopathy Treatment: A Network Pharmacology, Molecular Docking, and Molecular Dynamics Simulation Study

(1) Background: Diabetic retinopathy (DR) is a major complication of diabetes, marked by abnormal angiogenesis, microaneurysms, and retinal hemorrhages. Traditional Ayurvedic medicine advocates multi-target strategies for DR management. However, the mechanisms by which Solanum xanthocarpum (SX) and Plumbago zeylanica (PZ) exert therapeutic effects are not well understood; (2) Methods: To investigate these mechanisms, we employed network pharmacology and molecular docking techniques. Phytochemicals from SX and PZ were identified using the IMPPAT database and Swiss Target Prediction tool. DR-related protein targets were sourced from the GeneCards database, and common targets were identified through Venn diagram analysis. STRING and Cytoscape were used to construct and analyze protein-protein interaction networks. Pathway enrichment was performed with Gene Ontology and KEGG databases; (3) Results: We identified 28 active phytoconstituents, targeting proteins such as EGFR, SRC, STAT3, AKT1, and HSP90AA1. Molecular docking and dynamics simulations confirmed the strong binding affinities of these compounds to their targets; (4) Conclusions: The study highlights the multi-target activity of SX and PZ, particularly in pathways related to EGFR tyrosine kinase inhibitor resistance and PI3K-AKT signaling. These findings provide valuable insights into their therapeutic potential for DR, suggesting the effective modulation of key molecular pathways involved in the disease.

Benzenesulfonohydrazide-tethered non-fused and fused heterocycles as potential anti-mycobacterial agents targeting enoyl acyl carrier protein reductase (InhA) with antibiofilm activity

Because resistant variants of the disease are always emerging, tuberculosis is a global issue that affects economies. New antitubercular medications should be developed, and this can be done by inhibiting druggable targets. Enoyl acyl carrier protein (ACP) reductase (InhA) is a crucial enzyme for the survival of Mycobacterium tuberculosis (MTB). In this study, a series of small molecules based on non-fused and fused heterocycles (pyridine, coumarin, quinoline, and indole) tethered with benzenesulfonohydrazide were prepared via an aza-Michael reaction exploiting a one-pot synthesis approach. The synthesized molecules (2-7) were evaluated for their activity against tubercle bacilli. Three analogues showed efficacy against tuberculosis, with compound 7 demonstrating a MIC value as low as 8 μg mL-1. Consequently, compounds 3 and 7 successfully hindered the growth of mycobacteria in human monocyte-derived macrophages (MDMs), demonstrating their ability to penetrate human professional phagocytes. Furthermore, they restricted the ability of mycobacteria to produce biofilms. In addition, the inhibitory effects of compounds 3 and 7 against InhA were assessed. Compound 7 exhibited the best efficacy, with an IC50 value of 0.91 μM. The findings showed that the sulfonamide and methyl ester's carbonyl functionalities were engaged in hydrogen bonding with the essential Ile194 and Tyr158 residues, respectively.

β-Lactoglobulin Enhances Clay and Activated Carbon Binding and Protection Properties for Cadmium and Lead

The removal of heavy metals from wastewater remains a challenge due to the limitations of current remediation methods. This study aims to develop multicomponent composites as inexpensive and environmentally friendly sorbents with enhanced capture of cadmium (Cd) and lead (Pb). The composites are based on calcium montmorillonite (CM) and activated carbon (AC) because of their proven effectiveness as sorbents for diverse toxins in environmental settings. In this study, we used a combination of computational and experimental methods to delineate that β-lactoglobulin enhances CM and AC binding and protection properties for Cd and Pb. Modeling and molecular dynamics simulations investigated the formation of material systems formed by CM and AC in complex with β-lactoglobulin and predicted their capacity to bind heavy metal ions at neutral pH conditions. Our simulations suggest that the enhanced binding properties of the material systems are attributed to the presence of several binding pockets formed by β-lactoglobulin for the two heavy metal ions. At neutral pH conditions, divalent Cd and Pb shared comparable binding propensities in all material systems, with the former being consistently higher than the latter. To validate the interactions depicted in simulations, two ecotoxicological models (L. minor and H. vulgaris) were exposed to Cd, Pb, and a mixture of the two. The inclusion of CM-lactoglobulin (β-lactoglobulin amended CM) and AC-lactoglobulin (β-lactoglobulin amended AC) at 0.05-0.2% efficiently and dose-dependently reduced the severe toxicity of metals and increased the growth parameters. This high efficacy of protection shown in the ecotoxicological models may result from the numerous possible interaction pockets of the β-lactoglobulin-amended materials depicted in simulations. The ecotoxicological models support the agreement with computations. This study serves as a proof of concept on how computations in tandem with experiments can be used in the design of multicomponent clay- and carbon-based sorbent amended systems with augmented functionalities for particular toxins.

Cheminformatics-enhanced discovery of therapeutic agents targeting isocitrate lyase in Mycobacterium tuberculosis infections

Tuberculosis (TB) is a global health challenge; therefore, there is an urgent requirement to develop a novel and more effective anti-TB therapeutic. This study targeted the isocitrate lyase (ICL) protein due to its pivotal role in the pathogenicity of Mycobacterium tuberculosis (Mtb). Virtual screening of 8752 bioactive compounds used an ML-based QSAR model and molecular docking. ADMET testing was performed on the top three hits to identify the compound most closely mimicking a drug molecule. The top hits, 648 and 2785758, showed high binding affinity towards ICL with -7.3 and -7 kcal/mol, comparable to the control. These molecules also showed strong binding with the residue Asp108, which plays a vital role in ICL activity. Molecular dynamics simulations showed stability for 648 and 2785758, comparable to the control compound used in this study. It was found that 648 bound to the protein maintained the RMSD constant and consistent at 0.3 nm for a complete 100 ns simulation. 2785758 showed a comparable RMSD trend to the control. Both 648 and 2785758 showed high RMSF for critical residue Asp108. Further, PCA and FEL confirmed the formation of a stable complex. MM/GBSA estimations of binding free energy indicated that compounds 648 had an elevated level of stability (ΔGTOTAL = -28.11 kcal/mol) and 2785758 (ΔGTOTAL = -21.05 kcal/mol). This study suggests that compounds 648 and 2785758 can potentially affect the activity of ICL, leading to its inactivation and ultimately preventing the progression of tuberculosis.Communicated by Ramaswamy H. Sarma.

A Patent-Pending Ointment Containing Extracts of Five Different Plants Showed Antinociceptive and Anti-Inflammatory Mechanisms in Preclinical Studies

Background/Objectives: The antinociceptive and anti-inflammatory effects of a patent-pending ointment containing plant extracts from Eucalyptus globulus, Curcuma longa, Hamamelis virginiana, Echinacea purpurea, and Zingiber officinale were evaluated. Methods: Plant extracts were chemically characterized by gas chromatography-mass spectroscopy. The antinociceptive activity of the ointment was assessed using the hot plate, tail flick, and formalin tests, whereas the anti-inflammatory activity was measured using the acute and chronic TPA-induced ear edema tests. Mechanisms of action were evaluated using inhibitors from signaling pathways related to pain response and by using histological analysis and assessing the expression and activity of pro-inflammatory mediators. Results: The ointment showed antinociceptive and anti-inflammatory effects like those observed with diclofenac gel (1.16% v/v) and ketoprofen gel (2.5% v/v). The antinociceptive actions of the ointment are mediated by the possible participation of the opiodergic system and the nitric oxide pathway. The anti-inflammatory response was characterized by a decrease in myeloperoxidase (MPO) activity and by a reduction in ear swelling and monocyte infiltration in the acute inflammation model. In the chronic model, the mechanism of action relied on a decrease in pro-inflammatory mediators such as COX-2, IL-1β, TNF-α, and MPO. An in-silico study with myristic acid, one of the compounds identified in the ointment's plant mixture, corroborated the in vivo results. Conclusions: The ointment showed antinociceptive activities mediated by the decrease in COX-2 and NO levels, and anti-inflammatory activity due to the reduction in IL-1β and TNFα levels, a reduction in MPO activity, and a decrease in NF-κB and COX-2 expression.

Human cardiac β-myosin powerstroke energetics: Thin filament, Pi displacement, and mutation effects

The powerstroke of human cardiac β-myosin is an important stage of the cross-bridge cycle that generates force for muscle contraction. However, the starting structure of this process has never been resolved, and the relative timing of the powerstroke and inorganic phosphate (Pi) release is still controversial. In this study, we generated an atomistic model of myosin on the thin filament and utilized metadynamics simulations to predict the absent starting structure of the powerstroke. We demonstrated that the displacement of Pi from the active site during the powerstroke is likely necessary, reducing the energy barrier of the conformation change. The effects of the presence of the thin filament, the hypertrophic cardiomyopathy mutation R712L, and the binding of mavacamten on the powerstroke process were also investigated.

Peptide Bundlemer Networks or Lattices: Controlling Cross-Linking and Self-Assembly Using Protein-like Display of Chemistry

Coiled-coil 'bundlemer' peptides were selectively modified with allyloxycarbonyl (alloc)-protected lysine, a non-natural amino acid containing an alkene on its side chain. The specific display of this alkene from the coiled-coil surface with protein-like specificity enabled this residue to be used as a covalent linkage for creating peptide networks with controllable properties or as a physical linkage for the self-assembly of bundlemers into unexpected, intricate lattices driven by the hydrophobic nature of the side chain. For network formation, peptides were modified with both alloc-protected lysine and cysteine amino acids for solution assembly into solvent-swollen films and subsequent covalent cross-linking via thiol-ene photo click reactions. The degree of network cross-linking, as determined by rheometry, was finely tuned by varying the specific spatial display of reactive groups on the bundlemer building block particles, transitioning between intrabundle and interbundle cross-linking. The designed display of alloc groups from the center of the bundlemer building block also prompted particle self-assembly into an unexpected intricate lattice with a porous morphology. The lattices were studied in a variety of solution conditions using transmission electron microscopy, cryotransmission electron microscopy, and small-angle X-ray scattering. The approximate particle arrangement in the lattice was determined by using coarse-grained modeling and machine learning optimization techniques along with experimental methods. The proposed truss-like face-centered cubic packing of the alloc-functionalized bundlemers agrees well with the experimental results.

Integrative residue-intuitive machine learning and MD Approach to Unveil Allosteric Site and Mechanism for β2AR

Allosteric drugs offer a new avenue for modern drug design. However, the identification of cryptic allosteric sites presents a formidable challenge. Following the allostery nature of residue-driven conformation transition, we propose a state-of-the-art computational pipeline by developing a residue-intuitive hybrid machine learning (RHML) model coupled with molecular dynamics (MD) simulation, through which we can efficiently identify the allosteric site and allosteric modulator as well as reveal their regulation mechanism. For the clinical target β2-adrenoceptor (β2AR), we discover an additional allosteric site located around residues D792.50, F2826.44, N3187.45 and S3197.46 and one putative allosteric modulator ZINC5042. Using Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and protein structure network (PSN), the allosteric potency and regulation mechanism are probed to further improve identification accuracy. Benefiting from sufficient computational evidence, the experimental assays then validate our predicted allosteric site, negative allosteric potency and regulation pathway, showcasing the effectiveness of the identification pipeline in practice. We expect that it will be applicable to other target proteins.

How do Photoswitchable Lipids Influence the Intercalation of Anticancer Drug in Lipid Membrane? Investigation using Molecular Dynamics Simulation

Photoswitchable lipids, particularly azobenzene-derivatized phosphatidylcholine (azoPC) lipids, offer a unique mechanism for reversible modification of membrane properties upon exposure to ultraviolet (UV) radiation. Through all-atom molecular dynamics simulations, we explore how UV irradiation-induced trans-to-cis photoisomerization (TCPI) of AzoPC lipid influences the structure and dynamics of a lipid membrane, composed of dipalmitoylphosphatidylcholine (DPPC) and cholesterol with similar composition to that of the DOXIL®. Structural and dynamical analyses of two states of the membrane, 'dark' state (containing cis-azoPC lipid) and 'bright' state (containing 85 % cis-azoPC and 15 % trans-azoPC lipids) reveal that the TCPI reduces membrane packing density and increases diffusivity of lipids. We have demonstrated an enhanced intercalation of doxorubicin (DOX), an anticancer drug, in the 'bright' state of the membrane compared to that in the 'dark' state. This study - elucidating the complex interplay between lipid composition, photoswitching, and lipid-drug interactions - contributes to the design of lipid-based systems for targeted drug delivery and biomedical applications.

Species-specific metabolism of triphenyl phosphate and its mono-hydroxylated product by human and rat CYP2E1 and the carp ortholog

Organophosphorus flame retardants (PFRs) are a class of flame retardants and environmental pollutants with various biological effects. Recentstudies have evidenced activation of some PFRs by human CYP enzymes (including CYP2E1) for genotoxic effects. However, the activity of CYPs in fish species toward PFR metabolism remains unclear. This study was aimed on comparing the metabolism of triphenyl phosphate (TPHP) and 4-OH-TPHP in human, rat, and common carp, and the involvement of human CYP2E1 and its orthologs in the metabolism, by using fomepizole (4-MP, CYP2E1 inhibitor) as a modulator, in silico molecular docking and dynamics analyses. The rate of TPHP metabolism was apparently faster with human and rat, microsomes than with fish microsomes, the major metabolites were phosphodiester and hydroxylated phosphate, with 30-80 % of TPHP forming unidentified metabolites in the system of each species. 4-OH-TPHP was readily metabolized by both human and rat microsomes, whereas it was hardly metabolized in carp assays. Meanwhile, with 4-MP the transformation of TPHP to 4-OH-TPHP was enhanced in the human/rat systems while suppressed in the carp system. Moreover, the formation of unidentified metabolites in human and rat systems was mostly inhibited by 4-MP. Through molecular dynamics analysis TPHP and its primary metabolites showed high affinity for human and rat CYP2E1, as well as the carp ortholog (CYP2G1-like enzyme), however, the 4-OH-TPHP bond to the latter was too far from the heme to permit a biochemical reaction. This study suggests that the metabolism/activation of TPHP might be favored in mammals rather than carp, a fish species.

In-silico based discovery of potential Keap1 inhibitors using the strategies of pharmacophore screening, molecular docking, and MD simulation studies

Introduction

The main objective of this research is to identify potential leads for developing potent Keap1 inhibitors.

Methods

In the current research article, in-silico methods have been employed to discover potential Keap1 inhibitors. 3D-QSAR was generated using the ChemBL database of Keap1 inhibitors with IC50. The best pharmacophore was selected for the screening of three different libraries namely Asinex, MiniMaybridge, and Zinc. The molecules screened from the databases were filtered through druggability rules and molecular docking studies. The best binding molecules obtained after docking studies were subjected to physicochemical properties toxicity determination by in-silico methods. The best hits were studied for stability in the cavity of Keap1 by molecular dynamic simulations.

Results

The virtual screening of different databases was carried out separately and three leads, were obtained. These lead molecules ASINEX 508, MiniMaybridgeHTS_01719, and ZINC 0000952883 showed the best binding in the Keap1 cavity. The molecular dynamic simulations of the binding complexes of the leads support the docking analysis. The leads (ASINEX 508, MiniMaybridgeHTS_01719, and ZINC 0000952883) were stabilized in the Keap1 binding cavity throughout 100 ns simulation, with average RMSD values of 0.100, 0.114, and 0.106 nm, respectively.

Conclusion

This research proposes three lead molecules as potential Keap1 inhibitors based on high throughput screening, docking, and MD simulation studies. These hit molecules can be used for further design and development of Keap1 inhibitors. This research provides the preliminary data for discovering novel Keap1 inhibitors. It opens new avenues for medicinal chemists to explore antioxidant-stimulating molecules targeting the Keap1-Nrf2 pathway.

Efficient Capture and Traceless Release of Functional CD8+ T Cells with a Microfluidic Chip for Enhanced In Vitro and In Vivo CD4-CAR Transduction

The chimeric antigen receptor (CAR) T cells targeting CD4 expressed cells in acute lymphoblastic leukemia (T-ALL) and acute myeloid leukemia (AML) could reduce the risk of off target effects in normal tissues. However, the efficacy of adoptive cell therapy is predominantly attributed to CD8+ T cells, necessitating their purification before lentivirus transfection to enhance the production of CD4-CAR-T cells. In this study, we developed a microfluidic chip functionalized with an optimized CD8 aptamer, A3t-MU, to facilitate the enrichment and purification of CD8+ T cells. The presented chip showed efficient capture and seamless release of CD8+ T cells from cultured T cells and peripheral blood mononuclear cells (PBMCs). The purity of the released CD8+ T cells reached 98.1%, representing a 13% improvement over the conventional magnetic bead separation method. CD4-CAR was efficiently transduced into the purified CD8+ T cells to construct CAR-T cells. We evaluated the antitumor capability of the CD4-CAR transduced CD8+ T cells (anti-CD4 CD8-CAR T cells) both in vitro and in vivo. The anti-CD4 CD8-CAR T cells exhibited significant cancer-cell-killing capacity across multiple tumor cell lines, including CEM, Jurkat, and MV4-11. Meanwhile, anti-CD4 CD8-CAR T cells significantly inhibited tumor growth in vivo. In conclusion, the presented microfluidic chip offers a cost-effective and high-purity approach for CD8+ T cell separation, enhancing CD4-CAR transduction and achieving efficient antitumor capability both in vitro and in vivo.

Employing Metadynamics to Predict the Membrane Partitioning of Carboxy-2H-Azirine Natural Products

Natural products containing the carboxy-2H-azirine moiety are an exciting target for investigation due to their broad-spectrum antimicrobial activity and new chemical space they afford for novel therapeutic development. The carboxy-2H-azirine moiety, including those appended to well-characterized chemical scaffolds, is understudied, which creates a challenge for understanding potential modes of inhibition. In particular, some known natural product carboxy-2H-azirines have long hydrophobic tails, which could implicate them in membrane-associated processes. In this study, we examined a small set of carboxy-2H-azirine natural products with varied structural features that could alter membrane partitioning. We compared the predicted membrane partitioning and alignment of these compounds to those of established membrane embedders with similar chemical scaffolds. To accomplish this, we developed parameters within the framework of the CHARMM36 force field for the 2H-azirine functional group and performed metadynamics simulations of the partitioning into a model bacterial membrane from aqueous solution. We determined that the carboxy-2H-azirine functional group is strongly hydrophilic, imbuing the long-chain natural products with amphipathicity similar to the known membrane-embedding molecules to which they were compared. For the long-chain analogs, the carboxy-2H-azirine head group stays within 1 nm of the phosphate layer, while the hydrophobic tail sits within the membrane. The carboxy-2H-azirine lacking the long alkyl chain instead partitions completely into aqueous solution.

Metadynamics Study of Lipid-Mediated Antibacterial Toxin Binding to the EmrE Multiefflux Protein

EmrE is a bacterial efflux protein in the small multidrug-resistant (SMR) family present in Escherichia coli. Due to its small size, 110 residues in each dimer subunit, it is an ideal model system to study ligand-protein-membrane interactions. Here in our work, we have calculated the free energy landscape of benzyltrimetylammonium (BTMA) and tetraphenyl phosphonium (TPP) binding to EmrE using the enhanced sampling method-multiple walker metadynamics. We estimate that the free energy of BTMA binding to EmrE is -21.2 ± 3.3 kJ/mol and for TPP is -43.6 ± 3.8 kJ/mol. BTMA passes through two metastable states to reach the binding pocket, while TPP has a more complex binding landscape with four metastable states and one main binding site. Our simulations show that the ligands interact with the membrane lipids at a distance 1 nm away from the binding site which forms a broad local minimum, consistent for both BTMA and TPP. This site can be an alternate entry point for ligands to partition from the membrane into the protein, especially for bulky and/or branched ligands. We also observed the membrane lipid and C-terminal 110HisA form salt-bridge interactions with the helix-1 residue 22LysB. Our free energy estimates and clusters are in close agreement with experimental data and give us an atomistic view of the ligand-protein-lipid interactions. Understanding the binding pathway of these ligands can guide us in future design of ligands that can alter or halt the function of EmrE.

Functional Divergence in the Affinity and Stability of Non-Canonical Cysteines and Non-Canonical Disulfide Bonds: Insights from a VHH and VNAR Study

Single-domain antibodies, including variable domains of the heavy chains of heavy chain-only antibodies (VHHs) from camelids and variable domains of immunoglobulin new antigen receptors (VNARs) from cartilaginous fish, show the therapeutic potential of targeting antigens in a cytosol reducing environment. A large proportion of single-domain antibodies contain non-canonical cysteines and corresponding non-canonical disulfide bonds situated on the protein surface, rendering them vulnerable to environmental factors. Research on non-canonical disulfide bonds has been limited, with a focus solely on VHHs and utilizing only cysteine mutations rather than the reducing agent treatment. In this study, we examined an anti-lysozyme VNAR and an anti-BC2-tag VHH, including their non-canonical disulfide bond reduced counterparts and non-canonical cysteine mutants. Both the affinity and stability of the VNARs and VHHs decreased in the non-canonical cysteine mutants, whereas the reduced-state samples exhibited decreased thermal stability, with their affinity remaining almost unchanged regardless of the presence of reducing agents. Molecular dynamics simulations suggested that the decrease in affinity of the mutants resulted from increased flexibility of the CDRs, the disappearance of non-canonical cysteine-antigen interactions, and the perturbation of other antigen-interacting residues caused by mutations. These findings highlight the significance of non-canonical cysteines for the affinity of single-domain antibodies and demonstrate that the mutation of non-canonical cysteines is not equivalent to the disruption of non-canonical disulfide bonds with a reducing agent when assessing the function of non-canonical disulfide bonds.

Ion Transport at Polymer-Argyrodite Interfaces

Solid-state electrolytes, particularly polymer/ceramic composite electrolytes, are emerging as promising candidates for lithium-ion batteries due to their high ionic conductivity and mechanical flexibility. The interfaces that arise between the inorganic and organic materials in these composites play a crucial role in ion transport mechanisms. While lithium ions are proposed to diffuse across or parallel to the interface, few studies have directly examined the quantitative impact of these pathways on ion transport and little is known about how they affect the overall conductivity. Here, we present an atomistic study of lithium-ion (Li+) transport across well-defined polymer-argyrodite interfaces. We present a force field for polymer-argyrodite interfacial systems, and we carry out molecular dynamics and enhanced sampling simulations of several composite systems, including poly(ethylene oxide) (PEO)/Li6PS5Cl, hydrogenated nitrile butadiene rubber (HNBR)/Li6PS5Cl, and poly(vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP)/Li6PS5Cl. For the materials considered here, Li-ion exhibits a preference for the ceramic material, as revealed by free energy differences for Li-ion between the inorganic and the organic polymer phase in excess of 13 kBT. The relative free energy profiles of Li-ion for different polymeric materials exhibit similar shapes, but their magnitude depends on the strength of interaction between the polymers and Li-ion: the greater the interaction between the polymer and Li-ions, the smaller the free energy difference between the inorganic and organic materials. The influence of the interface is felt over a range of approximately 1.5 nm, after which the behavior of Li-ion in the polymer is comparable to that in the bulk. Near the interface, Li-ion transport primarily occurs parallel to the interfacial plane, and ion mobility is considerably slower near the interface itself, consistent with the reduced segmental mobility of the polymer in the vicinity of the ceramic material. These findings provide insights into ionic complexation and transport mechanisms in composite systems, and will help improve design of improved solid electrolyte systems.

Molecular-level design of alternative media for energy-saving pilot-scale fibrillation of nanocellulose

The outstanding mechanical properties, light weight, and biodegradability of cellulose nanofibrils (CNFs) make them promising components of renewable and sustainable next-generation reinforced composite biomaterials and bioplastics. Manufacturing CNFs at a pilot scale requires disc-refining fibrillation of dilute cellulose fibers in aqueous pulp suspensions to shear the fibers apart into their nanodimensional forms, which is, however, an energy-intensive process. Here, we used atomistic molecular dynamics (MD) simulation to examine media that might facilitate the reduction of interactions between cellulose fibers, thereby reducing energy consumption in fibrillation. The most suitable medium found by the simulations was an aqueous solution with 0.007:0.012 wt.% NaOH:urea, and indeed this was found in pilot-scale experiments to reduce the fibrillation energy by ~21% on average relative to water alone. The NaOH:urea-mediated CNFs have similar crystallinity, morphology, and mechanical strength to those formed in water. The NaOH and urea act synergistically on CNFs to aid fibrillation but at different length scales. NaOH deprotonates hydroxyl groups leading to mesoscale electrostatic repulsion between fibrils, whereas urea forms hydrogen bonds with protonated hydroxyl groups thus disrupting interfibril hydrogen bonds. This suggests a general mechanism in which an aqueous medium that contains a strong base and a small organic molecule acting as a hydrogen-bond acceptor and/or donor may be effectively employed in materials processes where dispersion of deprotonable polymers is required. The study demonstrates how atomic-detail computer simulation can be integrated with pilot-scale experiments in the rational design of materials processes for the circular bioeconomy.

A Coarse-Grained SPICA Makeover for Solvated and Bare Sodium and Chloride Ions

Aqueous ionic solutions are pivotal in various scientific domains due to their natural prevalence and vital roles in biological and chemical processes. Molecular dynamics has emerged as an effective methodology for studying the dynamic behavior of these systems. While all-atomistic models have made significant strides in accurately representing and simulating these ions, the challenge persists in achieving precise models for coarse-grained (CG) simulations. Our study introduces two optimized models for sodium and chloride ions within the nonpolarizable surface property fitting coarse-grained force field (SPICA-FF) framework. The two models represent solvated ions, such as the original FF model, and unsolvated or bare ions. The nonbonded Lennard-Jones interactions were reparameterized to faithfully reproduce bulk properties, including density and surface tension, in sodium chloride solutions at varying concentrations. Notably, these optimized models replicate experimental surface tensions at high ionic strengths, a property not well-captured by the ions of the original model in the SPICA-FF. The optimized unsolvated model also proved successful in reproducing experimental osmotic pressure. Additionally, the newly reparameterized ion models capture hydrophobic interactions within sodium chloride solutions and show qualitative agreement when modeling structural changes in phospholipid bilayers, aligning with experimental observations. For aqueous solutions, these optimized models promise a more precise representation of the ion behavior.

Molecular Basis of the Substrate Specificity of Phosphotriesterase from Pseudomonas diminuta: A Combined QM/MM MD and Electron Density Study

The occurrence of organophosphorus compounds, pesticides, and flame-retardants in wastes is an emerging ecological problem. Bacterial phosphotriesterases are capable of hydrolyzing some of them. We utilize modern molecular modeling tools to study the hydrolysis mechanism of organophosphorus compounds with good and poor leaving groups by phosphotriesterase from Pseudomonas diminuta (Pd-PTE). We compute Gibbs energy profiles for enzymes with different cations in the active site: native Zn2+cations and Co2+cations, which increase the steady-state rate constant. Hydrolysis occurs in two elementary steps via an associative mechanism and formation of the pentacoordinated intermediate. The first step, a nucleophilic attack, occurs with a low energy barrier independently of the substrate. The second step has a low energy barrier and considerable stabilization of products for substrates with good leaving groups. For substrates with poor leaving groups, the reaction products are destabilized relative to the ES complex that suppresses the reaction. The reaction proceeds with low energy barriers for substrates with good leaving groups with both Zn2+and Co2+cations in the active site; thus, the product release is likely to be a limiting step. Electron density and geometry analysis of the QM/MM MD trajectories of the intermediate states with all considered compounds allow us to discriminate substrates by their ability to be hydrolyzed by the Pd-PTE. For hydrolyzable substrates, the cleaving bond between a phosphorus atom and a leaving group is elongated, and electron density depletion is observed on the Laplacian of electron density maps.

Novel Autotaxin Inhibitor ATX-1d Significantly Enhances Potency of Paclitaxel-An In Silico and In Vitro Study

The development of drug resistance in cancer cells poses a significant challenge for treatment, with nearly 90% of cancer-related deaths attributed to it. Over 50% of ovarian cancer patients and 30-40% of breast cancer patients exhibit resistance to therapies such as Taxol. Previous literature has shown that cytotoxic cancer therapies and ionizing radiation damage tumors, prompting cancer cells to exploit the autotaxin (ATX)-lysophosphatidic acid (LPA)-lysophosphatidic acid receptor (LPAR) signaling axis to enhance survival pathways, thus reducing treatment efficacy. Therefore, targeting this signaling axis has become a crucial strategy to overcome some forms of cancer resistance. Addressing this challenge, we identified and assessed ATX-1d, a novel compound targeting ATX, through computational methods and in vitro assays. ATX-1d exhibited an IC50 of 1.8 ± 0.3 μM for ATX inhibition and demonstrated a significant binding affinity for ATX, as confirmed by MM-GBSA, QM/MM-GBSA, and SAPT in silico methods. ATX-1d significantly amplified the potency of paclitaxel, increasing its effectiveness tenfold in 4T1 murine breast carcinoma cells and fourfold in A375 human melanoma cells without inducing cytotoxic effects as a single agent.

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.

All-Atom Membrane Builder via Multiscale Simulation

I present an automated and flexible tool designed for constructing bilayer membranes at all-atom (AA) resolution. The builder initiates the construction and equilibration of bilayer membranes at Martini coarse-grained (CG) resolution, followed by resolution enhancement to the atomic level using the accompanying backmapping tool. Notably, this tool enables users to create bilayer membranes with user-defined lipid compositions and protein structures, while also offering the flexibility to accommodate new lipid types. To assess the simplicity and robustness of the tool, I demonstrate the construction of several membranes incorporating protein structures. The tool is freely available at github.com/ksy141/mstool.

Analysis of action of 1,4-naphthoquinone scaffold-derived compounds against acute myeloid leukemia based on network pharmacology, molecular docking and molecular dynamics simulation

1,4-Naphthoquinone scaffold-derived compounds has shown considerable pharmacological properties against cancer, including acute myeloid leukemia (AML) However, its impact and mechanisms in AML are uncertain. In this study, the mechanisms of 1,4-naphthoquinone scaffold-derived compounds against AML were investigated via network pharmacology, molecular docking and molecular dynamics simulation. ASINEX database was used to collect the 1,4-naphthoquinone scaffold-derived compounds, and compounds were extracted from the software to evaluate their drug similarity and toxicity. The potential targets of compounds were retrieved from the SwissTargetPrediction Database and the Similarity Ensemble Approach Database, while the potential targets of AML were obtained from the GeneCards databases and Gene Expression Omnibus. The STRING database was used to construct a protein-protein interaction (PPI) network, topologically and Cyto Hubb plugin of Cytoscape screen the central targets. After selecting the potential key targets, the gene ontology (GO) function annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed for the intersection targets, and a network map of "compounds-potential targets-pathway-disease" were constructed. Molecular docking of the compounds with the core target was performed, and core target with the strongest binding force and 1,4-naphthoquinone scaffold-derived compounds was selected for further molecular dynamics simulation and further molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) approach verification. In addition, the Bloodspot database was applied to perform the overall survival of core targets. A total of 19 1,4-naphthoquinone scaffold-derived compounds were chosen out, and then 836 targets of compounds, 96 intersection targets of AML were screened. Core targets include STAT3, TLR4, HSP90AA1, JUN, MMP9, PTPRC, JAK2, PTGS2, KIT and CSF1R. GO functional enrichment analysis revealed that 90 biological processes, 10 cell components and 12 molecular functions were enriched while KEGG pathway enrichment analysis revealed 34 enriched signaling pathways. Analysis of KEGG enrichment hinted that these 10 core genes were located in the pathways in cancer, suggesting that 1,4-naphthoquinone scaffold-derived compounds had potential activity against AML. Molecular docking analysis revealed that the binding energies between 1,4-naphthoquinone scaffold-derived compounds and the core proteins were all higher than - 6 kcal/mol, indicating that the 10 core targets all had strong binding ability with compounds. Moreover, a good binding capacity was inferred from molecular dynamics simulations between compound 7 and MMP9. The total binding free energy calculated using the MM/GBSA approach revealed values of - 6356.865 kcal/mol for the MMP9-7 complex. In addition, Bloodspot database results exhibited that HSP90AA1, MMP9 and PTPRC were associated with overall survival. The findings provide foundations for future studies into the interaction underlying the anti-AML potential of compounds with 1,4-naphthoquinone-based scaffold structures. Compounds with 1,4-naphthoquinone-based scaffold structures exhibits considerable potential in mitigating and treating AML through multiple targets and pathways.

Pharmacophore-Based Study: An In Silico Perspective for the Identification of Potential New Delhi Metallo-β-lactamase-1 (NDM-1) Inhibitors

In the ongoing battle against antibiotic-resistant bacteria, New Delhi metallo-β-lactamase-1 (NDM-1) has emerged as a significant therapeutic challenge due to its ability to confer resistance to a broad range of β-lactam antibiotics. This study presents a pharmacophore-based virtual screening, docking, and molecular dynamics simulation approach for the identification of potential inhibitors targeting NDM-1, a critical enzyme associated with antibiotic resistance. Through the generation of a pharmacophore model and subsequent virtual screening of compound libraries, candidate molecules (ZINC29142850 (Z1), ZINC78607001 (Z2), and ZINC94303138 (Z3)) were prioritized based on their similarity to known NDM-1 binder (hydrolyzed oxacillin (0WO)). Molecular docking studies further elucidated the binding modes and affinities of the selected compounds towards the active site of NDM-1. These compounds demonstrated superior binding affinities to the enzyme compared to a control compound (-7.30 kcal/mol), with binding scores of -7.13, -7.92, and -8.10 kcal/mol, respectively. Binding interactions within NDM-1's active site showed significant interactions with critical residues such as His250, Asn220, and Trp93 for these compounds. Subsequent molecular dynamics simulations were conducted to assess the stability of the ligand-enzyme complexes, showing low root mean square deviation (RMSD) values between 0.5 and 0.7 nm for Z1, Z2, which indicate high stability. Z2's compactness in principal component analysis (PCA) suggests that it can stabilize particular protein conformations more efficiently. Z2 displays a very cohesive landscape with a notable deep basin, suggesting a very persistent conformational state induced by the ligand, indicating robust binding and perhaps efficient inhibition. Z2 demonstrates the highest binding affinity among the examined compounds with a binding free energy of -25.68 kcal/mol, suggesting that it could offer effective inhibition of NDM-1. This study highlights the efficacy of computational tools in identifying novel antimicrobial agents against resistant bacteria, accelerating drug discovery processes.

Burkholderia cenocepacia epigenetic regulator M.BceJIV simultaneously engages two DNA recognition sequences for methylation

Burkholderia cenocepacia is an opportunistic and infective bacterium containing an orphan DNA methyltransferase called M.BceJIV with roles in regulating gene expression and motility of the bacterium. M.BceJIV recognizes a GTWWAC motif (where W can be an adenine or a thymine) and methylates N6 of the adenine at the fifth base position. Here, we present crystal structures of M.BceJIV/DNA/sinefungin ternary complex and allied biochemical, computational, and thermodynamic analyses. Remarkably, the structures show not one, but two DNA substrates bound to the M.BceJIV dimer, with each monomer contributing to the recognition of two recognition sequences. We also show that methylation at the two recognition sequences occurs independently, and that the GTWWAC motifs are enriched in intergenic regions in the genomes of B. cenocepacia strains. We further computationally assess the interactions underlying the affinities of different ligands (SAM, SAH, and sinefungin) for M.BceJIV, as a step towards developing selective inhibitors for limiting B. cenocepacia infection.

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.

Identification of a Clade-Specific HLA-C*03:02 CTL Epitope GY9 Derived from the HIV-1 p17 Matrix Protein

Efforts towards an effective HIV-1 vaccine have remained mainly unsuccessful. There is increasing evidence for a potential role of HLA-C-restricted CD8+ T cell responses in HIV-1 control, including our recent report of HLA-C*03:02 among African children. However, there are no documented optimal HIV-1 CD8+ T cell epitopes restricted by HLA-C*03:02; additionally, the structural influence of HLA-C*03:02 on epitope binding is undetermined. Immunoinformatics approaches provide a fast and inexpensive method to discover HLA-restricted epitopes. Here, we employed immunopeptidomics to identify HLA-C*03:02 CD8+ T cell epitopes. We identified a clade-specific Gag-derived GY9 (GTEELRSLY) HIV-1 p17 matrix epitope potentially restricted to HLA-C*03:02. Residues E62, T142, and E151 in the HLA-C*03:02 binding groove and positions p3, p6, and p9 on the GY9 epitope are crucial in shaping and stabilizing the epitope binding. Our findings support the growing evidence of the contribution of HLA-C molecules to HIV-1 control and provide a prospect for vaccine strategies.

Binding of steroid substrates reveals the key to the productive transition of the cytochrome P450 OleP

OleP is a bacterial cytochrome P450 involved in oleandomycin biosynthesis as it catalyzes regioselective epoxidation on macrolide intermediates. OleP has recently been reported to convert lithocholic acid (LCA) into murideoxycholic acid through a highly regioselective reaction and to unspecifically hydroxylate testosterone (TES). Since LCA and TES mainly differ by the substituent group at the C17, here we used X-ray crystallography, equilibrium binding assays, and molecular dynamics simulations to investigate the molecular basis of the diverse reactivity observed with the two steroids. We found that the differences in the structure of TES and LCA affect the capability of these molecules to directly form hydrogen bonds with N-terminal residues of OleP internal helix I. The establishment of these contacts, by promoting the bending of helix I, fosters an efficient trigger of the open-to-closed structural transition that occurs upon substrate binding to OleP and contributes to the selectivity of the subsequent monooxygenation reaction.

Marine-Derived Furanones Targeting Quorum-Sensing Receptors in Pseudomonas aeruginosa: Molecular Insights and Potential Mechanisms of Inhibition

The quorum-sensing (QS) machinery in disease-causing microorganisms is critical in developing antibiotic resistance. In Pseudomonas aeruginosa, QS is involved in biofilm formation, virulence factors production, and general tolerance to antimicrobials. Owing to the major role QS plays, interference in the process is probably a facile route to overcome antimicrobial resistance. Some furanone-derived compounds from marine sources have shown promising anti-QS activity. However, their protein targets and potential mechanisms of action have not been explored. To elucidate their potential protein targets in this study, marine metabolites with furanone backbones similar to their cognitive autoinducers (AIs) were screened against various QS receptors (LasR, RhlR, and PqsR) using molecular docking and molecular dynamics (MD) simulation techniques. The order by which the compounds bind to the receptors follows LasR > RhlR > PqsR. Compounds exhibited remarkable stability against LasR and RhlR, likely because the AIs of these receptors are structural analogs of furanones. Furanones with shorter alkyl side chains bound strongly against RhlR. The presence of halogens improved binding against various receptors. PqsR, with its hydrophobic-binding site and structurally different AIs, showed weaker binding. This study provides a molecular basis for the design of potent antagonists against QS receptors using marine-derived furanones.

Druglike Molecules Binding to Large Membrane Proteins: Absolute Binding Free Energy Computation

In this research, we employed the alchemical double-decoupling method alongside restraining potentials, coupled with the FEPMD method, to ascertain the standard binding free energy of a drug-like molecule termed BHQ and three analogous compounds engineered with progressive addition of bulky para-alkyl groups binding to SERCA (Ca2+-ATPase of skeletal muscle sarcoplasmic reticulum). Integral transmembrane proteins represent crucial drug targets in numerous therapeutic interventions, presenting computational challenges due to their considerable system sizes. Our approach integrated the generalized born potential method and the spherical solvent boundary potential method, allowing us to explicitly focus on the active binding site while treating the remainder of the system implicitly. We evaluated contributions to the standard binding free energy from distinct interaction potentials: electrostatic, repulsive, dispersive, and restraining potentials, computed separately. The resulting absolute binding free energy of BHQ (11.63 kcal/mol) closely aligns with experimental measurements (10.56 kcal/mol). Notably, an accurate estimation of the absolute binding free energy was achieved for the simplest analog, created with the addition of a single para-methyl group. However, the analog with two para-methyl groups exhibited the highest binding free energy, which disagreed with experimental results. Determining the binding free energy of the BHQ analog engineered with three para-methyl groups presented challenges in convergence and resulted in the lowest free energy among the three.

Bisphenol F and Bisphenol S in a Complex Biomembrane: Comparison with Bisphenol A

Bisphenols are a group of endocrine-disrupting chemicals used worldwide for the production of plastics and resins. Bisphenol A (BPA), the main bisphenol, exhibits many unwanted effects. BPA has, currently, been replaced with bisphenol F (BPF) and bisphenol S (BPS) in many applications in the hope that these molecules have a lesser effect on metabolism than BPA. Since bisphenols tend to partition into the lipid phase, their place of choice would be the cellular membrane. In this paper, I carried out molecular dynamics simulations to compare the localization and interactions of BPA, BPF, and BPS in a complex membrane. This study suggests that bisphenols tend to be placed at the membrane interface, they have no preferred orientation inside the membrane, they can be in the monomer or aggregated state, and they affect the biophysical properties of the membrane lipids. The properties of bisphenols can be attributed, at least in part, to their membranotropic effects and to the modulation of the biophysical membrane properties. The data support that both BPF and BPS, behaving in the same way in the membrane as BPA and with the same capacity to accumulate in the biological membrane, are not safe alternatives to BPA.

Kernel-Based Minimal Distributed Charges: A Conformationally Dependent ESP-Model for Molecular Simulations

A kernel-based method (kernelized minimal distributed charge model (kMDCM)) to represent the molecular electrostatic potential (ESP) in terms of off-center point charges is introduced. The positions of the charges adapt to the molecular geometry and allow the description of intramolecular charge flow. Using Gaussian kernels and atom-atom distances as the features, the ESPs for water and methanol are shown to improve by at least a factor of 2 compared with point charge models fit to an ensemble of structures. The conformationally fluctuating molecular dipole moment of water is reproduced almost twice as accurately using kMDCM compared with static PCs, despite not fitting to the dipole directly. The role of hyperparameters in the kernelization is investigated and their implication on model performance and simulation stability is discussed. Combining kMDCM for the electrostatics and reproducing kernels for the bonded terms allows energy-conserving simulations of 2000 water molecules with periodic boundary conditions on the nanosecond time scale. These MD simulations sample geometries outside the training set but remain stable, which demonstrates the robustness of the model and its implementation.

Machine learning, network pharmacology, and molecular dynamics reveal potent cyclopeptide inhibitors against dengue virus proteins

The dengue virus is a major global health hazard responsible for an estimated 390 million diseases yearly. This study focused on identifying cyclopeptide inhibitors for envelope structural proteins E, NS1, NS3, and NS5. Additionally, 5579 cyclopeptides were individually screened against the four target proteins using a machine learning-based quantitative structure-activity relationship model. Subsequently, the best 10 cyclopeptides from each protein were selected for molecular docking with their corresponding proteins. Moreover, the protein-peptide complexes with the highest affinity were subjected to a 100-ns molecular dynamics simulation. The protein-protein complexes exhibited superior structural stability and binding interactions. Based on the results of the MD simulation analyses, which included checking values for Root Mean Square Deviation, Root Mean Square Fluctuation, Principal Component Analysis (PCA), free energy landscapes, and energetic components, it was found that NS5-CP03714 complex is more stable and has stronger binding interactions than NS3-CP02054. PCA and free energy landscape plots have confirmed the higher conformational stability of NS5-CP03714. Analysis of the energetic components revealed that NS5-CP03714 (total binding energy = - 47.19 kcal/mol) exhibits more favorable interaction energies and overall binding energy compared to NS3-CP02054 (total binding energy = - 27.36 kcal/mol), suggesting a stronger and more stable formation of the complex. In addition, the drug-target network of two specific peptides (CP02950 and CP05582) and their associated target proteins were analyzed. This analysis revealed valuable information about their ability to target several proteins and their potential for broad-spectrum activity. Additional experimental investigations are necessary to validate these computational results and assess the efficacy of identified peptide inhibitors in biological systems.

Predicting 3D RNA structure from the nucleotide sequence using Euclidean neural networks

Fast and accurate 3D RNA structure prediction remains a major challenge in structural biology, mostly due to the size and flexibility of RNA molecules, as well as the lack of diverse experimentally determined structures of RNA molecules. Unlike DNA structure, RNA structure is far less constrained by basepair hydrogen bonding, resulting in an explosion of potential stable states. Here, we propose a convolutional neural network that predicts all pairwise distances between residues in an RNA, using a recently described smooth parametrization of Euclidean distance matrices. We achieve high-accuracy predictions on RNAs up to 100 nt in length in fractions of a second, a factor of 107 faster than existing molecular dynamics-based methods. We also convert our coarse-grained machine learning output into an all-atom model using discrete molecular dynamics with constraints. Our proposed computational pipeline predicts all-atom RNA models solely from the nucleotide sequence. However, this method suffers from the same limitation as nucleic acid molecular dynamics: the scarcity of available RNA crystal structures for training.

Binding Interaction of Coumarin Derivative Daphnetin with Ovalbumin: Molecular Insights into the Complexation Process and Effects of Metal Ions and pH in the Binding and Antifibrillation Studies

This study investigates the interaction between daphnetin and ovalbumin (OVA) as well as its potential to inhibit OVA fibrillation using both spectroscopic and computational analysis. A moderate binding affinity of 1 × 104 M-1 was observed between OVA and daphnetin, with a static quenched mechanism identified during the fluorescence quenching processes. Metal ions' (Cu2+ and Zn2+) presence led to an increase in the binding affinities of daphnetin toward OVA, mirroring a similar trend observed with the pH variation. Synchronous and 3D fluorescence studies indicated an increase in the polarity of the microenvironment surrounding the Trp residues during binding. Interestingly, circular dichroism and Fourier transform infrared studies showed a significant change in the secondary structure of OVA upon binding with daphnetin. The efficacy of daphnetin in inhibiting protein fibrillation was confirmed through thioflavin T and Congo Red binding assays along with fluorescence microscopic imaging analysis. The thermodynamic assessment showed positive ΔH° [+(29.34 ± 1.526) kJ mol-1] and ΔS° [+(181.726 ± 5.465) J mol-1] values, indicating the presence of the hydrophobic forces, while negative ΔG° signifies spontaneous binding interactions. These experimental findings were further correlated with computational analysis, revealing daphnetin dynamics within the binding site of OVA.

Adsorption of Staphylococcus aureus biofilm associated compounds on silica probed with molecular dynamics simulations

Staphylococcus aureus (S. aureus) is a potentially pathogenic bacterium that commonly colonizes surfaces through the formation of biofilms. Silica glass is a common material in the built environment, especially in laboratory and medical spaces. The chemical and physical mechanisms by which S. aureus initially adheres to surfaces are unclear. In this study, the adsorption of several S. aureus biofilm associated compounds on silica is probed using molecular dynamics simulations. Model compounds containing a phosphorylated backbone, N-acetylglucosamine (GlcNAc), or D-alanine (D-Ala) were simulated across a range of pH. GlcNAc adsorption is unfavorable and insensitive to pH. D-Ala adsorption is unfavorable across the range of tested pH. Phosphorylated backbone adsorption is unfavorable at low pH but favorable at high pH. Adsorbate titration and solution salt concentration were probed to establish effects of molecular charge and charge screening. Hydrogen bonding between compounds and the silica surface is a key factor for stronger adsorption. The findings of this study are important for the rational design of improved silica surfaces through chemical functionalization or through the application of optimal chemical disinfectants that discourage the initial stages of biofilm growth.

Location and interaction of idebenone and mitoquinone in a membrane similar to the inner mitochondrial membrane. Comparison with ubiquinone 10

Oxygen is essential for aerobic life on earth but it is also the origin of harmful reactive oxygen species (ROS). Ubiquinone is par excellence the endogenous cellular antioxidant, but a very hydrophobic one. Because of that, other molecules have been envisaged, such as idebenone (IDE) and mitoquinone (MTQ), molecules having the same redox active benzoquinone moiety but higher solubility. We have used molecular dynamics to determine the location and interaction of these molecules, both in their oxidized and reduced forms, with membrane lipids in a membrane similar to that of the mitochondria. Both IDE and reduced IDE (IDOL) are situated near the membrane interface, whereas both MTQ and reduced MTQ (MTQOL) locate in a position adjacent to the phospholipid hydrocarbon chains. The quinone moieties of both ubiquinone 10 (UQ10) and reduced UQ10 (UQOL10) in contraposition to the same moieties of IDE, IDOL, MTQ and MTQOL, located near the membrane interphase, whereas the isoprenoid chains remained at the middle of the hydrocarbon chains. These molecules do not aggregate and their functional quinone moieties are located in the membrane at different depths but near the hydrophobic phospholipid chains whereby protecting them from ROS harmful effects.

Identifying potential inhibitors of C-X-C motif chemokine ligand10 against vitiligo: structure-based virtual screening, molecular dynamics simulation, and principal component analysis

The research aims to envisage small molecule inhibitors targeting the C-X-C motif chemokine ligand 10 (CXCL10) of the JAK/STAT pathway. CXCL10 plays a significant role in inducing auto-immunity in vitiligo through JAK/STAT pathway. To accomplish the aim, structure-based virtual screening with fundamental search limits, e.g., molecular weight (MW ≤ 500 Da), hydrogen bond donor (HBD ≤ 5), hydrogen bond acceptor (HBA ≤ 10), and lipophilicity (logP ≤ 5) was used to screen investigational molecules from MCULE database. The SBVS-ligand hits were sifted through toxicity profiling followed by filtration through the Brain or IntestinaL EstimateD-Egg model to check the human intestinal abortion and blood-brain barrier permeation based on two physicochemical properties, including topological surface area and WLOGP. The BOILED-Egg filtered compounds were passed through drug-likeness features other than Pfizer's Lipinski rule of five, viz., Ghose filters, Muegge filters, Egan parameters, and Veber filters, followed by medicinal chemistry's pan assay interference structure and Brenk alert investigation. Chemical compounds that comply with the above ADME descriptors were docked with target protein CXCL10 via AutoDock Vina. The stability of the top two ligand hits was assessed through dynamics simulations of 100 ns and principal component analysis and compared with the reference drugs Baicalein and EGCG. Based on the findings of Gibbs free energy of binding, ADME profiling, medicinal chemistry attributes depiction, root-mean-square deviation, root-mean-square fluctuation, solvent accessible surface area, the free energy of solvation, the radius of gyration, and PCA, MCULE2726078782-0-2 was found better than potential reference drug Baicalein.Communicated by Ramaswamy H. Sarma.

Combined role of stearic acid and maleic anhydride in the development of thermoplastic starch-based materials with ultrahigh ductility and durability

The diverse properties reported for starch-based materials indicate their potential for use in the preparation of biodegradable flexible actuators. However, their natural brittleness and lack of durability after modification limit their practical application. Therefore, we propose a strategy for preparing flexible starch-based composites. The results of macro/micro property characterizations and molecular dynamics simulations indicated that using starch, maleic anhydride, and stearic acid (SA), the mobility of the starch chains was enhanced and retrogradation was inhibited through the synergistic effects induced by chain breaking, complex formation with SA, and esterification of the starch molecules. In addition, the elongation at break of the modified starch (MS) reached 2070 %, and considerable ductility (>1000 %) as well as well-complexed structure were maintained after six months. Furthermore, the MS was able to undergo self-healing after fracture or a temperature-controlled stiffness transition. Moreover, it underwent complete degradation in soil within 30 d. Finally, an actuator was prepared by doping the MS with nano-Fe3O4 particles to realize a dual magnetic and optical response. Dynamic monitoring was also achieved based on the electrical signal, thereby demonstrating the broad application scope of this material in the development of biodegradable flexible actuators.

Comprehensive views toward the biomolecular recognition of an anticancer drug, leflunomide with human serum albumin

Biomolecular association of an anticancer drug, leflunomide (LEF) with human serum albumin (HSA), the leading ligands carrier in human circulation was characterized using biophysical (i.e., fluorescence, absorption and voltammetric) methods and computational (i.e., molecular docking and molecular dynamics simulation) techniques. Evaluations of fluorescence, absorption and voltammetric findings endorsed the complex formation between LEF and HSA. An inverse relationship of Stern-Volmer constant-temperature and hyperchromic shift of the protein's absorption signal with addition of LEF confirmed the LEF quenched the HSA fluorescence through static process. Moderate nature of binding strength (binding constant = 2.76-4.77 × 104 M-1) was detected towards the LEF-HSA complexation, while the association process was naturally driven via hydrophobic interactions, van der Waals interactions and hydrogen bonds, as evident from changes in entropy (ΔS= + 19.91 J mol-1 K-1) and enthalpy (ΔH = - 20.09 kJ mol-1), and molecular docking assessments. Spectral analyses of synchronous and three-dimensional fluorescence validated microenvironmental fluctuations near Trp and Tyr residues upon LEF binding to the protein. LEF association with HSA significantly defended temperature-induced destabilization of the protein. Although LEF was found to attach to HSA at Sudlow's sites I and II, but exhibited greater preference toward its site I, as detected by the investigations of competitive site-marker displacement. Molecular dynamics simulation assessment revealed that the complex attained equilibrium throughout simulations, showing the LEF-HSA complex constancy.Communicated by Ramaswamy H. Sarma.

Deciphering the Potential Therapeutic Effects of Hydnocarpus wightianus Seed Extracts using in vitro and in silico approaches

Phytocompounds possess the potential to treat a broad spectrum of disorders due to their remarkable bioactivity. Naturally occurring compounds possess lower toxicity profiles, which making them attractive targets for drug development. Hydnocarpus wightianus seeds were extracted using ethanol, acetone, and hexane solvents. The extracts were evaluated for phytochemicals screening and other therapeutic characteristics, such as free radicals scavenging, anti α-amylase, anti α-glucosidase, and anti-bacterial activities. The ethanolic extract exhibited noteworthy antibacterial characteristics and demonstrated considerable antioxidant, and anti-diabetic effects. The IC50 value of the ethanolic extract for Dpph, α-amylase, and α-glucosidase were found to be 77.299 ± 3.381 μg/mL, 165.56 2.56 μg/mL, and 136.58 ± 5.82 μg/mL, respectively. The ethanolic extract was effective against Methicillin resistant Staphylococcus aureus (26 mm zone of inhibition at 100 μL concentration). Molecular docking investigations revealed the phytoconstituent's inhibitory mechanisms against diabetic, free radicals, and bacterial activity. Docking score for phytocompounds against targeted protein varies from -7.2 to -5.1 kcal/mol. The bioactive compounds present in the ethanolic extract were identified by Gas chromatography/Mass spectrometry analysis, followed by molecular docking and molecular dynamic simulation studies to further explore the phytoconstituent's inhibitory mechanism of α-glucosidase, ∝-amylase, radical scavenging, and bacterial activity. The electronic structure and possible pharmacological actions of the phytocompound were revealed through the use of Density Functional Theory (DFT) analysis. Computational and in vitro studies revealed that these identified compounds have anti-diabetic, anti-oxidant, and anti-bacterial activities against antibiotic-resistant strain of Staphylococcus aureus.

Discovery of Protease-activated receptor 2 antagonists derived from phenylalanine for the treatment of breast cancer

Protease-activated receptor 2 (PAR2) has garnered attention as a potential therapeutic target in breast cancer. PAR2 is implicated in the activation of extracellular signal-regulated kinase 1/2 (ERK 1/2) via G protein and beta-arrestin pathways, contributing to the proliferation and metastasis of breast cancer cells. Despite the recognized role of PAR2 in breast cancer progression, clinically effective PAR2 antagonists remain elusive. To address this unmet clinical need, we synthesized and evaluated a series of novel compounds that target the orthosteric site of PAR2. Using in silico docking simulations, we identified compound 9a, an optimized derivative of compound 1a ((S)-N-(1-(benzylamino)-1-oxo-3-phenylpropan-2-yl)benzamide), which exhibited enhanced PAR2 antagonistic activity. Subsequent molecular dynamics simulations comparing 9a with the partial agonist 9d revealed that variations in ligand-induced conformational changes and interactions dictated whether the compound acted as an antagonist or agonist of PAR2. The results of this study suggest that further development of 9a could contribute to the advancement of PAR2 antagonists as potential therapeutic agents for breast cancer.

Targeting COVID-19 and varicocele by blocking inflammasome: Ligand-based virtual screening

COVID-19 is a new generation of outbreaks that invade not only local emerging region, continental but also the whole globe. Varicocele on the other hand, is a testicular vascular disease that underlies 40 % of male infertility cases. Fortunately, the two diseases can be blocked through targeting one common target, NLRP3 inflammasome. Upon searching for similar drugs that gained FDA-approval in ChEMBL library along with examining their potential blockade of the receptor through docking using CB-DOCK-2, three potential approved drugs can be repurposed, ChEMBL 4297185, ChEMBL 1201749, ChEMBL 1200545 which had binding energy of -9.8 and -9.7 kcal/mol (stronger than the reference inhibitor, -9.3 kcal/mol). Also, ADME profile of the top 3 drugs showed better attributes. Also, the simulated proteins exhibited stable pattern with strong free binding energies. Among the potential inhibitor drugs ChEMBL 4297185 was found to remain inside the binding site of the protein during the 200 ns simulation time. Hence, it is anticipated to have the highest binding and thus inhibition potential against the protein. The suggested drugs, especially ChEMBL 4297185, are potentially repurposable toward treating COVID-19 and varicocele which deserve further experimental validation.

Viral methyltransferase inhibitors: berbamine, venetoclax, and ponatinib as efficacious antivirals against chikungunya virus

Chikungunya virus (CHIKV), transmitted by mosquitoes, poses a significant global health threat. Presently, no effective treatment options are available to reduce the disease burden. The lack of approved therapeutics against CHIKV and the complex spectrum of chronic musculoskeletal and neurological manifestations raise significant concerns, and repurposing drugs could offer swift avenues in the development of effective treatment strategies. RNA capping is a crucial step meditated by non-structural protein 1 (nsP1) in CHIKV replication. In this study, FDA-approved antivirals targeting CHIKV nsP1 methyltransferase (MTase) have been identified by structure-based virtual screening. Berbamine Hydrochloride (BH), ABT199/Venetoclax (ABT), and Ponatinib (PT) were the top-hits, which exhibited robust binding energies. Tryptophan fluorescence spectroscopy-based assay confirmed binding of BH-, ABT-, and PT to purified nsP1 with KD values ∼5.45 μM, ∼161.3 μM, and ∼3.83 μM, respectively. In a capillary electrophoresis-based assay, a decrease in CHIKV nsP1 MTase activity was observed in a dose-dependent manner. Treatment with BH, ABT, and PT lead to a dose-dependent reduction in the virus titer with IC50 < 100, ∼6.75, and <3.9 nM, respectively, and reduced viral mRNA levels. The nsP1 MTases are highly conserved among alphaviruses; therefore, BH, ABT, and PT, as expected, inhibited replication machinery in Sindbis virus (SINV) replicon assay with IC50 ∼1.94, ∼0.23, and >1.25 μM, respectively. These results highlight the potential of repurposing drugs as rapid and effective antiviral therapeutics against CHIKV.

Selective protease inhibitors from secondary metabolites of Philippine medicinal plants against porcine epidemic diarrhea virus: A computational veterinary drug discovery approach

Background

Porcine epidemic diarrhea virus (PEDV) is a recurring coronavirus that causes severe diarrhea in pigs with high mortality and morbidity rates, especially in neonatal pigs. Despite the availability of vaccines, their efficacy is limited owing to antigenic differences between the vaccine and field strains, which poses a challenge to infection control. Antiviral drugs targeting conserved PEDV proteins show promise for complementing vaccination strategies. PEDV Nsp3 (PL2Pro) and Nsp5 (3CLPro) are essential proteases vital for viral replication, making them attractive targets for drug development against PEDV.

Aim

To address the lack of therapeutics against recurring PEDV outbreaks and bridge the gap in the application of bioinformatics in veterinary drug discovery, this study aimed to discover compounds that inhibit PEDV proteases from Philippine medicinal plants by applying a modified virtual screening methodology that considers the physiology of swine hosts.

Methods

This study employed a library of 690 metabolites from Philippine medicinal plants to screen for potential protease inhibitors targeting PEDV PL2Pro and 3CLPro. This includes evaluating the binding affinity, pharmacokinetics, dynamic stability, and critical binding site residues. Compounds demonstrating high affinity underwent a modified ADMET analysis, considering the enteric localization of the virus and potential toxicity to swine hosts. Furthermore, molecular dynamics simulations assessed compound stability under physiological swine conditions.

Results

The study identified Bisandrographolide from Andrographis paniculata, CID 162866964 from Euphorbia neriifolia, and betulinic acid from Vitex negundo and Ocimum basilicum as metabolites that bind favorably and selectively to PEDV 3CLPro and have excellent pharmacokinetic properties and dynamic stability. In contrast, no selective inhibitor for PL2pro passed the same criteria.

Conclusion

Employing the modified virtual screening protocol tailored for swine host considerations, the compounds identified in this study are anticipated to exert inhibitory effects against PEDV without off-target binding to analogous swine proteases and receptors. CID 162866964, bisandrographolide, and betulinic acid show promise for developing potent antivirals against PEDV.

A novel step towards the heterologous biosynthesis of paclitaxel: Characterization of T1βOH taxane hydroxylase

In the quest for innovative cancer therapeutics, paclitaxel remains a cornerstone in clinical oncology. However, its complex biosynthetic pathway, particularly the intricate oxygenation steps, has remained a puzzle in the decades following the characterization of the last taxane hydroxylase. The high divergence and promiscuity of enzymes involved have posed significant challenges. In this study, we adopted an innovative approach, combining in silico methods and functional gene analysis, to shed light on this elusive pathway. Our molecular docking investigations using a library of potential ligands uncovered TB574 as a potential missing enzyme in the paclitaxel biosynthetic pathway, demonstrating auspicious interactions. Complementary in vivo assays utilizing engineered S. cerevisiae strains as novel microbial cell factory consortia not only validated TB574's critical role in forging the elusive paclitaxel intermediate, T5αAc-1β,10β-diol, but also achieved the biosynthesis of paclitaxel precursors at an unprecedented yield including T5αAc-1β,10β-diol with approximately 40 mg/L. This achievement is highly promising, offering a new direction for further exploration of a novel metabolic engineering approaches using microbial consortia. In conclusion, our study not only furthers study the roles of previously uncharacterized enzymes in paclitaxel biosynthesis but also forges a path for pioneering advancements in the complete understanding of paclitaxel biosynthesis and its heterologous production. The characterization of T1βOH underscores a significant leap forward for future advancements in paclitaxel production using heterologous systems to improve cancer treatment and pharmaceutical production, thereby holding immense promise for enhancing the efficacy of cancer therapies and the efficiency of pharmaceutical manufacturing.

Targeting MDM2-p53 interaction in Glioblastoma: Transcriptomic analysis and Peptide-Based inhibition strategy

MDM2 is a gene that encodes a protein involved in cell survival, growth, and DNA repair. It has been implicated in the development and progression of glioblastoma (GBM). Inhibition of the MDM2-p53 interaction has emerged as a promising strategy for treating GBM. In this study, we performed comprehensive transcriptomic expression analysis from diverse datasets and observed MDM2 overexpression in a subset of GBM cases. MDM2 negatively regulates the major onco-suppressor p53. The interaction between MDM2 and p53 is a promising target for cancer therapy, as it can trigger p53-mediated cell death in response to different stress conditions, such as oncogene activation or DNA damage. In this study, we have identified a peptide-based inhibition of MDM2 as a therapeutic strategy for GBM. We have further validated the stability of the MDM2-peptide interaction using a molecular structural dynamics approach. The major trajectories, including root mean square of deviation (RMSD), root mean square of fluctuation (RMSF), and radius of gyration (RoG), indicate that the candidate peptides have a more stable binding compared to the native ligand and control drug. The stability of the binding interaction was further estimated by MMGBSA analysis, which also suggests that MDM2 has a stable binding with both peptide molecules. Based on these results, peptides P-1843 and P-3837 could be tested further for experimental validation to confirm their targeted inhibition of MDM-2. This approach could provide a highly selective and efficient inhibitor with potentially fewer side effects and less toxicity compared to small drug-based molecules.

Inhibition of Monkeypox Virus DNA Polymerase Using Moringa oleifera Phytochemicals: Computational Studies of Drug-Likeness, Molecular Docking, Molecular Dynamics Simulation and Density Functional Theory

The emergence of zoonotic monkeypox (MPX) disease, caused by the double-stranded DNA monkeypox virus (MPXV), has become a global threat. Due to unavailability of a specific small molecule drug for MPX, this study investigated Moringa oleifera phytochemicals to find potent and safe inhibitors of DNA Polymerase (DNA Pol), a poxvirus drug target due to its role in the viral life cycle. For that, 146 phytochemicals were screened through drug-likeness and molecular docking analyses. Among these, 136 compounds exhibited drug-like properties, with Gossypetin showing the highest binding affinity (- 7.8 kcal/mol), followed by Riboflavin (- 7.6 kcal/mol) and Ellagic acid (- 7.6 kcal/mol). In comparison, the control drugs Cidofovir and Brincidofovir displayed lower binding affinities, with binding energies of - 6.0 kcal/mol and - 5.1 kcal/mol, respectively. Hydrogen bonds, electrostatic and hydrophobic interactions were the main non-bond interactions between inhibitors and protein active site. The identified compounds were further evaluated using molecular dynamics simulation, density functional theory analysis and ADMET analysis. Molecular dynamics simulations conducted over 200 ns revealed that DNA Pol-Gossypetin complex was not stable, however, Riboflavin and Ellagic acid complexes showed excellent stability indicating them as better DNA Pol inhibitors. The density functional theory analysis exhibited the chemical reactivity of these inhibitor compounds. The ADMET analysis suggested that the top phytochemicals were safe and showed no toxicity. In conclusion, this study has identified Riboflavin and Ellagic acid as potential DNA Pol inhibitors to control MPXV. Further experimental assays and clinical trials are needed to confirm their activity against the disease.

Graphical Abstract

Fatty acid conjugated EPI-X4 derivatives with increased activity and in vivo stability

Dysregulation of the CXCL12/CXCR4 axis is implicated in autoimmune, inflammatory, and oncogenic diseases, positioning CXCR4 as a pivotal therapeutic target. We evaluated optimized variants of the specific endogenous CXCR4 antagonist, EPI-X4, addressing existing challenges in stability and potency. Our structure-activity relationship study investigates the conjugation of EPI-X4 derivatives with long-chain fatty acids, enhancing serum albumin interaction and receptor affinity. Molecular dynamic simulations revealed that the lipid moieties stabilize the peptide-receptor interaction through hydrophobic contacts at the receptor's N-terminus, anchoring the lipopeptide within the CXCR4 binding pocket and maintaining essential receptor interactions. Accordingly, lipidation resulted in increased receptor affinities and antagonistic activities. Additionally, by interacting with human serum albumin lipidated EPI-X4 derivatives displayed sustained stability in human plasma and extended circulation times in vivo. Selected candidates showed significant therapeutic potential in human retinoblastoma cells in vitro and in ovo, with our lead derivative exhibiting higher efficacies compared to its non-lipidated counterpart. This study not only elucidates the optimization trajectory for EPI-X4 derivatives but also underscores the intricate interplay between stability and efficacy, crucial for delineating their translational potential in clinical applications.