Insights into the Pharmacological Effects of Flavonoids: The Systematic Review of Computer Modeling

Computer modeling is a method that is widely used in scientific investigations to predict the biological activity, toxicity, pharmacokinetics, and synthesis strategy of compounds based on the structure of the molecule. This work is a systematic review of articles performed in accordance with the recommendations of PRISMA and contains information on computer modeling of the interaction of classical flavonoids with different biological targets. The review of used computational approaches is presented. Furthermore, the affinities of flavonoids to different targets that are associated with the infection, cardiovascular, and oncological diseases are discussed. Additionally, the methodology of bias risks in molecular docking research based on principles of evidentiary medicine was suggested and discussed. Based on this data, the most active groups of flavonoids and lead compounds for different targets were determined. It was concluded that flavonoids are a promising object for drug development and further research of pharmacology by in vitro, ex vivo, and in vivo models is required

Molecular Dynamics Study of Citrullinated Proteins Associated with the Development of Rheumatoid Arthritis

Biological activity regulation by protein post-translational modification (PTM) is critical for cell function, development, differentiation, and survival. Dysregulation of PTM proteins is present in various pathological conditions, including rheumatoid arthritis (RA). RA is a systemic autoimmune disease that primarily affects joints, and there are three main types of protein PTMs associated with the development of this disease, namely, glycosylation, citrullination, and carbamylation. Glycosylation is important for the processing and presentation of antigen fragments on the cell surface and can modulate immunoglobulin activity. The citrullination of autoantigens is closely associated with RA, as evidenced by the presence of antibodies specific to citrullinated proteins in the serum of patients. Carbamylation and dysregulation have recently been associated with RA development in humans.In this study, we performed an overview analysis of proteins with post-translational modifications associated with the development of RA adverted in peer-reviewed scientific papers for the past 20 years. As a result of the search, a list of target proteins and corresponding amino acid sequences with PTM in RA was formed. Structural characteristics of the listed modified proteins were extracted from the Protein Data Bank. Then, molecular dynamics experiments of intact protein structures and corresponding structures with PTMs were performed regarding structures in the list announced in the ProtDB service. This study aimed to conduct a molecular dynamics study of intact proteins and proteins, including post-translational modification and protein citrullination, likely associated with RA development. We observed another exhibition of the fundamental physics concept, symmetry, at the submolecular level, unveiled as the autonomous repetitions of outside the protein structural motif performance globule corresponding to those in the whole protein molecule

Prediction of Monomeric and Dimeric Structures of CYP102A1 Using AlphaFold2 and AlphaFold Multimer and Assessment of Point Mutation Effect on the Efficiency of Intra- and Interprotein Electron Transfer

The three-dimensional structure of monomers and homodimers of CYP102A1/WT (wild-type) proteins and their A83F and A83I mutant forms was predicted using the AlphaFold2 (AF2) and AlphaFold Multimer (AFMultimer) programs, which were compared with the rate constants of hydroxylation reactions of these enzyme forms to determine the efficiency of intra- and interprotein electron transport in the CYP102A1 hydroxylase system. The electron transfer rate constants (ket), which determine the rate of indole hydroxylation by the CYP102A1 system, were calculated based on the distances (R) between donor-acceptor prosthetic groups (PG) FAD→FMN→HEME of these proteins using factor β, which describes an exponential decay from R the speed of electron transport (ET) according to the tunnelling mechanism. It was shown that the structure of monomers in the homodimer, calculated using the AlpfaFold Multimer program, is in good agreement with the experimental structures of globular domains (HEME-, FMN-, and FAD-domains) in CYP102A1/WT obtained by X-ray structural analysis, and the structure of isolated monomers predicted in AF2 does not coincide with the structure of monomers in the homodimer, although a high level of similarity in individual domains remains. The structures of monomers and homodimers of A83F and A83I mutants were also calculated, and their structures were compared with the wild-type protein. Significant differences in the structure of all isolated monomers with respect to the structures of monomers in homodimers were also found for them, and at the same time, insignificant differences were revealed for all homodimers. Comparative analysis for CYP102A1/WT between the calculated intra- and interprotein distances FAD→FMN→HEME and the rate constants of hydroxylation in these proteins showed that the distance between prosthetic groups both in the monomer and in the dimer allows the implementation of electron transfer between PGs, which is consistent with experimental literature data about kcat. For the mutant form of monomer A83I, an increase in the distance between PGs was obtained, which can restrict electron transportation compared to WT; however, for the dimer of this protein, a decrease in the distance between PGs was observed compared to the WT form, which can lead to an increase in the electron transfer rate constant and, accordingly, kcat. For the monomer and homodimer of the A83F mutant, the calculations showed an increase in the distance between the PGs compared to the WT form, which should have led to a decrease in the electron transfer rate, but at the same time, for the homodimer, the approach of the aromatic group F262 with heme can speed up transportation for this form and, accordingly, the rate of hydroxylation

Molecular Modeling Insights into Upadacitinib Selectivity upon Binding to JAK Protein Family

