Assembly of Biomolecular Gigastructures and Visualization with the Vulkan Graphics API

Enhanced Thermostability of Nattokinase by Computation-Based Rational Redesign of Flexible Regions

Nattokinase is a nutrient in healthy food natto that has the function of preventing and treating blood thrombus. However, its low thermostability and fibrinolytic activity limit its application in food and pharmaceuticals. In this study, we used bioinformatics analysis to identify two loops (loop10 and loop12) in the flexible region of nattokinase rAprY. Using this basis, we screened the G131S-S161T variant, which showed a 2.38-fold increase in half-life at 55 °C, and the M3 variant, which showed a 2.01-fold increase in activity, by using a thermostability prediction algorithm. Bioinformatics analysis revealed that the enhanced thermostability of the G131S-S161T variant was due to the increased rigidity and structural shrinkage of the overall structure. Additionally, the increased rigidity of the local region surrounding the active center and its mutated sites helps maintain its normal conformation in high-temperature environments. The increased catalytic activity of the M3 variant may be due to its more efficient substrate binding mechanism. We investigated strategies to improve the thermostability and fibrinolytic activity of nattokinase, and the resulting variants show promise for industrial production and application.

A refined picture of the native amine dehydrogenase family revealed by extensive biodiversity screening

Native amine dehydrogenases offer sustainable access to chiral amines, so the search for scaffolds capable of converting more diverse carbonyl compounds is required to reach the full potential of this alternative to conventional synthetic reductive aminations. Here we report a multidisciplinary strategy combining bioinformatics, chemoinformatics and biocatalysis to extensively screen billions of sequences in silico and to efficiently find native amine dehydrogenases features using computational approaches. In this way, we achieve a comprehensive overview of the initial native amine dehydrogenase family, extending it from 2,011 to 17,959 sequences, and identify native amine dehydrogenases with non-reported substrate spectra, including hindered carbonyls and ethyl ketones, and accepting methylamine and cyclopropylamine as amine donor. We also present preliminary model-based structural information to inform the design of potential (R)-selective amine dehydrogenases, as native amine dehydrogenases are mostly (S)-selective. This integrated strategy paves the way for expanding the resource of other enzyme families and in highlighting enzymes with original features.

Cystine-cored diphenylalanine appended peptide-based self-assembled fluorescent nanostructures direct redox-responsive drug delivery

Fabrication of stimuli-responsive superstructure capable of delivering chemotherapeutics directly to the cancer cell by sparing healthy cells is crucial. Herein, we developed redox-responsive hollow spherical assemblies through self-assembly of disulfide-linked cysteine-diphenylalanine (SN). These fluorescent hollow spheres display intrinsic green fluorescence, are proteolytically stable and biocompatible, and allow for real-time monitoring of their intracellular entry. The disulfide bond facilitates selective degradation in the presence of high glutathione (GSH) concentrations, prevalent in cancer cells. We achieved efficient encapsulation (68.72%) of the anticancer drug doxorubicin (Dox) and demonstrated GSH-dependent, redox-responsive drug release within cancerous cells. SN-Dox exhibited a 20-fold lower effective concentration (2.5 μM) for compromising breast cancer cell viability compared to non-malignant cells (50 μM). The ability of SN-Dox to initiate DNA damage signaling and trigger apoptosis was comparable to that of the unencapsulated drug. Our findings highlight the potential of SN for creating site-specific drug delivery vehicles for sustained therapeutic release.

From complex data to clear insights: visualizing molecular dynamics trajectories

Advances in simulations, combined with technological developments in high-performance computing, have made it possible to produce a physically accurate dynamic representation of complex biological systems involving millions to billions of atoms over increasingly long simulation times. The analysis of these computed simulations is crucial, involving the interpretation of structural and dynamic data to gain insights into the underlying biological processes. However, this analysis becomes increasingly challenging due to the complexity of the generated systems with a large number of individual runs, ranging from hundreds to thousands of trajectories. This massive increase in raw simulation data creates additional processing and visualization challenges. Effective visualization techniques play a vital role in facilitating the analysis and interpretation of molecular dynamics simulations. In this paper, we focus mainly on the techniques and tools that can be used for visualization of molecular dynamics simulations, among which we highlight the few approaches used specifically for this purpose, discussing their advantages and limitations, and addressing the future challenges of molecular dynamics visualization.

A Concise Review of Biomolecule Visualization

