Molecular recognition between volatile molecules and odorant binding proteins 7 by homology modeling, molecular docking and molecular dynamics simulation.

BACKGROUND: Odorant-binding proteins (OBPs) in insects are key to detection and recognition of external chemical signals associated with survival. OBP7 in Spodoptera frugiperda's larval stage (SfruOBP7) may search for host plants by sensing plant volatiles, which are important sources of pest attractants and repellents. However, the atomic-level basis of binding modes remains elusive. RESULTS: SfruOBP7 structure was constructed through homology modeling, and complex models of six plant volatiles ((E)-2-hexenol, α-pinene, (Z)-3-hexenyl acetate, lauric acid, O-cymene and 1-octanol) and SfruOBP7 were obtained through molecular docking. To study the detailed interactions between the six plant volatile molecules and SfruOBP7, we conducted three 300 ns molecular dynamics simulations for each study object. The correlation coefficients between binding free energy obtained by molecular mechanics/generalized Born surface area together with solvated interaction energy methods and experimental values are 0.90 and 0.88, respectively, showing a good correlation. By comparing binding free energy along with interaction patterns between SfruOBP7 and the six volatile molecules, hotspot residues of SfruOBP7 when binding with different volatile molecules were determined. Hydrophobic interactions stemming from van der Waals interactions play a significant role in SfruOBP7 and these plant volatile systems. CONCLUSION: The optimized three-dimensional structure of SfruOBP7 and its binding modes with six plant volatiles revealed their interactions, thus providing a means for estimating the binding energies of other plant volatiles. Our study will help to guide the rational design of effective and selective insect attractants. © 2024 Society of Chemical Industry.

Molecular basis of RNA recognition by TBP of HIV-1 from multiple molecular dynamics simulations and energy predictions.

The frequent outbreaks of the AIDS (Acquired Immune Deficiency Syndrome) pandemic and the limited availability of anti-Human Immunodeficiency Virus (HIV) drugs highlight the urgent need to develop new antiviral drugs. A detailed understanding of the interactions between TAR-Binding Proteins (TBP) and RNA will facilitate the discovery of new anti-AIDS drugs. In order to characterize and explore the key interactions between RNA and TBP, we focused on the wild type (WT) and three mutant TBPs (TBP6.9, TBP6.7, and TBP6.3) with RNA, multiple molecular dynamics simulation and energy computation were performed. The results showed that 12 key residues played a major role in the interaction between TBP and RNA. The mutated residues of TBP changed the interaction between their surrounding residues and RNA, thus affecting the binding of TBP to RNA. In addition, structural and energy analyses showed that in contrast with WT TBP-RNA complex, the mutated residues had little effect on the backbone structure of TBP, but changes in the van der Waals interactions and electrostatic interaction associated with the side chains are responsible for the altered the binding between three mutant TBPs and RNA complexes. The discovery of TBP-RNA recognition mechanism in our work provides some useful insights and new opportunities for the development of anti-aids drugs.

Electrochemical Enabled Cascade Phosphorylation of N-H/O-H/S-H Bonds with P-H Compounds: An Efficient Access to P(O)-X Bonds.

An electrochemical three component cascade phosphorylation reaction of various heteroatoms-containing nucleophiles including carbazoles, indoles, phenols, alcohols, and thiols with Ph2 PH has been established. Electricity is used as the "traceless" oxidant and water and air are utilized as the "green" oxygen source. All kinds of structurally diverse organophosphorus compounds with P(O)-N/P(O)-O/P(O)-S bonds are assembled in moderate to excellent yields (three categories of phosphorylation products, 50 examples, up to 97 % yield). A tentative free radical course is put forward to rationalize the reaction procedure.

Multiple Molecular Dynamics Simulations and Free-Energy Predictions Uncover the Susceptibility of Variants of HIV-1 Protease against Inhibitors Darunavir and KNI-1657.

