ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB

Phosphopeptide binding to the N-SH2 domain of tyrosine phosphatase SHP2 correlates with the unzipping of its central β-sheet

SHP2 is a tyrosine phosphatase that plays a regulatory role in multiple intracellular signaling cascades and is known to be oncogenic in certain contexts. In the absence of effectors, SHP2 adopts an autoinhibited conformation with its N-SH2 domain blocking the active site. Given the key role of N-SH2 in regulating SHP2, this domain has been extensively studied, often by X-ray crystallography. Using a combination of structural analyses and molecular dynamics (MD) simulations we show that the crystallographic environment can significantly influence the structure of the isolated N-SH2 domain, resulting in misleading interpretations. As an orthogonal method to X-ray crystallography, we use a combination of NMR spectroscopy and MD simulations to accurately determine the conformation of apo N-SH2 in solution. In contrast to earlier reports based on crystallographic data, our results indicate that apo N-SH2 in solution primarily adopts a conformation with a fully zipped central β-sheet, and that partial unzipping of this β-sheet is promoted by binding of either phosphopeptides or even phosphate/sulfate ions.

Cysteine hyperoxidation rewires communication pathways in the nucleosome and destabilizes the dyad

Gene activity is tightly controlled by reversible chemical modifications called epigenetic marks, which are of various types and modulate gene accessibility without affecting the DNA sequence. Despite an increasing body of evidence demonstrating the role of oxidative-type modifications of histones in gene expression regulation, there remains a complete absence of structural data at the atomistic level to understand the molecular mechanisms behind their regulatory action. Owing to μs time-scale MD simulations and protein communication networks analysis, we describe the impact of histone H3 hyperoxidation (i.e., S-sulfonylation) on the nucleosome core particle dynamics. Our results reveal the atomic-scale details of the intrinsic structural networks within the canonical histone core and their perturbation by hyperoxidation of the histone H3 C110. We show that this modification involves local rearrangements of the communication networks and destabilizes the dyad, and that one modification is enough to induce a maximal structural signature. Our results suggest that cysteine hyperoxidation in the nucleosome core particle might favor its disassembly.

Avian cryptochrome 4 binds superoxide

Flavin-binding cryptochromes are blue-light sensitive photoreceptors that have been implicated with magnetoreception in some species. The photocycle involves an intra-protein photo-reduction of the flavin cofactor, generating a magnetosensitive radical pair, and its subsequent re-oxidation. Superoxide (O2 • - ) is generated in the re-oxidation with molecular oxygen. The resulting O2 • - -containing radical pairs have also been hypothesised to underpin various magnetosensitive traits, but due to fast spin relaxation when tumbling in solution would require immobilisation. We here describe our insights in the binding of superoxide to cryptochrome 4 from C. livia based on extensive all-atom molecular dynamics studies and density-functional theory calculations. The positively charged "crypt" region that leads to the flavin binding pocket transiently binds O2 • - at 5 flexible binding sites centred on arginine residues. Typical binding times amounted to tens of nanoseconds, but exceptional binding events extended to several hundreds of nanoseconds and slowed the rotational diffusion, thereby realising rotational correlation times as large as 1 ns. The binding sites are particularly efficient in scavenging superoxide escaping from a putative generation site close to the flavin-cofactor, possibly implying a functional relevance. We discuss our findings in view of a potential magnetosensitivity of biological flavin semiquinone/superoxide radical pairs.

Conformational plasticity of SpyCas9 induced by AcrIIA4 and AcrIIA2: Insights from molecular dynamics simulation

CRISPR-Cas9 systems constitute bacterial adaptive immune systems that protect against phage infections. Bacteriophages encode anti-CRISPR proteins (Acrs) that mitigate the bacterial immune response. However, the structural basis for their inhibitory actions from a molecular perspective remains elusive. In this study, through microsecond atomistic molecular dynamics simulations, we demonstrated the remarkable flexibility of Streptococcus pyogenes Cas9 (SpyCas9) and its conformational adaptability during interactions with AcrIIA4 and AcrIIA2. Specifically, we demonstrated that the binding of AcrIIA4 and AcrIIA2 to SpyCas9 induces a conformational rearrangement that causes spatial separation between the nuclease and cleavage sites, thus making the endonuclease inactive. This separation disrupts the transmission of signals between the protospacer adjacent motif recognition and nuclease domains, thereby impeding the efficient processing of double-stranded DNA. The simulation also reveals that AcrIIA4 and AcrIIA2 cause different structural variations of SpyCas9. Our research illuminates the precise mechanisms underlying the suppression of SpyCas9 by AcrIIA4 and AcrIIA2, thus presenting new possibilities for controlling genome editing with higher accuracy.

Multiscale calculations reveal new insights into the reaction mechanism between KRASG12C and α, β-unsaturated carbonyl of covalent inhibitors

Utilizing α,β-unsaturated carbonyl group as Michael acceptors to react with thiols represents a successful strategy for developing KRASG12C inhibitors. Despite this, the precise reaction mechanism between KRASG12C and covalent inhibitors remains a subject of debate, primarily due to the absence of an appropriate residue capable of deprotonating the cysteine thiol as a base. To uncover this reaction mechanism, we first discussed the chemical reaction mechanism in solvent conditions via density functional theory (DFT) calculation. Based on this, we then proposed and validated the enzymatic reaction mechanism by employing quantum mechanics/molecular mechanics (QM/MM) calculation. Our QM/MM analysis suggests that, in biological conditions, proton transfer and nucleophilic addition may proceed through a concerted process to form an enolate intermediate, bypassing the need for a base catalyst. This proposed mechanism differs from previous findings. Following the formation of the enolate intermediate, solvent-assisted tautomerization results in the final product. Our calculations indicate that solvent-assisted tautomerization is the rate-limiting step in the catalytic cycle under biological conditions. On the basis of this reaction mechanism, the calculated kinact/ki for two inhibitors is consistent well with the experimental results. Our findings provide new insights into the reaction mechanism between the cysteine of KRASG12C and the covalent inhibitors and may provide valuable information for designing effective covalent inhibitors targeting KRASG12C and other similar targets.

Anilino-1,4-naphthoquinones as potent mushroom tyrosinase inhibitors: in vitro and in silico studies

Tyrosinase, a pivotal enzyme in melanin synthesis, is a primary target for the development of depigmenting agents. In this work, in vitro and in silico techniques were employed to identify novel tyrosinase inhibitors from a set of 12 anilino-1,4-naphthoquinone derivatives. Results from the mushroom tyrosinase activity assay indicated that, among the 12 derivatives, three compounds (1, 5, and 10) demonstrated the most significant inhibitory activity against mushroom tyrosinase, surpassing the effectiveness of the kojic acid. Molecular docking revealed that all studied derivatives interacted with copper ions and amino acid residues at the enzyme active site. Molecular dynamics simulations provided insights into the stability of enzyme-inhibitor complexes, in which compounds 1, 5, and particularly 10 displayed greater stability, atomic contacts, and structural compactness than kojic acid. Drug likeness prediction further strengthens the potential of anilino-1,4-naphthoquinones as promising candidates for the development of novel tyrosinase inhibitors for the treatment of hyperpigmentation disorders.

Impact of an irreversible β-galactosylceramidase inhibitor on the lipid profile of zebrafish embryos

Krabbe disease is a sphingolipidosis characterized by the genetic deficiency of the acid hydrolase β-galactosylceramidase (GALC). Most of the studies concerning the biological role of GALC performed on Krabbe patients and Galc-deficient twitcher mice (an authentic animal model of the disease) indicate that the pathogenesis of this disorder is the consequence of the accumulation of the neurotoxic GALC substrate β-galactosylsphingosine (psychosine), ignoring the possibility that this enzyme may exert a wider biological impact. Indeed, limited information is available about the effect of GALC downregulation on the cell lipidome in adult and developing organisms. The teleost zebrafish (Danio rerio) has emerged as a useful platform to model human genetic diseases, including sphingolipidoses, and two GALC co-orthologs have been identified in zebrafish (galca and galcb). Here, we investigated the effect of the competitive and irreversible GALC inhibitor β-galactose-cyclophellitol (GCP) on the lipid profile of zebrafish embryos. Molecular modelling indicates that GCP can be sequestered in the catalytic site of the enzyme and covalently binds human GALC, and the zebrafish Galca and Galcb proteins in a similar manner. Accordingly, GCP inhibits the β-galactosylceramide hydrolase activity of zebrafish in vitro and in vivo, leading to significant alterations of the lipidome of zebrafish embryos. These results indicate that the lack of GALC activity deeply affects the lipidome during the early stages of embryonic development, and thereby provide insights into the pathogenesis of Krabbe disease.

