An overlooked dispersal route of Cardueae (Asteraceae) from the Mediterranean to East Asia revealed by phylogenomic and biogeographical analyses of Atractylodes.

BACKGROUND AND AIMS: The East Asian-Tethyan disjunction pattern and its mechanisms of formation have long been of interest to researchers. Here, we studied the biogeographical history of Asteraceae tribe Cardueae, with a particular focus on the temperate East Asian genus Atractylodes DC., to understand the role of tectonic and climatic events in driving the diversification and disjunctions of the genus. METHODS: A total of 76 samples of Atractylodes from 36 locations were collected for RAD-sequencing. Three single nucleotide polymorphism (SNP) datasets based on different filtering strategies were used for phylogenetic analyses. Molecular dating and ancestral distribution reconstruction were performed using both chloroplast DNA sequences (127 Cardueae samples) and SNP (36 Atractylodes samples) datasets. KEY RESULTS: Six species of Atractylodes were well resolved as individually monophyletic, although some introgression was identified among accessions of A. chinensis, A. lancea and A. koreana. Dispersal of the subtribe Carlininae from the Mediterranean to East Asia occurred after divergence between Atractylodes and Carlina L. + Atractylis L. + Thevenotia DC. at ~31.57 Ma, resulting in an East Asian-Tethyan disjunction. Diversification of Atractylodes in East Asia mainly occurred from the Late Miocene to the Early Pleistocene. CONCLUSIONS: Aridification of Asia and the closure of the Turgai Strait in the Late Oligocene promoted the dispersal of Cardueae from the Mediterranean to East China. Subsequent uplift of the Qinghai-Tibet Plateau as well as changes in Asian monsoon systems resulted in an East Asian-Tethyan disjunction between Atractylodes and Carlina + Atractylis + Thevenotia. In addition, Late Miocene to Quaternary climates and sea level fluctuations played major roles in the diversification of Atractylodes. Through this study of different taxonomic levels using genomic data, we have revealed an overlooked dispersal route between the Mediterranean and far East Asia (Japan/Korea) via Central Asia and East China.

The complete chloroplast genome of Papaver setigerum and comparative analyses in Papaveraceae.

Papaver setigerum is an annual herb that is closely related to the opium poppy, P. somniferum. Genetic resources for P. setigerum are scarce. In the present study, we assembled the complete chloroplast (cp) genome of P. setigerum based on genome skimming data, and we conducted comparative cp genome analyses to study the evolutionary pattern in Papaveraceae. The cp genome of P. setigerum is 152,862 bp in length with a typical quadripartite structure. Comparative analyses revealed no gene rearrangement in the Papaveraceae family, although differences were evident in genome size, gene losses, as well as inverted repeats (IR) region expansion and contraction. The rps15 gene has been lost from the genomes of Meconopsis racemosa, Coreanomecon hylomeconoides, P. orientale, P. somniferum, and P. setigerum, and the ycf15 gene is found only in C. hylomeconoides. Moreover, 13 cpDNA markers, including psbA-trnH, rps16-trnQ, trnS-trnG, trnC-petN, trnE-trnT, trnL-trnF, trnF-ndhJ, petA-psbJ, ndhF-rpl32, rpl32-trnL, ccsA-ndhD, ndhE-ndhG, and rps15-ycf1, were identified with relatively high levels of variation within Papaver, which will be useful for species identification in this genus. Among those markers, psbA-trnH is the best one to distinguish P. somniferum and P. setigerum.

Chloroplast genome analyses and genomic resource development for epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae), using genome skimming data.

