RCSB Protein Data Bank: powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences

Review of machine learning methods for RNA secondary structure prediction

Secondary structure plays an important role in determining the function of noncoding RNAs. Hence, identifying RNA secondary structures is of great value to research. Computational prediction is a mainstream approach for predicting RNA secondary structure. Unfortunately, even though new methods have been proposed over the past 40 years, the performance of computational prediction methods has stagnated in the last decade. Recently, with the increasing availability of RNA structure data, new methods based on machine learning (ML) technologies, especially deep learning, have alleviated the issue. In this review, we provide a comprehensive overview of RNA secondary structure prediction methods based on ML technologies and a tabularized summary of the most important methods in this field. The current pending challenges in the field of RNA secondary structure prediction and future trends are also discussed.

Preliminary modelling as the first stage of targeted organic synthesis

The review aims to present a classification and applicability analysis of methods for preliminary molecular modelling for targeted organic, catalytic and biocatalytic synthesis. The following three main approaches are considered as a primary classification of the methods: modelling of the target – ligand coordination without structural information on both the target and the resulting complex; calculations based on experimentally obtained structural information about the target; and dynamic simulation of the target – ligand complex and the reaction mechanism with calculation of the free energy of the reaction. The review is meant for synthetic chemists to be used as a guide for building an algorithm for preliminary modelling and synthesis of structures with specified properties.

The bibliography includes 353 references.

Camptothecin shows better promise than Curcumin in the inhibition of the Human Telomerase: A computational study

OBJECTIVES

The Human Telomerase enzyme has become a drug target in the treatment of cancers and age-related disorders. This study aims to identify potential natural inhibitors of the Human Telomerase from compounds derived from edible African plants.

MATERIALS AND METHODS

A library of 1,126 natural compounds was molecularly docked against the Telomerase Reverse Transcriptase (PDB ID: 5ugw), the catalytic subunit of the target protein. Curcumin, a known Telomerase inhibitor was used as the standard. The front-runner compounds were screened for bioavailability, pharmacokinetic properties, and bioactivity using the SWISSADME, PKCSM, and Molinspiration webservers respectively. The molecular dynamic simulation and analyses of the apo and holo proteins were performed by the Galaxy supercomputing webserver.

RESULTS

The results of the molecular docking and virtual screening reveal Augustamine and Camptothecin as lead compounds. Augustamine has better drug-likeness and pharmacokinetic properties while Camptothecin showed better bioactivity and stronger binding affinity (-8.2 kcal/mol) with the target. The holo structure formed by Camptothecin showed greater inhibitory activity against the target with a total RMSF of 169.853, B-Factor of 20.164, and 108 anti-correlating residues.

CONCLUSION

Though they both act at the same binding site, Camptothecin induces greater Telomerase inhibition and better molecular stability than the standard, Curcumin. Further tests are required to investigate the inhibitory activities of the lead compounds.

Protein-protein interactions at a glance: Protocols for the visualization of biomolecular interactions

Protein-protein interactions (PPIs) play a key role in many biological processes and are intriguing targets for drug discovery campaigns. Advancements in experimental and computational techniques are leading to a growth of data accessibility, and, with it, an increased need for the analysis of PPIs. In this respect, visualization tools are essential instruments to represent and analyze biomolecular interactions. In this chapter, we reviewed some of the available tools, highlighting their features, and describing their functions with practical information on their usage.

Loss of the mitochondrial phosphate carrier SLC25A3 induces remodeling of the cardiac mitochondrial protein acylome

Mitochondria are recognized as signaling organelles because, under stress, mitochondria can trigger various signaling pathways to coordinate the cell's response. The specific pathway(s) engaged by mitochondria in response to mitochondrial energy defects in vivo and in high-energy tissues like the heart are not fully understood. Here, we investigated cardiac pathways activated in response to mitochondrial energy dysfunction by studying mice with cardiomyocyte-specific loss of the mitochondrial phosphate carrier (SLC25A3), an established model that develops cardiomyopathy as a result of defective mitochondrial ATP synthesis. Mitochondrial energy dysfunction induced a striking pattern of acylome remodeling, with significantly increased post-translational acetylation and malonylation. Mass spectrometry-based proteomics further revealed that energy dysfunction-induced remodeling of the acetylome and malonylome preferentially impacts mitochondrial proteins. Acetylation and malonylation modified a highly interconnected interactome of mitochondrial proteins, and both modifications were present on the enzyme isocitrate dehydrogenase 2 (IDH2). Intriguingly, IDH2 activity was enhanced in SLC25A3-deleted mitochondria, and further study of IDH2 sites targeted by both acetylation and malonylation revealed that these modifications can have site-specific and distinct functional effects. Finally, we uncovered a novel crosstalk between the two modifications, whereby mitochondrial energy dysfunction-induced acetylation of sirtuin 5 (SIRT5), inhibited its function. Because SIRT5 is a mitochondrial deacylase with demalonylase activity, this finding suggests that acetylation can modulate the malonylome. Together, our results position acylations as an arm of the mitochondrial response to energy dysfunction and suggest a mechanism by which focal disruption to the energy production machinery can have an expanded impact on global mitochondrial function.

Using Budding Yeast to Identify Molecules That Block Cancer Cell 'Mitotic Slippage' Only in the Presence of Mitotic Poisons

Research on the budding yeast Saccharomyces cerevisiae has yielded fundamental discoveries on highly conserved biological pathways and yeast remains the best-studied eukaryotic cell in the world. Studies on the mitotic cell cycle and the discovery of cell cycle checkpoints in budding yeast has led to a detailed, although incomplete, understanding of eukaryotic cell cycle progression. In multicellular eukaryotic organisms, uncontrolled aberrant cell division is the defining feature of cancer. Some of the most successful classes of anti-cancer chemotherapeutic agents are mitotic poisons. Mitotic poisons are thought to function by inducing a mitotic spindle checkpoint-dependent cell cycle arrest, via the assembly of the highly conserved mitotic checkpoint complex (MCC), leading to apoptosis. Even in the presence of mitotic poisons, some cancer cells continue cell division via 'mitotic slippage', which may correlate with a cancer becoming refractory to mitotic poison chemotherapeutic treatments. In this review, knowledge about budding yeast cell cycle control is explored to suggest novel potential drug targets, namely, specific regions in the highly conserved anaphase-promoting complex/cyclosome (APC/C) subunits Apc1 and/or Apc5, and in a specific N-terminal region in the APC/C co-factor cell division cycle 20 (Cdc20), which may yield molecules which block 'mitotic slippage' only in the presence of mitotic poisons.