HIV-1 protease (PR) is considered to be the main targets of anti-AIDS drug design because of its role in the proteolytic processing of viral polyproteins. However, the emergence of drug-resistant HIV has become a major problem in the therapy of HIV-1-infected patients. Focused on the complexes of wild type (WT) PR and two mutant PRs (V32I/L33F/I54M/V82I and V32I/L33F/I54M/I84 V) with inhibitors Darunavir (DRV) and KNI-1657 (KNI), respectively, we have conducted research on the conformational dynamics and the resistance mechanism caused by residue mutations through multiple molecular dynamics (MD) simulations combined with an energy (MM-PBSA and solvated interaction energy (SIE)) prediction. The results indicate that mutated residues of PR alter the distance between flap regions and catalytic sites, the volume of the inner catalytic site, and the curling degree of the flap tips, thereby affecting DRV and KNI inhibitor binding to PR. These mutated residues reduced the binding affinity of the two mutant PRs to DRV, resulting in drug resistance, whereas the two mutant PRs increase the binding affinity with KNI, indicating they enhance the sensitivity to KNI. Compared with the WT PR, the changes in van der Waals interaction and electrostatic interaction in the two variant PRs play a vital part in the binding of PR with DRV and KNI. These results may supply valuable guidance for the design of anti-AIDS drugs targeting PR.

Multiple Molecular Dynamics Simulations of the Inhibitor GRL-02031 Complex with Wild Type and Mutant HIV-1 Protease Reveal the Binding and Drug-Resistance Mechanism.

Human immunodeficiency virus type 1 (HIV-1) protease is regarded as a fascinating target for drug development against HIV infection. However, mutations causing drug resistance severely limit the efficiency of the recently marketed drugs in the treatment of HIV replication. To elucidate the binding mechanism of HIV-1 protease with promising inhibitor GRL-02031 and further to probe the resistance mechanism associated with mutations (I47V, L76V, V82A, and N88D) to the inhibitor, we applied multiple molecular dynamics (MMD) simulations along with energy analysis by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and solvated interaction energy (SIE) methodology on specific HIV-1 protease with GRL-0231 complexes. On the basis of detail analysis of the simulations, we revealed key characteristics that constitute the drug resistance of four mutation HIV-1 proteases toward GRL-02031: substitution of the side chain in these four mutation residues leads to a change in the distances between the flaps and catalytic sites, thereby reducing the affinity for GRL-02031 with these four mutation proteases, even though the L76V and N88D residues cannot directly contact GRL-02031. The results of energy analysis according to the MM-PBSA and SIE methods further indicated that hydrophobic interaction was considered to be the prime driving force for inhibitor GRL-02031 binding to protease and the decrease in van der Waals interactions between inhibitor GRL-02031 and mutant proteases as the primary cause of the drug resistance. Analyses of the hydrogen bonds and atomic interactions further provided detailed explanations for the resistance of these four mutation proteases toward inhibitor GRL-02031. The present study provides potential guidance on the structure-based inhibitors' design targeting HIV-1 protease.

Revealing the binding and drug resistance mechanism of amprenavir, indinavir, ritonavir, and nelfinavir complexed with HIV-1 protease due to double mutations G48T/L89M by molecular dynamics simulations and free energy analyses.

Infection by human immunodeficiency virus type 1 (HIV-1) not only destroys the immune system bringing about acquired immune deficiency syndrome (AIDS), but also induces serious neurological diseases including behavioral abnormalities, motor dysfunction, toxoplasmosis, and HIV-1 associated dementia. The emergence of HIV-1 multidrug-resistant mutants has become a major problem in the therapy of patients with HIV-1 infection. Focusing on the wild type (WT) and G48T/L89M mutated forms of HIV-1 protease (HIV-1 PR) in complex with amprenavir (APV), indinavir (IDV), ritonavir (RTV), and nelfinavir (NFV), we have investigated the conformational dynamics and the resistance mechanism due to the G48T/L89M mutations by conducting a series of molecular dynamics (MD) simulations and free energy (MM-PBSA and solvated interaction energy (SIE)) analyses. The simulation results indicate that alterations in the side-chains of G48T/L89M mutated residues cause the inner active site to increase in volume and induce more curling of the flap tips, which provide the main contributions to weaker binding of inhibitors to the HIV-1 PR. The results of energy analysis reveal that the decrease in van der Waals interactions of inhibitors with the mutated PR relative to the wild-type (WT) PR mostly drives the drug resistance of mutations toward these four inhibitors. The energy decomposition analysis further indicates that the drug resistance of mutations can be mainly attributed to the change in van der Waals and electrostatic energy of some key residues (around Ala28/Ala28' and Ile50/Ile50'). Our work can give significant guidance to design a new generation of anti-AIDS inhibitors targeting PR in the therapy of patients with HIV-1 infection.

