In silico and in vitro investigation of dual targeting Prima-1MET as precision therapeutic against lungs cancer

This study emphasizes the explorations of binding of Prima-1MET with two targets, p53 a tumor suppressor protein, and tyrosine kinase of epidermal growth factor receptor. In silico investigations reveal that Prima-1MET showed robust binding with both targets. Molecular docking simulations demonstrated the binding affinity of Prima-1MET with p53 and tyrosine kinase was found to be -38.601 kJ/mol and -38.976 kJ/mol. In addition, the stability of Prima-1MET was explored by molecular dynamics simulation. Prima-1MET attains stability in the binding site of the respective protein till the simulation period is over. Moreover, the free binding energy ΔGbind was calculated by the molecular mechanics Poisson Boltzmann surface area method. The ΔGbind of Prima-1MET with tyrosine kinase was found to be -58.585 ± 0.327 kJ/mol and with p53 it was -35.910 ± 0.335 kJ/mol. Next, cytotoxicity of the Prima-1MET was evaluated using multiple cancer cell lines and the IC50 value were ranging between 4.5 and 30 μM. The cell death was identified by apoptosis assay. Further, the p53 and tyrosine kinase expression was monitored using immunofluorescence techniques, it was found Prima-1MET induces the expression of p53 protein and mimics the level of tyrosine kinase oncogenic target. Also, reactive oxygen species (ROS) and membrane potential activity of Prima-1MET was evaluated by using a lung cancer cell line. A significant decrease in intracellular ROS was observed and resulted in disruption of mitochondrial transmembrane potential. This study uncovers the underlying mechanism of Prima-1MET and could be helpful to design further leads against lung cancers.Communicated by Ramaswamy H. Sarma.

Aggregation of graphene oxide upon the stripping of oxidized debris: An experimental and molecular dynamics simulation study

The most recent structural study of graphene oxide (GO) indicates that the oxidized debris (ODs) adhered to as-prepared GO will strip in certain aquatic settings. The impact of ODs stripping on the characteristics of GO has been widely reported, but its effects on GO aggregation have received less attention. Here, the influence of OD stripping on the GO aggregation property was identified, and the aggregation of as-prepared GO and GO upon OD stripping was compared. Upon ODs stripping, the pKa values of GO shifted from 3.91, 6.25, and 9.84 to 4.54, 6.65, and 10.21, respectively. Further analysis indicated the removal of ODs reduced the net negative charge and improved the hydrophobicity of GO, hence promoting the aggregation of GO. The acceleration of GO-Ca2+-OD aggregate formation was facilitated by the collective effects of ODs stripping, functional group deprotonation, double layer compression, OD bridging, and charge neutralization. The metal ions and stripped ODs attach to GO edges and link GO, which perform like bridges and contribute to further aggregation. In general, the existence of ODs adds complexity to the constructions and characteristics of GO, and it is important to take this into account while evaluating the aggregation characteristic of GO-based materials.

Identification of levomenthol derivatives as potential dipeptidyl peptidase-4 inhibitors: a comparative study with gliptins

Dipeptidyl peptidase-4 (DPP4) inhibitors are a potent therapeutic treatment for type 2 diabetes mellitus (T2DM). There is a family of compounds used as DPP4 inhibitors (DPP4Is) called gliptins. They bind tightly to DPP4 to form an inactive protein-ligand complex. However, there remains a need to identify novel DPP4Is that are more efficacious and safer due to the increasing prevalence of T2DM and the undesirable side effects of gliptins. To identify potential DPP4Is, we screened over 1800 novel compounds in a comparative study with gliptins. We performed dual-factor molecular docking to assess the binding affinity of the compounds to DPP4 and found four compounds with a higher binding affinity to DPP4 than currently used gliptins. The newly identified compounds interacted with the dyad glutamate (GLU205 and GLU206) and tyrosine (TYR662 and TYR666) residues in DPP4's active site. We performed molecular dynamics simulations to determine the stability of the protein-ligand complexes formed by the compounds and DPP4. Furthermore, we examined the toxicity and pharmacological profile of the compounds. The compounds are drug-like, easy to synthesize, and relatively less toxic than gliptins. Collectively, our results suggest that the novel compounds are potential DPP4Is and should be considered for further studies to develop novel antidiabetics.Communicated by Ramaswamy H. Sarma.

Pharmacoinformatics approach for the screening of Kovidra (Bauhinia variegata) phytoconstituents against tumor suppressor protein in triple negative breast cancer

Globally, 2.3 million women were diagnosed with breast cancer, with 6,85000 mortalities in year 2021; making it the world's most prevalent cancer. This growing global burden necessitates a new treatment option, and plant-based medicines offers a promising alternative to conventional cancer treatment. In this work, screening of phytoconstituents of an indigenous therapeutic plant, Bauhinia variegata carried out for potential regulator of tumor suppressor protein p53. Here, an in-silico analysis was employed to develop more effective, pharmaceutically potent small drug-like compounds that target tumor suppressor protein p53. The methanol and aqueous powdered extracts of Bauhinia variegata were prepared and phytochemically evaluated along with antioxidant property evaluation. The LC50 of methanol (325.33 µg/ml) and aqueous extract (361.15 µg/ml) showed their cytotoxic characteristics. Further, GCMS analysis of both the extracts reveals total 57 secondary metabolites. Among these, four lead compounds; compound 1, compound 2, compound 3 and compound 4 were found to have the highest binding ability (-8.15 to -5.40 kcal/mol) with p53. MD simulation and binding free energy validates these findings with highest binding free energy (-67.09 ± 4.87 kcal/mol) towards p53 by the lead phytocompound 2. Selected compounds exhibit excellent pharmacokinetic features and drug-like characteristics. The acute toxicity (LD50) values of the lead phytocompounds ranges from 670 mg/kg to 3100 mg/kg, with toxicity classes of IV and V. As a result, these druggable phytochemicals could serve as potential lead applicants for triple negative breast cancer treatment. However, more in vitro and in vivo research is planned to produce future breast cancer medicine. HIGHLIGHTSScreening of phytoconstituents of an indigenous therapeutic plant, Bauhinia variegata, for potential regulator of tumor suppressor protein p53.The LC50 of methanol (325.33µg/ml) and aqueous extract (361.15µg/ml) showed their cytotoxic characteristics.GCMS analysis of both the extracts reveals total 57 secondary metabolites. Among these, four lead compounds were found to have the highest binding affinity (-8.153 to -5.401 kcal/mol) with tumor suppressor protein p53.MD simulation along with the Prime MM/GBSA binding free energy validates this discovery with highest binding free energy (-67.09 ± 4.87 kcal/mol) towards p53 by the lead compound 2.The acute toxicity (LD50) values of the lead phytocompounds ranges from 670 mg/kg to 3100 mg/kg, with toxicity classes of IV and V.As a result, these druggable phytochemicals could serve as potential lead applicants for triple negative breast cancer treatment.Communicated by Ramaswamy H. Sarma.

Binding properties of selective inhibitors of P323L mutated RdRp of SARS-CoV-2: a combined molecular screening, docking and dynamics simulation study

Since 2019 the SARS-CoV-2 and its variants caused COVID-19, such incidents brought the world in pandemic situation. This happened due to furious mutations in SARS-CoV-2, in which some variants had high transmissibility and infective, this led the virus emerged as virulent and worsened the COVID-19 situation. Among the variants, P323L is one of the important mutants of RdRp in SARS-CoV-2. To inhibit the erroneous function of this mutated RdRp, we have screened 943 molecules against the P323L mutated RdRp with the criteria that the molecules with 90% similar to the structure of remdesivir (control drug) resulted nine molecules. Further, these molecules were evaluated by induced fit docking (IFD) identified two molecules (M2 & M4) which are forming strong intermolecular interactions with the key residues of mutated RdRp and has high binding affinity. Docking score of the M2 and M4 molecules with mutated RdRp are -9.24 and -11.87 kcal/mol, respectively. Further, to understand the intermolecular interactions, conformational stability, the molecular dynamics simulation and binding free energy calculations were performed. The binding free energy values of M2 and M4 molecules with the P323L mutated RdRp complexes are -81.60 and -83.07 kcal/mol, respectively. The results of this in silico study confirm that M4 is a potential molecule; hence, it may be considered as the potential inhibitor of P323L mutated RdRp to treat COVID-19 after clinical investigation.Communicated by Ramaswamy H. Sarma.