PDBeCIF: an open-source mmCIF/CIF parsing and processing package

Background

Biomacromolecular structural data outgrew the legacy Protein Data Bank (PDB) format which the scientific community relied on for decades, yet the use of its successor PDBx/Macromolecular Crystallographic Information File format (PDBx/mmCIF) is still not widespread. Perhaps one of the reasons is the availability of easy to use tools that only support the legacy format, but also the inherent difficulties of processing mmCIF files correctly, given the number of edge cases that make efficient parsing problematic. Nevertheless, to fully exploit macromolecular structure data and their associated annotations such as multiscale structures from integrative/hybrid methods or large macromolecular complexes determined using traditional methods, it is necessary to fully adopt the new format as soon as possible.

Results

To this end, we developed PDBeCIF, an open-source Python project for manipulating mmCIF and CIF files. It is part of the official list of mmCIF parsers recorded by the wwPDB and is heavily employed in the processes of the Protein Data Bank in Europe. The package is freely available both from the PyPI repository (http://pypi.org/project/pdbecif) and from GitHub (https://github.com/pdbeurope/pdbecif) along with rich documentation and many ready-to-use examples.

Conclusions

PDBeCIF is an efficient and lightweight Python 2.6+/3+ package with no external dependencies. It can be readily integrated with 3rd party libraries as well as adopted for broad scientific analyses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04271-9.

BENZ WS: the Bologna ENZyme Web Server for four-level EC number annotation

The Bologna ENZyme Web Server (BENZ WS) annotates four-level Enzyme Commission numbers (EC numbers) as defined by the International Union of Biochemistry and Molecular Biology (IUBMB). BENZ WS filters a target sequence with a combined system of Hidden Markov Models, modelling protein sequences annotated with the same molecular function, and Pfams, carrying along conserved protein domains. BENZ returns, when successful, for any enzyme target sequence an associated four-level EC number. Our system can annotate both monofunctional and polyfunctional enzymes, and it can be a valuable resource for sequence functional annotation.

A Brief History of De Novo Protein Design: Minimal, Rational, and Computational

Protein design has come of age, but how will it mature? In the 1980s and the 1990s, the primary motivation for de novo protein design was to test our understanding of the informational aspect of the protein-folding problem; i.e., how does protein sequence determine protein structure and function? This necessitated minimal and rational design approaches whereby the placement of each residue in a design was reasoned using chemical principles and/or biochemical knowledge. At that time, though with some notable exceptions, the use of computers to aid design was not widespread. Over the past two decades, the tables have turned and computational protein design is firmly established. Here, I illustrate this progress through a timeline of de novo protein structures that have been solved to atomic resolution and deposited in the Protein Data Bank. From this, it is clear that the impact of rational and computational design has been considerable: More-complex and more-sophisticated designs are being targeted with many being resolved to atomic resolution. Furthermore, our ability to generate and manipulate synthetic proteins has advanced to a point where they are providing realistic alternatives to natural protein functions for applications both in vitro and in cells. Also, and increasingly, computational protein design is becoming accessible to non-specialists. This all begs the questions: Is there still a place for minimal and rational design approaches? And, what challenges lie ahead for the burgeoning field of de novo protein design as a whole?

Antiproliferative Properties and G-Quadruplex-Binding of Symmetrical Naphtho[1,2-b:8,7-b']dithiophene Derivatives

Background: G-quadruplex (G4) forming sequences are recurrent in telomeres and promoter regions of several protooncogenes. In normal cells, the transient arrangements of DNA in G-tetrads may regulate replication, transcription, and translation processes. Tumors are characterized by uncontrolled cell growth and tissue invasiveness and some of them are possibly mediated by gene expression involving G-quadruplexes. The stabilization of G-quadruplex sequences with small molecules is considered a promising strategy in anticancer targeted therapy. Methods: Molecular virtual screening allowed us identifying novel symmetric bifunctionalized naphtho[1,2-b:8,7-b’]dithiophene ligands as interesting candidates targeting h-Telo and c-MYC G-quadruplexes. A set of unexplored naphtho-dithiophene derivatives has been synthesized and biologically tested through in vitro antiproliferative assays and spectroscopic experiments in solution. Results: The analysis of biological and spectroscopic data highlighted noteworthy cytotoxic effects on HeLa cancer cell line (GI₅₀ in the low μM range), but weak interactions with G-quadruplex c-MYC promoter. Conclusions: The new series of naphtho[1,2-b:8,7-b’]dithiophene derivatives, bearing the pharmacophoric assumptions necessary to stabilize G-quadruplexes, have been designed and successfully synthesized. The interesting antiproliferative results supported by computer aided rational approaches suggest that these studies are a significant starting point for a lead optimization process and the isolation of a more efficacious set of G-quadruplexes stabilizers.

Antithrombin Resistance Rescues Clotting Defect of Homozygous Prothrombin-Y510N Dysprothrombinemia

A patient with hematuria in our clinic was diagnosed with urolithiasis. Analysis of the patient's plasma clotting time indicated that both activated partial thromboplastin time (52.6 seconds) and prothrombin time (19.4 seconds) are prolonged and prothrombin activity is reduced to 12.4% of normal, though the patient exhibited no abnormal bleeding phenotype and a prothrombin antigen level of 87.9%. Genetic analysis revealed the patient is homozygous for prothrombin Y510N mutation. We expressed and characterized the prothrombin-Y510N variant in appropriate coagulation assays and found that the specificity constant for activation of the mutant zymogen by factor Xa is impaired approximately fivefold. Thrombin generation assay using patient's plasma and prothrombin-deficient plasma supplemented with either wild-type or prothrombin-Y510N revealed that both peak height and time to peak for the prothrombin mutant are decreased; however, the endogenous thrombin generation potential is increased. Further analysis indicated that the thrombin mutant exhibits resistance to antithrombin and is inhibited by the serpin with approximately 12-fold slower rate constant. Protein C activation by thrombin-Y510N was also decreased by approximately 10-fold; however, thrombomodulin overcame the catalytic defect. The Na⁺-concentration-dependence of the amidolytic activities revealed that the dissociation constant for the interaction of Na⁺ with the mutant has been elevated approximately 20-fold. These results suggest that Y510 (Y184a in chymotrypsin numbering) belongs to network of residues involved in binding Na⁺. A normal protein C activation by thrombin-Y510N suggests that thrombomodulin modulates the conformation of the Na⁺-binding loop of thrombin. The clotting defect of thrombin-Y510N appears to be compensated by its markedly lower reactivity with antithrombin, explaining patient's normal hemostatic phenotype.

Ligand Binding Site Comparison - LiBiSCo - a web-based tool for analyzing interactions between proteins and ligands to explore amino acid specificity within active sites

Interaction between protein and ligands are ubiquitous in a biological cell, and understanding these interactions at the atom level in protein-ligand complexes is crucial for structural bioinformatics and drug discovery. Here, we present a web-based protein-ligand interaction application named Ligand Binding Site Comparison (LiBiSCo) for comparing the amino acid residues interacting with atoms of a ligand molecule between different protein-ligand complexes available in the Protein Data Bank (PDB) database. The comparison is performed at the ligand atom level irrespectively of having binding site similarity or not between the protein structures of interest. The input used in LiBiSCo is one or several PDB IDs of protein-ligand complex(es) and the tool returns a list of identified interactions at ligand atom level including both bonded and non-bonded interactions. A sequence profile for the interaction for each ligand atoms is provided as a WebLogo. The LiBiSco is useful in understanding ligand binding specificity and structural promiscuity among families that are structurally unrelated. The LiBiSCo tool can be accessed through https://albiorix.bioenv.gu.se/LiBiSCo/HomePage.py.

Implications of the simple chemical structure of the odorant molecules interacting with the olfactory receptor 1A1

G protein‒coupled receptors (GPCRs), including olfactory receptors, account for the largest group of genes in the human genome and occupy a very important position in signaling systems. Although olfactory receptors, which belong to the broader category of GPCRs, play an important role in monitoring the organism’s surroundings, their actual three-dimensional structure has not yet been determined. Therefore, the specific details of the molecular interactions between the receptor and the ligand remain unclear. In this report, the interactions between human olfactory receptor 1A1 and its odorant molecules were simulated using computational methods, and we explored how the chemically simple odorant molecules activate the olfactory receptor.

Performance of Web tools for predicting changes in protein stability caused by mutations

BACKGROUND

Despite decades on developing dedicated Web tools, it is still difficult to predict correctly the changes of the thermodynamic stability of proteins caused by mutations. Here, we assessed the reliability of five recently developed Web tools, in order to evaluate the progresses in the field.

RESULTS

The results show that, although there are improvements in the field, the assessed predictors are still far from ideal. Prevailing problems include the bias towards destabilizing mutations, and, in general, the results are unreliable when the mutation causes a ΔΔG within the interval ± 0.5 kcal/mol. We found that using several predictors and combining their results into a consensus is a rough, but effective way to increase reliability of the predictions.

CONCLUSIONS

We suggest all developers to consider in their future tools the usage of balanced data sets for training of predictors, and all users to combine the results of multiple tools to increase the chances of having correct predictions about the effect of mutations on the thermodynamic stability of a protein.

IUPred3: prediction of protein disorder enhanced with unambiguous experimental annotation and visualization of evolutionary conservation

Intrinsically disordered proteins and protein regions (IDPs/IDRs) exist without a single well-defined conformation. They carry out important biological functions with multifaceted roles which is also reflected in their evolutionary behavior. Computational methods play important roles in the characterization of IDRs. One of the commonly used disorder prediction methods is IUPred, which relies on an energy estimation approach. The IUPred web server takes an amino acid sequence or a Uniprot ID/accession as an input and predicts the tendency for each amino acid to be in a disordered region with an option to also predict context-dependent disordered regions. In this new iteration of IUPred, we added multiple novel features to enhance the prediction capabilities of the server. First, learning from the latest evaluation of disorder prediction methods we introduced multiple new smoothing functions to the prediction that decreases noise and increases the performance of the predictions. We constructed a dataset consisting of experimentally verified ordered/disordered regions with unambiguous annotations which were added to the prediction. We also introduced a novel tool that enables the exploration of the evolutionary conservation of protein disorder coupled to sequence conservation in model organisms. The web server is freely available to users and accessible at https://iupred3.elte.hu.

Thermodynamic and kinetic insights into the interactions between functionalized CdTe quantum dots and human serum albumin: A surface plasmon resonance approach

To explore in vivo application of quantum dots (QDs), it is essential to understand the dynamics and energetics of interactions between QDs and proteins. Here, surface plasmon resonance (SPR) and molecular docking were employed to investigate the kinetics and thermodynamics of interactions between human serum albumin (HSA) and CdTe QDs (~3 nm) functionalized with mercaptopropionic acid (MPA) or thioglycolic acid (TGA). Kinetic analysis showed that HSA-QD interactions involved transition-complex formation. Despite the structural similarities between MPA and TGA, the [HSA-CdTe@TGA]^‡ formation by association of free HSA and QDs demanded 70% more energy and higher entropic gain (E_a-TGA^‡= 65.10 and T∆S_a-TGA^‡= 28.62 kJ mol^-1) than the formation of [HSA-CdTe@MPA]^‡ (E_a-MPA^‡ = 38.13 and T∆S_a-MPA^‡ = 0.53kJ mol^-1). While the [HSA-CdTe@MPA]^° dissociation required higher energy and lower entropy loss (E_d-MPA^‡ = 49.96 and T∆S_d-MPA^‡ = - 32.18kJ mol^-1) than the [HSA-CdTe@TGA]^° dissociation (E_d-TGA^‡= 30.78 and T∆S_d-TGA^‡= - 51.12 kJ mol^-1). The stability of [HSA-QDs]^° was independent of the temperature and functionalizing group. However, the enthalpic and entropic components were highly affected by the substitution of MPA (ΔH° = - 11.83 and TΔS° = 32.72 kJ mol^-1) with TGA (ΔH° = 34.31 and TΔS° = 79.73 kJ mol^-1). Furthermore, molecular docking results indicated that the metal site on the QDs contributes to the stabilization of [HSA-QDs]^°. Therefore, differences in QD functionalization and surface coverage densities can alter the HSA-QD interaction, thus their application.

Discovery of Selective Inhibitor Leads by Targeting an Allosteric Site in Insulin-Regulated Aminopeptidase

Insulin-Regulated aminopeptidase (IRAP) is a zinc-dependent aminopeptidase with several important biological functions and is an emerging pharmaceutical target for cognitive enhancement and immune system regulation. Aiming to discover lead-like IRAP inhibitors with enhanced selectivity versus homologous enzymes, we targeted an allosteric site at the C-terminal domain pocket of IRAP. We compiled a library of 2.5 million commercially available compounds from the ZINC database, and performed molecular docking at the target pocket of IRAP and the corresponding pocket of the homologous endoplasmic reticulum aminopeptidase 1 (ERAP1). Of the top compounds that showed high selectivity, 305 were further analyzed by molecular dynamics simulations and free energy calculations, leading to the selection of 33 compounds for in vitro evaluation. Two orthogonal functional assays were employed: one using a small fluorogenic substrate and one following the degradation of oxytocin, a natural peptidic substrate of IRAP. In vitro evaluation suggested that several of the compounds tested can inhibit IRAP, but the inhibition profile was dependent on substrate size, consistent with the allosteric nature of the targeted site. Overall, our results describe several novel leads as IRAP inhibitors and suggest that the C-terminal domain pocket of IRAP is a promising target for developing highly selective IRAP inhibitors.

Residue-Residue Interaction Prediction via Stacked Meta-Learning

Protein-protein interactions (PPIs) are the basis of most biological functions determined by residue-residue interactions (RRIs). Predicting residue pairs responsible for the interaction is crucial for understanding the cause of a disease and drug design. Computational approaches that considered inexpensive and faster solutions for RRI prediction have been widely used to predict protein interfaces for further analysis. This study presents RRI-Meta, an ensemble meta-learning-based method for RRI prediction. Its hierarchical learning structure comprises four base classifiers and one meta-classifier to integrate predictive strengths from different classifiers. It considers multiple feature types, including sequence-, structure-, and neighbor-based features, for characterizing other properties of a residue interaction environment to better distinguish between noninteracting and interacting residues. We conducted the same experiments using the same data as previously reported in the literature to demonstrate RRI-Meta's performance. Experimental results show that RRI-Meta is superior to several current prediction tools. Additionally, to analyze the factors that affect the performance of RRI-Meta, we conducted a comparative case study using different protein complexes.

Interplay among transacting factors around promoter in the initial phases of transcription

The initiation signals are raised around the promoter by one of the general transcription factors, triggering a sequence of events that lead to mRNA transcript formation from target genes. Both specific noncoding DNA regions and transacting, macromolecular assemblies are intricately involved and indispensable. The transition between the subsequent transcriptional stages is accompanied by stage-specific signals and structural changes in the macromolecular assemblies and facilitated by the recruitment/removal of other chromatin and transcription-associated elements. Here, we discuss the choreography of transacting factors around promoter in the establishment and effectuation of the initial phases of transcription such as NDR formation, +1 nucleosome positioning, promoter DNA opening, and RNAPII promoter escape from a structural viewpoint.

PU.1 Mimic Synthetic Peptides Selectively Bind with GATA-1 and Allow c-Jun PU.1 Binding to Enhance Myelopoiesis

BACKGROUND

Hematopoietic stem cells' commitment to myelopoiesis builds immunity to prevent infection. This process is controlled through transcription factor, especially Purine rich box 1 (PU.1) protein, which plays a central role in regulating myelopoiesis. The β3/β4 region of PU.1 accommodates a coactivator transcription factor, c-Jun, to activate myelopoiesis. However, an erythroid transcription factor, GATA-1, competes with c-Jun for the β3/β4 region, abolishing myelopoiesis and promoting erythropoiesis. This competitive regulation decides the hematopoietic stem cells' commitment towards either erythroid or myeloid lineage.

METHODS

Therefore, this study investigated the in vitro and in vivo effect of novel synthetic PU.1 β3/β4 mimic peptide analogs and peptide-loaded hydrophilic poly(D,L-lactide-co-glycolide) (PLGA) nanoparticles.

RESULTS

The designed peptides significantly increase the expression of corresponding myeloid markers, specifically CD33 in vitro. However, the in vivo delivery of peptide-loaded PLGA nanoparticles, which have sustained release effect of peptides, increases 10.8% of granulocytes as compared to control.

CONCLUSION

The observations showed that the fabricated nanoparticles protected the loaded peptides from the harsh intracellular environment for a longer duration without causing any toxicity. These findings highlight the possibility to use these peptides and peptide-loaded nanoparticles to increase hematopoietic stem cell commitment to myeloid cells in case of opportunistic infection.

Advances and perspectives in discovery and functional analysis of small secreted proteins in plants

Small secreted proteins (SSPs) are less than 250 amino acids in length and are actively transported out of cells through conventional protein secretion pathways or unconventional protein secretion pathways. In plants, SSPs have been found to play important roles in various processes, including plant growth and development, plant response to abiotic and biotic stresses, and beneficial plant-microbe interactions. Over the past 10 years, substantial progress has been made in the identification and functional characterization of SSPs in several plant species relevant to agriculture, bioenergy, and horticulture. Yet, there are potentially a lot of SSPs that have not been discovered in plant genomes, which is largely due to limitations of existing computational algorithms. Recent advances in genomics, transcriptomics, and proteomics research, as well as the development of new computational algorithms based on machine learning, provide unprecedented capabilities for genome-wide discovery of novel SSPs in plants. In this review, we summarize known SSPs and their functions in various plant species. Then we provide an update on the computational and experimental approaches that can be used to discover new SSPs. Finally, we discuss strategies for elucidating the biological functions of SSPs in plants.

New nimesulide derivatives with amide/sulfonamide moieties: Selective COX-2 inhibition and antitumor effects

Seventeen new amide/sulfonamide containing nimesulide derivatives were synthesized and characterized by several spectroscopic techniques and primarily investigated for their inhibitory potential on COX enzymes and other pro-inflammatory factors. Experimental analyses showed that among seventeen compounds, N8 and N10 have remarkable potency and selectivity for the COX-2 enzyme over COX-1 at very low doses as compared to nimesulide. Moreover, both N8 and N10 selectively reduced the Lipopolysaccharide (LPS)-stimulated COX-2 mRNA expression level while the COX-1 level remained stable. Both PGE₂ release and nitric oxide production in macrophage cells were significantly suppressed by the N8 and N10 treatment groups. In silico ADME/Tox, molecular docking and molecular dynamics (MD) simulations were also conducted. Additionally, all compounds were also screened in a panel of cancer cell lines for their antiproliferative properties by MTT and SRB assays. Compound N17 exhibited a considerable antiproliferative effect on the colon (IC₅₀: 9.24 μM) and breast (IC₅₀: 11.35 μM) cancer cell lines. N17 exposure for 48 h decreased expression of anti-apoptotic protein BCL-2 and increased the expression of apoptogenic BAX. Besides, the BAX/BCL-2 ratio was increased with visible ultrastructural changes and apoptotic bodies under scanning electron microscopy. In order to investigate the structural and dynamical properties of selected hits on the target structures, multiscale molecular modeling studies are also conducted. Our combined in silico and in vitro results suggest that N8 and N10 could be further developed as potential nonsteroidal anti-inflammatory drugs (NSAIDs), while cytotoxic N17 might be studied as a potential lead compound that could be developed as an anticancer agent.

Biomacromolecule-based photo-thermal agents for tumor treatment

Cancer treatment has become one of the biggest challenges in modern medicine. Recently, many efforts have been devoted to treat tumors by surgical resection, radiotherapy, or chemotherapy. In comparison to these methods, photo-thermal therapy (PTT) with noninvasive, controllable, direct, and precise characteristics has received tremendous attention in eliminating tumor cells over the past decades. In particular, PTT based on biomacromolecule-based photo-thermal agents (PTAs) outperforms other systems with high photo-thermal efficiency, simple coating, and low immunogenicity. Considering the unique advantages of biomacromolecule-based PTAs in tumor treatment, it is necessary to summarize the recent progress in the field of biomacromolecule-based PTAs for tumor treatment. Herein, this minireview outlines recent progress in the fabrication and applications of biomacromolecule-based PTAs. Within this framework, various types of biomacromolecule-based PTAs are highlighted, including cell-based agents, protein-based agents, nucleotide-based agents, and polysaccharide-based PTAs. In each section, the functional design, photo-thermal effects, and potential clinical applications of each type of PTA are discussed. Finally, a brief perspective for the development of biomacromolecule-based PTAs is presented.

Machine learning builds full-QM precision protein force fields in seconds

Full-quantum mechanics (QM) calculations are extraordinarily precise but difficult to apply to large systems, such as biomolecules. Motivated by the massive demand for efficient calculations for large systems at the full-QM level and by the significant advances in machine learning, we have designed a neural network-based two-body molecular fractionation with conjugate caps (NN-TMFCC) approach to accelerate the energy and atomic force calculations of proteins. The results show very high precision for the proposed NN potential energy surface models of residue-based fragments, with energy root-mean-squared errors (RMSEs) less than 1.0 kcal/mol and force RMSEs less than 1.3 kcal/mol/Å for both training and testing sets. The proposed NN-TMFCC method calculates the energies and atomic forces of 15 representative proteins with full-QM precision in 10-100 s, which is thousands of times faster than the full-QM calculations. The computational complexity of the NN-TMFCC method is independent of the protein size and only depends on the number of residue species, which makes this method particularly suitable for rapid prediction of large systems with tens of thousands or even hundreds of thousands of times acceleration. This highly precise and efficient NN-TMFCC approach exhibits considerable potential for performing energy and force calculations, structure predictions and molecular dynamics simulations of proteins with full-QM precision.

Visual pH Sensors: From a Chemical Perspective to New Bioengineered Materials

Many human activities and cellular functions depend upon precise pH values, and pH monitoring is considered a fundamental task. Colorimetric and fluorescence sensors for pH measurements are chemical and biochemical tools able to sense protons and produce a visible signal. These pH sensors are gaining widespread attention as non-destructive tools, visible to the human eye, that are capable of a real-time and in-situ response. Optical “visual” sensors are expanding researchers’ interests in many chemical contexts and are routinely used for biological, environmental, and medical applications. In this review we provide an overview of trending colorimetric, fluorescent, or dual-mode responsive visual pH sensors. These sensors include molecular synthetic organic sensors, metal organic frameworks (MOF), engineered sensing nanomaterials, and bioengineered sensors. We review different typological chemical entities of visual pH sensors, three-dimensional structures, and signaling mechanisms for pH sensing and applications; developed in the past five years. The progression of this review from simple organic molecules to biological macromolecules seeks to benefit beginners and scientists embarking on a project of pH sensing development, who needs background information and a quick update on advances in the field. Lessons learned from these tools will aid pH determination projects and provide new ways of thinking for cell bioimaging or other cutting-edge in vivo applications.

Open-access data: A cornerstone for artificial intelligence approaches to protein structure prediction

The Protein Data Bank (PDB) was established in 1971 to archive three-dimensional (3D) structures of biological macromolecules as a public good. Fifty years later, the PDB is providing millions of data consumers around the world with open access to more than 175,000 experimentally determined structures of proteins and nucleic acids (DNA, RNA) and their complexes with one another and small-molecule ligands. PDB data users are working, teaching, and learning in fundamental biology, biomedicine, bioengineering, biotechnology, and energy sciences. They also represent the fields of agriculture, chemistry, physics and materials science, mathematics, statistics, computer science, and zoology, and even the social sciences. The enormous wealth of 3D structure data stored in the PDB has underpinned significant advances in our understanding of protein architecture, culminating in recent breakthroughs in protein structure prediction accelerated by artificial intelligence approaches and deep or machine learning methods.

Discovery of Natural Products Targeting NQO1 via an Approach Combining Network-Based Inference and Identification of Privileged Substructures

NAD(P)H:quinone oxidoreductase 1 (NQO1) has been shown to be a potential therapeutic target for various human diseases, such as cancer and neurodegenerative disorders. Recent advances in computational methods, especially network-based methods, have made it possible to identify novel compounds for a target with high efficiency and low cost. In this study, we designed a workflow combining network-based methods and identification of privileged substructures to discover new compounds targeting NQO1 from a natural product library. According to the prediction results, we purchased 56 compounds for experimental validation. Without the assistance of privileged substructures, 31 compounds (31/56 = 55.4%) showed IC₅₀ < 100 μM, and 11 compounds (11/56 = 19.6%) showed IC₅₀ < 10 μM. With the assistance of privileged substructures, the two success rates were increased to 61.8 and 26.5%, respectively. Seven natural products further showed inhibitory activity on NQO1 at the cellular level, which may serve as lead compounds for further development. Moreover, network analysis revealed that osthole may exert anticancer effects against multiple cancer types by inhibiting not only carbonic anhydrases IX and XII but also NQO1. Our workflow and computational methods can be easily applied in other targets and become useful tools in drug discovery and development.

Structural biology in the time of COVID-19: perspectives on methods and milestones

The global COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has wreaked unprecedented havoc on global society, in terms of a huge loss of life and burden of morbidity, economic upheaval and social disruption. Yet the sheer magnitude and uniqueness of this event has also spawned a massive mobilization of effort in the scientific community to investigate the virus, to develop therapeutics and vaccines, and to understand the public health impacts. Structural biology has been at the center of these efforts, and so it is advantageous to take an opportunity to reflect on the status of structural science vis-à-vis its role in the fight against COVID-19, to register the unprecedented response and to contemplate the role of structural biology in addressing future outbreak threats. As the one-year anniversary of the World Health Organization declaration that COVID-19 is a pandemic has just passed, over 1000 structures of SARS-CoV-2 biomolecules have been deposited in the Worldwide Protein Data Bank (PDB). It is rare to obtain a snapshot of such intense effort in the structural biology arena and is of special interest as the 50th anniversary of the PDB is celebrated in 2021. It is additionally timely as it overlaps with a period that has been termed the 'resolution revolution' in cryoelectron microscopy (CryoEM). CryoEM has recently become capable of producing biomolecular structures at similar resolutions to those traditionally associated with macromolecular X-ray crystallo-graphy. Examining SARS-CoV-2 protein structures that have been deposited in the PDB since the virus was first identified allows a unique window into the power of structural biology and a snapshot of the advantages of the different techniques available, as well as insight into the complementarity of the structural methods.

Deep Learning in Virtual Screening: Recent Applications and Developments

Drug discovery is a cost and time-intensive process that is often assisted by computational methods, such as virtual screening, to speed up and guide the design of new compounds. For many years, machine learning methods have been successfully applied in the context of computer-aided drug discovery. Recently, thanks to the rise of novel technologies as well as the increasing amount of available chemical and bioactivity data, deep learning has gained a tremendous impact in rational active compound discovery. Herein, recent applications and developments of machine learning, with a focus on deep learning, in virtual screening for active compound design are reviewed. This includes introducing different compound and protein encodings, deep learning techniques as well as frequently used bioactivity and benchmark data sets for model training and testing. Finally, the present state-of-the-art, including the current challenges and emerging problems, are examined and discussed.

PyPLIF HIPPOS-Assisted Prediction of Molecular Determinants of Ligand Binding to Receptors

Identification of molecular determinants of receptor-ligand binding could significantly increase the quality of structure-based virtual screening protocols. In turn, drug design process, especially the fragment-based approaches, could benefit from the knowledge. Retrospective virtual screening campaigns by employing AutoDock Vina followed by protein-ligand interaction fingerprinting (PLIF) identification by using recently published PyPLIF HIPPOS were the main techniques used here. The ligands and decoys datasets from the enhanced version of the database of useful decoys (DUDE) targeting human G protein-coupled receptors (GPCRs) were employed in this research since the mutation data are available and could be used to retrospectively verify the prediction. The results show that the method presented in this article could pinpoint some retrospectively verified molecular determinants. The method is therefore suggested to be employed as a routine in drug design and discovery.

Experimental approaches for investigating ion atmospheres around nucleic acids and proteins

Ionic interactions are crucial to biological functions of DNA, RNA, and proteins. Experimental research on how ions behave around biological macromolecules has lagged behind corresponding theoretical and computational research. In the 21st century, quantitative experimental approaches for investigating ionic interactions of biomolecules have become available and greatly facilitated examinations of theoretical electrostatic models. These approaches utilize anomalous small-angle X-ray scattering, atomic emission spectroscopy, mass spectrometry, or nuclear magnetic resonance (NMR) spectroscopy. We provide an overview on the experimental methodologies that can quantify and characterize ions within the ion atmospheres around nucleic acids, proteins, and their complexes.

The fourth scientific discovery paradigm for precision medicine and healthcare: Challenges ahead

With the progression of modern information techniques, such as next generation sequencing (NGS), Internet of Everything (IoE) based smart sensors, and artificial intelligence algorithms, data-intensive research and applications are emerging as the fourth paradigm for scientific discovery. However, we face many challenges to practical application of this paradigm. In this article, 10 challenges to data-intensive discovery and applications in precision medicine and healthcare are summarized and the future perspectives on next generation medicine are discussed.

Computational evaluation of interactions between olfactory receptor OR2W1 and its ligands

Mammalian olfactory receptors are a family of G protein‒coupled receptors (GPCRs) that occupy a large part of the genome. In human genes, olfactory receptors account for more than 40% of all GPCRs. Several types of GPCR structures have been identified, but there is no single olfactory receptor whose structure has been determined experimentally to date. The aim of this study was to model the interactions between an olfactory receptor and its ligands at the molecular level to provide hints on the binding modes between the OR2W1 olfactory receptor and its agonists and inverse agonists. The results demonstrated the modes of ligand binding in a three-dimensional model of OR2W1 and showed a statistically significant difference in binding affinity to the olfactory receptor between agonists and inverse agonists.

A systematic approach to inserting split inteins for Boolean logic gate engineering and basal activity reduction

Split inteins are powerful tools for seamless ligation of synthetic split proteins. Yet, their use remains limited because the already intricate split site identification problem is often complicated by the requirement of extein junction sequences. To address this, we augment a mini-Mu transposon-based screening approach and devise the intein-assisted bisection mapping (IBM) method. IBM robustly reveals clusters of split sites on five proteins, converting them into AND or NAND logic gates. We further show that the use of inteins expands functional sequence space for splitting a protein. We also demonstrate the utility of our approach over rational inference of split sites from secondary structure alignment of homologous proteins, and that basal activities of highly active proteins can be mitigated by splitting them. Our work offers a generalizable and systematic route towards creating split protein-intein fusions for synthetic biology.

Network Pharmacology and Molecular Docking Suggest the Mechanism for Biological Activity of Rosmarinic Acid

Rosmarinic acid (RosA) is a natural phenolic acid compound, which is mainly extracted from Labiatae and Arnebia. At present, there is no systematic analysis of its mechanism. Therefore, we used the method of network pharmacology to analyze the mechanism of RosA. In our study, PubChem database was used to search for the chemical formula and the Chemical Abstracts Service (CAS) number of RosA. Then, the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) was used to evaluate the pharmacodynamics of RosA, and the Comparative Toxicogenomics Database (CTD) was used to identify the potential target genes of RosA. In addition, the Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of target genes were carried out by using the web-based gene set analysis toolkit (WebGestalt). At the same time, we uploaded the targets to the STRING database to obtain the protein interaction network. Then, we carried out a molecular docking about targets and RosA. Finally, we used Cytoscape to establish a visual protein-protein interaction network and drug-target-pathway network and analyze these networks. Our data showed that RosA has good biological activity and drug utilization. There are 55 target genes that have been identified. Then, the bioinformatics analysis and network analysis found that these target genes are closely related to inflammatory response, tumor occurrence and development, and other biological processes. These results demonstrated that RosA can act on a variety of proteins and pathways to form a systematic pharmacological network, which has good value in drug development and utilization.

Modeling Immunity with Rosetta: Methods for Antibody and Antigen Design

Structure-based antibody and antigen design has advanced greatly in recent years, due not only to the increasing availability of experimentally determined structures but also to improved computational methods for both prediction and design. Constant improvements in performance within the Rosetta software suite for biomolecular modeling have given rise to a greater breadth of structure prediction, including docking and design application cases for antibody and antigen modeling. Here, we present an overview of current protocols for antibody and antigen modeling using Rosetta and exemplify those by detailed tutorials originally developed for a Rosetta workshop at Vanderbilt University. These tutorials cover antibody structure prediction, docking, and design and antigen design strategies, including the addition of glycans in Rosetta. We expect that these materials will allow novice users to apply Rosetta in their own projects for modeling antibodies and antigens.

Integrating Multimeric Threading With High-throughput Experiments for Structural Interactome of Escherichia coli

Genome-wide protein-protein interaction (PPI) determination remains a significant unsolved problem in structural biology. The difficulty is twofold since high-throughput experiments (HTEs) have often a relatively high false-positive rate in assigning PPIs, and PPI quaternary structures are more difficult to solve than tertiary structures using traditional structural biology techniques. We proposed a uniform pipeline, Threpp, to address both problems. Starting from a pair of monomer sequences, Threpp first threads both sequences through a complex structure library, where the alignment score is combined with HTE data using a naïve Bayesian classifier model to predict the likelihood of two chains to interact with each other. Next, quaternary complex structures of the identified PPIs are constructed by reassembling monomeric alignments with dimeric threading frameworks through interface-specific structural alignments. The pipeline was applied to the Escherichia coli genome and created 35,125 confident PPIs which is 4.5-fold higher than HTE alone. Graphic analyses of the PPI networks show a scale-free cluster size distribution, consistent with previous studies, which was found critical to the robustness of genome evolution and the centrality of functionally important proteins that are essential to E. coli survival. Furthermore, complex structure models were constructed for all predicted E. coli PPIs based on the quaternary threading alignments, where 6771 of them were found to have a high confidence score that corresponds to the correct fold of the complexes with a TM-score >0.5, and 39 showed a close consistency with the later released experimental structures with an average TM-score = 0.73. These results demonstrated the significant usefulness of threading-based homologous modeling in both genome-wide PPI network detection and complex structural construction.

Molecular storytelling for online structural biology outreach and education

Knowledge about the structure and function of biomolecules continues to grow exponentially, enabling us to "see" structural snapshots of biomolecular interactions and functional assemblies. At PDB-101, the educational portal of the RCSB Protein Data Bank, we have taken a storytelling approach to make this body of knowledge accessible and comprehensible to a wide community of students, educators, and the general public. For over 20 years, the Molecule of the Month series has utilized a traditional illustrated storytelling approach that is regularly adapted for classroom instruction. Similar visual and interactive storytelling approaches are used to present topical subjects at PDB-101 and full curricular materials and case studies for building a detailed narrative around topics of particular interest. This emphasis on storytelling led to the Video Challenge for High School students, now in its 8th year. In this Article, we will present some of the lessons we have learned for teaching and communicating structural biology using the PDB archive of biomolecular structures.

Regulation of the conserved 3'-5' exoribonuclease EXOSC10/Rrp6 during cell division, development and cancer

The conserved 3'-5' exoribonuclease EXOSC10/Rrp6 processes and degrades RNA, regulates gene expression and participates in DNA double-strand break repair and control of telomere maintenance via degradation of the telomerase RNA component. EXOSC10/Rrp6 is part of the multimeric nuclear RNA exosome and interacts with numerous proteins. Previous clinical, genetic, biochemical and genomic studies revealed the protein's essential functions in cell division and differentiation, its RNA substrates and its relevance to autoimmune disorders and oncology. However, little is known about the regulatory mechanisms that control the transcription, translation and stability of EXOSC10/Rrp6 during cell growth, development and disease and how these mechanisms evolved from yeast to human. Herein, we provide an overview of the RNA- and protein expression profiles of EXOSC10/Rrp6 during cell division, development and nutritional stress, and we summarize interaction networks and post-translational modifications across species. Additionally, we discuss how known and predicted protein interactions and post-translational modifications influence the stability of EXOSC10/Rrp6. Finally, we explore the idea that different EXOSC10/Rrp6 alleles, which potentially alter cellular protein levels or affect protein function, might influence human development and disease progression. In this review we interpret information from the literature together with genomic data from knowledgebases to inspire future work on the regulation of this essential protein's stability in normal and malignant cells.

The 2021 Nucleic Acids Research database issue and the online molecular biology database collection

The 2021 Nucleic Acids Research database Issue contains 189 papers spanning a wide range of biological fields and investigation. It includes 89 papers reporting on new databases and 90 covering recent changes to resources previously published in the Issue. A further ten are updates on databases most recently published elsewhere. Seven new databases focus on COVID-19 and SARS-CoV-2 and many others offer resources for studying the virus. Major returning nucleic acid databases include NONCODE, Rfam and RNAcentral. Protein family and domain databases include COG, Pfam, SMART and Panther. Protein structures are covered by RCSB PDB and dispersed proteins by PED and MobiDB. In metabolism and signalling, STRING, KEGG and WikiPathways are featured, along with returning KLIFS and new DKK and KinaseMD, all focused on kinases. IMG/M and IMG/VR update in the microbial and viral genome resources section, while human and model organism genomics resources include Flybase, Ensembl and UCSC Genome Browser. Cancer studies are covered by updates from canSAR and PINA, as well as newcomers CNCdatabase and Oncovar for cancer drivers. Plant comparative genomics is catered for by updates from Gramene and GreenPhylDB. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). The NAR online Molecular Biology Database Collection has been substantially updated, revisiting nearly 1000 entries, adding 90 new resources and eliminating 86 obsolete databases, bringing the current total to 1641 databases. It is available at https://www.oxfordjournals.org/nar/database/c/.