BACKGROUND: Epilithic sister genera Oresitrophe and Mukdenia (Saxifragaceae) have an epilithic habitat (rocky slopes) and a parapatric distribution in East Asia, which makes them an ideal model for a more comprehensive understanding of the demographic and divergence history and the influence of climate changes in East Asia. However, the genetic background and resources for these two genera are scarce. RESULTS: The complete chloroplast (cp) genomes of two Oresitrophe rupifraga and one Mukdenia rossii individuals were reconstructed and comparative analyses were conducted to examine the evolutionary pattern of chloroplast genomes in Saxifragaceae. The cp genomes ranged from 156,738 bp to 156,960 bp in length and had a typical quadripartite structure with a conserved genome arrangement. Comparative analysis revealed the intron of rpl2 has been lost in Heuchera parviflora, Tiarella polyphylla, M. rossii and O. rupifraga but presents in the reference genome of Penthorum chinense. Seven cp hotspot regions (trnH-psbA, trnR-atpA, atpI-rps2, rps2-rpoC2, petN-psbM, rps4-trnT and rpl33-rps18) were identified between Oresitrophe and Mukdenia, while four hotspots (trnQ-psbK, trnR-atpA, trnS-psbZ and rpl33-rps18) were identified within Oresitrophe. In addition, 24 polymorphic cpSSR loci were found between Oresitrophe and Mukdenia. Most importantly, we successfully developed 126 intergeneric polymorphic gSSR markers between Oresitrophe and Mukdenia, as well as 452 intrageneric ones within Oresitrophe. Twelve randomly selected intergeneric gSSRs have shown that these two genera exhibit a significant genetic structure. CONCLUSIONS: In this study, we conducted genome skimming for Oresitrophe rupifraga and Mukdenia rossii. Using these data, we were able to not only assemble their complete chloroplast genomes, but also develop abundant genetic resources (cp hotspots, cpSSRs, polymorphic gSSRs). The genomic patterns and genetic resources presented here will contribute to further studies on population genetics, phylogeny and conservation biology in Saxifragaceae.

Spiroketones and a Biphenyl Analog from Stems and Leaves of Larrea nitida and Their Inhibitory Activity against IL-6 Production.

Bioactivity-guided fractionation for the stems of leaves of Larrea nitida Cav., using interleukin-6 (IL-6) inhibitory assay in human mast cells (HMC-1), led to the isolation of three new compounds with an unprecedented skeleton in nature (1-3) and three known compounds (4-6). Their structures were elucidated through extensive spectroscopic analysis. The three new compounds were elucidated as two new spiroketones, nitidaones A (1), and B (2) and one new biphenyl analog, nitidaol (3). The known compounds were identified as nordihydroguaiaretic acid (4), 7,3',4'-tri-O-methylquercetin (5) and ayanin (6). All the isolates were tested for their inhibitory activity against IL-6 production in HMC-1 cells. Of them, compounds 1, 3-6 showed potent anti-inflammatory activity, with IC50 values of 12.8, 17.5, 14.9, 22.9, and 17.8 µM, respectively.

Desmodianone H and uncinanone B, potential tyrosinase inhibitors obtained from Lespedeza maximowiczii by using bioactivity-guided isolation.

A new bioactive compound, namely desmodianone H(1), and another known compound uncinanone B(2) were first isolated using bioactivity-guided isolation from the leaves of Lespedeza maximowiczii and structures were elucidated by comprehensive analysis of their nuclear magnetic resonance and mass spectrometry data. Compounds 1 and 2 exhibited strong inhibitory effects on mushroom tyrosinase activity.

Comparative phylogeography of the Smilax hispida group (Smilacaceae) in eastern Asia and North America--implications for allopatric speciation, causes of diversity disparity, and origins of temperate elements in Mexico.

The Smilax hispida group (Smilacaceae) exhibits a discontinuous distribution in eastern Asia, eastern and western United States, and Mexico. A broad scale phylogeographic analysis was conducted for this group to evaluate the hypotheses of accelerated allopatric divergence in eastern Asia and a northern origin of the temperate elements in Mexico. Phylogeny was inferred using seven plastid and nuclear DNA sequences. Species delineation was assessed using genealogical sorting indices (GSI). Lineage divergence time, haplotype diversification rates, and ancestral distributions were estimated using Bayesian methods. Phylogeographic patterns in eastern Asia and North America were compared by analyzing 539 individuals from 64 populations to assess allopatric diversification. Results strongly supported delineation of six allopatric species, the origin of this group from a Mexican ancestor around 11.42mya, and Mexican origins of the temperate species in Mexico. Significant geographic structure of haplotypes was found in eastern Asia, and greater haplotype diversification rate was observed for the North American lineage. Our data support allopatric speciation in eastern Asia but do not find evidence of an elevated diversification rate. Greater species diversity of the study system in eastern Asia may be due to a longer evolutionary history. Our results do not support northern origins of the Mexican temperate species.