The structural characteristics of biomolecules are a major focus in the field of structural biology. Molecular visualization plays a crucial role in displaying structural information in an intuitive manner, aiding in the understanding of molecular properties. This paper provides a comprehensive overview of core concepts, key techniques, and tools in molecular visualization. Additionally, it presents the latest research findings to uncover emerging trends and highlights the challenges and potential directions for the development of the field.

Iron uptake pathway of Escherichia coli as an entry route for peptide nucleic acids conjugated with a siderophore mimic

Bacteria secrete various iron-chelators (siderophores), which scavenge Fe3+ from the environment, bind it with high affinity, and retrieve it inside the cell. After the Fe3+ uptake, bacteria extract the soluble iron(II) from the siderophore. Ferric siderophores are transported inside the cell via the TonB-dependent receptor system. Importantly, siderophore uptake paths have been also used by sideromycins, natural antibiotics. Our goal is to hijack the transport system for hydroxamate-type siderophores to deliver peptide nucleic acid oligomers into Escherichia coli cells. As siderophore mimics we designed and synthesized linear and cyclic Nδ-acetyl-Nδ-hydroxy-l-ornithine based peptides. Using circular dichroism spectroscopy, we found that iron(III) is coordinated by the linear trimer with hydroxamate groups but not by the cyclic peptide. The internal flexibility of the linear siderophore oxygen atoms and their interactions with Fe3+ were confirmed by all-atom molecular dynamics simulations. Using flow cytometry we found that the designed hydroxamate trimer transports PNA oligomers inside the E. coli cells. Growth recovery assays on various E. coli mutants suggest the pathway of this transport through the FhuE outer-membrane receptor, which is responsible for the uptake of the natural iron chelator, ferric-coprogen. This pathway also involves the FhuD periplasmic binding protein. Docking of the siderophores to the FhuE and FhuD receptor structures showed that binding of the hydroxamate trimer is energetically favorable corroborating the experimentally suggested uptake path. Therefore, this siderophore mimic, as well as its conjugate with PNA, is most probably internalized through the hydroxamate pathway.

Improvement of the fibrinolytic activity, acid resistance and thermostability of nattokinase by surface charge engineering

Nattokinase is a promising thrombolytic drug due to its powerful fibrinolytic effect and few side effects. However, the low fibrinolytic activity and stability of nattokinase have limited its industrial production and oral application. In this study, the basic and neutral amino acid residues on the surface of recombinant nattokinase AprY from Bacillus mojavensis LY-06 (rAprY) were mutated to acidic amino acid residues by surface charge engineering strategy, and two variants K12D and N109D with 92.6 % and 8.4 % increased fibrinolytic activity were obtained. The R45E variant with enhanced acid stability and thermostability was also screened, its acid stability at pH 4 and t1/2 at 55 °C were 3.7-fold and 1.8-fold higher than that of wild type rAprY, respectively. Bioinformatics analysis showed that the increased activities of K12D and N109D variants were related to the increased flexibility of the region around their active centers. The increased rigidity of 97-103 amino acid residues around the active center of R45E may be the reason for its enhanced stability and reduced catalytic activity. The multipoint mutation K12D-N109D (M2)'s catalytic activity did not increase cumulatively, but its pH stability did. The nattokinase variants generated in this study have potential for industrial production and application.

YASARA Model-Interactive Molecular Modeling from Two Dimensions to Virtual Realities

The industry's transition from three-dimensional (3D) glasses to virtual reality (VR) headsets has left modelers stranded without hardware supply, since walking around and waving arms in a virtual world is a great experience but also very tiring when doing time-intensive modeling work. We present a novel software implementation that uses a VR headset while sitting at a desk in front of the normal screen, which is beamed into the virtual reality together with keyboard, mouse, and chair using the headset's cameras and an extra tracker attached to the seat-back. Compared to 3D glasses, this yields a comparably relaxing but much more immersive workplace and provides additional possibilities such as taking molecules into one's hands, standing up, and walking or teleporting through the models. This VR functionality has been combined with a molecular graphics engine based on Vulkan, a next-generation cross-platform application programming interface (API) for GPUs and the successor of the widely used Open Graphics Library (OpenGL). It is built into the YASARA Model program, which includes many features like small and large molecule builders, electron densities, partial surfaces, contact analysis, coordinate manipulation, and animations. Interactive tutorials are provided to guide modelers into VR and familiarize them with the molecular modeling features. YASARA Model is available for Linux, Windows, Android, and MacOS (the latter without VR) with an introductory video at www.YASARA.org/vr.

