A Second-Generation Oral SARS-CoV-2 Main Protease Inhibitor Clinical Candidate for the Treatment of COVID-19

Despite the record-breaking discovery, development and approval of vaccines and antiviral therapeutics such as Paxlovid, coronavirus disease 2019 (COVID-19) remained the fourth leading cause of death in the world and third highest in the United States in 2022. Here, we report the discovery and characterization of PF-07817883, a second-generation, orally bioavailable, SARS-CoV-2 main protease inhibitor with improved metabolic stability versus nirmatrelvir, the antiviral component of the ritonavir-boosted therapy Paxlovid. We demonstrate the in vitro pan-human coronavirus antiviral activity and off-target selectivity profile of PF-07817883. PF-07817883 also demonstrated oral efficacy in a mouse-adapted SARS-CoV-2 model at plasma concentrations equivalent to nirmatrelvir. The preclinical in vivo pharmacokinetics and metabolism studies in human matrices are suggestive of improved oral pharmacokinetics for PF-07817883 in humans, relative to nirmatrelvir. In vitro inhibition/induction studies against major human drug metabolizing enzymes/transporters suggest a low potential for perpetrator drug-drug interactions upon single-agent use of PF-07817883.

Characterization and Applications of Permeabilized Hepatocytes in Drug Discovery

Hepatocytes are one of the most physiologically relevant in vitro liver systems for human translation of clearance and drug-drug interactions (DDI). However, the cell membranes of hepatocytes can limit the entry of certain compounds into the cells for metabolism and DDI. Passive permeability through hepatocytes can be different in vitro and in vivo, which complicates the human translation. Permeabilized hepatocytes offer a useful tool to probe mechanistic understanding of permeability-limited metabolism and DDI. Incubation with saponin of 0.01% at 0.5 million cells/mL and 0.05% at 5 million cells/mL for 5 min at 37°C completely permeabilized the plasma membrane of hepatocytes, while leaving the membranes of subcellular organelles intact. Permeabilized hepatocytes maintained similar enzymatic activity as intact unpermeabilized hepatocytes and can be stored at -80°C for at least 7 months. This approach reduces costs by preserving leftover hepatocytes. The relatively low levels of saponin in permeabilized hepatocytes had no significant impact on the enzymatic activity. As the cytosolic contents leak out from permeabilized hepatocytes, cofactors need to be added to enable metabolic reactions. Cytosolic enzymes will no longer be present if the media are removed after cells are permeabilized. Hence permeabilized hepatocytes with and without media removal may potentially enable reaction phenotyping of cytosolic enzymes. Although permeabilized hepatocytes work similarly as human liver microsomes and S9 fractions experimentally requiring addition of cofactors, they behave more like hepatocytes maintaining enzymatic activities for over 4 h. Permeabilized hepatocytes are a great addition to the drug metabolism toolbox to provide mechanistic insights.

Novel Multiplexed High Throughput Screening of Selective Inhibitors for Drug-Metabolizing Enzymes Using Human Hepatocytes

Selective chemical inhibitors are critical for reaction phenotyping to identify drug-metabolizing enzymes that are involved in the elimination of drug candidates. Although relatively selective inhibitors are available for the major cytochrome P450 enzymes (CYP), they are quite limited for the less common CYPs and non-CYPs. To address this gap, we developed a multiplexed high throughput screening (HTS) assay using 20 substrate reactions of multiple enzymes to simultaneously monitor the inhibition of enzymes in a 384-well format. Four 384-well assay plates can be run at the same time to maximize throughput. This is the first multiplexed HTS assay for drug-metabolizing enzymes reported. The HTS assay is technologically enabled with state-of-the-art robotic systems and highly sensitive modern LC-MS/MS instrumentation. Virtual screening is utilized to identify inhibitors for HTS based on known inhibitors and enzyme structures. Screening of ~4600 compounds generated many hits for many drug-metabolizing enzymes including the two time-dependent and selective aldehyde oxidase inhibitors, erlotinib and dibenzothiophene. The hit rate is much higher than that for the traditional HTS for biological targets due to the promiscuous nature of the drug-metabolizing enzymes and the biased compound selection process. Future efforts will focus on using this method to identify selective inhibitors for enzymes that do not currently have quality hits and thoroughly characterizing the newly identified selective inhibitors from our screen. We encourage colleagues from other organizations to explore their proprietary libraries using a similar approach to identify better inhibitors that can be used across the industry.

Non-coding CGG repeat expansion in LOC642361/NUTM2B-AS1 is associated with a phenotype of oculopharyngodistal myopathy

BACKGROUND

Oculopharyngodistal myopathy (OPDM) is a rare adult-onset neuromuscular disease, associated with CGG repeat expansions in the 5' untranslated region of LRP12, GIPC1, NOTCH2NLC and RILPL1. However, the genetic cause of a proportion of pathoclinically confirmed cases remains unknown.

METHODS

A total of 26 OPDM patients with unknown genetic cause(s) from 4 tertiary referral hospitals were included in this study. Clinical data and laboratory findings were collected. Muscle samples were observed by histological and immunofluorescent staining. Long-read sequencing was initially conducted in six patients with OPDM. Repeat-primed PCR was used to screen the CGG repeat expansions in LOC642361/NUTM2B-AS1 in all 26 patients.

RESULTS

We identified CGG repeat expansion in the non-coding transcripts of LOC642361/NUTM2B-AS1 in another two unrelated Chinese cases with typical pathoclinical features of OPDM. The repeat expansion was more than 70 times in the patients but less than 40 times in the normal controls. Both patients showed no leucoencephalopathy but one showed mild cognitive impairment detected by Montreal Cognitive Assessment. Rimmed vacuoles and p62-positive intranuclear inclusions (INIs) were identified in muscle pathology, and colocalisation of CGG RNA foci with p62 was also found in the INIs of patient-derived fibroblasts.

CONCLUSIONS

We identified another two unrelated cases with CGG repeat expansion in the long non-coding RNA of the LOC642361/NUTM2B-AS1 gene, presenting with a phenotype of OPDM. Our cases broadened the recognised phenotypic spectrum and pathogenesis in the disease associated with CGG repeat expansion in LOC642361/NUTM2B-AS1.

Comparison of Tumor Binding Across Tumor Types and Cell Lines to Support Free Drug Considerations for Oncology Drug Discovery

Tumor binding is an important parameter to derive unbound tumor concentration to explore pharmacokinetics (PK) and pharmacodynamics (PD) relationships for oncology disease targets. Tumor binding was evaluated using eleven matrices, including various commonly used ex vivo human and mouse xenograft and syngeneic tumors, tumor cell lines and liver as a surrogate tissue. The results showed that tumor binding is highly correlated among the different tumors and tumor cell lines except for the mouse melanoma (B16F10) tumor type. Liver fraction unbound (fu) has a good correlation with B16F10 tumor binding. Liver also demonstrates a two-fold equivalency, on average, with binding of other tumor types when a scaling factor is applied. Predictive models were developed for tumor binding, with correlations established with LogD (acids), predicted muscle fu (neutrals) and measured plasma protein binding (bases) to estimate tumor fu when experimental data are not available. Many approaches can be applied to obtain and estimate tumor binding values. One strategy proposed is to use a surrogate tumor tissue, such as mouse xenograft ovarian cancer (OVCAR3) tumor, as a surrogate for tumor binding (except for B16F10) to provide an early assessment of unbound tumor concentrations for development of PK/PD relationships.