Molecular docking analysis of HSV-1 proteins models with synthetic and plant derived compounds.

The atomic resolution model of US9, UL20, and gH protein of HSV is known. Hence, the ligand protein interaction of the US9, UL20, and gH protein models were carried out with synthetic drugs like acyclovir, bexarotene, vinorelbine, foscarnet, famciclovir, cidofovir and two plant derived natural drug acacetin and anthraquinone. Based on structure and docking study, it is predicted that protein US20 and gH binds with particular anti-HSV drug i.e. acyclovir, cidofovir, acacetin and famciclovir, acacetin respectively, while interaction of different protein is different with drugs.

Genome size estimation of false daisy, cheek weed, pot marigold and marigold.

We report the genome size estimated using flow cytometry for four closely related species, including false daisy (Eclipta prostrate), cheek weed (Ageratum conyzoides), pot marigold (Calendula officinalis), and marigold (Tagetes erecta) belonging to Asteraceae family. The detected genome size for false daisy, cheek weed, pot marigold, and marigold was, 2.435, 3.266, 3.413, and 1.897, Gbp, respectively, while their respective 2C DNA content was 2.5, 3.3, 3.5, and 1.9, pg. The information on genome size presented here will be useful for understanding genomic evolution and will also clear the way for additional genomic research in these species.

Neurological potency of native plants from sub-Himalayan West Bengal through reverse pharmacology.

Neurodegenerative diseases, such as Alzheimer's disease (AD), Parkinson's disease (PD), and epilepsy, pose a growing global health challenge due to an aging population. These conditions share common processes, including protein accumulation, oxidative stress, and neuro-inflammation, making their treatment complex and costly. Network pharmacology, an innovative approach integrating systems biology and computational biology, offers insights into multi-target formulations and the repurposing of existing medications for neurodegenerative diseases. We shortlisted 730 bioactive compounds from 25 traditional Himalayan plants, assessed their drug-like properties using ADME criteria, and predicted their potential target proteins through reverse docking and pharmacophore mapping. Our study identified 287 compounds with high gastrointestinal absorption and good blood-brain barrier permeability. These compounds were subjected to target prediction, yielding a list of 171 potential target proteins. Functional annotation and pathway enrichment analysis highlighted their involvement in steroid hormone-related pathways, MAPK signaling, FOXO signaling, TNF signaling, VEGF signaling, and neurotrophin signaling. Importantly, one plant, Valeriana jatamansi, exhibited an association with beta-amyloid binding activity, a potential therapeutic approach for AD. From our study we could understand how these plants modulate our body to manage these diseases. However, further in vitro and in vivo validation is needed before commercial and public use of this data.

Masitinib analogues with the N-methylpiperazine group replaced - A new hope for the development of anti-COVID-19 drugs.

Masitinib is an orally acceptable tyrosine kinase inhibitor that is currently investigated under clinical trials against cancer, asthma, Alzheimer's disease, multiple sclerosis and amyotrophic lateral sclerosis. A recent study confirmed the anti-severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) activity of masitinib through inhibition of the main protease (Mpro) enzyme, an important pharmacological drug target to block the replication of the coronavirus. However, due to the adverse effects and lower potency of the drug, there are opportunities to design better analogues of masitinib. Herein, we substituted the N-methylpiperazine group of Masitinib with different chemical moieties and evaluated their drug-likeness and toxicities. The filtered analogues were subjected to molecular docking studies which revealed that the analogues with substituents methylamine in M10 (CID10409602), morpholine in M23 (CID59789397) and 4-methylmorpholine in M32 (CID143003625) have a stronger affinity to the drug receptor compared to masitinib. The molecular dynamics (MD) simulation analysis reveals that the identified analogues alter the mobility, structural compactness, accessibility to solvent molecules, and the number of hydrogen bonds in the native target enzyme. These structural alterations can help explain the inhibitory mechanisms of these analogues against the target enzyme. Thus, our studies provide avenues for the design of new masitinib analogues as the SARS-CoV-2 Mpro inhibitors.

