Angiogenesis is an essential process in tumorigenesis, tumor invasion, and metastasis, and is an intriguing pathway for drug discovery. Targeting vascular endothelial growth factor receptor 2 (VEGFR2) to inhibit tumor angiogenic pathways has been widely explored and adopted in clinical practice. However, most drugs, such as the Food and Drug Administration -approved drug axitinib (ATC code: L01EK01), have considerable side effects and limited tolerability. Therefore, there is an urgent need for the development of novel VEGFR2 inhibitors. In this study, we propose a novel strategy to design potential candidates targeting VEGFR2 using three-dimensional (3D) deep learning and structural modeling methods. A geometric-enhanced molecular representation learning method (GEM) model employing a graph neural network (GNN) as its underlying predictive algorithm was used to predict the activity of the candidates. In the structural modeling method, flexible docking was performed to screen data with high affinity and explore the mechanism of the inhibitors. Small -molecule compounds with consistently improved properties were identified based on the intersection of the scores obtained from both methods. Candidates identified using the GEM-GNN model were selected for in silico modeling using molecular dynamics simulations to further validate their efficacy. The GEM-GNN model enabled the identification of candidate compounds with potentially more favorable properties than the existing drug, axitinib, while achieving higher efficacy.

In the current study, we utilized molecular modeling and simulation approaches to define putative potential molecular targets for Burdock Inulin, including inflammatory proteins such as iNOS, COX-2, TNF-alpha, IL-6, and IL-1β. Molecular docking results revealed potential interactions and good binding affinity for these targets; however, IL-1β, COX-2, and iNOS were identified as the best targets for Inulin. Molecular simulation-based stability assessment demonstrated that inulin could primarily target iNOS and may also supplementarily target COX-2 and IL-1β during DSS-induced colitis to reduce the role of these inflammatory mechanisms. Furthermore, residual flexibility, hydrogen bonding, and structural packing were reported with uniform trajectories, showing no significant perturbation throughout the simulation. The protein motions within the simulation trajectories were clustered using principal component analysis (PCA). The IL-1β-Inulin complex, approximately 70% of the total motion was attributed to the first three eigenvectors, while the remaining motion was contributed by the remaining eigenvectors. In contrast, for the COX2-Inulin complex, 75% of the total motion was attributed to the eigenvectors. Furthermore, in the iNOS-Inulin complex, the first three eigenvectors contributed to 60% of the total motion. Furthermore, the iNOS-Inulin complex contributed 60% to the total motion through the first three eigenvectors. To explore thermodynamically favorable changes upon mutation, motion mode analysis was carried out. The Free Energy Landscape (FEL) results demonstrated that the IL-1β-Inulin achieved a single conformation with the lowest energy, while COX2-Inulin and iNOS-Inulin exhibited two lowest-energy conformations each. IL-1β-Inulin and COX2-Inulin displayed total binding free energies of - 27.76 kcal/mol and - 37.78 kcal/mol, respectively, while iNOS-Inulin demonstrated the best binding free energy results at - 45.89 kcal/mol. This indicates a stronger pharmacological potential of iNOS than the other two complexes. Thus, further experiments are needed to use inulin to target iNOS and reduce DSS-induced colitis and other autoimmune diseases.

Small non-coding RNAs (sRNAs) play vital roles in gene expression regulation and RNA interference. To comprehend their molecular mechanisms and develop therapeutic approaches, determining the accurate three-dimensional structure of sRNAs is crucial. Although nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for structural biology, obtaining high-resolution structures of sRNAs using NMR data alone can be challenging. In such cases, structural modeling can provide additional details about RNA structures. In this context, we present a protocol for the structural modeling of sRNA using the SimRNA method based on sparse NMR constraints. To demonstrate the efficacy of our method, we provide selected examples of NMR spectra and RNA structures, specifically for the second stem-loop of DsrA sRNA.

Polyketide retrobiosynthesis, where the biosynthetic pathway of a given polyketide can be reversibly engineered due to the colinearity of the polyketide synthase (PKS) structure and function, has the potential to produce millions of organic molecules. Mixing and matching modules from natural PKSs is one of the routes to produce many of these molecules. Evolutionary analysis of PKSs suggests that traditionally used module boundaries may not lead to the most productive hybrid PKSs and that new boundaries around and within the ketosynthase domain may be more active when constructing hybrid PKSs. As this is still a nascent area of research, the generality of these design principles based on existing engineering efforts remains inconclusive. Recent advances in structural modeling and synthetic biology present an opportunity to accelerate PKS engineering by re-evaluating insights gained from previous engineering efforts with cutting edge tools.

