Drug Dissolution in Oral Drug Absorption: Workshop Report.

The in-person workshop "Drug Dissolution in Oral Drug Absorption" was held on May 23-24, 2023, in Baltimore, MD, USA. The workshop was organized into lectures and breakout sessions. Three common topics that were re-visited by various lecturers were amorphous solid dispersions (ASDs), dissolution/permeation interplay, and in vitro methods to predict in vivo biopharmaceutics performance and risk. Topics that repeatedly surfaced across breakout sessions were the following: (1) meaning and assessment of "dissolved drug," particularly of poorly water soluble drug in colloidal environments (e.g., fed conditions, ASDs); (2) potential limitations of a test that employs sink conditions for a poorly water soluble drug; (3) non-compendial methods (e.g., two-stage or multi-stage method, dissolution/permeation methods); (4) non-compendial conditions (e.g., apex vessels, non-sink conditions); and (5) potential benefit of having both a quality control method for batch release and a biopredictive/biorelevant method for biowaiver or bridging scenarios. An identified obstacle to non-compendial methods is the uncertainty of global regulatory acceptance of such methods.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa.

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

Interplay between gene expression and gene architecture as a consequence of gene and genome duplications: evidence from metabolic genes of Arabidopsis thaliana.

Gene and genome duplications have been widespread during the evolution of flowering plant which resulted in the increment of biological complexity as well as creation of plasticity of a genome helping the species to adapt to changing environments. Duplicated genes with higher evolutionary rates can act as a mechanism of generating novel functions in secondary metabolism. In this study, we explored duplication as a potential factor governing the expression heterogeneity and gene architecture of Primary Metabolic Genes (PMGs) and Secondary Metabolic Genes (SMGs) of Arabidopsis thaliana. It is remarkable that different types of duplication processes controlled gene expression and tissue specificity differently in PMGs and SMGs. A complex relationship exists between gene architecture and expression patterns of primary and secondary metabolic genes. Our study reflects, expression heterogeneity and gene structure variation of primary and secondary metabolism in Arabidopsis thaliana are partly results of duplication events of different origins. Our study suggests that duplication has differential effect on PMGs and SMGs regarding expression pattern by controlling gene structure, epigenetic modifications, multifunctionality and subcellular compartmentalization. This study provides an insight into the evolution of metabolism in plants in the light of gene and genome scale duplication. Supplementary Information: The online version contains supplementary material available at 10.1007/s12298-022-01188-2.

Mutations in membrane-fusion subunit of spike glycoprotein play crucial role in the recent outbreak of COVID-19.

Coronavirus disease-2019 (COVID-19), the ongoing pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a major threat to the entire human race. It is reported that SARS-CoV-2 seems to have relatively low pathogenicity and higher transmissibility than previously outbroke SARS-CoV. To explore the reason of the increased transmissibility of SARS-CoV-2 compared with SARS-CoV, we have performed a comparative analysis on the structural proteins (spike, envelope, membrane, and nucleoprotein) of two viruses. Our analysis revealed that extensive substitutions of hydrophobic to polar and charged amino acids in spike glycoproteins of SARS-CoV2 creates an intrinsically disordered region (IDR) at the beginning of membrane-fusion subunit and intrinsically disordered residues in fusion peptide. IDR provides a potential site for proteolysis by furin and enriched disordered residues facilitate prompt fusion of the SARS-CoV2 with host membrane by recruiting molecular recognition features. Here, we have hypothesized that mutation-driven accumulation of intrinsically disordered residues in spike glycoproteins play dual role in enhancing viral transmissibility than previous SARS-coronavirus. These analyses may help in epidemic surveillance and preventive measures against COVID-19.

Chaperone client proteins evolve slower than non-client proteins.

Chaperones are important molecular machinery that assists proteins to attain their native three-dimensional structure crucial for function. Earlier studies using experimental evolution showed that chaperones impose a relaxation of sequence constraints on their "client" proteins, which may lead to the fixation of slightly deleterious mutations on the latter. However, we hypothesized that such a phenomenon might be harmful to the organism in a natural physiological condition. In this study, we investigated the evolutionary rates of chaperone client and non-client proteins in five model organisms from both prokaryotic and eukaryotic lineages. Our study reveals a slower evolutionary rate of chaperone client proteins in all five organisms. Additionally, the slower folding rate and lower aggregation propensity of chaperone client proteins reveal that the chaperone may play an essential role in rescuing the slightly disadvantageous effects due to random mutations and subsequent protein misfolding. However, the fixation of such mutations is less likely to be selected in the natural population.