Detection of the peptidyl epitope for vaccine development against MPV.

Background: Monkeypox Virus (MPV) is the cause of zoonotic disease characterized by skin-eruption with pus cell formation and lymphadenopathy. This virus belongs to the Orthopoxvirus genus with DNA as its genetic material. Previously, this infection was reported from Africa and occasionally from USA and UK. However, recently there is a sudden surge of infection in non-epidemic countries and a new strain of MPVhas been discovered. Therefore it is important to revisit the phylogeny of MPV with the addition of new strains. Recently WHO also stressed the need of developing vaccines for new strains. In this scenario we have two objectives for this study -first, to reveal the exact phylogenetic position of the 2022 strain and second, to identify specific peptides which may be used for vaccine development in the future. Methods: The phylogenetic analysis was done with the help of Bayesian phylogeny. The dN/dS calculation was performed based on DNA polymerase genes of selected MPV strains. The peptidyl-epitope was searched in MPV2022/2 SLO strain with the help of several algorithms implemented in Allergen FP v.1.0, NetMHCII 2.3, MHCpred and Toxin Pred. The structure prediction of the proteins and peptides was performed through Hpepdock. The quality of the structures was validated through the Ramachandran plot. The molecular dynamics and simulation were performed through Gromacs software. The interaction between peptide and protein was assessed through Ligplot software. Results: The phylogenetic analysis revealed that the considered 2022 MPVstrains were close to the USA strains. The evolutionary analysis showed the volatile nature of the genome. Moreover, 9-mer peptide sequence was identified as an epitope for vaccine development. Conclusions: The emergence of more virulent strains in near future may not be ruled out. Immunocompromised patients are more susceptible to this virus hence sub-unit vaccine is a better choice than a recombinant or attenuated vaccine against monkeypox. We have identified a small stretch of specific peptide which may be used for developing a subunit vaccine against this virus.

Molecular docking analysis of Omt-A protein model from Aspergillus flavus with synthetic compounds.

Aflatoxin is a potent mycotoxin of Aspergillus flavus that has been classified as a Group I carcinogen. O-methyltransferase A (Omt-A) is a critical enzyme in the formation of aflatoxin. It catalyzes the methylation of norsalic acid to form the highly toxic intermediate averantin. The ligand-protein interaction of Omt-A was performed with piperlonguminin and blasticidins. The maximum affinity of -10.6 was found for the 5ICC_A piperlonguminine at site1 (X,Y,Z: -15.282, 21.785, 5.672). Compounds such as Blasticidin S, Neoeriocitrin, Blasticidin S - hydrochloric acid, 6,6''-Bigenkwanin, Pipernomaline, and Eriodictyol were found to have binding features to protein residues, as shown by computational interaction at the molecular level.

Anticancer potential of phytochemicals from Oroxylum indicum targeting Lactate Dehydrogenase A through bioinformatic approach.

In recent years, small molecule inhibition of LDHA (Lactate Dehydrogenase A) has evolved as an appealing option for anticancer therapy. LDHA catalyzes the interconversion of pyruvate and lactate in the glycolysis pathway to play a crucial role in aerobic glycolysis. Therefore, in the current investigation LDHA was targeted with bioactive phytochemicals of an ethnomedicinally important plant species Oroxylum indicum (L.) Kurz. A total of 52 phytochemicals were screened against LDHA protein through molecular docking, ADMET (Absorption, Distribution, Metabolism, Excretion and Toxicity) assay and molecular dynamics simulation to reveal three potential lead compounds such as Chrysin-7-O-glucuronide (-8.2 kcal/mol), Oroxindin (-8.1 kcal/mol) and Oroxin A (-8.0 kcal/mol). ADMET assay unveiled favorable pharmacokinetic, pharmacodynamic and toxicity properties for all the lead compounds. Molecular dynamics simulation exhibited significant conformational stability and compactness. MM/GBSA free binding energy calculations further corroborated the selection of top candidates where Oroxindin (-46.47 kcal/mol) was found to be better than Chrysin-7-O-glucuronide (-45.72 kcal/mol) and Oroxin A (-37.25 kcal/mol). Aldolase reductase and Xanthine dehydrogenase enzymes were found as potential drug targets and Esculin, the FDA approved drug was identified as structurally analogous to Oroxindin. These results could drive in establishing novel medications targeting LDHA to fight cancer.