Rheumatoid arthritis (RA) is a chronic disease characterized by bone joint damage and incapacitation. The mechanism underlying RA pathogenesis is autoimmunity in the connective tissue. Cytokines play an important role in the human immune system for signal transduction and in the development of inflammatory responses. Janus kinases (JAK) participate in the JAK/STAT pathway, which mediates cytokine effects, in particular interleukin 6 and IFNγ. The discovery of small molecule inhibitors of the JAK protein family has led to a revolution in RA therapy. The novel JAK inhibitor upadacitinib (RinvoqTM) has a higher selectivity for JAK1 compared to JAK2 and JAK3 in vivo. Currently, details on the molecular recognition of JAK1 by upadacitinib are not available. We found that characteristics of hydrogen bond formation with the glycine loop and hinge in JAKs define the selectivity. Our molecular modeling study could provide insight into the drug action mechanism and pharmacophore model differences in JAK isoforms

Changes in Protein Structural Motifs upon Post-Translational Modification in Kidney Cancer

Post-translational modification (PTM) leads to conformational changes in protein structure, modulates the biological function of proteins, and, consequently, changes the signature of metabolic transformations and the immune response in the body. Common PTMs are reversible and serve as a mechanism for modulating metabolic trans-formations in cells. It is likely that dysregulation of post-translational cellular signaling leads to abnormal proliferation and oncogenesis. We examined protein PTMs in the blood samples from patients with kidney cancer. Conformational changes in proteins after modification were analyzed. The proteins were analyzed using ultra-high resolution HPLC-MS/MS and structural analysis was performed with the AMBER and GROMACS software packages. Fifteen proteins containing PTMs were identified in blood samples from patients with kidney cancer. For proteins with PDB structures, a comparative analysis of the structural changes accompanying the modifications was performed. Results revealed that PTMs are localized in stable and compact space protein globule motifs that are exposed to a solvent. The phenomenon of modification is accompanied, as a rule, by an increase in the area available for the solvent of the modified amino acid residue and its active environment

Insights into the Cardiotoxic Effects of <i>Veratrum Lobelianum</i> Alkaloids: Pilot Study

Jervine, protoveratrine A (proA), and protoveratrine B (proB) are Veratrum alkaloids that are presented in some remedies obtained from Veratrum lobelianum, such as Veratrum aqua. This paper reports on a single-center pilot cardiotoxic mechanism study of jervine, proA, and proB in case series. The molecular aspects were studied via molecular dynamic simulation, molecular docking with cardiac sodium channel NaV1.5, and machine learning-based structure–activity relationship modeling. HPLC-MS/MS method in combination with clinical events were used to analyze Veratrum alkaloid cardiotoxicity in patients. Jervine demonstrates the highest docking score (−10.8 kcal/mol), logP value (4.188), and pKa value (9.64) compared with proA and proB. Also, this compound is characterized by the lowest calculated IC50. In general, all three analyzed alkaloids show the affinity to NaV1.5 that highly likely results in cardiotoxic action. The clinical data of seven cases of intoxication by Veratrum aqua confirms the results of molecular modeling. Patients exhibited nausea, muscle weakness, bradycardia, and arterial hypotension. The association between alkaloid concentrations in blood and urine and severity of patient condition is described. These experiments, while primary, confirmed that jervine, proA, and proB contribute to cardiotoxicity by NaV1.5 inhibition

Modification of Taxifolin Properties by Spray Drying

Taxifolin is known as an active pharmaceutical ingredient (API) and food supplement due to its high antioxidant activity, multiple pharmacological effects, and good safety profile. Previously, taxifolin spheres (TS) were obtained from industrially produced API taxifolin in Russia (RT). In our work, we perform a pharmaceutical analysis of this new taxifolin material versus RT. TS is an amorphous material; however, it is stable without the polymer carrier, as confirmed by Fourier transform infrared spectroscopy. Both RT and TS demonstrate high safety profiles and are assigned to Class 1 of the Biopharmaceutical Classification System based on the results of experiments with MDCK cells. The water solubility of the new taxifolin form was 2.225 times higher compared with RT. Hausner ratios for RT and TS were 1.421 and 1.219, respectively, while Carr indices were 29.63% and 19.00%, respectively. Additionally, TS demonstrated sustained release from tablets compared with RT: the half-life values of tablets were 14.56 min and 20.63 min for RT and TS, respectively. Thus, TS may be a promising object for developing oral antiseptics in the form of orally dispersed tablets with sustained release patterns because of its anti-inflammatory, -protozoal, and -viral activities

Modification of Taxifolin Properties by Spray Drying

Taxifolin is known as an active pharmaceutical ingredient (API) and food supplement due to its high antioxidant activity, multiple pharmacological effects, and good safety profile. Previously, taxifolin spheres (TS) were obtained from industrially produced API taxifolin in Russia (RT). In our work, we perform a pharmaceutical analysis of this new taxifolin material versus RT. TS is an amorphous material; however, it is stable without the polymer carrier, as confirmed by Fourier transform infrared spectroscopy. Both RT and TS demonstrate high safety profiles and are assigned to Class 1 of the Biopharmaceutical Classification System based on the results of experiments with MDCK cells. The water solubility of the new taxifolin form was 2.225 times higher compared with RT. Hausner ratios for RT and TS were 1.421 and 1.219, respectively, while Carr indices were 29.63% and 19.00%, respectively. Additionally, TS demonstrated sustained release from tablets compared with RT: the half-life values of tablets were 14.56 min and 20.63 min for RT and TS, respectively. Thus, TS may be a promising object for developing oral antiseptics in the form of orally dispersed tablets with sustained release patterns because of its anti-inflammatory, -protozoal, and -viral activities