The combined influence of codon composition and tRNA copy number regulates translational efficiency by influencing synonymous nucleotide substitution.

Codon usage bias is an important genomic phenomenon, where highly expressed genes use optimal codons for smoother translation with high yield, facilitated by the cognate tRNAs. Here, we presented the tRNA co-adaptation index (co-AI) by correlating tRNA gene copy number and codon composition in Saccharomyces cerevisiae. We observed that this co-AI is positively correlated with protein abundance and translation rate. Considering nucleotide substitutions, co-AI influences synonymous substitutions more than gene expression and protein abundance, the most important determinants of evolutionary rate. Co-AI correlates positively with mRNA secondary structure stability and mRNA half-life, which may lead to protein accumulation under high co-AI. However, the highly expressed proteins encoded by high co-AI genes are assisted by molecular chaperones to attain their proper functional conformation and prevent accumulation.

Deciphering the intrinsic properties of fungal proteases in optimizing phytopathogenic interaction.

Phytopathogenic fungi secrete a wide range of enzymes to penetrate and colonize host tissues. Of them protease activity is reported to increase disease aggressiveness in the plant. With the aim to explore the reason of the higher infection potential of proteases, we have compared several genomic and proteomic attributes among different hydrolytic enzymes coded by five pathogenic fungal species which are the potent infectious agents of plant. Categorizing the enzymes into four major groups, namely protease, lipase, amylase and cell-wall degraders, we observed that proteases are evolutionary more conserved, have higher expression levels, contain more hydrophobic buried residues, short linear motifs and post-translational modified (PTM) sites than the other three groups of enzymes. Again, comparing these features of protease between pathogenic and non-pathogenic Aspergillus sps, we have hypothesized that protein structural properties could play significant roles in imposing infection potency to the fungal proteases.

The optimization of mRNA expression level by its intrinsic properties-Insights from codon usage pattern and structural stability of mRNA.

Codon usage bias (CUB) and mRNA structural stability are important intrinsic features of mRNA that correlate positively with mRNA expression level. However, it remains unclear whether the mRNA expression level can be regulated by adjusting these two parameters, influencing the mRNAs' structure. Here we explored the influence of CUB and mRNA structural stability on mRNA expression levels in Saccharomyces cerevisiae, using both wild type and computationally mutated mRNAs. Although in wild type, both CUB and mRNA stability positively regulate the mRNA expression level, any deviation from natural situation breaks such equilibrium. The naturally occurring codon composition is responsible for optimizing the mRNA expression, and under such composition, the mRNA structure having highest stability is selected by nature.

Exploring the major cross-talking edges of competitive endogenous RNA networks in human Chronic and Acute Myeloid Leukemia.

BACKGROUND: Human Chronic and Acute Myeloid Leukemia are myeloproliferative disorders in myeloid lineage of blood cells characterized by accumulation of aberrant white blood cells. In cancer, the anomalous transcriptome includes deregulated expression of non-coding RNAs in conjunction with protein-coding mRNAs in human genome. The coding or non-coding RNA transcripts harboring miRNA-binding sites can converse with and regulate each other by explicitly contending for a limited pool of shared miRNAs and act as competitive endogenous RNAs (ceRNAs). An unifying hypothesis attributing 'modulation of expression of transcripts' in this fashion had been defined as 'competitive endogenous RNA hypothesis'. Network built with ceRNAs evidently offers a platform to elucidate complex regulatory interactions at post-transcriptional level in human cancers. METHODS: Contemplating cancers of human myeloid lineage we constructed ceRNA networks for CML and AML coding and non-coding repertoire utilizing patient sample data. Through functional enrichment analysis we selected the significant functional modules for transcripts being differentially expressed in Blastic phases of each cancer types with respect to Normal. After retrieving free energy of binding and duplex formation of shared miRNAs on ceRNAs, we performed statistical averaging of energy values over the ensemble of populations considering cellular system as in canonical (Iso-thermal) situation. RESULTS AND CONCLUSIONS: We aimed to shed light on 'Sibling Rivalry' in ceRNA partners from the perspective of statistical thermodynamics, identified major cross-talking tracks and ceRNAs influencing transcripts concerned in myeloid cancer systems. GENERAL SIGNIFICANCE: Insights into ceRNA-regulation will shed light on progression and prognosis of human Chronic and Acute Myeloid Leukemia.