Environmental factors affecting the risk of generalization for ocular-onset myasthenia gravis: a nationwide cohort study

BACKGROUND

The environmental effects on the prognosis of ocular myasthenia gravis (OMG) remain largely unexplored.

AIM

To investigate the association between specific environmental factors and the generalization of OMG.

DESIGN

The cohort study was conducted in China based on a nationwide multicenter database.

METHODS

Adult patients with OMG at onset, who were followed up for at least 2 years until May 2022, were included. We collected data on demographic and clinical factors, as well as environmental factors, including latitude, socioeconomic status (per capita disposable income [PDI] at provincial level and education) and smoking. The study outcome was the time to the development of generalized myasthenia gravis (GMG). Cox models were employed to examine the association between environmental exposures and generalization. Restricted cubic spline was used to model the association of latitude with generalization risk.

RESULTS

A total of 1396 participants were included. During a median follow-up of 5.15 (interquartile range [IQR] 3.37-9.03) years, 735 patients developed GMG within a median of 5.69 (IQR 1.10-15.66) years. Latitude of 20-50°N showed a U-shaped relation with generalization risk, with the lowest risk at around 30°N; both higher and lower latitudes were associated with the increased risk (P for non-linearity <0.001). Living in areas with lower PDI had 1.28-2.11 times higher risk of generalization. No significant association was observed with education or smoking.

CONCLUSIONS

Latitude and provincial-level PDI were associated with the generalization of OMG in China. Further studies are warranted to validate our findings and investigate their potential applications in clinical practice and health policy.

Characterization of CYP3A5 Selective Inhibitors for Reaction Phenotyping of Drug Candidates

CYP3A is one of the most important classes of enzymes and is involved in the metabolism of over 70% drugs. While several selective CYP3A4 inhibitors have been identified, the search for a selective CYP3A5 inhibitor has turned out to be rather challenging. Recently, several selective CYP3A5 inhibitors have been identified through high-throughput screening of ~ 11,000 compounds and hit expansion using human recombinant enzymes. We set forth to characterize the three most selective CYP3A5 inhibitors in a more physiologically relevant system of human liver microsomes to understand if these inhibitors can be used for reaction phenotyping studies in drug discovery settings. Gomisin A and T-5 were used as selective substrate reactions for CYP3A4 and CYP3A5 to determine IC50 values of the two enzymes. The results showed that clobetasol propionate and loteprednol etabonate were potent and selective CYP3A5 reversible inhibitors with selectivity of 24-fold against CYP3A4 and 39-fold or more against the other major CYPs. The selectivity of difluprednate in HLM is much weaker than that in the recombinant enzymes due to hydrolysis of the acetate group in HLM. Based on the selectivity data, loteprednol etabonate can be utilized as an orthogonal approach, when experimental fraction metabolized of CYP3A5 is greater than 0.5, to understand CYP3A5 contribution to drug metabolism and its clinical significance. Future endeavors to identify even more selective CYP3A5 inhibitors are warranted to enable accurate determination of CYP3A5 contribution to metabolism versus CYP3A4.

PBPK Modeling of PAXLOVIDTM: Incorporating Rotamer Conversion Kinetics to Advanced Dissolution and Absorption Model

PAXLOVID™ is a combination medicine of nirmatrelvir tablets co-packaged with ritonavir tablets. Nirmatrelvir is a peptidomimetic inhibitor of SARS-CoV2 main protease (Mpro), developed for the treatment of COVID-19. Ritonavir is co-administered as a pharmacokinetics (PK) enhancer to inhibit CYP3A mediated metabolism increasing exposures of nirmatrelvir. In the solid form, nirmatrelvir exists in a stable single conformational state (ANTI form). However, nirmatrelvir exhibits atropisomerism in solution whereby upon dissolution the ANTI rotational isomer reversibly converts to another conformation state (SYN form). Nirmatrelvir rotamer conversion follows pseudo first order kinetics with a conversion half-life of approximately 15 min in aqueous solutions, which is on a similar time scale of diffusion mediated dissolution from the solid form. In vitro dissolution studies further indicated that rotamer conversion is one of the processes controlling nirmatrelvir dissolution. It was hypothesized that rotamer conversion kinetics would affect oral absorption of nirmatrelvir in vivo. Consequently, a physiologically based pharmacokinetic (PBPK) model for Paxlovid was developed in Simcyp™ using the advanced dissolution, absorption, and metabolism model (ADAM) by incorporating rotamer conversion kinetics to achieve a more mechanistic description of nirmatrelvir oral absorption. The results demonstrate that the established absorption model with rotamer kinetics adequately described observed clinical data from various nirmatrelvir doses, dosage forms, and dosing regimens. The predicted vs. observed AUCinf and Cmax ratios were within 2-fold. The model has been internally used to inform clinical studies and dose recommendations for pediatrics.

The Comparison of Machine Learning and Mechanistic In Vitro-In Vivo Extrapolation Models for the Prediction of Human Intrinsic Clearance

Accurate prediction of human pharmacokinetics (PK) remains one of the key objectives of drug metabolism and PK (DMPK) scientists in drug discovery projects. This is typically performed by using in vitro-in vivo extrapolation (IVIVE) based on mechanistic PK models. In recent years, machine learning (ML), with its ability to harness patterns from previous outcomes to predict future events, has gained increased popularity in application to absorption, distribution, metabolism, and excretion (ADME) sciences. This study compares the performance of various ML and mechanistic models for the prediction of human IV clearance for a large (645) set of diverse compounds with literature human IV PK data, as well as measured relevant in vitro end points. ML models were built using multiple approaches for the descriptors: (1) calculated physical properties and structural descriptors based on chemical structure alone (classical QSAR/QSPR); (2) in vitro measured inputs only with no structure-based descriptors (ML IVIVE); and (3) in silico ML IVIVE using in silico model predictions for the in vitro inputs. For the mechanistic models, well-stirred and parallel-tube liver models were considered with and without the use of empirical scaling factors and with and without renal clearance. The best ML model for the prediction of in vivo human intrinsic clearance (CLint) was an in vitro ML IVIVE model using only six in vitro inputs with an average absolute fold error (AAFE) of 2.5. The best mechanistic model used the parallel-tube liver model, with empirical scaling factors resulting in an AAFE of 2.8. The corresponding mechanistic model with full in silico inputs achieved an AAFE of 3.3. These relative performances of the models were confirmed with the prediction of 16 Pfizer drug candidates that were not part of the original data set. Results show that ML IVIVE models are comparable to or superior to their best mechanistic counterparts. We also show that ML IVIVE models can be used to derive insights into factors for the improvement of mechanistic PK prediction.

Quantitative Translation of Substrate Intrinsic Clearance from Recombinant CYP1A1 to Humans