Design of a novel multi-epitope vaccine candidate against Chlamydia trachomatis using structural and nonstructural proteins: an immunoinformatics study

Chlamydia trachomatis (C. trachomatis) is an obligate intracellular bacterium which causes eye and sexually transmitted infections. During pregnancy, the bacterium is associated with preterm complications, low weight of neonates, fetal demise and endometritis leading to infertility. The aim of our study was design of a multi-epitope vaccine (MEV) candidate against C. trachomatis. After protein sequence adoption from the NCBI, potential epitopes toxicity, antigenicity, allergenicity, MHC-I and MHC-II binding, cytotoxic T lymphocytes (CTLs), Helper T lymphocytes (HTLs) and interferon-γ (IFN-γ)- induction were predicted. The adopted epitopes were fused together using appropriate linkers. In the next step, the MEV structural mapping and characterization, three-dimensional (3D) structure homology modeling and refinement were also performed. The MEV candidate interaction with the toll-like receptor 4 (TLR4) was also docked. The immune responses simulation was assessed using the C-IMMSIM server. Molecular dynamic (MD) simulation verified the structural stability of the TLR4-MEV complex. The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach demonstrated the MEV high affinity of binding to the TLR4, MHC-I and MHC-II. The MEV construct was also stable and water soluble and had enough antigenicity and lacked allergenicity with stimulation of T cells and B cells and INF-γ release. The immune simulation confirmed acceptable responses of both the humoral and cellular arms. It is proposed that in vitro and in vivo studies are needed to evaluate the findings of this study.Communicated by Ramaswamy H. Sarma.

Interaction analysis of SARS-CoV-2 omicron BA1 and BA2 of RBD with fifty monoclonal antibodies: Molecular dynamics approach

This report provides detailed insights into the interaction of fifty monoclonal antibodies with two recent Omicron variants, BA1 and BA2. It has been observed that numerous mutations in the receptor binding domain (RBD) result in significant structural changes in Omicron, enhancing its ability to mediate viral infections compared to other variants of concern. The following antibodies, namely JX3S304, 7KMG, 7CH4, 7BELCOVOX45, 7CDJ, 7C01, 7JX3S2H14, 6XCA, 7CDI, 7JMO, 7B3O, 6ZER, 6XC7CR3022, JX3S309, 6XC7CC123, 7CM4, 7KMI, 7L7EAZD8895, exhibit a superior binding affinity towards the Spike when compared to the reference CR3022. Four best-docked systems were subjected to further testing through molecular dynamics (MD) simulations. The MM/GBSA free energy for the top-scored complexes of BA1 variant are BA1_JX3S3O4, BA1_7KMI, BA1_7CH4, and BA1_7KMG, with respective values of -56.120 kcal/mol, -41.30 kcal/mol, -17.546 kcal/mol, and -8.527 kcal/mol; and of BA2 variant are BA2_JX3S3O4, BA2_7CM4, BA2_KMG, and BA2_7CH4, with respective values of -40.903 kcal/mol, -23.416 kcal/mol, -17.350 kcal/mol, and -5.460 kcal/mol. Detailed structural/energetic parameters, principal component analysis, and free energy landscape (FEL) studies reveal a significant decrease in antibody resistance due to the disappearance of numerous hydrogen bond interactions and various metastable states. We believe that these crucial mechanistic insights will contribute to breakthroughs in SARS-CoV-2 research.

Understanding cellulose pyrolysis via ab initio deep learning potential field

Comprehensive and dynamic studies of cellulose pyrolysis reaction mechanisms are crucial in designing experiments and processes with enhanced safety, efficiency, and sustainability. The details of the pyrolysis mechanism are not readily available from experiments but can be better described via molecular dynamics (MD) simulations. However, the large size of cellulose molecules challenges accurate ab initio MD simulations, while existing reactive force field parameters lack precision. In this work, precise ab initio deep learning potentials field (DPLF) are developed and applied in MD simulations to facilitate the study of cellulose pyrolysis mechanisms. The formation mechanism and production rate of both valuable and greenhouse products from cellulose at temperatures larger than 1073 K are comprehensively described. This study underscores the critical role of advanced simulation techniques, particularly DLPF, in achieving efficient and accurate understanding of cellulose pyrolysis mechanisms, thus promoting wider industrial applications.

Identification of potential anti-mucor agents by targeting endothelial cell receptor glucose-regulated protein-78 using in silico approach

Mucormycosis is a fungal infection of the sinuses, brain and lungs that is the cause of approximately 50% mortality rate despite the available first-line therapy. Glucose-Regulated Protein 78 (GRP78) is already reported to be a novel host receptor that mediates invasion and damage of human endothelial cells by Rhizopus oryzae and Rhizopus delemar, the most common etiologic species of Mucorales. The expression of GRP78 is also regulated by the levels of iron and glucose in the blood. There are several antifungal drugs in the market but they pose a serious side effect to the vital organs of the body. Therefore, there is an immediate need to discover effective drug molecules having increased efficacy with no side effects. With the help of various computational tools, the current study was attempted to determine potential antimucor agents against GRP78. The receptor molecule GRP78 was screened against 8820 known drugs deposited in DrugBank library using high-throughput virtual screening method. Total top 10 compounds were selected based on the binding energies greater than the reference co-crystal molecule. Furthermore, molecular dynamic (MD) simulations using AMBER were performed to calculate the stability of the top-ranked compounds in the active site of GRP78. After extensive computational studies, we propose that two compounds (CID439153 and CID5289104) have inhibitory potency against mucormycosis and can serve as potential drugs that can form the basis of treating mucormycosis disease.Communicated by Ramaswamy H. Sarma.

Explosive weld joint characteristics of Copper-Tantalum via simulation

This report aims to examine the effects of impact velocity, impact depth, and impact orientation on the Cu-Ta weld joint of the explosive welding process via MD simulation. The findings indicate that the residual shear stress in the welded block mostly increases as the impact velocity rises. The bottom Ta block is more severely distorted than the higher Cu block due to the impact direction. During the tensile test, three stress zones can be identified including the low-stress Cu block, the high-stress Ta block, and the medium-stress weld joint in the middle of the samples. The weld joint position is lower than the median line of the welded block. The Cu-Ta welded block with 500 m/s impact velocities had the highest ultimate tensile strength (UTS) value of 6.49 GPa. With increasing impact depth, the atomic strain level, residual shear stress, and weld joint dimensions all noticeably increase. The Cu-Ta welded block with an impact depth of 7.5 Å has the greatest UTS values, measuring 11.65 GPa, because of its well-crystal structure. Changing the impact orientation does not result in a dramatic change in atomic strain. Orientation (001) vs (001) has the highest strain and stress rates. With an impact orientation of (110) vs. (111), the Cu-Ta welded block gets the highest UTS value of 8.03 GPa compared to other orientations.

Synthesis, in vitro, and in silico studies of morpholine-based thiosemicarbazones as ectonucleotide pyrophosphatase/phosphodiesterase-1 and -3 inhibitors

An extensive range of new biologically active morpholine based thiosemicarbazones derivatives 3a-r were synthesized, characterized by spectral techniques and evaluated as inhibitors of ENPP isozymes. Most of the novel thiosemicarbazones exhibit potent inhibition towards NPP1 and NPP3 isozymes. Compound 3 h was potent inhibitor of NPP1 with IC50 value of 0.55 ± 0.02. However, the most powerful inhibitor of NPP3 was 3e with an IC50 value of 0.24 ± 0.02. Furthermore, Lineweaver-Burk plot for compound 3 h against NPP1 and for compound 3e against NPP3 was devised through enzymes kinetics studies. Molecular docking and in silico studies was also done for analysis of interaction pattern of all newly synthesized compounds. The results were further validated by molecular dynamic (MD) simulation where the stability of conformational transformation of the best protein-ligand complex (3e) were justified on the basis of RMSD and RMSF analysis.