β-Lactam TRPM8 Antagonists Derived from Phe-Phenylalaninol Conjugates: Structure-Activity Relationships and Antiallodynic Activity

The protein transient receptor potential melastatin type 8 (TRPM8), a non-selective, calcium (Ca2+)-permeable ion channel is implicated in several pathological conditions, including neuropathic pain states. In our previous research endeavors, we have identified β-lactam derivatives with high hydrophobic character that exhibit potent and selective TRPM8 antagonist activity. This work describes the synthesis of novel derivatives featuring C-terminal amides and diversely substituted N'-terminal monobenzyl groups in an attempt to increase the total polar surface area (TPSA) in this family of compounds. The primary goal was to assess the influence of these substituents on the inhibition of menthol-induced cellular Ca2+ entry, thereby establishing critical structure-activity relationships. While the substitution of the tert-butyl ester by isobutyl amide moieties improved the antagonist activity, none of the N'-monobencyl derivatives, regardless of the substituent on the phenyl ring, achieved the activity of the model dibenzyl compound. The antagonist potency of the most effective compounds was subsequently verified using Patch-Clamp electrophysiology experiments. Furthermore, we evaluated the selectivity of one of these compounds against other members of the transient receptor potential (TRP) ion channel family and some receptors connected to peripheral pain pathways. This compound demonstrated specificity for TRPM8 channels. To better comprehend the potential mode of interaction, we conducted docking experiments to uncover plausible binding sites on the functionally active tetrameric protein. While the four main populated poses are located by the pore zone, a similar location to that described for the N-(3-aminopropyl)-2-[(3-methylphenyl)methoxy]-N-(2-thienylmethyl)-benzamide (AMTB) antagonist cannot be discarded. Finally, in vivo experiments, involving a couple of selected compounds, revealed significant antinociceptive activity within a mice model of cold allodynia induced by oxaliplatin (OXA).

In Silico Analysis, Anticonvulsant Activity, and Toxicity Evaluation of Schisandrin B in Zebrafish Larvae and Mice

The aim of this study is to evaluate the anticonvulsant potential of schisandrin B, a main ingredient of Schisandra chinensis extracts. Schisandrin B showed anticonvulsant activity in the zebrafish larva pentylenetetrazole acute seizure assay but did not alter seizure thresholds in the intravenous pentylenetetrazole test in mice. Schisandrin B crosses the blood-brain barrier, which we confirmed in our in silico and in vivo analyses; however, the low level of its unbound fraction in the mouse brain tissue may explain the observed lack of anticonvulsant activity. Molecular docking revealed that the anticonvulsant activity of the compound in larval zebrafish might have been due to its binding to a benzodiazepine site within the GABAA receptor and/or the inhibition of the glutamate NMDA receptor. Although schisandrin B showed a beneficial anticonvulsant effect, toxicological studies revealed that it caused serious developmental impairment in zebrafish larvae, underscoring its teratogenic properties. Further detailed studies are needed to precisely identify the properties, pharmacological effects, and safety of schisandrin B.

Combined Computer-Aided Predictors to Improve the Thermostability of Nattokinase

Food-derived nattokinase has strong thrombolytic activity and few side effects. In the field of medicine, nattokinase has been developed as an adjuvant drug for the treatment of thrombosis, and nattokinase-rich beverages and health foods have also shown great potential in the field of food development. At present, the poor thermostability of nattokinase limits its industrial production and application. In this study, we used several thermostability-prediction algorithms to predict nattokinase from Bacillus mojavensis LY-06 (AprY), and screened two variants S33T and T174V with increased thermostability and fibrinolytic activity. The t1/2 of S33T and T174V were 8.87-fold and 2.51-fold those of the wild type AprY, respectively, and their enzyme activities were also increased (1.17-fold and 1.28-fold, respectively). Although the thermostability of N218L was increased by 2.7 times, the fibrinolytic activity of N218L was only 73.3% of that of wild type AprY. The multiple-point mutation results showed that S33T-N218L and S33T-T174V-N218L variants lost their activity, and the T174V-N218L variant did not show any significant change in catalytic performance, while S33T-T174V increased its thermostability and activity by 21.3% and 24.8%, respectively. Although the S33T-T174V variant did not show the additive effect of thermostability, it combined the excellent transient thermostability of S33T with the better thrombolytic activity of T174V. Bioinformatics analysis showed that the overall structure of S33T and T174V variants tended to be stable, while the structure of S33T-T174V variant was more flexible. Local structure analysis showed that the increased rigidity of the active center region (positions 64-75) and the key loop region (positions 129-130, 155-163, 187-192, 237-241, and 268-270) determined the increased thermostability of all variants. In addition, the enhanced flexibility of S33T-T174V variant in the Ca1 binding region (positions 1-4, 75-82) and the peripheral region of the catalytic pocket (positions 210-216) may account for the inability to superpose its thermostability. We explored the effective strategy to enhance the thermostability of nattokinase, and the resulting variants have potential industrial production and application.