CYP1A1 is a cytochrome P450 family 1 enzyme that is mostly expressed in the extrahepatic tissues. To understand the CYP1A1 contribution to drug clearance in humans, we examined the in vitro-in vivo extrapolation (IVIVE) of intrinsic clearance (CLint) for a set of drugs that are in vitro CYP1A1 substrates. Despite being strong in vitro CYP1A1 substrates, 82% of drugs gave good IVIVE with predicted CLint within 2-3-fold of the observed values using human liver microsomes and hepatocytes, suggesting they were not in vivo CYP1A1 substrates due to the lack of extrahepatic contribution to CLint. Only three drugs (riluzole, melatonin and ramelteon) that are CYP1A2 substrates yielded significant underprediction of in vivo CLint up to 11-fold. The fold of CLint underprediction was linearly proportional to human recombinant CYP1A1 (rCYP1A1) CLint, indicating they were likely to be in vivo CYP1A1 substrates. Using these three substrates, a calibration curve can be developed to enable direct translation from in vitro rCYP1A1 CLint to in vivo extrahepatic contributions in humans. In vivo CYP1A1 substrates are planar and small, which is consistent with the structure of the active site. This is in contrast to the in vitro substrates, which include large and nonplanar molecules, suggesting rCYP1A1 is more accessible than what is in vivo. The impact of CYP1A1 on first-pass intestinal metabolism was also evaluated and shown to be minimal. This is the first study providing new insights on in vivo translation of CYP1A1 contributions to human clearance using in vitro rCYP1A1 data.

Clinical features of MOGAD with brainstem involvement in the initial attack versus NMOSD and MS

OBJECTIVE

To assess the characteristics of Myelin oligodendrocyte glycoprotein (MOG) antibody-associated disorder (MOGAD) with brainstem involvement in the first event (BSIFE) and make comparisons with aquaporin-4-IgG seropositive neuromyelitis optica spectrum disorder (AQP4-IgG-NMOSD) and multiple sclerosis (MS).

METHODS

From 2017 to 2022, this study identified MOG-IgG-positive patients with brainstem or both brainstem and cerebellum lesions in the first episode. As a comparison group, AQP4-IgG-NMOSD (n = 30) and MS (n = 30) patients with BSIFE were enroled.

RESULTS

Thirty-five patients (35/146, 24.0%) were the BSIFE of MOGAD. Isolated brainstem episodes occurred in 9 of the 35 (25.7%) MOGAD patients, which was similar to MS (7/30, 23.3%) but was lower than AQP4-IgG-NMOSD (17/30, 56.7%, P = 0.011). Pons (21/35, 60.0%), medulla oblongata (20/35, 57.1%) and middle cerebellar peduncle (MCP, 19/35, 54.3%) were the most frequently affected areas. Intractable nausea (n = 7), vomiting (n = 8) and hiccups (n = 2) happened in MOGAD patients, but EDSS of MOGAD was lower than AQP4-IgG-NMOSD (P = 0.001) at the last follow-up. MOGAD patients with or without BSIFE did not significantly differ in terms of the ARR (P = 0.102), mRS (P = 0.823), or EDSS (P = 0.598) at the most recent follow-up. Specific oligoclonal bands appeared in MOGAD (13/33, 39.4%) and AQP4-IgG-NMOSD (7/24, 29.2%) in addition to MS (20/30, 66.7%). Fourteen MOGAD patients (40.0%) experienced relapse in this study. When the brainstem was involved in the first attack, there was an increased likelihood of a second attack occurring at the same location (OR=12.22, 95%CI 2.79 to 53.59, P = 0.001). If the first and second events were both in the brainstem, the third event was likely to occur at the same location (OR=66.00, 95%CI 3.47 to 1254.57, P = 0.005). Four patients experienced relapses after the MOG-IgG turned negative.

CONCLUSION

BSIFE occurred in 24.0% of MOGAD. Pons, medulla oblongata and MCP were the most frequently involved regions. Intractable nausea, vomiting and hiccups occurred in MOGAD and AQP4-IgG-NMOSD, but not MS. The prognosis of MOGAD was better than AQP4-IgG-NMOSD. In contrast to MS, BSIFE may not indicate a worse prognosis for MOGAD. When patients with BSIFE, MOGAD tent to reoccur in the brainstem. Four of the 14 recurring MOGAD patients relapsed after the MOG-IgG test turned negative.

Physiologically-Based Pharmacokinetic Modeling of PAXLOVID™ with First-Order Absorption Kinetics

PURPOSE

PAXLOVID™ is nirmatrelvir tablets co-packaged with ritonavir tablets. Ritonavir is used as a pharmacokinetics (PK) enhancer to reduce metabolism and increase exposure of nirmatrelvir. This is the first disclosure of Paxlovid physiologically-based pharmacokinetic (PBPK) model.

METHODS

Nirmatrelvir PBPK model with first-order absorption kinetics was developed using in vitro, preclinical, and clinical data of nirmatrelvir in the presence and absence of ritonavir. Clearance and volume of distribution were derived from nirmatrelvir PK obtained using a spray-dried dispersion (SDD) formulation where it is considered to be dosed as an oral solution, and absorption is near complete. The fraction of nirmatrelvir metabolized by CYP3A was estimated based on in vitro and clinical ritonavir drug-drug interaction (DDI) data. First-order absorption parameters were established for both SDD and tablet formulation using clinical data. Nirmatrelvir PBPK model was verified with both single and multiple dose human PK data, as well as DDI studies. Simcyp® first-order ritonavir compound file was also verified with additional clinical data.

RESULTS

The nirmatrelvir PBPK model described the observed PK profiles of nirmatrelvir well with predicted AUC and Cmax values within ± 20% of the observed. The ritonavir model performed well resulting in predicted values within twofold of observed.

CONCLUSIONS

Paxlovid PBPK model developed in this study can be applied to predict PK changes in special populations, as well as model the effect of victim and perpetrator DDI. PBPK modeling continues to play a critical role in accelerating drug discovery and development of potential treatments for devastating diseases such as COVID-19. NCT05263895, NCT05129475, NCT05032950 and NCT05064800.

Recent advances in measurement of metabolic clearance, metabolite profile and reaction phenotyping of low clearance compounds

INTRODUCTION

Low metabolic clearance is usually a highly desirable property of drug candidates in order to reduce dose and dosing frequency. However, measurement of low clearance can be challenging in drug discovery. A number of new tools have recently been developed to address the gaps in the measurement of intrinsic clearance, identification of metabolites, and reaction phenotyping of low clearance compounds.

AREAS COVERED

The new methodologies of low clearance measurements are discussed, including the hepatocyte relay, HepatoPac®, HμREL®, and spheroid systems. In addition, metabolite formation rate determination and in vivo allometric scaling approaches are covered as alternative methods for low clearance measurements. With these new methods, measurement of ~ 20-fold lower limit of intrinsic clearance can be achieved. The advantages and limitations of each approach are highlighted.

EXPERT OPINION

Although several novel methods have been developed in recent years to address the challenges of low clearance, these assays tend to be time and labor intensive and costly. Future innovations focusing on developing systems with high enzymatic activities, ultra-sensitive universal quantifiable detectors, and artificial intelligence will further enhance our ability to explore the low clearance space.

A synonymous codon variant altering splicing of RBCK1 expands the phenotype and genotype spectra of polyglucosan body myopathy 1