LC-MS-based chemotaxonomic classification of wild-type Lespedeza sp. and its correlation with genotype.

In this study, 39 specimens belonging to Lespedeza species (Lespedeza cyrtobotrya, L. bicolor, L. maximowiczii, and Lespedeza cuneata) (Leguminosae) were classified phenotypically and genotypically. We constructed a phylogenetic tree based on the combined nrDNA (internal transcribed spacer; ITS) and cpDNA (trnL-trnF) sequences with the aim of classifying the genotypes. Samples were mainly divided into three genotypes. Samples of L. cyrtobotrya and L. bicolor were mixed in a single branch, whereas samples of L. maximowiczii and L. cuneata were clustered within species, respectively. We performed a liquid chromatography-electrospray ionization-mass spectrometry-based metabolite profiling analysis to classify the phenotypes. Multivariate statistical analyses such as principal component analysis (PCA) and hierarchical clustering analysis (HCA) were used for the clustering pattern analysis and distance analysis between species, respectively. According to the PCA and HCA results, leaves were classified into four phenotypes according to species. In both the genetic and chemotaxonomic classification methods, the distance between L. cyrtobotrya and L. bicolor was the closest between species, and L. cuneata was the farthest away from the other three species. Additionally, orthogonal partial least squares-discriminant analysis was employed to identify significantly different phytochemicals between species. We classified L. cyrtobotrya and L. bicolor by identifying significantly different phytochemicals. Interestingly, leaves and stems showed different phenotypic classifications based on the chemotaxonomic classification. Stem samples of the other three species were mixed regardless of species, whereas L. cyrtobotrya stem samples were clustered within species. The phenotypic classification of leaves coincided more with the genotypic classification than that of stems. Key message We classified four wild-type Lespedeza sp. by analyzing the combined nrDNA (ITS) and cpDNA (trnL-trnF) sequences. We also classified leaves and stems of Lespedeza sp. by applying liquid chromatography-mass spectroscopy-based metabolite profiling.

Macroevolutionary patterns of defense and pollination in Dalechampia vines: adaptation, exaptation, and evolutionary novelty.

We conducted phylogenetically informed comparative analyses of 81 taxa of Dalechampia (Euphorbiaceae) vines and shrubs to assess the roles of historical contingency and trait interaction in the evolution of plant-defense and pollinator-attraction systems. We asked whether defenses can originate by exaptation from preexisting pollinator attractants, or vice versa, whether plant defenses show escalation, and if so, whether by enhancing one line of defense or by adding new lines of defense. Two major patterns emerged: (i) correlated evolution of several complementary lines of defense of flowers, seeds, and leaves, and (ii) 5 to 6 losses of the resin reward, followed by redeployment of resin for defense of male flowers in 3 to 4 lineages, apparently in response to herbivore-mediated selection for defense of staminate flowers upon relaxation of pollinator-mediated selection on resin. In all cases, redeployment of resin involved reversion to the inferred ancestral arrangement of flowers and resiniferous bractlets. Triterpene resin has also been deployed for defense of leaves and developing seeds. Other unique defenses against florivores include nocturnal closure of large involucral bracts around receptive flowers and permanent closure around developing fruits (until opening again upon dehiscence). Escalation in one major clade occurred through an early dramatic increase in the number of lines of defense and in the other major clade by more limited increases throughout the group's evolution. We conclude that preaptations played important roles in the evolution of unique defense and attraction systems, and that the evolution of interactions with herbivores can be influenced by adaptations for pollination, and vice versa.

In silico analyses of major active constituents of fingerroot (Boesenbergia rotunda) unveils inhibitory activities against SARS-CoV-2 main protease enzyme.