Development of multi-epitope mRNA vaccine against Clostridioides difficile using reverse vaccinology and immunoinformatics approaches

Clostridioides difficile (C. difficile), as the major pathogen of diarrhea in healthcare settings, has become increasingly prevalent within community populations, resulting in significant morbidity and mortality. However, the therapeutic options for Clostridioides difficile infection (CDI) remain limited, and as of now, no authorized vaccine is available to combat this disease. Therefore, the development of a novel vaccine against C. difficile is of paramount importance. In our study, the complete proteome sequences of 118 strains of C. difficile were downloaded and analyzed. We found four antigenic proteins that were highly conserved and can be used for epitope identification. We designed two vaccines, WLcd1 and WLcd2, that contain the ideal T-cell and B-cell epitopes, adjuvants, and the pan HLA DR-binding epitope (PADRE) sequences. The biophysical and chemical assessments of these vaccine candidates indicated that they were suitable for immunogenic applications. Molecular docking analyses revealed that WLcd1 bonded with higher affinity to Toll-like receptors (TLRs) than WLcd2. Furthermore, molecular dynamics (MD) simulations, performed using Gmx_MMPBSA v1.56, confirmed the binding stability of WLcd1 with TLR2 and TLR4. The preliminary findings suggested that this multi-epitope vaccine could be a promising candidate for protection against CDI; however, experimental studies are necessary to confirm these predictions.

Electrochemical biosensor based on the integration of maple leaf-like gold nanocrystal and truncated aptamer for detection of α-amanitin with high sensitivity, selectivity and rapidity

This study reports the fabrication of three-dimensional gold nanocrystals as sensing material in the presence of l-glutathion and high-performance aptamer with 20 bases of α-amanitin via truncation and optimization of along aptamer. The resulting maple leaf-like gold nanocrystal (ML-Au) exhibits an improved catalytic activity due to more exposed high-index facets. The use of truncated aptamer increases the sensitivity by 15 times and reduces the reaction time by two times compared with those of original aptamer. An α-amanitin electrochemical biosensor constructed by integrating ML-Au nanocrystals with truncated aptamer exhibits high sensitivity, selectivity and rapidity. An increase of the α-amanitin concentration in the range of 1 × 10-14-1 × 10-9 M causes a linear decrease in the amperometric current with a limit of detection of 2.9 × 10-15 M (S/N = 3). The proposed analytical method is satisfactorily used for electrochemical sensing of α-amanitin in urine and wild mushroom samples.

Colorimetric aptasensor utilizing MOF-235 with exceptional peroxidase-like activity for the detection of oxytetracycline residues in raw milk

In this work, a simple, convenient and cost-effective colorimetric aptasensor was successfully constructed for the detection of antibiotic residues in raw milk based on the property that aptamer (Apt) synergistically enhances the catalase-like activity of MOF-235. Under optimised conditions, the proposed colorimetric aptasensor exhibited a wide detection range (15-1500 nM) with a low detection limit (6.92 nM). Furthermore, the proposed aptasensor demonstrated high selectivity, good resistance to interference and storage stability. The proposed aptasensor was validated by spiking recovery in camel milk, cow milk and goat milk with satisfactory recoveries, which demonstrated the great potential of the aptasensor for further application in real food samples, and also suggested that MOF-235 can be used as a potential universal platform to build a sensitive detection platform for other targets.

A potential therapeutic agent for the treatment of hyperuricemia and gout: 3,4-Dihydroxy-5-nitrobenzaldehyde phenylthiosemicarbazide

Uric acid, the metabolic product of purines, relies on xanthine oxidase (XOD) for production. XOD is a target for the development of drugs for hyperuricemia (HUA) and gout. Currently, treatment options remain limited for gout patients. 3, 4-Dihydroxy-5-nitrobenzaldehyde (DHNB) is a derivative of the natural product protocatechualdehyde with good biological activity. In this work, we identify a DHNB thiosemicarbazide class of compounds that targets XOD. 3,4-Dihydroxy-5-nitrobenzaldehyde phenylthiosemicarbazone can effectively inhibit XOD activity (IC50 value: 0.0437 μM) and exhibits a mixed inhibitory effect. In a mouse model of acute hyperuricemia, a moderate dose (10 mg/kg.w) of 3,4-dihydroxy-5-nitrobenzaldehyde phenylthiosemicarbazide effectively controlled the serum uric acid content and significantly inhibited serum XOD activity. In addition, 3,4-Dihydroxy-5-nitrobenzaldehyde phenylthiosemicarbazide showed favorable safety profiles, and mice treated with the target compound did not show any symptoms of general toxicity following a single dose of 500 mg/kg. In the allopurinol group, 50 % of the mice died. These results provide a structural framework and mechanism of XOD inhibition that may facilitate the design of hyperuricemia and gout treatments.

In silico-based screening of natural products as potential inhibitors of SARS-CoV-2 macrodomain 1

The rapid spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) worldwide has led to over 600 million cases of coronavirus disease 2019 (COVID-19). Identifying effective molecules that can counteract the virus is imperative. SARS-CoV-2 macrodomain 1 (Mac1) represents a promising antiviral drug target. In this study, we predicted potential inhibitors of SARS-CoV-2 Mac1 from natural products using in silico-based screening. Based on the high-resolution crystal structure of Mac1 bound to its endogenous ligand ADP-ribose (ADPr), we first performed a docking-based virtual screening of Mac1 inhibitors against a natural product library and obtained five representative compounds (MC1-MC5) by clustering analysis. All five compounds were stably bound to Mac1 during 500 ns long molecular dynamics simulations. The binding free energy of these compounds to Mac1 was calculated using molecular mechanics generalized Born surface area and further refined with localized volume-based metadynamics. The results demonstrated that both MC1 (-9.8 ± 0.3 kcal/mol) and MC5 (-9.6 ± 0.3 kcal/mol) displayed more favorable affinities to Mac1 with respect to ADPr (-8.9 ± 0.3 kcal/mol), highlighting their potential as potent SARS-CoV-2 Mac1 inhibitors. Overall, this study provides potential SARS-CoV-2 Mac1 inhibitors, which may pave the way for developing effective therapeutics for COVID-19.Communicated by Ramaswamy H. Sarma.

Study on multiscale structures and digestibility of cassava starch and medium-chain fatty acids complexes using molecular simulation techniques

Effect of complexation of three medium-chain fatty acids (octanoic, decylic and lauric acid, OA, DA and LA, respectively) on structural characteristics, physicochemical properties and digestion behaviors of cassava starch (CS) was investigated. Current study indicated that LA was more easily to combine with CS (complex index 88.9%), followed by DA (80.9%), which was also consistent with their corresponding complexed lipids content. Following the investigation of morphology, short-range ordered structure, helical structure, crystalline/amorphous region and fractal dimension of the various complexes, all cassava starch-fatty acids complexes (CS-FAs) were characterized with a flaked morphology rather than a round morphology in native starch (control CS). X-ray diffraction demonstrated that all CS-FAs had a V-type crystalline structure, and nuclear magnetic resonance spectroscopy confirmed that the complexes made from different fatty acids displayed similar V6 or V7 type polymorphs. Interestingly, small-angle X-ray scattering analysis revealed that α value became greater following increased carbon chain length of fatty acids, indicating the formation of a more ordered fractal structure in the aggregates. Changes in rheological parameters G' and G'' indicated that starch complexed with fatty acids was more likely to form a gel network, but difference among three CS-FAs complexes was significant, which might be contributed to their corresponding hydrophobicity and hydrophilicity raised from individual fatty acids. Importantly, digestion indicated that CS-LA complexes had the lowest hydrolysis degree, followed by the greatest RS content, indicating the importance of chain length of fatty acids for manipulating the fine structure and functionality of the complexes.

Effect of quercetin on the protein-substrate interactions in SIRT6: Insight from MD simulations