Polyglucosan body myopathy type 1 (PGBM1, OMIM #615895.) is a rare autosomal recessive disorder caused by RBCK1 mutations. The patients displayed polyglucosan accumulation in skeletal and cardiac muscles, giving rise to loss of ambulation and heart failure with or without immune system dysregulation. So far, only 24 patients have been reported, all of whom exhibited symptoms before adulthood. Here, we reported the first case of an adult-onset PGBM1 patient with a novel compound heterozygous RBCK1 gene mutation consisting of a nonsense and synonymous variant affecting splicing.

In Vitro-In Vivo Extrapolation and Scaling Factors for Clearance of Human and Preclinical Species with Liver Microsomes and Hepatocytes

In vitro-in vivo extrapolation ((IVIVE) and empirical scaling factors (SF) of human intrinsic clearance (CL_int) were developed using one of the largest dataset of 455 compounds with data from human liver microsomes (HLM) and human hepatocytes (HHEP). For extended clearance classification system (ECCS) class 2/4 compounds, linear SFs (SF_lin) are approximately 1, suggesting enzyme activities in HLM and HHEP are similar to those in vivo under physiological conditions. For ECCS class 1A/1B compounds, a unified set of SFs was developed for CL_int. These SFs contain both SF_lin and an exponential SF (SF_β) of fraction unbound in plasma (f_u,p). The unified SFs for class 1A/1B eliminate the need to identify the transporters involved prior to clearance prediction. The underlying mechanisms of these SFs are not entirely clear at this point, but they serve practical purposes to reduce biases and increase prediction accuracy. Similar SFs have also been developed for preclinical species. For HLM-HHEP disconnect (HLM > HHEP) ECCS class 2/4 compounds that are mainly metabolized by cytochrome P450s/FMO, HLM significantly overpredicted in vivo CL_int, while HHEP slightly underpredicted and geometric mean of HLM and HHEP slightly overpredicted in vivo CL_int. This observation is different than in rats, where rat liver microsomal CL_int correlates well with in vivo CL_int for compounds demonstrating permeability-limited metabolism. The good CL_int IVIVE developed using HLM and HHEP helps build confidence for prospective predictions of human clearance and supports the continued utilization of these assays to guide structure-activity relationships to improve metabolic stability.

Comparison between mono-tacrolimus and mono-glucocorticoid in the treatment of myasthenia gravis

OBJECTIVE

Use of tacrolimus in mild to moderate myasthenia gravis (MG) is generally limited to glucocorticoid-refractory cases; the advantage of mono-tacrolimus over mono-glucocorticoids is unknown.

METHODS

We included mild to moderate MG patients treated with mono-tacrolimus (mono-TAC) or mono-glucocorticoids (mono-GC). The correlation between the immunotherapy options and the treatment efficacy and side effects were examined in 1:1 propensity-score matching. The main outcome was time to minimal manifestations status or better (MMS or better). Secondary outcomes include time to relapse, the mean changes in Myasthenia Gravis-specific Activities of Daily Living (MG-ADL) scores and the rate of adverse events.

RESULTS

Baseline characteristics showed no difference between matched groups (49 matched pairs). There were no differences in median time to MMS or better between the mono-TAC group and mono-GC group (5.1 vs. 2.8 months: unadjusted hazard ratio [HR], 0.73; 95% CI, 0.46-1.16; p = 0.180), as well as in median time to relapse (data unavailable for the mono-TAC group since 44 of 49 [89.8%] participants remained in MMS or better; 39.7 months in mono-GC group: unadjusted HR, 0.67; 95% CI, 0.23-1.97; p = 0.464). Changes in MG-ADL scores between the two groups were similar (mean differences, 0.3; 95% CI, -0.4 to 1.0; p = 0.462). The rate of adverse events was lower in the mono-TAC group compared to the mono-GC group (24.5% vs. 55.1%, p = 0.002).

INTERPRETATION

Mono-tacrolimus performs superior tolerability with non-inferior efficacy compared to mono-glucocorticoids in mild to moderate myasthenia gravis patients who refuse or have a contraindication to glucocorticoids.

Species differences in plasma protein binding of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease inhibitor nirmatrelvir

Plasma protein binding (PPB) studies on the SARS-CoV-2 main protease inhibitor nirmatrelvir revealed considerable species differences primarily in dog and rabbit, which prompted further investigations into the biochemical basis for these differences.The unbound fraction (f_u) of nirmatrelvir in dog and rabbit plasma was concentration (2-200 µM)-dependent (dog f_u,p 0.024-0.69, rabbit f_u,p 0.010-0.82). Concentration (0.1-100 µM)-dependent binding in serum albumin (SA) (f_u,SA 0.040-0.82) and alpha-1-acid glycoprotein (AAG) (f_u,AAG 0.050-0.64) was observed in dogs. Nirmatrelvir showed minimal binding to rabbit SA (1-100 µM: f_u,SA 0.70-0.79), while binding to rabbit AAG was concentration-dependent (0.1-100 µM: f_u,AAG 0.024-0.66). In contrast, nirmatrelvir (2 µM) revealed minimal binding (f_u,AAG 0.79-0.88) to AAG from rat and monkeys. Nirmatrelvir showed minimal-to-moderate binding to SA (1-100 µM; f_u,SA 0.70-1.0) and AAG (0.1-100 µM; f_u,AAG 0.48-0.58) from humans across tested concentrations.Nirmatrelvir molecular docking studies using published crystal structures and homology models of human and preclinical species SA and AAG were used to rationalise the species differences to plasma proteins. This suggested that species differences in PPB are primarily driven by molecular differences in albumin and AAG resulting in differences in binding affinity.

Addressing the Accuracy of Plasma Protein Binding Measurement for Highly Bound Compounds Using the Dilution Method

Currently, regulatory guidelines recommend using 0.01 as the lower limit of plasma fraction unbound (f_u) for prediction of drug-drug interactions (DDI) to err on the conservative side. One way to increase experimental f_u of highly bound compounds is to dilute the plasma. With the dilution method, a diluted f_u, or f_u,d, of ≥ 0.01 can be achieved by adjusting the dilution factor. The undiluted f_u can be calculated from f_u,d and be used for DDI prediction. In this study, the dilution method was evaluated, and the results showed that it gave similar f_u values as those determined using the pre-saturation method without plasma dilution. The dilution method enables generation of accurate f_u values and alignment with the regulatory recommendation of reportable f_u values of ≥ 0.01 for DDI prediction. We recommend using the dilution method to bridge the regulatory recommended f_u limit of 0.01 for DDI prediction and the pre-saturation or equivalent methods for definitive plasma protein binding studies. As the pharmaceutical industry continues to generate high quality PPB data, regulatory agencies will gain confidence in the accuracy of f_u measurements for highly bound compounds, and the f_u lower limit may no longer be needed in the future.

Quantum Mechanics Modeling of Oxetanes as Epoxide Hydrolase Substrates

Background:

Epoxide hydrolases comprise an important class of enzymes that have critical functions in the detoxification of xenobiotics and regulation of signaling molecules. In addition to epoxides, oxetanes have recently been identified as novel substrates of microsomal epoxide hydrolase (mEH). Oxetanes are common scaffolds used in medicinal chemistry design to improve potency and drug-like properties. Metabolism of oxetanes by mEH can result in high uncertainties in the prediction of human clearance due to extrahepatic contribution and large interindividual variability. Therefore, reducing mEH-mediated oxetane metabolism is highly desirable to minimize its contribution to clearance.

Objective:

The aim of the study is to evaluate whether quantum mechanical parameters are able to predict the hydrolytic rate of mEH-mediated oxetane metabolism in order to guide medicinal chemistry design in order to minimize mEH contribution to clearance.

Methods:

Quantum mechanics modeling was used to evaluate the hydrolytic rate of twenty-three oxetanes by mEH. All modeling studies were performed with the Maestro software package.

Results:

The results show that LUMO energy is highly correlated with the diol formation rate of oxetane hydrolysis by mEH for compounds that are structurally similar, while other quantum mechanical parameters are less predictive. The data suggest that the intrinsic reactivity determines the hydrolytic rate of oxetanes. This occurs when the orientations of the molecules in the mEH active site are similar. Predictions of mEH substrate metabolic rates using LUMO are most accurate when comparing subtle structural changes without drastic changes in MW and chemotype.

Conclusion:

The study suggests that LUMO energy can be used to rank-order oxetanes for their hydrolytic rate by mEH for structurally similar compounds. This finding enables the medicinal chemistry design to reduce mEH-mediated oxetane metabolism based on the calculated LUMO energy.

Unbound Brain-to-Plasma Partition Coefficient, Kp,uu,brain-a Game Changing Parameter for CNS Drug Discovery and Development

PURPOSE

More than 15 years have passed since the first description of the unbound brain-to-plasma partition coefficient (Kp,uu,brain) by Prof. Margareta Hammarlund-Udenaes, which was enabled by advancements in experimental methodologies including cerebral microdialysis. Since then, growing knowledge and data continue to support the notion that the unbound (free) concentration of a drug at the site of action, such as the brain, is the driving force for pharmacological responses. Towards this end, Kp,uu,brain is the key parameter to obtain unbound brain concentrations from unbound plasma concentrations.

METHODS

To understand the importance and impact of the Kp,uu,brain concept in contemporary drug discovery and development, a survey has been conducted amongst major pharmaceutical companies based in Europe and the USA. Here, we present the results from this survey which consisted of 47 questions addressing: 1) Background information of the companies, 2) Implementation, 3) Application areas, 4) Methodology, 5) Impact and 6) Future perspectives.