The coat protein of tobacco mosaic virus as an anti-tobacco mosaic virus: a molecular dynamics simulation

The Coat Protein (CP) of the Tobacco Mosaic Virus (TMV) executes an important duty in the protection of virus RNA. The interaction between the virus CP and host plant proteins induces infection in the host and creates dark and light green mosaics on crops, which disturb the growth and function of the plant. The interaction between the virus CP and the modified CP, expressed in transgenic plants, causes Coat Protein-Mediated Resistance (CP-MR), which reduces virus infection in transgenic plants. In this study, a model is suggested for resistance as "stop assembly of CP" in the virus. It is based on the fact that the CP, when mutated, acts as a dead-end in virus assembly. For evaluation of the model, we investigated the effect of four mutants including CBT28I, ABT42W, ABD77R, and ABT89W complexes on plant resistance against TMV infection by molecular dynamics simulation. Previous studies had shown the influence of such mutations on the CP-MR. The MD results of in the present study further confirmed the mentioned effect and demonstrated how the mutations could be the cause of CP-MR. The results are calculated by the RMSD, Rg, H-bond, and g-MMPBSA scripts. The change in binding energy between two chains is consistent with CP-MR such that with increase in binding energy, the affinity between two chains was reduced and the CP-MR increased. Based on this model, it is possible to design mutants with a high level of efficiency.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 proteases Mpro and PLpro: Design of inhibitors with predicted high potency and low mammalian toxicity using artificial neural networks, ligand-protein docking, molecular dynamics simulations, and ADMET calculations

The main (Mpro) and papain-like (PLpro) proteases are highly conserved viral proteins essential for replication of the COVID-19 virus, SARS-COV-2. Therefore, a logical plan for producing new drugs against this pathogen is to discover inhibitors of these enzymes. Accordingly, the goal of the present work was to devise a computational approach to design, characterize, and select compounds predicted to be potent dual inhibitors - effective against both Mpro and PLpro. The first step employed LigDream, an artificial neural network, to create a virtual ligand library. Ligands with computed ADMET profiles indicating drug-like properties and low mammalian toxicity were selected for further study. Initial docking of these ligands into the active sites of Mpro and PLpro was done with GOLD, and the highest-scoring ligands were redocked with AutoDock Vina to determine binding free energies (ΔG). Compounds 89-00, 89-07, 89-32, and 89-38 exhibited favorable ΔG values for Mpro (-7.6 to -8.7 kcal/mol) and PLpro (-9.1 to -9.7 kcal/mol). Global docking of selected compounds with the Mpro dimer identified prospective allosteric inhibitors 89-00, 89-27, and 89-40 (ΔG -8.2 to -8.9 kcal/mol). Molecular dynamics simulations performed on Mpro and PLpro active site complexes with the four top-scoring ligands from Vina demonstrated that the most stable complexes were formed with compounds 89-32 and 89-38. Overall, the present computational strategy generated new compounds with predicted drug-like characteristics, low mammalian toxicity, and high inhibitory potencies against both target proteases to form stable complexes. Further preclinical studies will be required to validate the in silico findings before the lead compounds could be considered for clinical trials.

Synthesis, Structural and Behavioral Studies of Indole Derivatives D2AAK5, D2AAK6 and D2AAK7 as Serotonin 5-HT1A and 5-HT2A Receptor Ligands

Serotonin receptors are involved in a number of physiological functions and regulate aggression, anxiety, appetite, cognition, learning, memory, mood, nausea, sleep, and thermoregulation. Here we report synthesis and detailed structural and behavioral studies of three indole derivatives: D2AAK5, D2AAK6, and D2AAK7 as serotonin 5-HT_1A and 5-HT_2A receptor ligands. X-ray studies revealed that the D2AAK5 compound crystallizes in centrosymmetric triclinic space group with one molecule in the asymmetric unit. The main interaction between the ligands and the receptors is the salt bridge between the protonatable nitrogen atom of the ligands and the conserved Asp (3.32) of the receptors. The complexes were stable in the molecular dynamic simulations. MD revealed that the studied ligands are relatively stable in their binding sites, with the exception of D2AAK7 in the serotonin 5-HT_1A receptor. D2AAK7 exerts anxiolytic activity in the EPM test, while D2AAK5 has a beneficial effect on the memory processes in the PA test.