Molecular docking and molecular dynamics simulation decoding molecular mechanism of EDCs binding to hERRγ.

CONTEXT: Human estrogen-related receptor γ (hERRγ) is a key protein involved in various endocrines and metabolic signaling. Numerous environmental endocrine-disrupting chemicals (EDCs) can impact related physiological activities through receptor signaling pathways. Focused on hERRγ with 4-isopropylphenol, bisphenol-F (BPF), and BP(2,2)(Un) complexes, we executed molecular docking and multiple molecular dynamics (MD) simulations along with molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) and solvation interaction energy (SIE) calculation to study the detailed dynamical structural characteristics and interactions between them. Molecular docking showed that hydrogen bonds and hydrophobic interactions were the prime interactions to keep the stability of BPF-hERRγ and hERRγ-BP(2,2)(Un) complexes. Through MD simulations, we observed that all complexes reach equilibrium during the initial 50 ns of simulation, but these three EDCs lead to local structure changes in hERRγ. Energy results further identified key residues L268, V313, L345, and F435 around the binding pockets through CH-π, π-π, and hydrogen bonds interactions play an important stabilizing role in the recognition with EDCs. And most noticeable of all, hydrophobic methoxide groups in BP(2,2)(Un) is useful for decreasing the binding ability between EDCs and hERRγ. These results may contribute to evaluate latent diseases associated with EDCs exposure at the micro level and find potential substitutes. METHOD: Autodock4.2 was used to conduct the molecular docking, sietraj program was performed to calculate the energy, and VMD software was used to visualize the structure. Amber18 was conducted to perform the MD simulation and other analyses.

2D MoB MBene: An Efficient Co-Catalyst for Photocatalytic Hydrogen Production under Visible Light.

Highly active and low-cost co-catalysts have a positive effect on the enhancement of solar H2 production. Here, we employ two-dimensional (2D) MBene as a noble-metal-free co-catalyst to boost semiconductor for photocatalytic H2 production. MoB MBene is a 2D nanoboride, which is directly made from MoAlB by a facile hydrothermal etching and manual scraping off process. The as-synthesized MoB MBene with purity >95 wt % is treated by ultrasonic cell pulverization to obtain ultrathin 2D MoB MBene nanosheets (∼0.61 nm) and integrated with CdS via an electrostatic interaction strategy. The CdS/MoB composites exhibit an ultrahigh photocatalytic H2 production activity of 16,892 µmol g-1 h-1 under visible light, surpassing that of pure CdS by an exciting factor of ≈1135%. Theoretical calculations and various measurements account for the high performance in terms of Gibbs free energy, work functions, and photoelectrochemical properties. This work discovers the huge potential of these promising 2D MBene family materials as high-efficiency and low-cost co-catalysts for photocatalytic H2 production.

Molecular basis of ssDNA recognition by RBM45 protein of neurodegenerative disease from multiple molecular dynamics simulations and energy predictions.

Neurodegenerative diseases (NDD) are a group of cognitive and behavioral disorders characterized by progressive loss of neuronal structure and function. As the population ages, the incidence is getting higher and higher, but there is currently no effective treatment. The details of RNA/DNA recognition by the RNA-binding protein RBM45 closely related to neurodegenerative diseases through its two tandem RNA-recognition domains at its N-terminus have important implications for structure-based drug discovery against degenerative diseases. To explore the key characteristics of interaction between ssDNA and RBM45, we performed multiple molecular dynamics (MD) simulations along with MM-PBSA energy prediction on the complexes of wild type (WT) and three mutant RBM45s (K100A, F124A/Y165A, and F29A/F70A/F124A/Y165A) with ssDNA, respectively. The findings suggest that these mutated residues of RBM45 modify the interaction of their surrounding residues with ssDNA, thereby affecting RBM45 protein binding to ssDNA. In contrast with WT RBM45 protein, variations in van der Waals and electrostatic interactions with ssDNA caused by these three RBM45 mutants are critical to affect binding between them. In addition, energy analysis showed that RBM45 is a specific ssDNA-binding protein. The results of our work provide valuable theoretical guidance for the design effective drugs of NDD.