Computation-Guided Rational Design of Cysteine-Less Protein Variants in Engineered hCGL

The engineered human cystathionine-γ-lyase (hCGL) resulting in enhanced activity toward both cysteine and cystine unveils a potential robust antitumor activity. However, the presence of cysteine residues has the potential to induce oligomerization or incorrect disulfide bonding, which may decrease the bioavailability of biopharmaceuticals. Through a meticulous design process targeting the cysteine residues within engineered hCGL, a set of potential beneficial mutants were obtained by virtual screening employing Rosetta and ABACUS. Experimental measurements have revealed that most of the mutants showed increased activity toward both substrates l-Cys and CSSC. Furthermore, mutants C109V and C229D demonstrated Tm value increases of 8.2 and 1.8 °C, respectively. After an 80 min incubation at 60 °C, mutant C229D still maintained high residual activity. Unexpectedly, mutant C109V, displaying activity approximately 2-fold higher than the activity of wild type (WT) for both substrates, showed disappointing instability in plasma, which suggests that computational design still requires further consideration. Analysis of their structure and molecular dynamics (MD) simulation revealed the impact of hydrophobic interaction, hydrogen bonds, and near-attack conformation (NAC) stability on activity and stability. This study acquired information about mutants that exhibit enhanced activity or thermal resistance and serve as valuable guidance for subsequent specific cysteine modifications.

Collision induced unfolding and molecular dynamics simulations of norovirus capsid dimers reveal strain-specific stability profiles

Collision induced unfolding (CIU) is a method used with ion mobility mass spectrometry to examine protein structures and their stability. Such experiments yield information about higher order protein structures, yet are unable to provide details about the underlying processes. That information can however be provided using molecular dynamics simulations. Here, we investigate the gas-phase unfolding of norovirus capsid dimers from the Norwalk and Kawasaki strains by employing molecular dynamics simulations over a range of temperatures, representing different levels of activation, together with CIU experiments. The dimers have highly similar structures, but their CIU reveals different stability that can be explained by the different dynamics that arises in response to the activation seen in the simulations, including a part of the sequence with previously observed strain-specific dynamics in solution. Our findings show how similar protein variants can be examined using mass spectrometric techniques in conjunction with atomistic molecular dynamics simulations to reveal differences in stability as well as differences in how and where unfolding takes place upon activation.

Identification of the potassium-binding site in serotonin transporter

Clearance of serotonin (5-hydroxytryptamine, 5-HT) from the synaptic cleft after neuronal signaling is mediated by serotonin transporter (SERT), which couples this process to the movement of a Na+ ion down its chemical gradient. After release of 5-HT and Na+ into the cytoplasm, the transporter faces a rate-limiting challenge of resetting its conformation to be primed again for 5-HT and Na+ binding. Early studies of vesicles containing native SERT revealed that K+ gradients can provide an additional driving force, via K+ antiport. Moreover, under appropriate conditions, a H+ ion can replace K+. Intracellular K+ accelerates the resetting step. Structural studies of SERT have identified two binding sites for Na+ ions, but the K+ site remains enigmatic. Here, we show that K+ antiport can drive substrate accumulation into vesicles containing SERT extracted from a heterologous expression system, allowing us to study the residues responsible for K+ binding. To identify candidate binding residues, we examine many cation binding configurations using molecular dynamics simulations, predicting that K+ binds to the so-called Na2 site. Site-directed mutagenesis of residues in this site can eliminate the ability of both K+ and H+ to drive 5-HT accumulation into vesicles and, in patch clamp recordings, prevent the acceleration of turnover rates and the formation of a channel-like state by K+ or H+. In conclusion, the Na2 site plays a pivotal role in orchestrating the sequential binding of Na+ and then K+ (or H+) ions to facilitate 5-HT uptake in SERT.

Four-dimensional imaging for cryo-electron microscopy experiments using molecular simulations and manifold learning

Elucidating protein conformational changes is essential because conformational changes are closely related to the functions of proteins. Cryo-electron microscopy (cryo-EM) experiment can be used to reconstruct protein conformational changes via a method that involves using the experimental data (two-dimensional protein images). In this study, a reconstruction method, referred to as the "four-dimensional imaging," was proposed. In our four-dimensional imaging technique, the protein conformational change was obtained using the two-dimensional protein images (the three-dimensional electron density maps used in previously proposed techniques were not used). The protein conformation for each two-dimensional protein image was obtained using our original protocol with molecular dynamics simulations. Using a manifold-learning technique and two-dimensional protein images, the protein conformations were arranged according to the conformational change of the protein. By arranging the protein conformations according to the arrangement of the protein images, four-dimensional imaging is constructed. A simulation for a cryo-EM experiment demonstrated the validity of our four-dimensional imaging technique.

Oxyresveratrol-β-cyclodextrin mitigates streptozotocin-induced Alzheimer's model cognitive impairment, histone deacetylase activity in rats: in silico & in vivo studies

Alzheimer's disease (AD) is associated with cognitive deficits and epigenetic deacetylation that can be modulated by natural products. The role of natural oxyresveratrol-β-cyclodextrin (ORV) on cognition and histone deacetylase activity in AD is unclear. Herein, in-silico docking and molecular dynamics simulation analysis determined that oxyresveratrol potentially targets histone deacetylase-2 (HDAC2). We therefore evaluated the in vivo ameliorative effect of ORV against cognitive deficit, cerebral and hippocampal expression of HDAC in experimental AD rats. Intracerebroventricular injection of STZ (3 mg/kg) induced experimental AD and the rats were treated with low dose (200 mg/kg), high dose (400 mg/kg) of ORV and donepezil (10 mg/kg) for 21 days. The STZ-induced AD caused cognitive and behavioural deficits demonstrated by considerable increases in acetylcholinesterase activity and escape latency compared to sham control. The levels of malondialdehyde (MDA) and HDAC activity were significantly increased in AD disease group comparison to the sham. Interestingly, the ORV reversed the cognitive-behavioural deficit and prominently reduced the MDA and HDAC levels comparable to the effect of the standard drug, donepezil. The findings suggest anti-AD role of ORV via antioxidant effect and inhibition of HDAC in the hippocampal and frontal cortical area of rats for AD.

Iridium(III)-based minor groove binding complexes as DNA photocleavage agents