RESULTS AND CONCLUSIONS

From the responses, it is clear that the majority of the companies (93%) has established a common understanding across disciplines of the concept and utility of Kp,uu,brain as compared to other parameters related to brain exposure. Adoption of the Kp,uu,brain concept has been mainly driven by individual scientists advocating its application in the various companies rather than by a top-down approach. Remarkably, 79% of all responders describe the portfolio impact of Kp,uu,brain implementation in their companies as 'game-changing'. Although most companies (74%) consider the current toolbox for Kp,uu,brain assessment and its validation satisfactory for drug discovery and early development, areas of improvement and future research to better understand human brain pharmacokinetics/pharmacodynamics translation have been identified.

Early Drug-Induced Liver Injury (DILI) Risk Screening: "Free", as good as it gets

For all the promise of and need for clinical drug-induced liver injury (DILI) risk screening systems, demonstrating the predictive value of these systems versus readily available physicochemical properties and inherent dosing information has not been thoroughly evaluated. Therefore, we utilized a systematic approach to evaluate the predictive value of in vitro safety assays including Bile Salt Export Pump (BSEP) transporter inhibition and cytotoxicity in HepG2 and transformed human liver epithelial (THLE) along with physicochemical properties. We also evaluated the predictive value of in vitro ADME assays including hepatic partition coefficient (Kp) and its unbound counterpart since they provide insight on hepatic accumulation potential. The datasets comprised of 569 marketed drugs with FDA DILIrank annotation (Most vs Less/None), dose and physicochemical information, 384 drugs with Kp and plasma protein binding data, and 279 drugs with safety assay data. For each dataset and combination of input parameters, we developed random forest machine learning models and measured model performance using the receiver operator characteristic area-under-the-curve (ROC AUC). The median ROC AUC across the various data and parameters sets ranged from 0.67 to 0.77 with little evidence of additive predictivity when including safety or ADME assay data. Subsequent machine learning models consistently demonstrated daily dose, fraction sp3 or ionization, and cLogP/D inputs produced the best, simplest model for predicting clinical DILI risk with an ROC AUC of 0.75. This systematic framework should be used for future assay predictive value assessments and highlights the need for continued improvements to clinical DILI risk annotation.

A Randomized Open-Labeled Trial of Methotrexate as a Steroid-Sparing Agent for Patients With Generalized Myasthenia Gravis

Background and Purpose

Two clinical trials assessing the steroid-sparing effect of methotrexate (MTX) yielded conflicting results. Our objective was to investigate whether MTX would show a steroid-sparing effect in the treatment of generalized myasthenia gravis (MG) patients who fitted Myasthenia Gravis Foundation of America (MGFA) Class II and Class III.

Methods

We performed an 18-month prospective, randomized, open-labeled trial of prednisone combined with MTX 10 mg orally every week versus prednisone alone in 40 recently diagnosed MG patients of MGFA Class II and Class III between July 2014 and July 2018. The primary endpoint was the prednisone area under the dose–time curve (AUDTC) from months 3 to 18. Secondary endpoints included changes of the Quantitative Myasthenia Gravis Score (QMG), the Myasthenia Gravis Activity of Daily Living Score (MG-ADL), initial time of prednisone reduction, the median prednisone daily dose in each month, adverse events, and treatment failures in each group.

Results

Forty participants were included; among those, 5 individuals withdrew. A total of 35 participants completed 18 months of follow-up (18 in prednisone+MTX, 17 in prednisone group). Combined use of MTX reduced the month 3–18 prednisone AUDTC (prednisone+MTX 5,663.44 ± 1,678.08 mg, prednisone 6,683.94 ± 678.08 mg, p = 0.03, 95% confidence interval -1916.01 to -124.98). The initial times of prednisone reduction were 4.34 ± 1.44 months in the prednisone+MTX group and 5.56 ± 2.05 months in the prednisone group (p = 0.04, 95% CI -2.41 to -0.03). The median daily prednisone dose was significantly lower in the prednisone+MTX group at month 6 and months 9–18. No significant differences were found in QMG and MG-ADL scores between the two groups. No serious drug-related adverse events were observed in both groups.

Conclusions

This study provides evidence that MTX has the steroid-sparing ability in generalized MG patients of MGFA Class II and Class III.

Clinical Trial Registration

http://www.chictr.org.cn/showproj.aspx?proj=10563 identifier ChiCTR-IPR-15006081.

The CGG repeat expansion in RILPL1 is associated with oculopharyngodistal myopathy type 4

Recent studies indicate that CGG repeat expansions in LRP12, GIPC1, and NOTCH2NLC are associated with oculopharyngodistal myopathy (OPDM) types 1, 2, and 3, respectively. However, some clinicopathologically confirmed OPDM cases continue to have unknown genetic causes. Here, through a combination of long-read whole-genome sequencing (LRS), repeat-primed polymerase chain reaction (RP-PCR), and fluorescence amplicon length analysis PCR (AL-PCR), we found that a CGG repeat expansion in the 5' UTR of RILPL1 is associated with familial and simplex OPDM type 4 (OPDM4). The number of repeats ranged from 139 to 197. Methylation analysis indicates that the methylation levels in RILPL1 were unaltered in OPDM4 individuals. Analyses of muscle biopsies suggested that the expanded CGG repeat might be translated into a toxic poly-glycine protein that co-localizes with p62 in intranuclear inclusions. Moreover, analyses suggest that the toxic RNA gain-of-function effects also contributed to the pathogenesis of this disease. Intriguingly, all four types of OPDM have been found to be associated with the CGG repeat expansions located in 5' UTRs. This finding suggests that a common pathogenic mechanism, driven by the CGG repeat expansion, might underlie all cases of OPDM.