Boesenbergia rotunda (L.) Mansf., commonly known as fingerroot is a perennial herb in the Zingiberaceae family with anticancer, anti-leptospiral, anti-inflammatory, antioxidant, antiulcer, and anti-herpes viral activities. While the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) inhibitory activity of B. rotunda extract has been recently found, the active compounds contributing to this activity are yet unknown. The main protease (Mpro) enzyme is one of the most well established therapeutic targets among coronaviruses which plays a vital role in the maturation and cleavage of polyproteins during viral replication. The current work aims to identify active phytochemical substances from B. rotunda extract that can inhibit the replication of SARS-CoV-2 by using a combined molecular docking and dynamic simulation approaches. The virtual screening experiment revealed that fifteen molecules out of twenty-three major active compounds in the plant extract have acceptable drug-like characteristics. Alpinetin, Pinocembrin, and Pinostrobin have binding energies of -7.51 kcal/mol, -7.21 kcal/mol, and -7.18 kcal/mol, respectively, and can suppress Mpro activity. The stability of the simulated complexes of the lead compounds with the drug-receptor is demonstrated by 100-ns MD simulations. The binding free energies study utilizing molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and molecular mechanics generalized Born surface area (MM-GBSA) show that the compounds and Mpro enzyme have favourable thermodynamic interactions, which are majorly driven by van der Waals forces. Thus, the selected bioactive phytochemicals from B. rotunda might be used as anti-SARS-CoV-2 candidates that target the Mpro enzyme.

An in silico approach unveils the potential of antiviral compounds in preclinical and clinical trials as SARS-CoV-2 omicron inhibitors.

The increased transmissibility and highly infectious nature of the new variant of concern (VOC) that is severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron and lack of effective therapy need the rapid discovery of therapeutic antivirals against it. The present investigation aimed to identify antiviral compounds that would be effective against SARS-CoV-2 Omicron. In this study, molecular docking experiments were carried out using the recently reported experimental structure of omicron spike protein in complex with human angiotensin-converting enzyme 2 (ACE2) and various antivirals in preclinical and clinical trial studies. Out of 36 tested compounds, Abemaciclib, Dasatinib and Spiperone are the three top-ranked molecules which scored binding energies of -10.08 kcal/mol, -10.06 kcal/mol and -9.54 kcal/mol respectively. Phe338, Asp339, and Asp364 are crucial omicron receptor residues involved in hydrogen bond interactions, while other residues were mostly involved in hydrophobic interactions with the lead molecules. The identified lead compounds also scored well in terms of drug-likeness. Molecular dynamics (MD) simulation, essential dynamics (ED) and entropic analysis indicate the ability of these molecules to modulate the activity of omicron spike protein. Therefore, Abemaciclib, Dasatinib and Spiperone are likely to be viable drug-candidate molecules that can block the interaction between the omicron spike protein and the host cellular receptor ACE2. Though our findings are compelling, more research into these molecules is needed before they can be employed as drugs to treat SARS-CoV-2 omicron infections.

Identification of Novel Drug Candidates for the Inhibition of Catalytic Cleavage Activity of Coronavirus 3CL-like Protease Enzyme.