Evolutionary rate heterogeneity between multi- and single-interface hubs across human housekeeping and tissue-specific protein interaction network: Insights from proteins' and its partners' properties.

Integrating gene expression into protein-protein interaction network (PPIN) leads to the construction of tissue-specific (TS) and housekeeping (HK) sub-networks, with distinctive TS- and HK-hubs. All such hub proteins are divided into multi-interface (MI) hubs and single-interface (SI) hubs, where MI hubs evolve slower than SI hubs. Here we explored the evolutionary rate difference between MI and SI proteins within TS- and HK-PPIN and observed that this difference is present only in TS, but not in HK-class. Next, we explored whether proteins' own properties or its partners' properties are more influential in such evolutionary discrepancy. Statistical analyses revealed that this evolutionary rate correlates negatively with protein's own properties like expression level, miRNA count, conformational diversity and functional properties and with its partners' properties like protein disorder and tissue expression similarity. Moreover, partial correlation and regression analysis revealed that both proteins' and its partners' properties have independent effects on protein evolutionary rate.

The role of introns in the conservation of the metabolic genes of Arabidopsis thaliana.

In Arabidopsis thaliana, primary metabolic genes (PMGs) are more evolutionarily conserved and intron-rich than secondary metabolic genes. We observed that PMGs are more primitive and pan-taxonomically persistent as compared to secondary (SMGs) and non-metabolic genes (NMGs). This difference in primitiveness and persistence is primarily correlated with intron number and is independent of gene expression level. We propose a twofold explanation behind higher intron enrichment in PMGs. Firstly, introns might increase protein versatility amongst PMGs through alternative splicing, providing selective advantage of PMGs and making them more persistent across diverse plant taxa. Also, multifunctional PMGs may acquire functional domains by increasing the intronic burden. Additionally, single nucleotide polymorphisms (SNPs) accumulate at a higher rate in introns as compared to exons. Moreover, a strong negative correlation between cumulative exonic SNPs density and intron number indicates that introns may protect the exonic regions against the deleterious effect of these mutations, making them more conserved.

Impact of the US FDA "Biopharmaceutics Classification System" (BCS) Guidance on Global Drug Development.

The FDA guidance on application of the biopharmaceutics classification system (BCS) for waiver of in vivo bioequivalence (BE) studies was issued in August 2000. Since then, this guidance has created worldwide interest among biopharmaceutical scientists in regulatory agencies, academia, and industry toward its implementation and further expansion. This article describes how the review implementation of this guidance was undertaken at the FDA and results of these efforts over last dozen years or so across the new, and the generic, drug domains are provided. Results show that greater than 160 applications were approved, or tentatively approved, based on the BCS approach across multiple therapeutic areas; an additional significant finding was that at least 50% of these approvals were in the central nervous system (CNS) area. These findings indicate a robust utilization of the BCS approach toward reducing unnecessary in vivo BE studies and speeding up availability of high quality pharmaceutical products. The article concludes with a look at the adoption of this framework by regulatory and health policy organizations across the globe, and FDA's current thinking on areas of improvement of this guidance.

Insights into human intrinsically disordered proteins from their gene expression profile.

Expression level provides important clues about gene function. Previously, various efforts have been undertaken to profile human genes according to their expression level. Intrinsically disordered proteins (IDPs) do not adopt any rigid conformation under physiological conditions, however, are considered as an important functional class in all domains of life. Based on a human tissue-averaged gene expression level, previous studies showed that IDPs are expressed at a lower level than ordered globular proteins. Here, we examined the gene expression pattern of human ordered and disordered proteins in 32 normal tissues. We noticed that in most of the tissues, ordered and disordered proteins are expressed at a similar level. Moreover, in a number of tissues IDPs were found to be expressed at a higher level than ordered proteins. Rigorous statistical analyses suggested that the lower tissue-averaged gene expression level of IDPs (reported earlier) may be the consequence of their biased gene expression in some specific tissues and higher protein length. When we considered the gene repertory of each tissue we noticed that a number of human tissues (brain, testes, etc.) selectively express a higher fraction of disordered proteins, which help them to maintain higher protein connectivity by forming disordered binding motifs and to sustain their functional specificities. Our results demonstrated that the disordered proteins are indispensable in these tissues for their functional advantages.