Tacrolimus treatment for relapsing-remitting chronic inflammatory demyelinating polyradiculoneuropathy: Two case reports

BACKGROUND

This study describes the efficacy of a tacrolimus treatment regimen used to treat two patients with relapsing-remitting chronic inflammatory demyelinating polyradiculoneuropathy (CIDP).

CASE SUMMARY

Two patients (17-year-old female and 27-year-old male) were enrolled in the current study and were followed up for 12 mo. The first patient was administered tacrolimus (2 mg/d) for 12 mo and prednisolone (40 mg/d) for six months. The second patient was administered tacrolimus (3 mg/d) for six months. Both patients were followed up for 12 mo and the degree of recurrent weakness or normalized motor function was monitored. In addition, nerve conduction studies and tacrolimus levels were recorded. Following tacrolimus treatment, both patients showed marked improvement in clinical outcomes. In the first patient, prednisolone treatment was successfully withdrawn after six months. Sensory as well as motor nerve conduction velocities showed evident recovery following treatment. However, conduction velocities did not completely return to normal, suggesting that electrophysiological recovery can be slower than clinical recovery.

CONCLUSION

Neither patient exhibited any adverse effects due to the tacrolimus therapy. Therefore, tacrolimus can be effective for the treatment of patients with steroid-resistant CIDP.

Clinical Predictors of Relapse in a Cohort of Steroid-Treated Patients With Well-Controlled Myasthenia Gravis

Objective

Despite the high efficiency of glucocorticoids (GCs), ~18–34% patients with myasthenia gravis (MG) may experience relapses of the disease. Here, we aim to identify clinical factors related to relapses during steroid tapering or after withdrawal in MG patients who were well-managed on steroid monotherapy.

Methods

We conducted a retrospective study on 125 MG patients from the Xuanwu Hospital MG Trial Database. Patients were treated with corticosteroids and achieved minimal manifestation status (MMS) or better. Patients were divided into steroid reduction subset (N = 74) and steroid withdrawal subset (N = 51). Clinical characteristics and therapeutic data were compared between patients with disease relapse and those who maintained clinical remission at the last follow-ups. Cox proportional hazards regression models were used to identify risk factors of relapse in each subset.

Results

Thirty-seven (29.6%) patients experienced relapses during the follow-up periods. Relapse during the steroid reduction was significantly associated with drug reducing duration (HR = 0.81, 95%CI 0.74–0.89, P &lt; 0.001). Risk of relapse was augmented if the drug reducing duration was &lt;11.5 months (HR 27.80, 95%CI 5.88–131.57, P &lt; 0.001). Among patients who discontinued the steroids, those with onset symptoms of bulbar weakness (adjusted HR 3.59, 95%CI 1.19–10.81, P = 0.023) were more likely to experience relapse.

Conclusion

Our study demonstrated that patients could benefit from prolonged steroid-reducing duration to prevent disease relapse. Patients with bulbar weakness at disease onset should be proposed to take long-term steroids or other immunosuppressants.

Missense Mutations of Codon 116 in the SOD1 Gene Cause Rapid Progressive Familial ALS and Predict Short Viability With PMA Phenotype

Amyotrophic lateral sclerosis (ALS) is the most common form of motor neuron disease, characterized by a great variety of both clinical presentations and genetic causes. Previous studies had identified two different missense mutations in SOD1 (p.R116C and p.R116G) causing familial ALS. In this study, we report a novel heterozygous missense mutation in the SOD1 gene (p.R116S) in a family with inherited ALS manifested as fast-deteriorating pure lower motor neuron symptoms. The patient displayed similar clinical picture and prognostic value to previous reported cases with different R116 substitution mutations. Modeling of all R116 substitutions in the resolved SOD1 protein structure revealed a shared mechanism with destroyed hydrogen bonds between R116 and other two residues, which might lead to protein unfolding and oligomer formation, ultimately conferring neurotoxicity.

Delayed Respiratory Insufficiency and Extramuscular Abnormalities in Selenoprotein N-Related Myopathies

Background: Selenoprotein N-related myopathies (SEPN1-RMs) are a subset of congenital myopathies caused by mutations of Selenoprotein N gene (SELENON or SEPN1). Clinical phenotype is considered as highly consistent and little attention has been given to the extramuscular abnormalities. Methods: We reported clinical, histopathological, and genetic features of four Chinese patients with SEPN1-RM and performed literature review on delayed respiratory insufficiency and extramuscular involvement. Results: A total of four patients exhibited both the typical and atypical clinical features of SEPN1-RM. The classical manifestations included axial and limb girdle weakness, spinal rigidity, scoliosis, respiratory insufficiency, and multiminicore morphological lesions. However, high interindividual variability was noticed on disease severity, especially the onset of respiratory involvement. Two adult patients postponed respiratory insufficiency to the third decade of life, while two juvenile patients manifested early hypoventilation with puberty exacerbation. As atypical features, extramuscular involvement of weight gain, subcutaneous adipose tissue accumulation, intellectual disability, and mild cardiac changes were observed. Molecular findings revealed three novel mutations of SELENON such as c.1286_1288 del CCT, c.1078_1086dupGGCTACATA, and c.785 G>C. Ten cases with delayed respiratory insufficiency were identified from previous publications. A total of 18 studies described extramuscular abnormalities including joint contractures, alterations of body mass index (BMI), mild cardiac changes, and insulin resistance. Intellectual impairment was extremely rare. Conclusion: SEPN1-RM should be considered as a differential diagnosis in adult patients with delayed respiratory involvement. Extramuscular involvement such as body composition alterations deserves more clinical attention. The novel mutations of SELENON widened the genetic spectrum of patients with SEPN1-RM.

An oral SARS-CoV-2 Mpro inhibitor clinical candidate for the treatment of COVID-19

[Figure: see text].

An update on the importance of plasma protein binding in drug discovery and development

Introduction: Plasma protein binding (PPB) remains a controversial topic in drug discovery and development. Fraction unbound (f_u) is a critical parameter that needs to be measured accurately, because it has significant impacts on the predictions of drug-drug interactions (DDI), estimations of therapeutic indices (TI), and developments of PK/PD relationships.  However, it is generally not advisable to change PPB through structural modifications, because PPB on its own has little relevance for in vivo efficacy.Areas covered: PPB fundamentals are discussed including the three main classes of drug binding proteins (i.e., albumin, alpha1-acid glycoprotein, and lipoproteins) and their physicochemical properties, in vivo half-life, and synthesis rate.  State-of-the-art methodologies for PPB are highlighted. Applications of PPB in drug discovery and development are presented.Expert opinion: PPB is an old topic in pharmacokinetics, but there are still many misconceptions. Improving the accuracy of PPB for highly bound compounds is an ongoing effort in the field with high priority. As the field continues to generate high quality data, the regulatory agencies will increase their confidence in our ability to accurately measure PPB of highly bound compounds, and experimental fu values below 0.01 will more likely be used for DDI predictions in the future.

Evolving approaches on measurements and applications of intracellular free drug concentration and Kpuu in drug discovery