Anhydroicaritin Inhibits EMT in Breast Cancer by Enhancing GPX1 Expression: A Research Based on Sequencing Technologies and Bioinformatics Analysis

Background: Breast cancer (BC) is the leading cause of cancer-related deaths among women worldwide. The application of advanced technology has promoted accurate diagnosis and treatment of cancer. Anhydroicaritin (AHI) is a flavonoid with therapeutic potential in BC treatment. The current study aimed to determine AHI's mechanism in BC treatment via RNA sequencing, comprehensive bioinformatics analysis, and experimental verification. Methods: Network pharmacology and MTT (3-(4,5)-dimethylthiazolyl-3,5- diphenyltetrazolium bromide) experiments were conducted to first confirm AHI's anti-BC effect. RNA sequencing was performed to identify the genes affected by AHI. Differential expression analysis, survival analysis, gene set enrichment analysis, and immune infiltration analysis were performed via bioinformatics analysis. Western blot analysis, reverse transcription-polymerase chain reaction (RT-PCR) experiment, molecular docking, and drug affinity responsive target stability (DARTS) experiments were also performed to confirm AHI's direct effect on glutathione peroxidase 1 (GPX1) expression. Confocal immunofluorescence analysis was conducted to verify AHI's effect on the occurrence and development of epithelial-mesenchymal transition (EMT). Finally, BC nude mouse xenografts were established, and AHI's molecular mechanism on BC was explored. Results: Network pharmacology results demonstrated that AHI's therapeutic targets on BC were related to the proliferation, invasion, and metastasis of BC cells. AHI significantly inhibited the proliferation of 4T1 and MDA-MB-231 BC cells in the MTT experiments. RNA sequencing results showed that AHI upregulated the GPX1 expression in the 4T1 and MDA-MB-231 BC cells. Next, bioinformatics analysis revealed that GPX1 is less expressed in BC than in normal breast tissues. Patients with high GPX1 expression levels tended to have prolonged overall survival and disease-free survival than patients with low GPX1 expression levels in BC. Western blot and RT-PCR experiments revealed that AHI increased the protein and mRNA levels of GPX1. Molecular docking and DARTS experiments confirmed the direct binding combination between AHI and GPX1. After the evaluation of the EMT scores of 1,078 patients with BC, we found a potential anti-BC role of GPX1 possibly via suppression of the malignant EMT. The confocal immunofluorescence analysis showed that AHI increased E-cadherin expression levels and reduced vimentin expression levels in BC cells. Animal experiments showed that AHI significantly inhibited tumor growth. AHI also inhibited EMT by enhancing GPX1 and caspase3 cleavage, hence inhibiting EMT markers (i.e., N-cadherin and vimentin) and Ki-67. Conclusion: GPX1 plays a critical role in BC, which may be a biomarker for the prognosis. In addition, AHI suppressed EMT by increasing GPX1 expression, which may serve as a potential therapy for BC treatment.