BACKGROUND: There has been tremendous pressure on healthcare facilities globally due to the recent emergence of novel coronavirus infection known as COVID-19 and its rapid spread across the continents. The lack of effective therapeutics for the management of the pandemic calls for the discovery of new drugs and vaccines. OBJECTIVE: In the present study, a chemical library was screened for molecules against three coronavirus 3CL-like protease enzymes (SARS-CoV-2 3CLpro, SARS-CoV 3CLpro and MERS-CoV 3CLpro), which are a key player in the viral replication cycle. METHODS: Extensive computational methods such as virtual screening and molecular docking were employed in this study. RESULTS: Two lead molecules, ZINC08825480 (4-bromo-N'-{(E)-[1-phenyl-3-(pyridin-3-yl)-1H-pyrazol- 4-yl]methylidene}benzene-1-sulfonohydrazide) and ZINC72009942 (N-[[2-[[(3S)-3-methyl-1-piperidyl] methyl]phenyl]methyl]-6-oxo-1-(p-tolyl)-4,5-dihydro-1,2,4-triazine-3-carboxamide), were identified with better affinity with the three target enzymes as compared to the approved antiviral drugs. Both the lead molecules possessed favorable drug-like properties, fit well into the active site pocket close to His- Cys dyad and showed a good number of hydrogen bonds with the backbone as well as side chains of key amino acid residues. CONCLUSION: Thus, the present study offers two novel chemical entities against coronavirus infections which can be validated through various biological assays.

Exploring the phytochemicals of Platycodon grandiflorus for TMPRSS2 inhibition in the search for SARS-CoV-2 entry inhibitors.

Platycodon grandiflorus (Jacq.) A. DC. (Campanulaceae) is commonly known as a balloon flower whose rhizomes have been widely utilized in traditional Chinese medicine (TCM) and in various Japanese prescriptions for the treatment of respiratory diseases, diabetes, and inflammatory disorders. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of coronavirus disease 2019 (COVID-19) global pandemic requires priming of the virus's spike (S) protein by cleavage of the S proteins by a multi-domain type II transmembrane serine protease, transmembrane protease serine 2 (TMPRSS2) to gain entry into the host cell. The current research aims at the screening of active phytocompounds of P. grandiflorus as potential inhibitors of cellular TMPRSS2 using molecular docking and molecular dynamics simulations approach. In silico toxicity analyses show that out of a total of 34 phytocompounds selected for the study, 12 compounds obey Lipinski's rule of five and have favourable pharmacokinetic properties. The top three lead molecules identified here were Apigenin, Luteolin and Ferulic acid which exhibited binding energies of -7.47 kcal/mol, -6.8 kcal/mol and -6.62 kcal/mol respectively with corresponding inhibition constants of 3.33 µM, 10.39 µM and 13.95 µM. The complexes between the lead molecules and the receptor were held by hydrogen bond interactions with key residues such as Gly383, Gly385, Glu389, Lys390, Asp435, Ser436, Ser441, Cys465 and Lys467, and hydrophobic interactions with surrounding residues. The stability of the protein-ligand complexes was evaluated during 100 ns molecular dynamics (MD) simulation by analysing key geometric properties such as RMSD, RMSF, radius of gyration, total solvent accessible surface area and the number of hydrogen bonds. The binding free energies analysis using MD simulations revealed that the compounds and TMPRSS2 have favourable thermodynamic interactions, which are primarily driven by van der Waals forces. As a result, the selected bioactive phytochemicals from P. grandiflorus that target the cellular TMPRSS2 could offer an alternative treatment option against SARS-CoV-2 infections.

Structural and functional insights into the major mutations of SARS-CoV-2 Spike RBD and its interaction with human ACE2 receptor.

Coronavirus Disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has rapidly spread around the world jeopardizing the global economy and health. The rapid proliferation and infectivity of the virus can be attributed to many accumulating mutations in the spike protein leading to continuous generation of variants. The spike protein is a glycoprotein that recognizes and binds to cell surface receptor known as angiotensin-converting enzyme 2 (ACE2) leading to the fusion of the viral and host cell membranes and entry into the host cells. These circulating variants in the population have greatly impacted the virulence, transmissibility, and immunological evasion of the host. The present study is aimed at understanding the impact of the major mutations (L452R, T478K and N501Y) in the receptor-binding domain (RBD) of spike protein and their consequences on the binding affinity to human ACE2 through protein-protein docking and molecular dynamics simulation approaches. Protein-protein docking and Molecular mechanics with generalised Born and surface area solvation (MM/GBSA) binding free energy analysis reveal that the spike mutants-L452R, T478K and N501Y have a higher binding affinity to human ACE2 as compared to the native spike protein. The increase in the number of interface residues, interface area and intermolecular forces such as hydrogen bonds, salt bridges and non-bonded contacts corroborated with the increase in the binding affinity of the spike mutants to ACE2. Further, 75 ns all-atom molecular dynamics simulation investigations show variations in the geometric properties such as root mean square deviation (RMSD), radius of gyration (Rg), total solvent accessible surface area (SASA) and number of hydrogen bonds (NHBs) in the mutant spike:ACE2 complexes with respect to the native spike:ACE2 complex. Therefore, the findings of this study unravel plausible molecular mechanisms of increase in binding affinity of spike mutants (L452R, T478K and N501Y) to human ACE2 leading to higher virulence and infectivity of emerging SARS-CoV-2 variants. The study will further aid in designing novel therapeutics targeting the interface residues between spike protein and ACE2 receptor.