Introduction: Intracellular-free drug concentration (C_u,cell) and unbound partition coefficient (Kp_uu) are two important parameters to develop pharmacokinetic and pharmacodynamic relationships, predict drug-drug interaction potentials and estimate therapeutic indices.Area covered: Methods on measurements of C_u,cell, Kp_uu, partition coefficient (Kp) and fraction unbound of cells (f_uc) are discussed. Advantages and limitations of several f_uc methods are reviewed. Applications highlighted here are bridging the potency gaps between biochemical and cell-based assays, in vitro hepatocyte assay to predict in vivo liver-to-plasma Kp_uu, the role of Kp_uu in prediction of hepatic clearance for enzyme- and transporter-mediated mechanisms using extended clearance equation, and structural attributes governing tissue Kp_uu.Expert opinion: C_u,cell and Kp_uu are of growing applications in drug discovery. Methods for measurements of these properties continue to evolve in order to achieve higher precision/accuracy and obtain more detailed information at the subcellular levels. Future directions of the field include the development of in vitro and in silico models to predict tissue Kp_uu, direct measurement of free drug concentration in subcellular organelles, and further investigations into the critical elements governing cell and tissue Kp_uu. Significant innovation is needed to advance this complex, but highly impactful and exciting area of science.

Framework and performance analysis of college English testing system based on data mining technology

At present, the college English testing system has become an indispensable system in many universities. However, the English test system is not highly humanized due to problems such as unreasonable framework structure. This paper combines data mining technology to build a college English test framework. The college English test system software based on data mining mainly realizes the computer program to automatically generate test papers, set the test time to automatically judge the test takers’ test results, and give out results on the spot. The test takers log in to complete the test through the test system software. The examination system software solves the functions of printing test papers, arranging invigilation classrooms, invigilating teachers, invigilating process, collecting test papers, scoring and analyzing test papers in traditional examinations. Finally, this paper analyzes the performance of this paper through experimental research. The research results show that the system constructed in this paper has certain practical effects.

CYP3A5*3 polymorphism and age affect tacrolimus blood trough concentration in myasthenia gravis patients