SIRT6 is of interest for its promising effect in the treatment of aging-related diseases. Studies have shown quercetin (QUE) and its derivatives have varying degrees of effect on the catalytic effect of SIRT6. In the research, the effect of QUE on the protein-substrate interaction in the SIRT6-mediated mono-ADP ribosylation system was investigated by conventional molecular dynamics (MD) simulations combined with MM/PBSA binding free energy calculations. The results show that QUE can bind stably to SIRT6 with the binding energy of -22.8 kcal/mol and further affect the atomic interaction between SIRT6 and NAD+ (or H3K9), resulting in an increased affinity between SIRT6-NAD+ and decreased SIRT6-H3K9 binding capacity. At the same time, the binding of QUE can also alter some structural characteristics of the protein, with large shifts occurring in the residue regions involving the N-terminal (residues 1-27), Rossmann fold regions (residues 55-92), and ZBD (residues 164-179). Thus, QUE shows great potential as a scaffold for the design of novel potent SIRT6 modulators.

In silico identification of cyclosporin derivatives as potential inhibitors for RdRp of rotavirus by molecular docking and molecular dynamic studies

Rotavirus is one of the most common gastrointestinal viral diseases. Till date, there are only two vaccines available in the markets, which are specifically to be administered to young babies. In this study, VP1 RdRp is selected as potential target to carry out inhibitory activities. Cyclosporin A (Cys A) derivatives were designed via FBDD, pharmacokinetics, molecular docking, molecular dynamics (MD) simulation and molecular mechanics generalized born surface area was applied on these compounds. The results from these investigations were analyzed and it was found that the considered derivatives in this study were nontoxic and docking results revealed that the derivatives made some important bonds inside the active site of the receptors within a catalytic triad (Serine-Histidine-Aspartate). After analyzing the mean values of root mean square density (RMSD), root mean square fluctuation (RMSF), radius of gyration (RoG) and solvent accessible surface area (SASA) at 100 ns MD simulation of the selected compounds, it was found that compound 1 exhibits RMSD of 0.74 ± 0.10 Å, RMSF of 0.85 ± 0.15 Å, RoG of 16.45 ± 0.40 Å, SASA of 66.55 ± 0.35 nm2 and ΔGbind of -32.76 ± 0.02 kcal/mol. Therefore, the study revealed that amongst the designed and reported compounds, compound 1 was more stable within the active region of the RdRp and also this compound possesses lower binding free energy as compared to other selected compounds and Cys A as well.Communicated by Ramaswamy H. Sarma.

Effects of theacrine on the astringency of EGCG by affecting salivary protein - EGCG interactions through different molecular mechanisms

(-)-Epigallocatechin gallate (EGCG) and theacrine are involved in imparting tea with its astringent and bitter tastes. This study investigated the effect of theacrine on the astringency of EGCG and its molecular mechanism. Sensory evaluation was used to study the astringent intensities of EGCG solutions in the presence and absence of various concentrations of theacrine. The results indicated a considerable increase in the astringency values of EGCG-theacrine solutions compared with those of EGCG solutions alone. Furthermore, dynamic light scattering (DLS) and molecular dynamics simulations (MD) were to explore the interaction mechanisms. The results revealed that theacrine increased the particle size of EGCG-proline-rich proteins (PRPs) aggregates with that of EGCG and PRPs alone. MD revealed that theacrine potentially acted as a bridge between EGCG and PRPs, promoting their interaction and intensifying the EGCG's astringency. However, theacrine could not bridge two or more mucins owing to the substantial spatial structure of mucin.

Improved physicochemical properties and in vitro digestion of walnut oil microcapsules with soy protein isolate and highly oxidized konjac glucomannan as wall materials

This study investigated the effect of the oxidation degrees of oxidized konjac glucomannan (OKGM) on the encapsulation efficiency (EE), physicochemical and in vitro digestive properties of soy protein isolate (SPI)-based microcapsules walnut oil using experimental and computational approaches. Microcapsules had the highest EE when the ratio of OKGM and SPI to oil was 2.5:1. With increasing the oxidation degree of OKGM, the EE of microcapsules was increased and the hygroscopicity was decreased. Molecular dynamics simulation results showed that SPI/oil/highly OKGM had relatively low binding energy (-4.03 × 106 kJ/mol) and strong electrostatic interactions, which may contribute to a higher EE and lower hygroscopicity of microcapsules, respectively. The oxidative stability of the oil was markedly improved by SPI and OKGM, and microcapsules prepared with SPI and highly OKGM had the highest in vitro digestion. This study provided theoretical support for broadening the application of microcapsules prepared with SPI and OKGM.

Insights into the ATP / GTP selectivity of a GTPase, adenylosuccinate synthetase from Leishmania donovani

Many GTPases have been shown to utilize ATP too as the phosphoryl donor. Both GTP and ATP are important molecules in the cellular environments and play multiple and discrete functional role within the cells. In our present study, we showed that one of the purine metabolic enzymes Adenylosuccinate synthetase from Leishmania donovani (LdAdSS) which belongs to the BioD-superfamily of GTPases can also carry out the catalysis by hydrolysing ATP instead of its cognate substrate GTP albeit with less efficiency. Biochemical and biophysical studies indicated its ability to bind to ATP too but at a higher concentration of ATP compared to that of GTP. Sequence analysis and molecular dynamic simulations suggested that residues of the switch loop and the G4-G5 (593SAXD596) connected motif of LdAdSS plays a role in determining the nucleotide specificity. Though the crucial interaction between Asp596 and the nucleotide is broken when ATP is bound, interactions between the Ala594 and the adenine ring of ATP could still hold ATP in the GTP binding site. The results of the present study suggested that though LdAdSS is GTP specific, it still shows ATP hydrolysing activity.

Assessment of Trachyspermum ammi essential oil against Aspergillus flavus, aflatoxin B1 contamination, and post-harvest quality of Sorghum bicolor

The present investigation explored the antifungal effectiveness of Trachyspermum ammi essential oil (TAEO) against Aspergillus flavus, aflatoxin B1 (AFB1) contamination, and its mechanism of action using biochemical and computational approaches. The GC-MS result revealed the chemical diversity of TAEO with the highest percentage of γ-terpinene (39 %). The TAEO exhibited minimum inhibitory concentration against A. flavus growth (0.5 µL/mL) and AFB1 (0.4 µL/mL) with radical scavenging activity (IC50 = 2.13 µL/mL). The mechanism of action of TAEO was associated with the alteration in plasma membrane functioning, antioxidative defense, and carbon source catabolism. The molecular dynamic result shows the multi-regime binding of γ-terpinene with the target proteins (Nor1, Omt1, and Vbs) of AFB1 biosynthesis. Furthermore, TAEO exhibited remarkable in-situ protection of Sorghum bicolor seed samples against A. flavus and AFB1 contamination and protected the nutritional deterioration. Hence, the study recommends TAEO as a natural antifungal agent for food protection against A. flavus mediated biodeterioration.

Computation-guided engineering of distal mutations in an artificial enzyme

Artificial enzymes are valuable biocatalysts able to perform new-to-nature transformations with the precision and (enantio-)selectivity of natural enzymes. Although they are highly engineered biocatalysts, they often cannot reach catalytic rates akin those of their natural counterparts, slowing down their application in real-world industrial processes. Typically, their designs only optimise the chemistry inside the active site, while overlooking the role of protein dynamics on catalysis. In this work, we show how the catalytic performance of an already engineered artificial enzyme can be further improved by distal mutations that affect the conformational equilibrium of the protein. To this end, we subjected a specialised artificial enzyme based on the lactococcal multidrug resistance regulator (LmrR) to an innovative algorithm that quickly inspects the whole protein sequence space for hotpots which affect the protein dynamics. From an initial predicted selection of 73 variants, two variants with mutations distant by more than 11 Å from the catalytic pAF residue showed increased catalytic activity towards the new-to-nature hydrazone formation reaction. Their recombination displayed a 66% higher turnover number and 14 °C higher thermostability. Microsecond time scale molecular dynamics simulations evidenced a shift in the distribution of productive enzyme conformations, which are the result of a cascade of interactions initiated by the introduced mutations.

Peptide Drug Design Using Alchemical Free Energy Calculation: An Application and Validation on Agonists of Ghrelin Receptor