Structural and energetic basis of interaction between human estrogen-related receptor γ and environmental endocrine disruptors from multiple molecular dynamics simulations and free energy predictions.

Environmental endocrine disruptors (EEDs), a class of molecules that are widespread in our environment, may adversely affect the endocrine system. Exploring the interactions between these compounds and their potential targets is important for assessing their role in the organism. Focused on the human estrogen-related receptor γ (hERRγ) with BPA, BPB, HPTE, BPE, BP(2,2)(Et), and BP(2,2)(MeO) complexes, respectively, we groped for the mechanisms of conformational changes and interactions of hERRγ when binding to these six EEDs by combining multiple molecular dynamics (MD) simulations with energy prediction (MM-PBSA and solvated interaction energy (SIE)). Dynamics analysis results revealed these six EEDs have different effects on the internal dynamics of hERRγ, resulting in significant changes in the interaction of key residues around Leu268, Val313, Leu345, and Phe435 with EEDs, and thus affected its binding energy with these EEDs. The energy calculations further demonstrated that van der Waals interactions are critical for these EEDs binding to hERRγ. These results present detailed molecular insight into the interaction features between EEDs and hERRγ and help guide the search for safer alternatives to BPA.

Multiple Molecular Dynamics Simulations and Energy Analysis Unravel the Dynamic Properties and Binding Mechanism of Mutants HIV-1 Protease with DRV and CA-p2.

PRS17, a variant of human immunodeficiency virus type I protease (HIV-1 PR), has 17 mutated residues showing high levels of multidrug resistance. To describe the effects of these mutated residues on the dynamic properties and the binding mechanism of PR with substrate and inhibitor, focused on six systems (two complexes of WT PR and PRS17 with inhibitor Darunavir (DRV), two complexes of WT PR and PRS17 with substrate analogue CA-p2, two unligand WT PR and PRS17), we performed multiple molecular dynamics (MD) simulations combined with MM-PBSA and solvated interaction energy (SIE) methods. For both the unligand PRs and ligand-PR complexes, the results from simulations revealed 17 mutated residues alter the flap-flap distance, the distance from flap regions to catalytic sites, and the curling degree of the flap tips. These mutated residues changed the flexibility of the flap region in PR, and thus affected its binding energy with DRV and CA-p2, resulting in differences in sensitivity. Hydrophobic cavity makes an important contribution to the binding of PR and ligands. And most noticeable of all, the binding of the guanidine group in CA-p2 and Arg8' of PRS17 is useful for increasing their binding ability. These results have important guidance for the further design of drugs against multidrug resistant PR. IMPORTANCE Developing effective anti-HIV inhibitors is the current requirement to cope with the emergence of the resistance of mutants. Compared with the experiments, MD simulations along with energy calculations help reduce the time and cost of designing new inhibitors. Based on our simulation results, we propose two factors that may help design effective inhibitors against HIV-1 PR: (i) importance of hydrophobic cavity, and (ii) introduction of polar groups similar to the guanidine group.

Electrochemical Oxidative Dehydrogenative Phosphorylation of N-Heterocycles with P(O)-H Compounds in Imidazolium-Based Ionic Liquid.

We report a direct and green electrochemical oxidative cross-dehydrogenative coupling reaction of N-heterocycles with hydrogen phosphoryl compounds under external oxidant-free conditions. Various phosphorylation products of substituted carbazoles and indoles are assembled in modest to excellent yields. A hydrogen release process is preliminarily demonstrated and H2 is the sole byproduct. An imidazolium based ionic liquid is selected as the optimal electrolyte.

CCL18 promotes the invasion and metastasis of breast cancer through Annexin A2.