Transition metal complexes containing the qtpy ligand (2':4,4'':4',4'''-quaterpyridyl) are known to be DNA intercalators or minor groove binders. In this study, new tricationic iridium(III) complexes of qtpy are reported. Both [Ir(bpy)2(qtpy)]3+1 and [Ir(phen)2(qtpy)]3+2 display good water solubility as chloride salts. The complexes possess high-energy excited states, which are quenched in the presence of duplex DNA and even by the mononucleotides guanosine monophosphate and adenosine monophosphate. Further studies reveal that although the complexes bind to quadruplex DNA, they display a preference for duplex structures, which are bound with an order of magnitude higher affinities than their isostructural dicationic RuII-analogues. Detailed molecular dynamics simulations confirm that the complexes are groove binders through the insertion of, predominantly, the qtpy ligand into the minor groove. Photoirradiation of 1 in the presence of plasmid DNA confirms that this class of complexes can function as synthetic photonucleases by cleaving DNA.

Network pharmacology combined with molecular docking and experimental verification to elucidate the effect of flavan-3-ols and aromatic resin on anxiety

This study investigated the potential anxiolytic properties of flavan-3-ols and aromatic resins through a combined computational and experimental approach. Network pharmacology techniques were utilized to identify potential anxiolytic targets and compounds by analyzing protein-protein interactions and KEGG pathway data. Molecular docking and simulation studies were conducted to evaluate the binding interactions and stability of the identified targets. Behavioral tests, including the elevated plus maze test, open field test, light-dark test, actophotometer, and holeboard test, were used to assess anxiolytic activity. The compound-target network analysis revealed complex interactions involving 306 nodes and 526 edges, with significant interactions observed and an average node degree of 1.94. KEGG pathway analysis highlighted pathways such as neuroactive ligand-receptor interactions, dopaminergic synapses, and serotonergic synapses as being involved in anxiety modulation. Docking studies on EGCG (Epigallocatechin gallate) showed binding energies of -9.5 kcal/mol for MAOA, -9.2 kcal/mol for SLC6A4, and -7.4 kcal/mol for COMT. Molecular dynamic simulations indicated minimal fluctuations, suggesting the formation of stable complexes between small molecules and proteins. Behavioral tests demonstrated a significant reduction in anxiety-like behavior, as evidenced by an increased number of entries into and time spent in the open arm of the elevated plus maze test, light-dark test, open field center activity, hole board head dips, and actophotometer beam interruptions (p < 0.05 or p < 0.01). This research provides a comprehensive understanding of the multi-component, multi-target, and multi-pathway intervention mechanisms of flavan-3-ols and aromatic resins in anxiety treatment. Integrated network and behavioral analyses collectively support the anxiolytic potential of these compounds and offer valuable insights for future research in this area.

Cheminformatics approach to identify andrographolide derivatives as dual inhibitors of methyltransferases (nsp14 and nsp16) of SARS-CoV-2

The Covid-19 pandemic outbreak has accelerated tremendous efforts to discover a therapeutic strategy that targets severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to control viral infection. Various viral proteins have been identified as potential drug targets, however, to date, no specific therapeutic cure is available against the SARS-CoV-2. To address this issue, the present work reports a systematic cheminformatic approach to identify the potent andrographolide derivatives that can target methyltransferases of SARS-CoV-2, i.e. nsp14 and nsp16 which are crucial for the replication of the virus and host immune evasion. A consensus of cheminformatics methodologies including virtual screening, molecular docking, ADMET profiling, molecular dynamics simulations, free-energy landscape analysis, molecular mechanics generalized born surface area (MM-GBSA), and density functional theory (DFT) was utilized. Our study reveals two new andrographolide derivatives (PubChem CID: 2734589 and 138968421) as natural bioactive molecules that can form stable complexes with both proteins via hydrophobic interactions, hydrogen bonds and electrostatic interactions. The toxicity analysis predicts class four toxicity for both compounds with LD50 value in the range of 500-700 mg/kg. MD simulation reveals the stable formation of the complex for both the compounds and their average trajectory values were found to be lower than the control inhibitor and protein alone. MMGBSA analysis corroborates the MD simulation result and showed the lowest energy for the compounds 2734589 and 138968421. The DFT and MEP analysis also predicts the better reactivity and stability of both the hit compounds. Overall, both andrographolide derivatives exhibit good potential as potent inhibitors for both nsp14 and nsp16 proteins, however, in-vitro and in vivo assessment would be required to prove their efficacy and safety in clinical settings. Moreover, the drug discovery strategy aiming at the dual target approach might serve as a useful model for inventing novel drug molecules for various other diseases.

Ligand-based drug design against Herpes Simplex Virus-1 capsid protein by modification of limonene through in silico approaches

The pharmacological effects of limonene, especially their derivatives, are currently at the forefront of research for drug development and discovery as well and structure-based drug design using huge chemical libraries are already widespread in the early stages of therapeutic and drug development. Here, various limonene derivatives are studied computationally for their potential utilization against the capsid protein of Herpes Simplex Virus-1. Firstly, limonene derivatives were designed by structural modification followed by conducting a molecular docking experiment against the capsid protein of Herpes Simplex Virus-1. In this research, the obtained molecular docking score exhibited better efficiency against the capsid protein of Herpes Simplex Virus-1 and hence we conducted further in silico investigation including molecular dynamic simulation, quantum calculation, and ADMET analysis. Molecular docking experiment has documented that Ligands 02 and 03 had much better binding affinities (- 7.4 kcal/mol and - 7.1 kcal/mol) to capsid protein of Herpes Simplex Virus-1 than Standard Acyclovir (- 6.5 kcal/mol). Upon further investigation, the binding affinities of primary limonene were observed to be slightly poor. But including the various functional groups also increases the affinities and capacity to prevent viral infection of the capsid protein of Herpes Simplex Virus-1. Then, the molecular dynamic simulation confirmed that the mentioned ligands might be stable during the formation of drug-protein complexes. Finally, the analysis of ADMET was essential in establishing them as safe and human-useable prospective chemicals. According to the present findings, limonene derivatives might be a promising candidate against the capsid protein of Herpes Simplex Virus-1 which ultimately inhibits Herpes Simplex Virus-induced encephalitis that causes interventions in brain inflammation. Our findings suggested further experimental screening to determine their practical value and utility.

Vesicle protrusion induced by antimicrobial peptides suggests common carpet mechanism for short antimicrobial peptides

Short-cationic alpha-helical antimicrobial peptides (SCHAMPs) are promising candidates to combat the growing global threat of antimicrobial resistance. They are short-sequenced, selective against bacteria, and have rapid action by destroying membranes. A full understanding of their mechanism of action will provide key information to design more potent and selective SCHAMPs. Molecular Dynamics (MD) simulations are invaluable tools that provide detailed insights into the peptide-membrane interaction at the atomic- and meso-scale level. We use atomistic and coarse-grained MD to look into the exact steps that four promising SCHAMPs-BP100, Decoralin, Neurokinin-1, and Temporin L-take when they interact with membranes. Following experimental set-ups, we explored the effects of SCHAMPs on anionic membranes and vesicles at multiple peptide concentrations. Our results showed all four peptides shared similar binding steps, initially binding to the membrane through electrostatic interactions and then flipping on their axes, dehydrating, and inserting their hydrophobic moieties into the membrane core. At higher concentrations, fully alpha-helical peptides induced membrane budding and protrusions. Our results suggest the carpet mode of action is fit for the description of SCHAMPs lysis activity and discuss the importance of large hydrophobic residues in SCHAMPs design and activity.

Anti-inflammatory and anti-diabetic properties of indanone derivative isolated from Fernandoa adenophylla in vitro and in silico studies

Fernandoa adenophylla, due to the presence of phytochemicals, has various beneficial properties and is used in folk medicine to treat many conditions. This study aimed to isolate indanone derivative from F. adenophylla root heartwood and assess in-vitro anti-inflammatory and anti-diabetic characteristics at varying concentrations. Heat-induced hemolysis and glucose uptake by yeast cells assays were conducted to evaluate these properties. Besides, docking analyses were performed on four molecular targets. These studies were combined with molecular dynamics simulations to elucidate the time-evolving inhibitory effect of selected inhibitors within the active pockets of the target proteins (COX-1 and COX-2). Indanone derivative (10-100 µM) inhibited the lysis of human red blood cells from 9.12 ± 0.75 to 72.82 ± 4.36% and, at 5-100 µM concentrations, it significantly increased the yeast cells' glucose uptake (5.16 ± 1.28% to 76.59 ± 1.62%). Concluding, the isolated indanone might act as an anti-diabetic agent by interacting with critical amino acid residues of 5' adenosine monophosphate-activated protein kinase (AMPK), and it showed a binding affinity with anti-inflammatory targets COX-1, COX-2, and TNF-α. Besides, the obtained results may help to consider the indanone derivative isolated from F. adenophylla as a promising candidate for drug delivery, subject to outcomes of further in vivo and clinical studies.

High-level production of chitinase by multi-strategy combination optimization in Bacillus licheniformis

In view of the extensive potential applications of chitinase (ChiA) in various fields such as agriculture, environmental protection, medicine, and biotechnology, the development of a high-yielding strain capable of producing chitinase with enhanced activity holds significant importance. The objective of this study was to utilize the extracellular chitinase from Bacillus thuringiensis as the target, and Bacillus licheniformis as the expression host to achieve heterologous expression of ChiA with enhanced activity. Initially, through structural analysis and molecular dynamics simulation, we identified key amino acids to improve the enzymatic performance of chitinase, and the specific activity of chitinase mutant D116N/E118N was 48% higher than that of the natural enzyme, with concomitant enhancements in thermostability and pH stability. Subsequently, the expression elements of ChiA(D116N/E118N) were screened and modified in Bacillus licheniformis, resulting in extracellular ChiA activity reached 89.31 U/mL. Further efforts involved the successful knockout of extracellular protease genes aprE, bprA and epr, along with the gene clusters involved in the synthesis of by-products such as bacitracin and lichenin from Bacillus licheniformis. This led to the development of a recombinant strain, DW2△abelA, which exhibited a remarkable improvement in chitinase activity, reaching 145.56 U/mL. To further improve chitinase activity, a chitinase expression frame was integrated into the genome of DW2△abelA, resulting in a significant increas to 180.26 U/mL. Optimization of fermentation conditions and medium components further boosted shake flask enzyme activity shake flask enzyme activity, achieving 200.28 U/mL, while scale-up fermentation experiments yielded an impressive enzyme activity of 338.79 U/mL. Through host genetic modification, expression optimization and fermentation optimization, a high-yielding ChiA strain was successfully constructed, which will provide a solid foundation for the extracellular production of ChiA.

Elucidating the Mechanism of Metabolism of Cannabichromene by Human Cytochrome P450s

Cannabichromene (CBC) is a nonpsychoactive phytocannabinoid well-known for its wide-ranging health advantages. However, there is limited knowledge regarding its human metabolism following CBC consumption. This research aimed to explore the metabolic pathways of CBC by various human liver cytochrome P450 (CYP) enzymes and support the outcomes using in vivo data from mice. The results unveiled two principal CBC metabolites generated by CYPs: 8'-hydroxy-CBC and 6',7'-epoxy-CBC, along with a minor quantity of 1″-hydroxy-CBC. Notably, among the examined CYPs, CYP2C9 demonstrated the highest efficiency in producing these metabolites. Moreover, through a molecular dynamics simulation spanning 1 μs, it was observed that CBC attains stability at the active site of CYP2J2 by forming hydrogen bonds with I487 and N379, facilitated by water molecules, which specifically promotes the hydroxy metabolite's formation. Additionally, the presence of cytochrome P450 reductase (CPR) amplified CBC's binding affinity to CYPs, particularly with CYP2C8 and CYP3A4. Furthermore, the metabolites derived from CBC reduced cytokine levels, such as IL6 and NO, by approximately 50% in microglia cells. This investigation offers valuable insights into the biotransformation of CBC, underscoring the physiological importance and the potential significance of these metabolites.

Computationally Designed Epitope-Mediated Imprinted Polymers versus Conventional Epitope Imprints for the Detection of Human Adenovirus in Water and Human Serum Samples

Detection of pathogenic viruses for point-of-care applications has attracted great attention since the COVID-19 pandemic. Current virus diagnostic tools are laborious and expensive, while requiring medically trained staff. Although user-friendly and cost-effective biosensors are utilized for virus detection, many of them rely on recognition elements that suffer major drawbacks. Herein, computationally designed epitope-imprinted polymers (eIPs) are conjugated with a portable piezoelectric sensing platform to establish a sensitive and robust biosensor for the human pathogenic adenovirus (HAdV). The template epitope is selected from the knob part of the HAdV capsid, ensuring surface accessibility. Computational simulations are performed to evaluate the conformational stability of the selected epitope. Further, molecular dynamics simulations are executed to investigate the interactions between the epitope and the different functional monomers for the smart design of eIPs. The HAdV epitope is imprinted via the solid-phase synthesis method to produce eIPs using in silico-selected ingredients. The synthetic receptors show a remarkable detection sensitivity (LOD: 102 pfu mL-1) and affinity (dissociation constant (Kd): 6.48 × 10-12 M) for HAdV. Moreover, the computational eIPs lead to around twofold improved binding behavior than the eIPs synthesized with a well-established conventional recipe. The proposed computational strategy holds enormous potential for the intelligent design of ultrasensitive imprinted polymer binders.

Allosteric Modulation of the YAP/TAZ-TEAD Interaction by Palmitoylation and Small-Molecule Inhibitors

The Hippo signaling pathway is a highly conserved signaling network that plays a central role in regulating cellular growth, proliferation, and organ size. This pathway consists of a kinase cascade that integrates various upstream signals to control the activation or inactivation of YAP/TAZ proteins. Phosphorylated YAP/TAZ is sequestered in the cytoplasm; however, when the Hippo pathway is deactivated, it translocates into the nucleus, where it associates with TEAD transcription factors. This partnership is instrumental in regulating the transcription of progrowth and antiapoptotic genes. Thus, in many cancers, aberrantly hyperactivated YAP/TAZ promotes oncogenesis by contributing to cancer cell proliferation, metastasis, and therapy resistance. Because YAP and TAZ exert their oncogenic effects by binding with TEAD, it is critical to understand this key interaction to develop cancer therapeutics. Previous research has indicated that TEAD undergoes autopalmitoylation at a conserved cysteine, and small molecules that inhibit TEAD palmitoylation disrupt effective YAP/TAZ binding. However, how exactly palmitoylation contributes to YAP/TAZ-TEAD interactions and how the TEAD palmitoylation inhibitors disrupt this interaction remains unknown. Utilizing molecular dynamics simulations, our investigation not only provides detailed atomistic insight into the YAP/TAZ-TEAD dynamics but also unveils that the inhibitor studied influences the binding of YAP and TAZ to TEAD in distinct manners. This discovery has significant implications for the design and deployment of future molecular interventions targeting this interaction.

Interaction between Yersinia pestis Ail Outer Membrane Protein and the C-Terminal Domain of Human Vitronectin

Yersinia pestis, the causative agent of plague, is capable of evading the human immune system response by recruiting the plasma circulating vitronectin proteins, which act as a shield and avoid its lysis. Vitronectin recruitment is mediated by its interaction with the bacterial transmembrane protein Ail, protruding from the Y. pestis outer membrane. By using all-atom long-scale molecular dynamic simulations of Ail embedded in a realistic model of the bacterial membrane, we have shown that vitronectin forms a stable complex, mediated by interactions between the disordered moieties of the two proteins. The main amino acids driving the complexation have also been evidenced, thus favoring the possible rational design of specific peptides which, by inhibiting vitronectin recruitment, could act as original antibacterial agents.

Structural Rearrangements of Pigeon Cryptochrome 4 Undergoing a Complete Redox Cycle

Cryptochrome is currently the major contender of a protein to underpin magnetoreception, the ability to sense the Earth's magnetic field. Among various types of cryptochromes, cryptochrome 4 has been identified as the likely magnetoreceptor in migratory birds. All-atom molecular dynamics (MD) studies have offered first insights into the structural dynamics of cryptochrome but are limited to a short time scale due to large computational demands. Here, we employ coarse-grained MD simulations to investigate the emergence of long-lived states and conformational changes in pigeon cryptochrome 4. Our coarse-grained simulations complete the picture by permitting observation on a significantly longer time scale. We observe conformational transitions in the phosphate-binding loop of pigeon cryptochrome 4 upon activation and identify prominent motions in residues 440-460, suggesting a possible role as a signaling state of the protein or as a gated interaction site for forming protein complexes that might facilitate downstream processes. The findings highlight the importance of considering longer time scales in studying cryptochrome dynamics and magnetoreception. Coarse-grained MD simulations offer a valuable tool to unravel the complex behavior of cryptochrome proteins and shed new light on the mechanisms underlying their role in magnetoreception. Further exploration of these conformational changes and their functional implications may contribute to a deeper understanding of the molecular mechanisms of magnetoreception in birds.

Dehydration of Lipid Membranes Drives Redistribution of Cholesterol Between Lateral Domains

Cholesterol-rich lipid rafts are found to facilitate membrane fusion, central to processes like viral entry, fertilization, and neurotransmitter release. While the fusion process involves local, transient membrane dehydration, the impact of reduced hydration on cholesterol's structural organization in biological membranes remains unclear. Here, we employ confocal fluorescence microscopy and atomistic molecular dynamics simulations to investigate cholesterol behavior in phase-separated lipid bilayers under controlled hydration. We unveiled that dehydration prompts cholesterol release from raft-like domains into the surrounding fluid phase. Unsaturated phospholipids undergo more significant dehydration-induced structural changes and lose more hydrogen bonds with water than sphingomyelin. The results suggest that cholesterol redistribution is driven by the equalization of biophysical properties between phases and the need to satisfy lipid hydrogen bonds. This underscores the role of cholesterol-phospholipid-water interplay in governing cholesterol affinity for a specific lipid type, providing a new perspective on the regulatory role of cell membrane heterogeneity during membrane fusion.

Mechanisms of the Drug Penetration Enhancer Propylene Glycol Interacting with Skin Lipid Membranes

Very few drugs have the necessary physicochemical properties to cross the skin's main permeability barrier, the stratum corneum (SC), in sufficient amounts. Propylene glycol (PG) is a chemical penetration enhancer that could be included in topical formulations in order to overcome the barrier properties of the skin and facilitate the transport of drugs across it. Experiments have demonstrated that PG increases the mobility and disorder of SC lipids and may extract cholesterol from the SC, but little is known about the molecular mechanisms of drug permeation enhancement by PG. In this work, we have performed molecular dynamics (MD) simulations to investigate the molecular-level effects of PG on the structure and properties of model SC lipid bilayers. The model bilayers were simulated in the presence of PG concentrations over the range of 0-100% w/w PG, using both an all-atom and a united atom force field. PG was found to localize in the hydrophilic headgroup regions at the bilayer interface, to occupy the lipid-water hydrogen-bonding sites, and to slightly increase lipid tail disorder in a concentration-dependent manner. We showed with MD simulation that PG enhances the permeation of small molecules such as water by interacting with the bilayer interface; the results of our study may be used to guide the design of formulations for transdermal drug delivery with enhanced skin permeation, as well as topical formulations and cosmetic products.

Design and Synthesis of Novel Macrocyclic Derivatives as Potent and Selective Cyclin-Dependent Kinase 7 Inhibitors

The duality of function (cell cycle regulation and gene transcription) of cyclin-dependent kinase 7 (CDK7) makes it an attractive oncology target and the discovery of CDK7 inhibitors has been a long-term pursuit by academia and pharmaceutical companies. However, achieving selective leading compounds is still difficult owing to the similarities among the ATP binding pocket. Herein, we detail the design and synthesis of a series of macrocyclic derivatives with pyrazolo[1,5-a]-1,3,5-triazine core structure as potent and selective CDK7 inhibitors. The diverse manners of macrocyclization led to distinguished selectivity profiles of the CDK family. Molecular dynamics (MD) simulation explained the binding difference between 15- and 16-membered macrocyclic compounds. Further optimization generated compound 37 exhibiting good CDK7 inhibitory activity and high selectivity over other CDKs. This work clearly demonstrated macrocyclization is a versatile method to finely tune the selectivity profile of small molecules and MD simulation can be a valuable tool in prioritizing designs of the macrocycle.

Hydrophobic Clusters Regulate Surface Hydration Dynamics of Bacillus subtilis Lipase A

The surface hydration diffusivity of Bacillus subtilis Lipase A (BSLA) has been characterized by low-field Overhauser dynamic nuclear polarization (ODNP) relaxometry using a series of spin-labeled constructs. Sites for spin-label incorporation were previously designed via an atomistic computational approach that screened for surface exposure, reflective of the surface hydration comparable to other proteins studied by this method, as well as minimal impact on protein function, dynamics, and structure of BSLA by excluding any surface site that participated in greater than 30% occupancy of a hydrogen bonding network within BSLA. Experimental ODNP relaxometry coupling factor results verify the overall surface hydration behavior for these BSLA spin-labeled sites similar to other globular proteins. Here, by plotting the ODNP parameters of relative diffusive water versus the relative bound water, we introduce an effective "phase-space" analysis, which provides a facile visual comparison of the ODNP parameters of various biomolecular systems studied to date. We find notable differences when comparing BSLA to other systems, as well as when comparing different clusters on the surface of BSLA. Specifically, we find a grouping of sites that correspond to the spin-label surface location within the two main hydrophobic core clusters of the branched aliphatic amino acids isoleucine, leucine, and valine cores observed in the BSLA crystal structure. The results imply that hydrophobic clustering may dictate local surface hydration properties, perhaps through modulation of protein conformations and samplings of the unfolded states, providing insights into how the dynamics of the hydration shell is coupled to protein motion and fluctuations.

Origin of Enantioselectivity in Engineered Cytochrome c-Catalyzed Carbon-Radical FePP Hydrolysis Revealed Using QM/MM (ABEEM Polarizable Force Field) and MD Simulations

The origin of highly efficient asymmetric aminohydroxylation of styrene catalyzed by engineered cytochrome c is investigated by the developed Atom-Bond Electronegativity Equalization Method polarizable force field (ABEEM PFF), which is a combined outcome of electronic and steric effects. Model molecules were used to establish the charge parameters of the ABEEM PFF, for which the bond-stretching and angle-bending parameters were obtained by using a combination of modified Seminario and scan methods. The interactions between carbon-radical Fe-porphyrin (FePP) and waters are simulated by molecular dynamics, which shows a clear preference for the pre-R over the pre-S. This preference is attributed to the hydrogen-bond between the mutated 100S and 101P residues as well as van der Waals interactions, enforcing a specific conformation of the carbon-radical FePP complex within the binding pocket. Meanwhile, the hydrogen-bond between water and the nitrogen atom in the active intermediate dictates the stereochemical outcome. Quantum mechanics/molecular mechanics (QM/MM (ABEEM PFF)) and free-energy perturbation calculations elucidate that the 3RTS is characterized by sandwich-like structure among adjacent amino acid residues, which exhibits greater stability than crowed arrangement in 3STS and enables the R enantiomer to form more favorably. Thus, this study provides mechanistic insight into the catalytic reaction of hemoproteins.

Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations

The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.

Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences

Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.

Effects of Proline on Internal Friction in Simulated Folding Dynamics of Several Alanine-Based α-Helical Peptides

We have studied in silico the effect of proline, a model cosolvent, on local and global friction coefficients in (un)folding of several typical alanine-based α-helical peptides. Local friction is related to dwell times of a single, ensemble-averaged hydrogen bond (HB) within each peptide. Global friction is related to energy dissipated in a series of configurational changes of each peptide experienced by increasing the number of HBs during folding. Both of these approaches are important in relation to future atomic force microscopic-based measurements of internal friction via force-clamp single-molecule force spectroscopy. Molecular dynamics (MD) simulations for six peptides, namely, ALA5, ALA8, ALA15, ALA21, (AAQAA)3, and H2N-GN(AAQAA)2G-COONH2, have been conducted at 2 and 5 M proline solutions in water. Using previously obtained MD data for these peptides in pure water as well as upgraded theoretical models, we obtained variations of local and global internal friction coefficients as a function of solution viscosity. The results showed the substantial role of proline in stabilizing the folded state and slowing the overall folding dynamics. Consequently, larger friction coefficients were obtained at larger viscosities. The local and global internal friction, i.e., respective, friction coefficients approximated to zero viscosity, was also obtained. The evolution of friction coefficients with viscosity was weakly dependent on the number of concurrent folding pathways but was rather dominated by a stabilizing effect of proline on the folded states. Obtained values of local and global internal friction showed qualitatively similar results and a clear dependency on the structure of the studied peptide.

Models of the Three-Component Bilayer of Stratum Corneum: A Molecular Simulation Study

The construction of the stratum corneum (SC) is crucial to the problems of transdermal drug delivery. SC consists of the keratinocyte layers and the lipid matrix surrounding it. Among them, the lipid matrix is the barrier for many exogenous molecules, mainly composed of ceramides (CERs), free fatty acids (FFA), and cholesterol (CHOL). In this work, we developed single-component (CERs, CER-NS, and CER-EOS) and six three-component models, and each model was simulated by using the GROMOS-54A7 force field. Short-period phase (SPP) and long-period phase (LPP) systems were established separately, and area per lipid (APL), thickness, order of carbon chain (SCD), and density distribution were analyzed. The transition of CER-NS and CER-EOS in LPP was observed. The results of hydrogen bonds in the lipid systems indicated that a strong hydrogen-bond network was formed between the skin-lipid bilayers. Umbrella sampling method simulations were performed to calculate the free energy change of ethanol moving into the skin-lipid bilayer. The results revealed that ethanol molecules pulled some water molecules into the membrane when they passed through SPP-1. Our findings provided some insights and models of the stratum corneum that could be used for the subsequent mechanism of macromolecule permeation through membranes in drugs, cosmetics, and so on.

[Computer-aided prediction and molecular mechanism investigation of active components in compound Kushen injection inhibiting p21-activated kinase 1]

Targeting p21-activated kinase 1 (PAK1) is a novel strategy for pancreatic cancer treatment. Compound Kushen injection contains many anti-pancreatic cancer components, but the specific targets are unknown. In this study, 14α-hydroxymatrine, an active component of Kushen injection, was found to possess high binding free energy with the allosteric site of PAK1 by molecular docking based virtual screening. Molecular dynamics simulations suggested that 14α-hydroxymatrine caused the α1 and α2 helices of the allosteric site of PAK1 to extend outward to form a deep allosteric regulatory pocket. Meanwhile, 14α-hydroxymatrine induced the β-folding region at the adenosine triphosphate (ATP)-binding pocket of PAK1 to close inward, resulting in the ATP-binding pocket in a "semi-closed" state which caused the inactivation of PAK1. After removal of 14α-hydroxymatrine, PAK1 showed a tendency to change from the inactive conformation to the active conformation. We supposed that 14α-hydroxymatrine of compound Kushen injection might be a reversible allosteric inhibitor of PAK1. This study used modern technologies and methods to study the active components of traditional Chinese medicine, which laid a foundation for the development and utilization of natural products and the search for new treatments for pancreatic cancer.

The search for novel proline analogs for viscosity reduction and stabilization of highly concentrated monoclonal antibody solutions

Administration of monoclonal antibodies (mAbs) is currently focused on subcutaneous injection associated with increased patient adherence and reduced treatment cost, leading to sustainable healthcare. The main bottleneck is low volume that can be injected, requiring highly concentrated mAb solutions. The latter results in increased solution viscosity with pronounced mAb aggregation propensity because of intensive protein-protein interactions. Small molecule excipients have been proposed to restrict the protein-protein interactions, contributing to reduced viscosity. The aim of the study was to discover novel compounds that reduce the viscosity of highly concentrated mAb solution. First, the chemical space of proline analogs was explored and 35 compounds were determined. Viscosity measurements revealed that 18 proline analogs reduced the mAb solution viscosity similar to or more than proline. The compounds forming both electrostatic and hydrophobic interactions with mAb reduced the viscosity of the formulation more efficiently without detrimentally effecting mAb physical stability. A correlation between the level of interaction and viscosity-reducing effect was confirmed with molecular dynamic simulations. Structure rigidity of the compounds and aromaticity contributed to their viscosity-reducing effect, dependent on molecule size. The study results highlight the novel proline analogs as an effective approach in viscosity reduction in development of biopharmaceuticals for subcutaneous administration.

The antimicrobial fibupeptide lugdunin forms water-filled channel structures in lipid membranes

Recently, a novel cyclo-heptapeptide composed of alternating D,L-amino acids and a unique thiazolidine heterocycle, called lugdunin, was discovered, which is produced by the nasal and skin commensal Staphylococcus lugdunensis. Lugdunin displays potent antimicrobial activity against a broad spectrum of Gram-positive bacteria, including challenging-to-treat methicillin-resistant Staphylococcus aureus (MRSA). Lugdunin specifically inhibits target bacteria by dissipating their membrane potential. However, the precise mode of action of this new class of fibupeptides remains largely elusive. Here, we disclose the mechanism by which lugdunin rapidly destabilizes the bacterial membrane potential using an in vitro approach. The peptide strongly partitions into lipid compositions resembling Gram-positive bacterial membranes but less in those harboring the eukaryotic membrane component cholesterol. Upon insertion, lugdunin forms hydrogen-bonded antiparallel β-sheets by the formation of peptide nanotubes, as demonstrated by ATR-FTIR spectroscopy and molecular dynamics simulations. These hydrophilic nanotubes filled with a water wire facilitate not only the translocation of protons but also of monovalent cations as demonstrated by voltage-clamp experiments on black lipid membranes. Collectively, our results provide evidence that the natural fibupeptide lugdunin acts as a peptidic channel that is spontaneously formed by an intricate stacking mechanism, leading to the dissipation of a bacterial cell's membrane potential.

PURA syndrome-causing mutations impair PUR-domain integrity and affect P-body association

Mutations in the human PURA gene cause the neurodevelopmental PURA syndrome. In contrast to several other monogenetic disorders, almost all reported mutations in this nucleic acid-binding protein result in the full disease penetrance. In this study, we observed that patient mutations across PURA impair its previously reported co-localization with processing bodies. These mutations either destroyed the folding integrity, RNA binding, or dimerization of PURA. We also solved the crystal structures of the N- and C-terminal PUR domains of human PURA and combined them with molecular dynamics simulations and nuclear magnetic resonance measurements. The observed unusually high dynamics and structural promiscuity of PURA indicated that this protein is particularly susceptible to mutations impairing its structural integrity. It offers an explanation why even conservative mutations across PURA result in the full penetrance of symptoms in patients with PURA syndrome.

Exceptional stability of ultrasmall cubic copper metal nanoclusters - a molecular dynamics study

The fabrication of shape-selective coinage metal nanoclusters (MNCs) has promising applications due to their exceptional physical and chemical molecule-like properties. However, the stability of the specific geometry of the nanoclusters, such as their cubic shapes, is unclear and has been unraveled by assessing the nanoclusters' interactions with different environments. In this work, we investigate the morphological stability of cubic structured, coinage metal nanoclusters of varying sizes ranging from 14 to 1099 atoms. The impact of solvent environments like water and the presence of ionic liquids (IL) on the stabilization of the MNCs were assessed using molecular dynamics (MD) simulations. In general, smaller MNCs composed of less than 256 atoms encountered structural distortion easily compared to the larger ones, which preserved their cubic morphology with minimal surface aberrations in water. However, in the presence of 4M 1-butyl-1,1,1-trimethyl ammonium methane sulfonate [N1114][C1SO3] IL solution, the overall cubic shape of the MNCs was successfully preserved. Strikingly, it is observed that in contrast to the noble MNCs like Au and Ag, the cubic morphology for Cu MNCs with sizes less than 256 atoms exhibited significant stability even in the absence of IL.

Oxidative Defect Detection Within Free and Packed DNA Systems: A Quantum Mechanical/Molecular Mechanics (QM/ MM) Approach

Base excision repair enzymes (BERs) detect and repair oxidative DNA damage with efficacy despite the small size of the defects and their often only minor structural impact. A charge transfer (CT) model for rapid scanning of DNA stretches has been evoked to explain the high detection rate in the face of numerous, small lesions. The viability of CT DNA defect detection is explored via hybrid QM/MM computational studies that leverage the accuracy of quantum mechanics (QM) for a region of interest and the descriptive power of molecularmechanics (MM) for the remainder of the system. We find that the presence of an oxidative lesion lowers theredox free energy of oxidation by approximately 1.0 eV regardless of DNA compaction (free DNA versus packed DNA in nucleosome core particles) and damage location indicating the high feasibility of a CT-based process for defect detection in DNA.

Bacillus velezensis iturins inhibit the hemolytic activity of Staphylococcus aureus

Bovine mastitis caused by S. aureus has a major economic impact on the dairy sector. With the crucial need for new therapies, anti-virulence strategies have gained attention as alternatives to antibiotics. Here we aimed to identify novel compounds that inhibit the production/activity of hemolysins, a virulence factor of S. aureus associated with mastitis severity. We screened Bacillus strains obtained from diverse sources for compounds showing anti-hemolytic activity. Our results demonstrate that lipopeptides produced by Bacillus spp. completely prevented the hemolytic activity of S. aureus at certain concentrations. Following purification, both iturins, fengycins, and surfactins were able to reduce hemolysis caused by S. aureus, with iturins showing the highest anti-hemolytic activity (up to 76% reduction). The lipopeptides showed an effect at the post-translational level. Molecular docking simulations demonstrated that these compounds can bind to hemolysin, possibly interfering with enzyme action. Lastly, molecular dynamics analysis indicated general stability of important residues for hemolysin activity as well as the presence of hydrogen bonds between iturins and these residues, with longevous interactions. Our data reveals, for the first time, an anti-hemolytic activity of lipopeptides and highlights the potential application of iturins as an anti-virulence therapy to control bovine mastitis caused by S. aureus.

Effect of GM1 concentration change on plasma membrane: molecular dynamics simulation and analysis

Ganglioside GM1 is a class of glycolipids predominantly located in the nervous system. Comprising a ceramide anchor and an oligosaccharide chain containing sialic acid, GM1 plays a pivotal role in various cellular processes, including signal transduction, cell adhesion, and membrane organization. Moreover, GM1 has been implicated in the pathogenesis of several neurological disorders, such as Parkinson's disease, Alzheimer's disease, and stroke. In this study, by creating a neural cell model membrane simulation system and employing rigorous molecular models, we utilize a coarse-grained molecular dynamics approach to explore the structural and dynamic characteristics of multi-component neuronal plasma membranes at varying GM1 ganglioside concentrations. The simulation results reveal that as GM1 concentration increases, a greater number of hydrogen bonds form between GM1 molecules, resulting in the formation of larger clusters, which leads to reduced membrane fluidity, increased lipid ordering, decreased membrane thickness and surface area and higher levels of GM1 dissociation. Through a meticulous analysis, while considering GM1's structural attributes, we offer valuable insights into the structural and dynamic traits of the cell membrane. This study provides a robust methodology for exploring membrane characteristics and enhances our comprehension of GM1 molecules, serving as a resource for both experimental and computational researchers in this field.

Symmetry-breaking malachite green as a near-infrared light-activated fluorogenic photosensitizer for RNA proximity labeling

Cellular RNA is asymmetrically distributed in cells and the regulation of RNA localization is crucial for proper cellular functions. However, limited chemical tools are available to capture dynamic RNA localization in complex biological systems with high spatiotemporal resolution. Here, we developed a new method for RNA proximity labeling activated by near-infrared (NIR) light, which holds the potential for deep penetration. Our method, termed FAP-seq, utilizes a genetically encoded fluorogen activating protein (FAP) that selectively binds to a set of substrates known as malachite green (MG). FAP binding restricts the rotation of MG and rapidly activates its fluorescence in a wash-free manner. By introducing a monoiodo modification to MG, we created a photosensitizer (MG-HI) with the highest singlet oxygen generation ability among various MG derivatives, enabling both protein and RNA proximity labeling in live cells. New insights are provided in the transcriptome analysis with FAP-seq, while a deeper understanding of the symmetry-breaking structural arrangement of FAP-MG-HI was obtained through molecular dynamics simulations. Overall, our wash-free and NIR light-inducible RNA proximity labeling method (FAP-seq) offers a powerful and versatile approach for investigating complex mechanisms underlying RNA-related biological processes.

The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study

Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.

Computational insights into the binding modes, keto-enol tautomerization and stereo-electronically controlled decarboxylation of oxaloacetate in the active site of macrophomate synthase

Oxaloacetic acid (OAA) is a β-ketocarboxylic acid, which plays an important role as an intermediate in some metabolic pathways, including the tricarboxylic acid cycle, gluconeogenesis and fatty acid biosynthesis. Animal studies have indicated that supplementing oxaloacetic acid shows an increase of lifespan and other substantial health benefits including mitochondrial DNA protection, and protection of retinal, neural and pancreatic tissues. Most of the chemical transformations of OAA in the metabolic pathways have been extensively studied; however, the understanding of decarboxylation of OAA at the atomic level is relatively lacking. Here, we carried out MD simulations and combined quantum mechanical/molecular mechanical (QM/MM) calculations as an example to systematically elucidate the binding modes, keto-enol tautomerization and decarboxylation of OAA in the active site of macrophomate synthase (MPS), which is a Mg(II)-dependent bifunctional enzyme that catalyzes both the decarboxylation of OAA and [4+2] cycloaddition of 2-pyrone with the decarboxylated intermediate of OAA (pyruvate enolate). On the basis of our calculations, it was found that the Mg2+-coordinated oxaloacetate may exist in enol forms and keto forms. The four keto forms can be transformed into each other by simply rotating the C2-C3 single bond, nevertheless, the keto-enol tautomerization strictly requires the assistance of pocket water molecules. In addition, the decarboxylation is stereo-electronically controlled, i.e., it is the relative orientation of the terminal carboxyl anion that determines the rate of decarboxylation. As such, the chemistry of oxaloacetate in the active site of MPS is complex. On one hand, the most stable binding mode (K-I) may undergo enol-keto tautomerization to isomerize to the enol form, which may further react with the second substrate; on the other hand, K-I may isomerize to another binding mode K-II to proceed decarboxylation to generate pyruvate enolate and CO2. Starting from K-I, the enol-keto tautomerization corresponds to a barrier of 16.2 kcal mol-1, whereas the decarboxylation is associated with an overall barrier of 19.7 kcal mol-1. These findings may provide useful information for understanding the chemistry of OAA and the catalysis of related enzymes, and they are basically in agreement with the available experimental kinetic data.

Predicting FFAR4 agonists using structure-based machine learning approach based on molecular fingerprints

Free Fatty Acid Receptor 4 (FFAR4), a G-protein-coupled receptor, is responsible for triggering intracellular signaling pathways that regulate various physiological processes. FFAR4 agonists are associated with enhancing insulin release and mitigating the atherogenic, obesogenic, pro-carcinogenic, and pro-diabetogenic effects, normally associated with the free fatty acids bound to FFAR4. In this research, molecular structure-based machine-learning techniques were employed to evaluate compounds as potential agonists for FFAR4. Molecular structures were encoded into bit arrays, serving as molecular fingerprints, which were subsequently analyzed using the Bayesian network algorithm to identify patterns for screening the data. The shortlisted hits obtained via machine learning protocols were further validated by Molecular Docking and via ADME and Toxicity predictions. The shortlisted compounds were then subjected to MD Simulations of the membrane-bound FFAR4-ligand complexes for 100 ns each. Molecular analyses, encompassing binding interactions, RMSD, RMSF, RoG, PCA, and FEL, were conducted to scrutinize the protein-ligand complexes at the inter-atomic level. The analyses revealed significant interactions of the shortlisted compounds with the crucial residues of FFAR4 previously documented. FFAR4 as part of the complexes demonstrated consistent RMSDs, ranging from 3.57 to 3.64, with minimal residue fluctuations 5.27 to 6.03 nm, suggesting stable complexes. The gyration values fluctuated between 22.8 to 23.5 nm, indicating structural compactness and orderliness across the studied systems. Additionally, distinct conformational motions were observed in each complex, with energy contours shifting to broader energy basins throughout the simulation, suggesting thermodynamically stable protein-ligand complexes. The two compounds CHEMBL2012662 and CHEMBL64616 are presented as potential FFAR4 agonists, based on these insights and in-depth analyses. Collectively, these findings advance our comprehension of FFAR4's functions and mechanisms, highlighting these compounds as potential FFAR4 agonists worthy of further exploration as innovative treatments for metabolic and immune-related conditions.

Gaussian accelerated molecular dynamics simulations facilitate prediction of the permeability of cyclic peptides

Despite their widespread use as therapeutics, clinical development of small molecule drugs remains challenging. Among the many parameters that undergo optimization during the drug development process, increasing passive cell permeability (i.e., log(P)) can have some of the largest impact on potency. Cyclic peptides (CPs) have emerged as a viable alternative to small molecules, as they retain many of the advantages of small molecules (oral availability, target specificity) while being highly effective at traversing the plasma membrane. However, the relationship between the dominant conformations that typify CPs in an aqueous versus a membrane environment and cell permeability remain poorly characterized. In this study, we have used Gaussian accelerated molecular dynamics (GaMD) simulations to characterize the effect of solvent on the free energy landscape of lariat peptides, a subset of CPs that have recently shown potential for drug development (Kelly et al., JACS 2021). Differences in the free energy of lariat peptides as a function of solvent can be used to predict permeability of these molecules, and our results show that permeability is most greatly influenced by N-methylation and exposure to solvent. Our approach lays the groundwork for using GaMD as a way to virtually screen large libraries of CPs and drive forward development of CP-based therapeutics.