Synthesis, Docking Studies and Pharmacological Evaluation of Serotoninergic Ligands Containing a 5-Norbornene-2-Carboxamide Nucleus

A new series of 5-norbornene-2-carboxamide derivatives was prepared and their affinities to the 5-HT_1A, 5-HT_2A, and 5-HT_2C receptors were evaluated and compared to a previously synthesized series of derivatives characterized by exo-N-hydroxy-5-norbornene-2,3-dicarboximidenucleus, in order to identify selective ligands for the above-mentioned subtype receptors. Arylpiperazines represents one of the most important classes of 5-HT_1AR ligands, and recent research concerning new derivatives has been focused on the modification of one or more portions of such pharmacophore. The combination of structural elements (heterocyclic nucleus, propyl chain and 4-substituted piperazine), known to be critical to the affinity to 5-HT_1A receptors, and the proper selection of substituents led to compounds with high specificity and affinity towards serotoninergic receptors. The most active compounds were selected for further in vivo assays to determine their functional activity. Finally, to rationalize the obtained results, molecular docking studies were performed. The results of the pharmacological studies showed that Norbo-4 and Norbo-18 were the most active and promising derivatives for the serotonin receptor considered in this study.

Mask-Point: Automatic 3D Surface Defects Detection Network for Fiber-Reinforced Resin Matrix Composites

Surface defects of fiber-reinforced resin matrix composites (FRRMCs) adversely affect their appearance and performance. To accurately and efficiently detect the three-dimensional (3D) surface defects of FRRMCs, a novel lightweight and two-stage semantic segmentation network, i.e., Mask-Point, is proposed. Stage 1 of Mask-Point is the multi-head 3D region proposal extractors (RPEs), generating several 3D regions of interest (ROIs). Stage 2 is the 3D aggregation stage composed of the shared classifier, shared filter, and non-maximum suppression (NMS). The two stages work together to detect the surface defects. To evaluate the performance of Mask-Point, a new 3D surface defects dataset of FRRMCs containing about 120 million points is produced. Training and test experiments show that the accuracy and the mean intersection of union (mIoU) increase as the number of different 3D RPEs increases in Stage 1, but the inference speed becomes slower when the number of different 3D RPEs increases. The best accuracy, mIoU, and inference speed of the Mask-Point model could reach 0.9997, 0.9402, and 320,000 points/s, respectively. Moreover, comparison experiments also show that Mask-Point offers relatively the best segmentation performance compared with several other typical 3D semantic segmentation networks. The mIoU of Mask-Point is about 30% ahead of the sub-optimal 3D semantic segmentation network PointNet. In addition, a distributed surface defects detection system based on Mask-Point is developed. The system is applied to scan real FRRMC products and detect their surface defects, and it achieves the relatively best detection performance in competition with skilled human workers. The above experiments demonstrate that the proposed Mask-Point could accurately and efficiently detect 3D surface defects of FRRMCs, and the Mask-Point also provides a new potential solution for the 3D surface defects detection of other similar materials

The Role of Lipids in Allosteric Modulation of Dopamine D2 Receptor-In Silico Study

The dopamine D₂ receptor, belonging to the class A G protein-coupled receptors (GPCRs), is an important drug target for several diseases, including schizophrenia and Parkinson’s disease. The D₂ receptor can be activated by the natural neurotransmitter dopamine or by synthetic ligands, which in both cases leads to the receptor coupling with a G protein. In addition to receptor modulation by orthosteric or allosteric ligands, it has been shown that lipids may affect the behaviour of membrane proteins. We constructed a model of a D₂ receptor with a long intracellular loop (ICL3) coupled with G_iα1 or G_iα2 proteins, embedded in a complex asymmetric membrane, and simulated it in complex with positive, negative or neutral allosteric ligands. In this study, we focused on the influence of ligand binding and G protein coupling on the membrane–receptor interactions. We show that there is a noticeable interplay between the cell membrane, G proteins, D₂ receptor and its modulators.