With recent large-scale applications and validations, the relative binding free energy (RBFE) calculated using alchemical free energy methods has been proven to be an accurate measure to probe the binding of small-molecule drug candidates. On the other hand, given the flexibility of peptides, it is of great interest to find out whether sufficient sampling could be achieved within the typical time scale of such calculation, and a similar level of accuracy could be reached for peptide drugs. However, the systematic evaluation of such calculations on protein-peptide systems has been less reported. Most reported studies of peptides were restricted to a limited number of data points or lacking experimental support. To demonstrate the applicability of the alchemical free energy method for protein-peptide systems in a typical real-world drug discovery project, we report an application of the thermodynamic integration (TI) method to the RBFE calculation of ghrelin receptor and its peptide agonists. Along with the calculation, the synthesis and in vitro EC50 activity of relamorelin and 17 new peptide derivatives were also reported. A cost-effective criterion to determine the data collection time was proposed for peptides in the TI simulation. The average of three TI repeats yielded a mean absolute error of 0.98 kcal/mol and Pearson's correlation coefficient (R) of 0.77 against the experimental free energy derived from the in vitro EC50 activity, showing good repeatability of the proposed method and a slightly better agreement than the results obtained from the arbitrary time frames up to 20 ns. Although it is limited by having one target and a deduced binding pose, we hope that this study can add some insights into alchemical free energy calculation of protein-peptide systems, providing theoretical assistance to the development of peptide drugs.

Molecular Insights into Converting Hydroxide Adenosyltransferase into Halogenase

Halogenation plays a unique role in the design of agrochemicals. Enzymatic halogenation reactions have attracted great attention due to their excellent specificity and mild reaction conditions. S-adenosyl-l-methionine (SAM)-dependent halogenases mediate the nucleophilic attack of halide ions (X-) to SAM to produce 5'-XDA. However, only 11 SAM-dependent fluorinases and 3 chlorinases have been reported, highlighting the desire for additional halogenases. SAM-dependent hydroxide adenosyltransferase (HATase) has a similar reaction mechanism as halogenases but uses water as a substrate instead of halide ions. Here, we explored a HATase from the thermophile Thermotoga maritima MSB8 and transformed it into a halogenase. We identified a key dyad W8L/V71T for the halogenation reaction. We also obtained the best performing mutants for each halogenation reaction: M1, M2 and M4 for Cl-, Br- and I-, respectively. The M4 mutant retained the thermostability of HATase in the iodination reaction at 80 °C, which surpasses the natural halogenase SalL. QM/MM revealed that these mutants bind halide ions with more suitable angles for nucleophilic attack of C5' of SAM, thus conferring halogenation capabilities. Our work achieved the halide ion specificity of halogenases and generated thermostable halogenases for the first time, which provides new opportunities to expand the halogenase repertoire from hydroxylase.

New Mechanism for the Apoptosis of Human Neuroblastoma Cells by the Interaction between Fluorene-9-Bisphenol and the G Protein-Coupled Estrogen Receptor 1

Fluorene-9-bisphenol (BHPF) is an emerging contaminant. Presently, there is no report on its interaction with G protein-coupled estrogen receptor 1 (GPER). By using an integrated toxicity research scenario that combined theoretical study with experimental methods, BHPF was found to inhibit the GPER-mediated effect via direct receptor binding. Molecular dynamics simulations found that Trp2726.48 and Glu2756.51 be the key amino acids of BHPF binding with GPER. Moreover, the calculation indicated that BHPF was a suspected GPER inhibitor, which neither can activate GPER nor is able to form water channels of GPER. The role of two residues was successfully verified by following gene knockout and site-directed mutagenesis assays. Further in vitro assays showed that BHPF could attenuate the increase in intracellular concentration of free Ca2+ induced by G1-activated GPER. Besides, BHPF showed an enhanced cytotoxicity compared with G15, indicating that BHPF might be a more potent GPER inhibitor than G15. In addition, a statistically significant effect on the mRNA level of GPER was observed for BHPF. In brief, the present study proposes that BHPF be a GPER inhibitor, and its GPER molecular recognition mechanism has been revealed, which is of great significance for the health risk and assessment of BHPF.

Exploring the potential mechanisms of Danshen against COVID-19 via network pharmacology analysis and molecular docking

Danshen, a prominent herb in traditional Chinese medicine (TCM), is known for its potential to enhance physiological functions such as blood circulation, immune response, and resolve blood stasis. Despite the effectiveness of COVID-19 vaccination efforts, some individuals still face severe complications post-infection, including pulmonary fibrosis, myocarditis arrhythmias and stroke. This study employs a network pharmacology and molecular docking approach to investigate the potential mechanisms underlying the therapeutic effects of candidate components and targets from Danshen in the treatment of complications in COVID-19. Candidate components and targets from Danshen were extracted from the TCMSP Database, while COVID-19-related targets were obtained from Genecards. Venn diagram analysis identified common targets. A Protein-Protein interaction (PPI) network and gene enrichment analysis elucidated potential therapeutic mechanisms. Molecular docking evaluated interactions between core targets and candidate components, followed by molecular dynamics simulations to assess stability. We identified 59 potential candidate components and 123 targets in Danshen for COVID-19 treatment. PPI analysis revealed 12 core targets, and gene enrichment analysis highlighted modulated pathways. Molecular docking showed favorable interactions, with molecular dynamics simulations indicating high stability of key complexes. Receiver operating characteristic (ROC) curves validated the docking protocol. Our study unveils candidate compounds, core targets, and molecular mechanisms of Danshen in COVID-19 treatment. These findings provide a scientific foundation for further research and potential development of therapeutic drugs.

Discordant results between Xpert MTB/RIF assay and Bactec MGIT 960 culture system regarding the detection of rifampin-resistant Mycobacterium tuberculosis isolates in Wenzhou, China

This study aimed to assess the possible causes of discordant results between Xpert MTB/RIF (Xpert) and Bactec MGIT 960 Culture System (MGIT960) regarding rifampicin (RIF) susceptibility in Mycobacterium tuberculosis. Patients with previous RIF-resistant tuberculosis who were admitted to Wenzhou Central Hospital from January 2020 to December 2022 were enrolled. The isolates obtained from these patients were subjected to RIF susceptibility tests using Xpert and MGIT960, and the minimum inhibitory concentration (MIC) of RIF was determined by the MYCOTB MIC plate test. Additionally, molecular docking and molecular dynamics (MD) simulations were performed to evaluate the binding efficacy of rpoB and RIF based on rpoB mutations detected in the isolates with discordant RIF susceptibility results. A total of 28 isolates with discordant RIF susceptibility test results were detected, 15 of them were RIF susceptible with MICs ≤ 0.5 µg/mL. Twelve out of 15 isolates contained borderline RIF resistance-associated mutations [L430P (n = 6), H445N (n = 6)], 1 isolate had D435Y and Q429H double mutation, and the remaining 2 isolates had a silent (Q432Q) mutation. Compared with the affinity of RIF toward the wild type (WT) (-45.83 kcal/mol) by MD, its affinity toward L452P (-55.52 kcal/mol), D435Y (-47.39 kcal/mol), L430P (approximately -69.72 kcal/mol), H445N (-49.53 kcal/mol), and Q429H (-55.67 kcal/mol) increased. Borderline RIF resistance-associated mutations were the main cause for the discordant RIF susceptibility results between Xpert and MGIT960, and the mechanisms of the resistance need further investigated.IMPORTANCEThis study is aimed at assessing discordant results between Xpert MTB/RIF (Xpert) assay and Bactec MGIT 960 Culture System (MGIT960) regarding the detection of rifampicin (RIF)-resistant Mycobacterium tuberculosis isolates in Wenzhou, China. The discordant results of RIF between these two assays were mainly caused by borderline RIF resistance-associated mutations, subsequently by silent mutations of rpoB. Borderline RIF resistance- associated mutations detected in our study were demonstrated to not be affected by the affinity of rpoB and RIF by molecular dynamics, and the mechanism of resistance was needed to be clarified. For the discordant results of RIF by Xpert and MGIT960 that occurred, rpoB DNA sequencing was recommended to investigate its association with resistance to RIF.

Novel benzimidazole angiotensin receptor blockers with anti-SARS-CoV-2 activity equipotent to that of nirmatrelvir: computational and enzymatic studies

BACKGROUND

Pre-existing hypertensive pathologies are linked to worsened outcomes in patients infected with severe-acute respiratory syndrome (SARS-CoV-2). Sartans, a family of anti-hypertensive angiotensin receptor blockers (ARBs), reduce morbidity and mortality in coronavirus 2019 patients by aeffecting angiotensin-converting enzyme-2 (ACE2). The purpose of this study was to determine the anti-SARS-CoV-2viral and antihypertensive abilities of nirmatrelvir, and commercially available (e.g. candesartan, losartan, and losartan carboxylic (Exp3174))) and newly synthesized sartans (e.g. benzimidazole-N-biphenyl carboxyl (ACC519C), and benzimidazole-N-biphenyl tetrazole (ACC519T)), compared to nirmatrelvir, the antiviral component of paxlovid.