Systematic Analyses and Prediction of Human Drug Side Effect Associated Proteins from the Perspective of Protein Evolution.

Identification of various factors involved in adverse drug reactions in target proteins to develop therapeutic drugs with minimal/no side effect is very important. In this context, we have performed a comparative evolutionary rate analyses between the genes exhibiting drug side-effect(s) (SET) and genes showing no side effect (NSET) with an aim to increase the prediction accuracy of SET/NSET proteins using evolutionary rate determinants. We found that SET proteins are more conserved than the NSET proteins. The rates of evolution between SET and NSET protein primarily depend upon their noncomplex (protein complex association number = 0) forming nature, phylogenetic age, multifunctionality, membrane localization, and transmembrane helix content irrespective of their essentiality, total druggability (total number of drugs/target), m-RNA expression level, and tissue expression breadth. We also introduced two novel terms-killer druggability (number of drugs with killing side effect(s)/target), essential druggability (number of drugs targeting essential proteins/target) to explain the evolutionary rate variation between SET and NSET proteins. Interestingly, we noticed that SET proteins are younger than NSET proteins and multifunctional younger SET proteins are candidates of acquiring killing side effects. We provide evidence that higher killer druggability, multifunctionality, and transmembrane helices support the conservation of SET proteins over NSET proteins in spite of their recent origin. By employing all these entities, our Support Vector Machine model predicts human SET/NSET proteins to a high degree of accuracy (∼86%).