The study aims to identify clinical factors affecting tacrolimus blood trough concentration (C0) in myasthenia gravis (MG) patients and to optimize the initial dose of tacrolimus in MG treatment. A total of 103 MG patients participated in this study, and their clinical factors, medication regimens, C0 values and CYP3A5*3 polymorphisms were collected in detail. We used a linear mixed model to analyze the effect of multiple factors on the dosage-weighted C0 (C0:D) and performed subgroup analyses to investigate the consistency of correlations between influencing factors and the C0:D ratios. Among all factors, CYP3A5*3 polymorphism and age showed a strong positive correlation with C0:D ratios. The C0:D ratios (ng/ml·mg^-1) were higher for CYP3A5*3/*3 than for CYP3A5*1 (mean difference: 1.038, 95% confidence interval [CI]: 0.820-1.256, P-value <0.001), and for age in the range of 45-64 and ≥ 65 years than for age < 45 years (mean difference [95% CI] and P-value: 0.531[0.257-0.805] and P-value <0.001, 0.703 [0.377-1.029] and P-value <0.001, respectively). The C0:D ratios were not related to corticosteroid dosage, body weight, sex, hematocrit or the concomitant use of calcium channel blockers. The consistencies of the correlations between C0:D ratios and CYP3A5*3 polymorphism or age were confirmed by subgroup analyses. Thus, CYP3A5*3 polymorphism and age should be considered in optimizing the initial dose of tacrolimus for MG treatment.

VNTR2/VNTR3 genotype in the FCGRT gene is associated with reduced effectiveness of intravenous immunoglobulin in patients with myasthenia gravis

Background:

Intravenous immunoglobulin (IVIG) has been commonly used to treat myasthenia gravis exacerbation, but is still ineffective in nearly 30% of patients. A variable number of tandem repeat (VNTR) polymorphism in the FCGRT gene has been found to reduce the efficiency of IgG biologics. However, whether the polymorphism influences the efficacy of IVIG in generalized myasthenia gravis (MG) patients with exacerbations remains unknown.

Methods:

The distribution of VNTR genotypes was analyzed in 334 patients with MG. Varied VNTR alleles were determined by capillary electrophoresis and confirmed by Sanger sequencing. Information of endogenous IgG levels were collected in patients without previous immunotherapy (n = 26). Medical records of patients who received IVIG therapy were retrospectively analyzed for therapeutic outcomes of IVIG treatment (n = 61). Patients whose Activities of Daily Living scores decreased by 2 or more points on day 14 were considered responders to the treatment.

Results:

The VNTR3/3 and VNTR2/3 genotypes were detected in 96.7% (323/334) and 3.4% (11/334) patients, respectively. Patients with VNTR2/3 heterozygosity had lower endogenous IgG levels than those with VNTR3/3 homozygosity (9.81 ± 2.61 g/L versus 12.41 ± 2.45g/L, p = 0.016). The response rate of IVIG therapy was 78.7% (48/61). All responders and nine non-responders were VNTR3/3 homozygotes, whereas all the patients with VNTR2/3 genotypes were non-responders (n = 4). In patients who took IVIG treatments, endogenous IgG levels were significantly lower in non-responders compared with responders (12.93 ± 2.24 g/L versus 8.85 ± 2.69 g/L, p = 0.006), especially in VNTR2/3 heterozygotes (7.86 ± 1.78 g/L, p = 0.001).

Conclusion:

The VNTR2/3 genotype could influence endogenous IgG levels and serve as a predictive marker for poor responses to IVIG in MG patients.

Consideration of the Unbound Drug Concentration in Enzyme Kinetics

The study of enzyme kinetics in drug metabolism involves assessment of rates of metabolism and inhibitory potencies over a suitable concentration range. In all but the very simplest in vitro system, these drug concentrations can be influenced by a variety of nonspecific binding reservoirs that can reduce the available concentration to the enzyme system(s) under investigation. As a consequence, the apparent kinetic parameters, such as K_m or K_i, that are derived can deviate from the true values. There are a number of sources of these nonspecific binding depots or barriers, including membrane permeation and partitioning, plasma or serum protein binding, and incubational binding. In the latter case, this includes binding to the assay apparatus as well as biological depots, depending on the characteristics of the in vitro matrix being used. Given the wide array of subcellular, cellular, and recombinant enzyme systems utilized in drug metabolism, each of these has different components which can influence the free drug concentration. The physicochemical properties of the test compound are also paramount in determining the influential factors in any deviation between true and apparent kinetic behavior. This chapter describes the underlying mechanisms determining the free drug concentration in vitro and how these factors can be accounted for in drug metabolism studies, illustrated with case studies from the literature.

Evaluation of Prediction Accuracy for Volume of Distribution in Rat and Human Using In Vitro, In Vivo, PBPK and QSAR Methods

Volume of distribution at steady state (V_ss) is an important pharmacokinetic parameter of a drug candidate. In this study, V_ss prediction accuracy was evaluated by using: (1) seven methods for rat with 56 compounds, (2) four methods for human with 1276 compounds, and (3) four in vivo methods and three Kp (partition coefficient) scalar methods from scaling of three preclinical species with 125 compounds. The results showed that the global QSAR models outperformed the PBPK methods. Tissue fraction unbound (f_u,t) method with adipose and muscle also provided high V_ss prediction accuracy. Overall, the high performing methods for human V_ss prediction are the global QSAR models, Øie-Tozer and equivalency methods from scaling of preclinical species, as well as PBPK methods with Kp scalar from preclinical species. Certain input parameter ranges rendered PBPK models inaccurate due to mass balance issues. These were addressed using appropriate theoretical limit checks. Prediction accuracy of tissue Kp were also examined. The f_u,t method predicted Kp values more accurately than the PBPK methods for adipose, heart and muscle. All the methods overpredicted brain Kp and underpredicted liver Kp due to transporter effects. Successful V_ss prediction involves strategic integration of in silico, in vitro and in vivo approaches.

Favorable Effects of Tacrolimus Monotherapy on Myasthenia Gravis Patients

Background and Purpose: Tacrolimus (TAC) has been proven to be a rapid-acting, steroid-sparing agent for myasthenia gravis (MG) therapy. However, evidence related to the effectiveness of TAC alone is rare. Therefore, this study was performed to investigate the effect of TAC monotherapy in MG patients.

Methods: Forty-four MG patients who received TAC monotherapy were retrospectively analyzed. A mixed effect model was used to analyze improvements in MG-specific activities of daily living scale (MG-ADL), quantitative MG score (QMG) and MG-ADL subscores. Kaplan-Meier analysis was used to estimate the cumulative probability of minimal manifestations (MM) or better. Adverse events (AEs) were recorded for safety analyses.

Results: Of the patients receiving TAC monotherapy, MG-ADL scores were remarkably improved at 3, 6 and 12 months compared with scores at baseline (mean difference and 95% CIs: −3.29 [−4.94, −1.64], −3.97 [−5.67, −2.27], and −4.67 [−6.48, −2.85], respectively). QMG scores significantly decreased at 6 and 12 months, with mean differences and 95% CIs of −4.67(−6.88, −2.45) and −5.77 (−7.55, −4.00), respectively. Estimated median period to achieve “MM or better” was 5.0 (95% CIs, 2.8, 7.2) months. Ocular MG (OMG) and generalized MG (GMG) showed similar therapeutic effects in cumulative probabilities of “MM or better” (P-value = 0.764). A better response was observed in MG-ADL subscores for ptosis and bulbar symptoms. AEs occurred in 37.5% of patients and were generally mild and reversible.

Conclusions: TAC monotherapy is a promising option to rapidly alleviate all symptoms of MG, especially for ptosis and bulbar symptoms.

The Critical Role of Passive Permeability in Designing Successful Drugs

Passive permeability is a key property in drug disposition and delivery. It is critical for gastrointestinal absorption, brain penetration, renal reabsorption, defining clearance mechanisms and drug-drug interactions. Passive diffusion rate is translatable across tissues and animal species, while the extent of absorption is dependent on drug properties, as well as in vivo physiology/pathophysiology. Design principles have been developed to guide medicinal chemistry to enhance absorption, which combine the balance of aqueous solubility, permeability and the sometimes unfavorable compound characteristic demanded by the target. Permeability assays have been implemented that enable rapid development of structure-permeability relationships for absorption improvement. Future advances in assay development to reduce nonspecific binding and improve mass balance will enable more accurately measurement of passive permeability. Design principles that integrate potency, selectivity, passive permeability and other ADMET properties facilitate rapid advancement of successful drug candidates to patients.

Determinants of the repetitive-CMAP occurrence and therapy efficacy in slow-channel myasthenia

OBJECTIVE

To find determinants of the occurrence of repetitive compound muscle action potential (R-CMAP) and to assess the efficacy of channel blocker therapy in slow-channel congenital myasthenic syndrome (SCCMS).

METHODS

Neurologic examination, EMG study, laboratory test, muscle biopsy, and next-generation and Sanger sequencing; literature review of reported patients with SCCMS, including EMG, kinetics of mutant acetylcholine receptors (AChRs), and response to therapy; and simulation of the decay phase of endplate potential (EPP) were performed.

RESULTS

Three newly characterized and 57 reported patients with SCCMS with mutations of AChR subunits were included. In patients with R-CMAP, the length of channel opening bursts of mutant AChR was increased 8.68 ± 2.82 (mean ± SD)-fold compared to wild-type; in patients without R-CMAP, the length was increased 3.84 ± 0.65-fold (95% confidence interval 3.18-6.50, p = 0.000014). The EPP amplitude after refractory period of action potential in muscle fiber is above the threshold in patients with R-CMAP but below the threshold in patients without R-CMAP. In patients with good results from channel blocker therapy, treatment was initiated 11.60 ± 5.17 years after onset of symptoms; in patients with no to moderate benefit from channel blocker therapy, treatment was initiated 30.70 ± 12.72 years after onset (95% confidence interval -28.57 to -9.63, p = 0.00089).

CONCLUSIONS

In SCCMS, the R-CMAP occurrence is related to the extent of prolongation of the opening episodes of mutant AChR channel. Channel blocker treatment is more effective the sooner it is started after the onset of symptoms.

CLASSIFICATION OF EVIDENCE

This study provides Class IV evidence that channel blocker therapy in patients with SCCMS improves symptoms.

A novel fast-channel myasthenia caused by mutation in β subunit of AChR reveals subunit-specific contribution of the intracellular M1-M2 linker to channel gating

Genetic variants causing the fast-channel congenital myasthenic syndrome (CMS) have been identified in the α, δ, and ε but not the β subunit of acetylcholine receptor (AChR). A 16-year-old girl with severe myasthenia had low-amplitude and fast-decaying miniature endplate potentials. Mutation analysis revealed two heteroallelic variants in CHRNB1 encoding the AChR β subunit: a novel c.812C>T (p.P248L) variant in M1-M2 linker (p.P271L in HGVS nomenclature), and a ~430 bp deletion causing loss of exon 8 leading to frame-shift and a premature stop codon (p.G251Dfs*21). P248 is conserved in all β subunits of different species, but not in other AChR subunits. Measurements of radio-labeled α-bungarotoxin binding show that βP248L reduces AChR expression to 60% of wild-type. Patch clamp recordings of ACh-elicited single channel currents demonstrate that βP248L shortens channel opening bursts from 3.3 ms to 1.2 ms, and kinetic analyses predict that the decay of the synaptic response is accelerated 2.4-fold due to reduced probability of channel reopening. Substituting βP248 with threonine, alanine or glycine reduces the burst duration to 2.3, 1.7, and 1.5 ms, respectively. In non-β subunits, substituting leucine for residues corresponding to βP248 prolongs the burst duration to 4.5 ms in the α subunit, shortens it to 2.2 ms in the δ subunit, and has no effect in the ε subunit. Conversely, substituting proline for residues corresponding to βP248 prolongs the burst duration to 8.7 ms in the α subunit, to 4.6 ms in the δ subunit, but has no effect in the ε subunit. Thus, this fast channel CMS is caused by the dual defects of βP248L in reducing expression of the mutant receptor and accelerating the decay of the synaptic response. The results also reveal subunit-specific contributions of the M1-M2 linker to the durations of channel opening bursts.

Dose Predictions for Drug Design

The efficacious dose of a drug is perhaps the most holistic metric reflecting its therapeutic potential. Dose is predicted at many stages in drug discovery and development. Prior to the 1990s, dose prediction was limited to the drug "working" at a reasonable dose and dose regimen in an animal model. Through the early 2000s, dose predictions were generated at candidate nomination and then refined during clinical development. Currently, dose predictions can be made early in drug discovery to enable drug design. Dose predictions at this stage can identify critical drug properties for a viable dose regimen and provide clinically relevant context to lead optimization. In this paper, we give an overview of the opportunities and challenges associated with dose prediction for drug design. A number of general considerations, approaches, and case examples are discussed.