RESEARCH DESIGN AND METHODS

Surface plasmon resonance-binding (SPR) and enzymatic studies were used to determine drug effect on ACE2, antiviral abilities were evaluated using Vero E6 cells infected with SARS-CoV-2 and antihypertensive effects were investigated using rabbit iliac arteries contracted to angiotensin II-dose response.

RESULTS

The antiviral activity of benzimidazole-based candesartan and ACC519C were equipotent with nirmatrelvir (95% inhibition), while the imidazole-based losartan and E×p3174and ACC519T were less potent (75-80 and 50%, respectively%) and E×p3174was the least effective (50%). Furthermore, SPR analysis revealed a high binding affinity of sartans for ACE2 protein. Moreover, candesartan and nirmatrelvir exhibited the greatest inhibitory and cytopathic effect (3.96%) when given in combination rather than individually (6.10% and 5.08%, respectively). In addition,The ACE2 enzymatic activity assays demonstrated enhancement or inhibition properties of novel sartans toward ACE2 enzyme. Lastly, ACC519T potently reduced angiotensin II-mediated contraction while nirmatrelvir and ACC519T(2) had no effect.

CONCLUSION

This study describes the discovery of a new class of benzimidazole-based sartans that drastically inhibit SARS-CoV-2, which could be attributed, at least in part, to their interaction with ACE2.

Energy Landscapes and Structural Ensembles of Glucagon-like Peptide-1 Monomers

While GLP-1 and its analogues are important pharmaceutical agents in the treatment of type 2 diabetes and obesity, their susceptibility to aggregate into amyloid fibrils poses a significant safety issue. Many factors may contribute to the aggregation propensity, including pH. While it is known that the monomeric structure of GLP-1 has a strong impact on primary nucleation, probing its diverse structural ensemble is challenging. Here, we investigated the monomer structural ensembles at pH 3, 4, and 7.5 using state-of-the-art computational methods in combination with experimental data. We found significant stabilization of β-strand structures and destabilization of helical structures at lower pH, correlating with observed aggregation lag times, which are lower under these conditions. We further identified helical defects at pH 4, which led to the fastest observed aggregation, in agreement with our far-UV circular dichroism data. The detailed atomistic structures that result from the computational studies help to rationalize the experimental results on the aggregation propensity of GLP-1. This work provides a new insight into the pH-dependence of monomeric structural ensembles of GLP-1 and connects them to experimental observations.

De novo design of proteins housing excitonically coupled chlorophyll special pairs

Natural photosystems couple light harvesting to charge separation using a 'special pair' of chlorophyll molecules that accepts excitation energy from the antenna and initiates an electron-transfer cascade. To investigate the photophysics of special pairs independently of the complexities of native photosynthetic proteins, and as a first step toward creating synthetic photosystems for new energy conversion technologies, we designed C2-symmetric proteins that hold two chlorophyll molecules in closely juxtaposed arrangements. X-ray crystallography confirmed that one designed protein binds two chlorophylls in the same orientation as native special pairs, whereas a second designed protein positions them in a previously unseen geometry. Spectroscopy revealed that the chlorophylls are excitonically coupled, and fluorescence lifetime imaging demonstrated energy transfer. The cryo-electron microscopy structure of a designed 24-chlorophyll octahedral nanocage with a special pair on each edge closely matched the design model. The results suggest that the de novo design of artificial photosynthetic systems is within reach of current computational methods.

The archaeal highly thermostable GH35 family β-galactosidase DaβGal has a unique seven domain protein fold

The most extensively studied β-d-galactosidases (EC3.2.1.23) belonging to four glycoside hydrolase (GH) families 1, 2, 35, and 42 are widely distributed among Bacteria, Archaea and Eukaryotes. Here, we report a novel GH35 family β-galactosidase from the hyperthermophilic Thermoprotei archaeon Desulfurococcus amylolyticus (DaβGal). Unlike fungal monomeric six-domain β-galactosidases, the DaβGal enzyme is a dimer; it has an extra jelly roll domain D7 and three composite domains (D4, D5, and D6) that are formed by the distantly located polypeptide chain regions. The enzyme possesses a high specificity for β-d-galactopyranosides, and its distinguishing feature is the ability to cleave pNP-β-d-fucopyranoside. DaβGal efficiently catalyzes the hydrolysis of lactose at high temperatures, remains stable and active at 65 °С, and retains activity at 95 °С with a half-life time value equal to 73 min. These properties make archaeal DaβGal a more attractive candidate for biotechnology than the widely used fungal β-galactosidases.

Design, synthesis, and analysis of macrobicyclic peptides for targeting the Gαi protein

Bicyclic peptides are important chemical tools that can function, for example, as bioactive ligands switching on/off signaling pathways mediated by guanine nucleotide-binding proteins as bicycles are more broadly applicable. Despite their relevance in medicinal chemistry, the synthesis of such peptides is challenging, and the final yield is highly dependent on the chemical composition and physicochemical properties of the scaffold. We recently discovered novel, state-specific peptide modulators targeting the Gαi protein, namely, GPM-2/GPM-3, by screening a one-bead-two-compound combinatorial library. A more detailed analysis, including sequence alignments and computer-assisted conformational studies based on the hit compounds, revealed the new peptide 10 as a potential macrobicyclic Gαi ligand sharing high sequence similarity to the known Gαi modulators. The Gαs protein was included in this study for comparison and to unravel the criteria for the specificity of modulator binding to Gαi versus Gαs. This work provides in-depth computer-assisted experimental studies for the analysis of novel macrobicyclic, library-derived Gαi protein ligands. The sequence and structural comparison of 10 with the lead compounds GPM-2 and GPM-3 reveals the importance of the size and amino acid composition of one ring of the bicyclic system and suggests features enhancing the binding affinity of the peptides to the Gαi protein.

Molecular dynamics of structural effects of reactive carbonyl species derivate of lipid peroxidation on bovine serum albumin

BACKGROUND

Serum albumin is the most abundant protein in the Mammalia blood plasma at where plays a decisive role in the transport wide variety of hydrophobic ligands. BSA undergoes oxidative modifications like the carbonylation by the reactive carbonyl species (RCSs) 4-hydroxy-2-nonenal (HNE), 4 hydroxy-2-hexenal (HHE), malondialdehyde (MDA) and 4-oxo-2-nonenal (ONE), among others. The structural and functional changes induced by protein carbonylation have been associated with the advancement of neurodegenerative, cardiovascular, metabolic and cancer diseases.

METHODS

To elucidate structural effects of protein carbonylation with RCSs on BSA, parameters for six new non-standard amino acids were designated and molecular dynamics simulations of its mono‑carbonylated-BSA systems were conducted in the AMBER force field. Trajectories were evaluated by RMSD, RMSF, PCA, RoG and SASA analysis.

RESULTS

An increase in the conformational instability for all proteins modified with local changes were observed, without significant changes on the BSA global three-dimensional folding. A more relaxed compaction level and major solvent accessible surface area for modified systems was found. Four regions of high molecular fluctuation were identified in all modified systems, being the subdomains IA and IIIB those with the most remarkable local conformational changes. Regarding essential modes of domain movements, it was evidenced that the most representatives were those related to IA subdomain, while IIIB subdomain presented discrete changes.

CONCLUSIONS

RCSs induces local structural changes on mono‑carbonylated BSA. Also, this study extends our knowledge on how carbonylation by RCSs induce structural effects on proteins.

Discovering potential WRN inhibitors from natural product database through computational methods

Microsatellite instability (MSI) is a relatively common feature associated with multiple cancers, and Werner syndrome (WRN) ATP-dependent helicase has been recognized as a novel target for treating MSI cancers, such as colorectal cancer. A small-molecule inhibitor targeting WRN would be a promising strategy for treating colorectal cancer with high MSI expression. In this study, we employed a computer-assisted drug discovery strategy to screen over 30,000 natural product molecules. By using a combination of docking, ligand efficiency, Molecular Mechanics/Generalized Born Surface Area (MM/GBSA), and thermodynamic integration (TI) calculations, we identified MOL008980, MOL010740, MOL011832, T4743, TN1166, and TNP-002173 as potential WRN inhibitors. Subsequent molecular dynamics simulation revealed that these screened natural products possessed better binding dynamic characteristics than ATP substrate and were capable of inhibiting the dynamic process of WRN, making them potential strong ATP competitive inhibitors. In conclusion, our computational approach revealed potential WRN inhibitors from a natural product database, providing a theoretical basis for future research.