Essential oils (EOs) and their volatile secondary metabolites have been proved to be effective on storage pests control, while restricted on the application due to unclear mechanism. Molecular dynamics (MD) simulations and binding free energies analysis provided an effective approach to reveal mechanism on conformational calculation. In this work, the insecticidal and repellent capacities of Praxelis clematidea and Ageratum houstonianum oils and their main components identified by gas chromatography-mass spectrometry (GC-MS) were scientifically measured. Interestingly, P. clematidea oil exhibited strong fumigant toxicity against Tribolium castaneum (LC50 = 7.07 mg/L air). Moreover, two EOs exhibited over 80% repellent rate against T. castaneum at the highest concentration of 78.63 nL/cm2. Furthermore, hundreds of enzymes related to the regulation of biological processes of T. castaneum were screened to explore the underlying molecular mechanism and develop promising insecticides. Besides, top hits were subjected to MD simulations and binding free energies analysis to elucidate complex inter-molecular stability and affinity over simulated time. The results demonstrated that isolongifolene, δ-cadinene, β-caryophyllene and caryophyllene oxide were prioritized as they were establishing conserved and stable interactions with residues of nuclear hormone receptor 3 (TcHR3) of T. castaneum, which suggested that the four sesquiterpenes have potential to be promising insecticides on storage pests control.

α-Ketoglutarate decarboxylase is a crucial enzyme in the tricarboxylic acid cycle of cyanobacteria, catalyzing the non-oxidative decarboxylation of α-ketoglutarate to produce succinate semialdehyde and CO2. The decarboxylation process is reliant on the cofactor of thiamine diphosphate. However, this enzyme\'s biochemical and structural properties have not been well characterized. In this work, two α-ketoglutarate decarboxylases encoded by MAE_06010 and MiAbw_01735 genes from Microcystis aeruginosa NIES-843 (MaKGD) and NIES-4325 (MiKGD), respectively, were overexpressed and purified by using an Escherichia coli expression system. It was found that MaKGD exhibited 9.2-fold higher catalytic efficiency than MiKGD, which may be attributed to the absence of glutamate decarboxylase in Microcystis aeruginosa NIES-843. Further biochemical investigation of MaKGD demonstrated that it displayed optimum activity at pH 6.5-7.0 and was most activated by Mg2+. Additionally, MaKGD showed substrate specificity towards α-ketoglutarate. Structural modeling and autodocking results revealed that the active site of MaKGD contained a distinct binding pocket where α-ketoglutarate and thiamine diphosphate interacted with specific amino acid residues via hydrophobic interactions, hydrogen bonds and salt bridges. Furthermore, the mutagenesis study provided strong evidence supporting the importance of certain residues in the catalysis of MaKGD. These findings provide new insights into the structure-function relationships of α-ketoglutarate decarboxylases from cyanobacteria.

Rift valley fever virus (RVFV) is the causative agent of a viral zoonosis that causes a significant clinical burden in domestic and wild ruminants. Major outbreaks of the virus occur in livestock, and contaminated animal products or arthropod vectors can transmit the virus to humans. The viral RNA-dependent RNA polymerase (RdRp; L protein) of the RVFV is responsible for viral replication and is thus an appealing drug target because no effective and specific vaccine against this virus is available. The current study reported the structural elucidation of the RVFV-L protein by in-depth homology modeling since no crystal structure is available yet. The inhibitory binding modes of known potent L protein inhibitors were analyzed. Based on the results, further molecular docking-based virtual screening of Selleckchem Nucleoside Analogue Library (156 compounds) was performed to find potential new inhibitors against the RVFV L protein. ADME (Absorption, Distribution, Metabolism, and Excretion) and toxicity analysis of these compounds was also performed. Besides, the binding mechanism and stability of identified compounds were confirmed by a 50 ns molecular dynamic (MD) simulation followed by MM/PBSA binding free energy calculations. Homology modeling determined a stable multi-domain structure of L protein. An analysis of known L protein inhibitors, including Monensin, Mycophenolic acid, and Ribavirin, provide insights into the binding mechanism and reveals key residues of the L protein binding pocket. The screening results revealed that the top three compounds, A-317491, Khasianine, and VER155008, exhibited a high affinity at the L protein binding pocket. ADME analysis revealed good pharmacodynamics and pharmacokinetic profiles of these compounds. Furthermore, MD simulation and binding free energy analysis endorsed the binding stability of potential compounds with L protein. In a nutshell, the present study determined potential compounds that may aid in the rational design of novel inhibitors of the RVFV L protein as anti-RVFV drugs.