Codon degeneracy and amino acid abundance influence the measures of codon usage bias: improved Nc (N̂c ) and ENCprime (N̂'c ) measures.

Effective number of codons (N^c) and its variant N^'c (effective number of codons prime) are the two widely used methods for measuring unequal usage of synonymous codons in coding sequences, known as the codon usage bias (CUB). The mathematical formula used in calculating N^c and N^'c values is giving inappropriate measures of CUB in case of low abundance of amino acids. In addition, the magnitude of error also varies according to codon degeneracy. In this study, a modified formula for N^c and N^'c has been developed to measure the CUB more accurately. Online implementations of the modified formula are available in the web portal at http://agnigarh.tezu.ernet.in/~ssankar/cub.php.

Overlapping Regions in HIV-1 Genome Act as Potential Sites for Host-Virus Interaction.

More than a decade, overlapping genes in RNA viruses became a subject of research which has explored various effect of gene overlapping on the evolution and function of viral genomes like genome size compaction. Additionally, overlapping regions (OVRs) are also reported to encode elevated degree of protein intrinsic disorder (PID) in unspliced RNA viruses. With the aim to explore the roles of OVRs in HIV-1 pathogenesis, we have carried out an in-depth analysis on the association of gene overlapping with PID in 35 HIV1- M subtypes. Our study reveals an over representation of PID in OVR of HIV-1 genomes. These disordered residues endure several vital, structural features like short linear motifs (SLiMs) and protein phosphorylation (PP) sites which are previously shown to be involved in massive host-virus interaction. Moreover, SLiMs in OVRs are noticed to be more functionally potential as compared to that of non-overlapping region. Although, density of experimentally verified SLiMs, resided in 9 HIV-1 genes, involved in host-virus interaction do not show any bias toward clustering into OVR, tat and rev two important proteins mediates host-pathogen interaction by their experimentally verified SLiMs, which are mostly localized in OVR. Finally, our analysis suggests that the acquisition of SLiMs in OVR is mutually exclusive of the occurrence of disordered residues, while the enrichment of PPs in OVR is solely dependent on PID and not on overlapping coding frames. Thus, OVRs of HIV-1 genomes could be demarcated as potential molecular recognition sites during host-virus interaction.

Deciphering the cross-talking of human competitive endogenous RNAs in K562 chronic myelogenous leukemia cell line.

Chronic myelogenous leukemia (CML) is a myeloproliferative disorder characterized by increased proliferation or abnormal accumulation of granulocytic cell line without the depletion of their capacity to differentiate. A reciprocal chromosomal translocation proceeding to the 'Philadelphia chromosome', involving the ABL proto-oncogene and BCR gene residing on Chromosome 9 and 22 respectively, is observed to be attributed to CML pathogenesis. Recent studies have been unraveling the crucial role of genomic 'dark matter' or the non-coding repertoire in cancer initiation and progression. The intricate cross-talk between competitive endogenous RNAs (ceRNAs) provides a scaffold to systematically functionalize the miRNA response element harboring non-coding RNAs and incorporate them with the protein-coding RNA dimension in complex ceRNA networks. This network of coding and non-coding transcriptome linked by shared miRNAs evidently offers a platform to elucidate the complex regulatory interactions at the post-transcriptional level in human cancers. In this context, analyzing CML, from the perspective of the ceRNA hypothesis, surely craves intensive attention and a comprehensive discussion. Here, we performed RNA-seq data analysis to retrieve Lymphoblastoid and CML coding as well as non-coding repertoire and constructed a ceRNA network for the CML cell line, considering the non-cancer lymphoblastoid cell line as the control. We investigated if any alteration exists in the ceRNA landscape of the transcripts which are exhibiting differential expression across the two cell lines and observed that the major ceRNA regulators vary in cancer network when compared with the Lymphoblastoid network. The top ranked significant functional modules in the ceRNA network display cancer associated attributes and reveal putative regulators in CML pathogenesis.

Navigating sticky areas in transdermal product development.

The benefits of transdermal delivery over the oral route to combat such issues of low bioavailability and limited controlled release opportunities are well known and have been previously discussed by many in the field (Prausnitz et al. (2004) [1]; Hadgraft and Lane (2006) [2]). However, significant challenges faced by developers as a product moves from the purely theoretical to commercial production have hampered full capitalization of the dosage forms vast benefits. While different technical aspects of transdermal system development have been discussed at various industry meetings and scientific workshops, uncertainties have persisted regarding the pharmaceutical industry's conventionally accepted approach for the development and manufacturing of transdermal systems. This review provides an overview of the challenges frequently faced and the industry's best practices for assuring the quality and performance of transdermal delivery systems and topical patches (collectively, TDS). The topics discussed are broadly divided into the evaluation of product quality and the evaluation of product performance; with the overall goal of the discussion to improve, advance and accelerate commercial development in the area of this complex controlled release dosage form.

Insights into the dN/dS ratio heterogeneity between brain specific genes and widely expressed genes in species of different complexity.

In mammals, it has long been suggested that brain-specific genes (BSGs) and widely expressed genes (WEGs) have seemingly lower dN/dS ratio than any other gene sets. However, to what extent these genes differ in their dN/dS ratio has still remained controversial. Here, we have revealed lower dN/dS ratio of BSGs than WEGs in human-mouse, human-orangutan, human-chimpanzee and mouse-rat orthologous pair. The significance level of dN/dS ratio difference indicates a trend of decreasing difference as complexity of compared pairs increases. Further studies with the human-mouse pair revealed that, removal of the duplicated genes from both the dataset has nullified this difference which dictates a vital role of duplicated genes in governing the selection pressure. Conclusively, higher paralog number, expression level, and longer regulatory region length of BSGs allow fewer nucleotide substitutions within them. Our results show for the first time to our knowledge lower dN/dS ratio of BSGs than WEGs.

Regulatory Perspectives on Strength-Dependent Dissolution Profiles and Biowaiver Approaches for Immediate Release (IR) Oral Tablets in New Drug Applications.

Dissolution profile comparisons are used by the pharmaceutical industry to assess the similarity in the dissolution characteristics of two formulations to decide whether the implemented changes, usually minor/moderate in nature, will have an impact on the in vitro/in vivo performance of the drug product. When similarity testing is applied to support the approval of lower strengths of the same formulation, the traditional approach for dissolution profile comparison is not always applicable for drug products exhibiting strength-dependent dissolution and may lead to incorrect conclusions about product performance. The objective of this article is to describe reasonable biopharmaceutic approaches for developing a biowaiver strategy for low solubility, proportionally similar/non-proportionally similar in composition immediate release drug products that exhibit strength-dependent dissolution profiles. The paths highlighted in the article include (1) approaches to address biowaiver requests, such as the use of multi-unit dissolution testing to account for sink condition differences between the higher and lower strengths; (2) the use of a single- vs. strength-dependent dissolution method; and (3) the use of single- vs. strength-dependent dissolution acceptance criteria. These approaches are cost- and time-effective and can avoid unnecessary bioequivalence studies.