Identifying and characterization of novel broad-spectrum bacteriocins from the Shanxi aged vinegar microbiome: Machine learning, molecular simulation, and activity validation

Shanxi aged vinegar microbiome encodes a wide variety of bacteriocins. The aim of this study was to mine, screen and characterize novel broad-spectrum bacteriocins from the large-scale microbiome data of Shanxi aged vinegar through machine learning, molecular simulation and activity validation. A total of 158 potential bacteriocins were innovatively mined from 117,552 representative genes based on metatranscriptomic information from the Shanxi aged vinegar microbiome using machine learning techniques and 12 microorganisms were identified to secrete bacteriocins at the genus level. Subsequently, employing AlphaFold2 structure prediction and molecular dynamics simulations, eight bacteriocins with high stability were further screened, and all of them were confirmed to have bacteriostatic activity by the Escherichia coli BL21 expression system. Then, gene_386319 (named LAB-3) and gene_403047 (named LAB-4) with the strongest antibacterial activities were purified by two-step methods and analyzed by mass spectrometry. The two bacteriocins have broad-spectrum antimicrobial activity with minimum inhibitory concentration values of 6.79 μg/mL-15.31 μg/mL against Staphylococcus aureus and Escherichia coli. Furthermore, molecular docking analysis indicated that LAB-3 and LAB-4 could interact with dihydrofolate reductase through hydrogen bonds, salt-bridge forces and hydrophobic forces. These findings suggested that the two bacteriocins could be considered as promising broad-spectrum antimicrobial agents.

Structural and functional analysis of Cyanovirin-N homologs: Carbohydrate binding affinities and antiviral potential of cyanobacterial peptides

Cyanobacteria, a group of photosynthetic prokaryotes, can sinthesize several substances due to their secondary metabolism, with notable properties, such as Cyanovirin-N(CVN), a carbohydrate-binding lectin, that exhibits antiviral activity against several pathogens, due to its ability to bind viral surface carbohydrates such as mannose, thus interfering with the viral entry on the cell. CVN has been described in several cyanobacterial strains and shows biotechnological potential for the development of drugs of pharmaceutical interest. This study focuses on the genomic exploration and characterization of Cyanovirin-N homologs to assess the conservation of carbohydrate-binding affinity within the group. The analysis of their antiviral properties was carried out using bioinformatics tools to study protein models through an in silico pipeline, following the steps of genomic prospection on public databases, homology modeling, docking, molecular dynamics and energetic analysis. Mannose served as the reference ligand, and the lectins' binding affinity with mannose was assessed across Cyanovirin-N homologs. Genomic mining identified 33 cyanobacterial lectin sequences, which underwent structural and functional characterization. The results obtained from this work indicate strong carbohydrate affinity on several homologs, pointing to the conservation of antiviral properties alongside the group. However, this affinity was not uniformly distributed among sequences, exhibiting significant heterogeneity in binding site residues, suggesting potential multi-ligand binding capabilities on the Cyanovirin-N homologs group. Studies focused on the properties involved in these molecules and the investigation of the genetic diversity of Cyanovirin-N homologs could provide valuable insights into the discovery of new drug candidates, harvesting the potential of bioinformatics for large-scale functional and structural analysis.

The role of the active site lysine residue on FAD reduction by NADPH in glutathione reductase

Glutathione reductase (GR) is a two dinucleotide binding domain flavoprotein (tDBDF) that catalyzes the reduction of glutathione disulfide to glutathione coupled to the oxidation of NADPH to NADP+. An interesting feature of GR and other tDBDFs is the presence of a lysine residue (Lys-66 in human GR) at the active site, which interacts with the flavin group, but has an unknown function. To better understand the role of this residue, the dynamics of GR was studied using molecular dynamics simulations, and the reaction mechanism of FAD reduction by NADPH was studied using QM/MM molecular modeling. The two possible protonation states of Lys-66 were considered: neutral and protonated. Molecular dynamics results suggest that the active site is more structured for neutral Lys-66 than for protonated Lys-66. QM/MM modeling results suggest that Lys-66 should be in its neutral state for a thermodynamically favorable reduction of FAD by NADPH. Since the reaction is unfavorable with protonated Lys-66, the reverse reaction (the reduction of NADP+ by FADH-) is expected to take place. A phylogenetic analysis of various tDBDFs was performed, finding that an active site lysine is present in different the tDBDFs enzymes, suggesting that it has a conserved biological role. Overall, these results suggest that the protonation state of the active site lysine determines the energetics of the reaction, controlling its reversibility.

Nature-inspired substituted 3-(imidazol-2-yl) morpholines targeting human topoisomerase IIα: Dynophore-derived discovery

The molecular nanomachine, human DNA topoisomerase IIα, plays a crucial role in replication, transcription, and recombination by catalyzing topological changes in the DNA, rendering it an optimal target for cancer chemotherapy. Current clinical topoisomerase II poisons often cause secondary tumors as side effects due to the accumulation of double-strand breaks in the DNA, spurring the development of catalytic inhibitors. Here, we used a dynamic pharmacophore approach to develop catalytic inhibitors targeting the ATP binding site of human DNA topoisomerase IIα. Our screening of a library of nature-inspired compounds led to the discovery of a class of 3-(imidazol-2-yl) morpholines as potent catalytic inhibitors that bind to the ATPase domain. Further experimental and computational studies identified hit compound 17, which exhibited selectivity against the human DNA topoisomerase IIα versus human protein kinases, cytotoxicity against several human cancer cells, and did not induce DNA double-strand breaks, making it distinct from clinical topoisomerase II poisons. This study integrates an innovative natural product-inspired chemistry and successful implementation of a molecular design strategy that incorporates a dynamic component of ligand-target molecular recognition, with comprehensive experimental characterization leading to hit compounds with potential impact on the development of more efficient chemotherapies.

Advanced computational approaches to understand protein aggregation

Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.

Evaluation of surfactant-aided polycyclic aromatic hydrocarbon biodegradation by molecular docking and molecular dynamic simulation in the marine environment

Marine oil spills directly cause polycyclic aromatic hydrocarbons (PAHs) pollution and affect marine organisms due to their toxic property. Chemical and bio-based dispersants composed of surfactants and solvents are considered effective oil spill-treating agents. Dispersants enhance oil biodegradation in the marine environment by rapidly increasing their solubility in the water column. However, the effect of dispersants, especially surfactants, on PAHs degradation by enzymes produced by microorganisms has not been studied at the molecular level. The role of the cytochrome P450 (CYP) enzyme in converting contaminants into reactive metabolites during the biodegradation process has been evidenced, but the activity in the presence of surfactants is still ambiguous. Thus, this study focused on the evaluation of the impact of chemical and bio-surfactants (i.e., Tween 80 (TWE) and Surfactin (SUC)) on the biodegradation of naphthalene (NAP), chrysene (CHR), and pyrene (PYR), the representative components of PAHs, with CYP enzyme from microalgae Parachlorella kessleri using molecular docking and molecular dynamics (MD) simulation. The molecular docking analysis revealed that PAHs bound to residues at the CYP active site through hydrophobic interactions for biodegradation. The MD simulation showed that the surfactant addition changed the enzyme conformation in the CYP-PAH complexes to provide more interactions between the enzyme and PAHs. This led to an increase in the enzyme's capability to degrade PAHs. Binding free energy (ΔG‌‌Bind) calculations confirmed that surfactant treatment could enhance PAHs degradation by the enzyme. The SUC gave a better result on NAP and PYR biodegradation based on ΔG‌‌Bind, while TWE facilitated the biodegradation of CHR. The research outputs could greatly facilitate evaluating the behaviors of oil spill-treating agents and oil spill response operations in the marine environment.

Binding selectivity analysis of new delhi metallo-beta-lactamase-1 inhibitors using molecular dynamics simulations: Exploring possibilities for decoding antimicrobial drug resistance

BACKGROUND

New Delhi metallo-beta-lactamase-1 (NDM1) confers resistance to several bacterial species against a broad range of beta-lactam antibiotics and turning them into superbugs that pose a significant threat to healthcare systems worldwide. As such, it is a potentially relevant biological target for counteracting bacterial infections. Given the lack of effective treatment options against NDM1 producing bacteria, finding a reliable inhibitor for the NDM1 enzyme is crucial.