Impact of structural biologists and the Protein Data Bank on small-molecule drug discovery and development

The Protein Data Bank (PDB) is an international core data resource central to fundamental biology, biomedicine, bioenergy, and biotechnology/bioengineering. Now celebrating its 50th anniversary, the PDB houses >175,000 experimentally determined atomic structures of proteins, nucleic acids, and their complexes with one another and small molecules and drugs. The importance of three-dimensional (3D) biostructure information for research and education obtains from the intimate link between molecular form and function evident throughout biology. Among the most prolific consumers of PDB data are biomedical researchers, who rely on the open access resource as the authoritative source of well-validated, expertly curated biostructures. This review recounts how the PDB grew from just seven protein structures to contain more than 49,000 structures of human proteins that have proven critical for understanding their roles in human health and disease. It then describes how these structures are used in academe and industry to validate drug targets, assess target druggability, characterize how tool compounds and other small-molecules bind to drug targets, guide medicinal chemistry optimization of binding affinity and selectivity, and overcome challenges during preclinical drug development. Three case studies drawn from oncology exemplify how structural biologists and open access to PDB structures impacted recent regulatory approvals of antineoplastic drugs.

Evolution of the SARS-CoV-2 proteome in three dimensions (3D) during the first six months of the COVID-19 pandemic

Three-dimensional structures of SARS-CoV-2 and other coronaviral proteins archived in the Protein Data Bank were used to analyze viral proteome evolution during the first six months of the COVID-19 pandemic. Analyses of spatial locations, chemical properties, and structural and energetic impacts of the observed amino acid changes in >48,000 viral proteome sequences showed how each one of the 29 viral study proteins have undergone amino acid changes. Structural models computed for every unique sequence variant revealed that most substitutions map to protein surfaces and boundary layers with a minority affecting hydrophobic cores. Conservative changes were observed more frequently in cores versus boundary layers/surfaces. Active sites and protein-protein interfaces showed modest numbers of substitutions. Energetics calculations showed that the impact of substitutions on the thermodynamic stability of the proteome follows a universal bi-Gaussian distribution. Detailed results are presented for six drug discovery targets and four structural proteins comprising the virion, highlighting substitutions with the potential to impact protein structure, enzyme activity, and functional interfaces. Characterizing the evolution of the virus in three dimensions provides testable insights into viral protein function and should aid in structure-based drug discovery efforts as well as the prospective identification of amino acid substitutions with potential for drug resistance.