Chemokine (C­C motif) ligand 18 (CCL18) is derived from breast tumor­associated macrophages (TAMs), which are primarily a macrophage subpopulation with an M2 phenotype. CCL18 binds to its receptor, PYK2 N­terminal domain interacting receptor 1 (Nir1), and promotes tumor progression and metastasis by inducing epithelial­mesenchymal transition (EMT) via the PI3K/Akt/GSK3ß/Snail signaling pathway in breast cancer cells. Recent research shows that Annexin A2 (AnxA2) plays a significant role in the invasion, metastasis, angiogenesis, proliferation, F­actin polymerization and multidrug resistance to chemotherapy of breast cancer. The present study aimed to elucidate the molecular mechanisms by which CCL18 promotes breast cancer progression through AnxA2 which are not fully understood. Western blot analysis showed that the expression of AnxA2 was upregulated in highly invasive breast cancer cell lines and invasive ductal carcinoma. Furthermore, through chemotaxis, scratch, Matrigel invasion, and spontaneous metastasis assays, it was demonstrated that AnxA2 enhanced the invasion of breast cancer cells and the metastasis of human breast cancer cells to lungs of SCID mice with CCL18 stimulation. Cellular F­actin measurement assay showed that reduction of AnxA2 suppressed CCL18­induced F­actin polymerization though phosphorylation of integrin ß1 in breast cancer cells. Immunofluorescence and western blot analysis revealed that AnxA2 promoted CCL18­induced EMT via the PI3K/Akt/GSK3ß/Snail signaling pathway, and LY294002 inhibited the phosphorylation of AnxA2 in vitro. In brief, AnxA2, as a downstream molecule of Nir 1 binding to CCL18, promotes invasion and metastasis by EMT through the PI3K/Akt/GSK3ß/Snail signaling pathway in breast cancer. This study suggests that AnxA2 is a potential anti­invasion/metastasis target for therapeutic intervention in breast cancer.

Influence of the R823W mutation on the interaction of the ANKS6-ANKS3: insights from molecular dynamics simulation and free energy analysis.

The sterile alpha motif (SAM) domain of the protein ANKS6, a protein-protein interaction domain, is responsible for autosomal dominant polycystic kidney disease. Although the disease is the result of the R823W point mutation in the SAM domain of the protein ANKS6, the molecular details are still unclear. We applied molecular dynamics simulations, the principal component analysis, and the molecular mechanics Poisson-Boltzmann surface area binding free energy calculation to explore the structural and dynamic effects of the R823W point mutation on the complex ANKS6-ANKS3 (PDB ID: 4NL9) in comparison to the wild proteins. The energetic analysis presents that the wild type has a more stable structure than the mutant. The R823W point mutation not only disrupts the structure of the ANKS6 SAM domain but also negatively affects the interaction of the ANKS6-ANKS3. These results further clarify the previous experiments to understand the ANKS6-ANKS3 interaction comprehensively. In summary, this study would provide useful suggestions to understand the interaction of these proteins and their fatal action on mediating kidney function.

Thermal stability and unfolding pathways of hyperthermophilic and mesophilic periplasmic binding proteins studied by molecular dynamics simulation.

The ribose binding protein (RBP), a sugar-binding periplasmic protein, is involved in the transport and signaling processes in both prokaryotes and eukaryotes. Although several cellular and structural studies have been reported, a description of the thermostability of RBP at the molecular level remains elusive. Focused on the hyperthermophilic Thermoytoga maritima RBP (tmRBP) and mesophilic Escherichia coli homolog (ecRBP), we applied molecular dynamics simulations at four different temperatures (300, 380, 450, and 500 K) to obtain a deeper insight into the structural features responsible for the reduced thermostability of the ecRBP. The simulations results indicate that there are distinct structural differences in the unfolding pathway between the two homologs and the ecRBP unfolds faster than the hyperthermophilic homologs at certain temperatures in accordance with the lower thermal stability found experimentally. Essential dynamics analysis uncovers that the essential subspaces of ecRBP and tmRBP are non-overlapping and these two proteins show different directions of motion within the simulations trajectories. Such an understanding is required for designing efficient proteins with characteristics for a particular application.