METHODS

Using molecular dynamics simulations, the binding selectivities and affinities of three ligands, viz. PNK, 3S0, and N1G were investigated against NDM1.

RESULTS

The results indicate that N1G binds with more affinity to NDM1 than PNK and 3S0. The binding energy decomposition analysis revealed that residues I35, W93, H189, K211, and N220 showed significant binding energies with PNK, 3S0, and N1G, and hence are crucially involved in the binding of the ligands to NDM1. Molecular dynamics trajectory analysis further elicited that the ligands influence dynamic flexibility of NDM1 morphology, which contributes to the partial selectivities of PNK, 3S0, and N1G.

CONCLUSIONS

This in silico study offers a vital information for developing potential NDM1 inhibitors with high selectivity. Nevertheless, in vitro and in vivo experimental validation is mandated to extend the possible applications of these ligands as NDM1 inhibitors that succor in combating antimicrobial resistance.

The reactivation kinetic analysis, molecular docking, and dynamics of oximes against three V-type nerve agents inhibited four human cholinesterases

Nerve agents pose significant threats to civilian and military populations. The reactivation of acetylcholinesterase (AChE) is critical in treating acute poisoning, but there is still lacking broad-spectrum reactivators, which presents a big challenge. Therefore, insights gained from the reactivation kinetic analysis and molecular docking are essential for understanding the behavior of reactivators towards intoxicated AChE. In this research, we present a systematic determination of the reactivation kinetics of three V agents-inhibited four human ChEs [(AChE and butyrylcholinesterase (BChE)) from either native or recombinant resources, namely, red blood cell (RBC) AChE, rhAChE, hBChE, rhBChE) reactivated by five standard oximes. We unveiled the effect of native and recombinant ChEs on the reactivation kinetics of V agents ex vitro, where the reactivation kinetics characteristic of Vs-inhibited BChE was reported for the first time. In terms of the inhibition type, all of the five oxime reactivators exhibited noncompetitive inhibition. The inhibition potency of these reactivators would not lead to the difference in the reactivation kinetics between native and recombinant ChE. Despite the significant differences between the native and recombinant ChEs observed in the inhibition, aging, and spontaneous reactivation kinetics, the reactivation kinetics of V agent-inhibited ChEs by oximes were less differentiated, which were supported by the ligand docking results. We also found differences in the reactivation efficiency between five reactivators and the phosphorylated enzyme, and molecular dynamic simulations can further explain from the perspectives of conformational stability, hydrogen bonding, binding free energies, and amino acid contributions. By Poisson-Boltzmann surface area (MM-PBSA) calculations, the total binding free energy trends aligned well with the experimental kr2 values.

Engineering and finetuning expression of SerC for balanced metabolic flux in vitamin B6 production

Vitamin B6 plays a crucial role in cellular metabolism and stress response, making it an essential component for growth in all known organisms. However, achieving efficient biosynthesis of vitamin B6 faces the challenge of maintaining a balanced distribution of metabolic flux between growth and production. In this study, our focus is on addressing this challenge through the engineering of phosphoserine aminotransferase (SerC) to resolve its redundancy and promiscuity. The enzyme SerC was semi-designed and screened based on sequences and predicted kcat values, respectively. Mutants and heterologous proteins showing potential were then fine-tuned to optimize the production of vitamin B6. The resulting strain enhances the production of vitamin B6, indicating that different fluxes are distributed to the biosynthesis pathway of serine and vitamin B6. This study presents a promising strategy to address the challenge posed by multifunctional enzymes, with significant implications for enhancing biochemical production through engineering processes.

Helical reorganization in the context of membrane protein folding: Insights from simulations with bacteriorhodopsin (BR) fragments

Membrane protein folding is distinct from folding of soluble proteins. Conformational acquisition in major membrane protein subclasses can be delineated into insertion and folding processes. An exception to the "two stage" folding, later developed to "three stage" folding, is observed within the last two helices in bacteriorhodopsin (BR), a system that serves as a model membrane protein. We employ a reductionist approach to understand interplay of molecular factors underlying the apparent defiance. Leveraging available solution NMR structures, we construct, sample in silico, and analyze partially (PIn) and fully inserted (FIn) BR membrane states. The membrane lateral C-terminal helix (CH) in PIn is markedly prone to transient structural distortions over microsecond timescales; a disorder prone region (DPR) is thereby identified. While clear transmembrane propensities are not acquired, the distortions induce alterations in local membrane curvature and area per lipid. Importantly, energetic decompositions reveal that overall, the N-terminal helix (NH) is thermodynamically more stable in the PIn. Higher overall stability of the FIn arises from favorable interactions between the NH and the CH. Our results establish lack of spontaneous transition of the PIn to the FIn, and attributes their partitioning to barriers that exceed those accessible with thermal fluctuations. This work paves the way for further detailed studies aimed at determining the thermo-kinetic roles of the initial five helices, or complementary external factors, in complete helical folding and insertion in BR. We comment that complementing such efforts with the growing field of machine learning assisted energy landscape searches may offer unprecedented insights.

On the interactions of peptides with gold nanoparticles: effects of sequence and size

Peptide-based self-assembly and synthesis techniques have emerged as a viable approach to designing active and stable inorganic nanostructures in aqueous media. In the present study, we use all-atom molecular dynamic (MD) simulations to study the interactions of ten short peptides (namely A3, AgBP1, AgBP2, AuBP1, AuBP2, GBP1, Midas2, Pd4, Z1, and Z2) with different gold nanoparticles (of different diameters ranging from 2 to 8 nm). Our MD simulation results imply that the gold nanoparticles have a remarkable effect on the stability and conformational properties of peptides. Moreover, the size of the gold nanoparticles and the type of peptide amino acid sequences play important roles in the stability of the peptide-AuNP complexes. Our results reveal that some amino acids such as Tyr, Phe, Met, Lys, Arg, and Gln have direct contact with the metal surface in comparison with Gly, Ala, Pro, Thr, and Val residues. The peptide adsorption on the surface of the gold nanoparticles is favorable from the energetic viewpoint, in which the van der Waals (vdW) interactions between the peptides and metal surface can be considered as one of the driving forces for the complexation process. The calculated Gibbs binding energies indicate that AuNPs have more sensitivity against the GBP1 peptide in the presence of different peptides. Overall, the results of this study can provide new insight into the peptide interaction with the gold nanoparticles from the molecular viewpoint, which can be important for designing new biomaterials based on the peptides and gold nanoparticles.Communicated by Ramaswamy H. Sarma.

Phosphate uptake in PhoX: Molecular mechanisms

PhoX is a high-affinity phosphate binding protein, present in Xanthomonas citri, a phytopathogen responsible for the citrus canker disease. Performing molecular dynamics simulations and different types of computational analyses, we study the molecular mechanisms at play in relation to phosphate binding, revealing the global functioning of the protein: PhoX naturally oscillates along its global normal modes, which allow it to explore both bound and unbound conformations, eventually attracting a nearby negative phosphate ion to the highly positive electrostatic potential on its surface, particularly close to the binding pocket. There, several hydrogen bonds are formed with the two main domains of the structure. Phosphate creates, in this way, a strong bridge that connects the domains, keeping itself between them, in a tight closed conformation, explaining its high binding affinity.

Molecular choreography: Unveiling the dynamic landscape of type IIA DNA topoisomerases before T-segment passage through all-atom simulations

Type IIA DNA topoisomerases are molecular nanomachines responsible for controlling topological states of DNA molecules. Here, we explore the dynamic landscape of yeast topoisomerase IIA during key stages of its catalytic cycle, focusing in particular on the events preceding the passage of the T-segment. To this end, we generated six configurations of fully catalytic yeast topo IIA, strategically inserted a T-segment into the N-gate in relevant configurations, and performed all-atom simulations. The essential motion of topo IIA protein dimer was characterized by rotational gyrating-like movement together with sliding motion within the DNA-gate. Both appear to be inherent properties of the enzyme and an inbuilt feature that allows passage of the T-segment through the cleaved G-segment. Coupled dynamics of the N-gate and DNA-gate residues may be particularly important for controlled and smooth passage of the T-segment and consequently the prevention of DNA double-strand breaks. QTK loop residue Lys367, which interacts with ATP and ADP molecules, is involved in regulating the size and stability of the N-gate. The unveiled features of the simulated configurations provide insights into the catalytic cycle of type IIA topoisomerases and elucidate the molecular choreography governing their ability to modulate the topological states of DNA topology.