A simple heat treatment method was used to optimize the three-dimensional network structure of the hydrophobic aerogel, and during the heat treatment process at 200-1000 °C, the thermal conductivity of the aerogel reached the lowest to 0.02240 W/m·K between 250 °C and 300 °C, which was mainly due to the optimization of microstructure and pyrolysis of surface groups. Further Fluent heat-transfer simulation also confirmed the above results. Synchrotron vacuum ultraviolet photoionization mass spectrometry (SVUV-PIMS) was used to finely measure the pyrolysis process of aerogels, and the pyrolysis process of aerogel was divided into four stages. (I) Until 419 °C, as the temperature continued to rise, surface methyl groups were oxidized to form hydroxyl. (II) As the temperature reached to 232 °C, the oxidation proceeded. In addition, inside the aerogel, because of lacking oxygen, the reaction produced CH4 and C-Si bonds would form. (III) After 283 °C, Si-OH groups began to condense to form Si-O-Si, which optimized the three-dimensional network structures to be beneficial to improve the thermal insulation performance of silica aerogel. (IV) When it reached 547 °C, the chemical reaction was terminated, and all the primary particles gradually fused into secondary particles and sintered to form clusters.

The inter-subunit interaction at the protein interfaces plays a key role in protein self-assembly, through which enabling protein self-assembly controllable is of great importance for preparing the novel nanoscale protein materials with unexplored properties. Different from normal 24-meric ferritin, archaeal ferritin, Thermotoga maritima ferritin (TmFtn) naturally occurs as a dimer, which can assemble into a 24-mer nanocage induced by salts. However, the regulation mechanism of protein self-assembly underlying this phenomenon remains unclear. Here, a combination of the computational energy simulation and key interface reconstruction revealed that a short helix involved interactions at the C4 interface are mainly responsible for the existence of such dimer. Agreeing with this idea, deletion of such short helix of each subunit triggers it to be a stable dimer, which losses the ability to reassemble into 24-meric ferritin in the presence of salts in solution. Further support for this idea comes from the observation that grafting a small helix from human H ferritin onto archaeal subunit resulted in a stable 24-mer protein nanocage even in the absence of salts. Thus, these findings demonstrate that adjusting the interactions at the protein interfaces appears to be a facile, effective approach to control subunit assembly into different protein architectures.

The human zinc finger NFX1-type containing 1 (ZNFX1) is an interferon-stimulated protein associated to the outer mitochondrial membrane, able to bind dsRNAs and interact with MAVS proteins, promoting type I IFN response in the early stage of viral infection. An N-terminal Armadillo (ARM)-type fold and a large helicase core (P-loop) and zinc fingers confer RNA-binding and ATPase activities to ZNFX1. We studied the phylogenetic distribution of metazoan ZNFX1s, ZNFX1 gene expression trends and genomic and protein signatures during viral infection of invertebrates. Based on 221 ZNFX1 sequences, we obtained a polyphyletic tree with a taxonomy-consistent branching at the phylum-level only. In metazoan genomes, ZNFX1 genes were found either in single copy, with up to some tens of exons in vertebrates, or in multiple copies, with one or a few exons and one of them sometimes encompassing most of the coding sequence, in invertebrates like sponges, sea urchins and mollusks. Structural analyses of selected ZNFX1 proteins showed high conservation of the helicase region (P-loop), an overall conserved region and domain architecture, an ARM-fold mostly traceable, and the presence of intrinsically disordered regions of varying length and position. The remarkable over-expression of ZNFX1 in bivalve and gastropod mollusks infected with dsDNA viruses underscores the antiviral role of ZNFX1, whereas nothing similar was found in virus-infected nematodes and corals. Whether the functional diversification reported in the C. elegans ZNFX1 occurs in other metazoan proteins remains to be established.