Network pharmacology-based screening of active constituents of Avicennia marina and their clinical biochemistry related mechanism against breast cancer

Breast cancer is the second major cause of cancer death in women globally. Avicennia marina is a medicinal plant that belongs to the family Acanthaceae and is known as grey or white mangrove. It has antioxidant, antiviral, anticancer, anti-inflammatory, and antibacterial activity in the treatment of various diseases including cancer. The goal of the study is to use a network pharmacology method to identify the potential phenomena of bioactive compounds of A. marina in the treatment of breast cancer and explore clinical biochemistry related aspects. A total of 74 active compounds of A. marina were retrieved from various databases as well as a literature review and collectively 429 targets of these compounds were identified by STITCH and Swiss Target Prediction databases. Breast cancer related 15606 potential targets were retrieved from the GeneCards database. A Venn diagram was drawn to find common key targets. To check the biological functions, the GO enrichment and KEGG pathways analysis of 171 key targets were performed through the DAVID database. To understand the interactions among key targets, Protein-protein interaction (PPI) studies were completed using the STRING database, and the Protein-Protein Interaction (PPI) network, as well as the compound-target-pathway network, was constructed using Cytoscape 3.9.0. Finally, molecular docking analysis of 5 hub genes named tumor protein 53 (TP53), catenin beta 1 (CTNNB1), interleukin 6 (IL6), tumor necrosis factor (TNF), and RAC-alpha serine/threonine protein kinases 1 (AKT1) with the active constituent of A. marina against breast cancer were performed. Additionally, a molecular docking study demonstrates that active drugs have a higher affinity for the target that may be used to decrease breast cancer. The molecular dynamic simulation analysis predicted the very stable behavior of docked complexes with no global structure deviations seen. The MMGBSA further supported strong intermolecular interactions with net energy values as; AKT1_Betulinic_acid (-20.97 kcal/mol), AKT1_Stigmasterol (-44.56 kcal/mol), TNF_Betulinic_acid (-28.68 kcal/mol) and TNF_Stigmastero (-29.47 kcal/mol).Communicated by Ramaswamy H. Sarma.

Decoding Inhibitor Egression from Wild-Type and G2019S Mutant LRRK2 Kinase: Insights into Unbinding Mechanisms for Precision Drug Design in Parkinson's Disease

Leucine-rich repeat kinase 2 (LRRK2) remains a viable target for drug development since the discovery of the association of its mutations with Parkinson's disease (PD). G2019S (in the kinase domain) is the most common mutation for LRRK2-based PD. Though various types of inhibitors have been developed for the kinase domain to reduce the effect of the mutation, understanding the working of these inhibitors at the molecular level is still ongoing. This study focused on the exploration of the dissociation mechanism (pathways) of inhibitors from (WT and G2019S) LRRK2 kinase (using homology model CHK1 kinase), which is one of the crucial aspects in drug discovery. Here, two ATP-competitive type I inhibitors, PF-06447475 and MLi-2 (Comp1 and Comp2 ), and one non-ATP-competitive type II inhibitor, rebastinib (Comp3), were considered for this investigation. To study the unbinding process, random accelerated molecular dynamics simulations were performed. The binding free energies of the three inhibitors for different egression paths were determined using umbrella sampling. This work found four major egression pathways that were adopted by the inhibitors Comp1 (path1, path2, and path3), Comp2 (path1, path2 and path3), and Comp3 (path3 and path4). Also, the mechanism of unbinding for each path and key residues involved in unbinding were explored. Mutation was not observed to impact the preference of the particular egression pathways for both LRRK2-Comp1 and -Comp2 systems. However, the findings suggested that the size of the inhibitor molecules might have an effect on the preference of the egression pathways. The binding energy and residence time of the inhibitors followed a similar trend to experimental observations. The findings of this work might provide insight into designing more potent inhibitors for the G2019S LRRK2 kinase.

Identification of FTY720 and COH29 as novel topoisomerase I catalytic inhibitors by experimental and computational studies

The development of novel topoisomerase I (TOP1) inhibitors is crucial for overcoming the drawbacks and limitations of current TOP1 poisons. Here, we identified two potential TOP1 inhibitors, namely, FTY720 (a sphingosine 1-phosphate antagonist) and COH29 (a ribonucleotide reductase inhibitor), through experimental screening of known active compounds. Biological experiments verified that FTY720 and COH29 were nonintercalative TOP1 catalytic inhibitors that did not induce the formation of DNA-TOP1 covalent complexes. Molecular docking revealed that FTY720 and COH29 interacted favorably with TOP1. Molecular dynamics simulations revealed that FTY720 and COH29 could affect the catalytic domain of TOP1, thus resulting in altered DNA-binding cavity size. The alanine scanning and interaction entropy identified Arg536 as a hotspot residue. In addition, the bioinformatics analysis predicted that FTY720 and COH29 could be effective in treating malignant breast tumors. Biological experiments verified their antitumor activities using MCF-7 breast cancer cells. Their combinatory effects with TOP1 poisons were also investigated. Further, FTY720 and COH29 were found to cause less DNA damage compared with TOP1 poisons. The findings provide reliable lead compounds for the development of novel TOP1 catalytic inhibitors and offer new insights into the potential clinical applications of FTY720 and COH29 in targeting TOP1.

TMV-CP based rational design and discovery of α-Amide phosphate derivatives as anti plant viral agents

The tobacco mosaic virus coat protein (TMV-CP) is indispensable for the virus's replication, movement and transmission, as well as for the host plant's immune system to recognize it. It constitutes the outermost layer of the virus particle, and serves as an essential component of the virus structure. TMV-CP is essential for initiating and extending viral assembly, playing a crucial role in the self-assembly process of Tobacco Mosaic Virus (TMV). This research employed TMV-CP as a primary target for virtual screening, from which a library of 43,417 compounds was sourced and SH-05 was chosen as the lead compound. Consequently, a series of α-amide phosphate derivatives were designed and synthesized, exhibiting remarkable anti-TMV efficacy. The synthesized compounds were found to be beneficial in treating TMV, with compound 3g displaying a slightly better curative effect than Ningnanmycin (NNM) (EC50 = 304.54 µg/mL) at an EC50 of 291.9 µg/mL. Additionally, 3g exhibited comparable inactivation activity (EC50 = 63.2 µg/mL) to NNM (EC50 = 67.5 µg/mL) and similar protective activity (EC50 = 228.9 µg/mL) to NNM (EC50 = 219.7 µg/mL). Microscale thermal analysis revealed that the binding of 3g (Kd = 4.5 ± 1.9 µM) to TMV-CP showed the same level with NNM (Kd = 5.5 ± 2.6 µM). Results from transmission electron microscopy indicated that 3g could disrupt the structure of TMV virus particles. The toxicity prediction indicated that 3g was low toxicity. Molecular docking showed that 3g interacted with TMV-CP through hydrogen bond, attractive charge interaction and π-Cation interaction. This research provided a novel α-amide phosphate structure target TMV-CP, which may help the discovery of new anti-TMV agents in the future.

A conserved polar residue plays a critical role in mismatch detection in A-family DNA polymerases

In A-family DNA polymerases (dPols), a functional 3'-5' exonuclease activity is known to proofread newly synthesized DNA. The identification of a mismatch in substrate DNA leads to transfer of the primer strand from the polymerase active site to the exonuclease active site. To shed more light regarding the mechanism responsible for the detection of mismatches, we have utilized DNA polymerase 1 from Aquifex pyrophilus (ApPol1). The enzyme synthesized DNA with high fidelity and exhibited maximal exonuclease activity with DNA substrates bearing mismatches at the -2 and - 3 positions. The crystal structure of apo-ApPol1 was utilized to generate a computational model of the functional ternary complex of this enzyme. The analysis of the model showed that N332 forms interactions with minor groove atoms of the base pairs at the -2 and - 3 positions. The majority of known A-family dPols show the presence of Asn at a position equivalent to N332. The N332L mutation led to a decrease in the exonuclease activity for representative purine-pyrimidine, and pyrimidine-pyrimidine mismatches at -2 and - 3 positions, respectively. Overall, our findings suggest that conserved polar residues located towards the minor groove may facilitate the detection of position-specific mismatches to enhance the fidelity of DNA synthesis.