Investigations of albumin-insulin detemir complexes using molecular dynamics simulations and free energy calculations

Insulin detemir is a lipidated insulin analogue that obtains a half-life extension by oligomerization and reversible binding to human serum albumin. In the present study, the complex between a detemir hexamer and albumin is investigated by an integrative approach combining molecular dynamics (MD) simulations, molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) free energy calculations, and dynamic light scattering (DLS) experiments. Recent reported small-angle X-ray scattering data could not unambiguously resolve the exact binding site of detemir on albumin. We therefore applied MD simulations to deduce the binding site and key protein–protein interactions. MD simulations were started from initial complex structures based on the SAXS models, and free energies of binding were estimated from the simulations by using the MM-PBSA approach for the different binding positions. The results suggest that the overlapping FA3–FA4 binding site (named FA4) is the most favorable site with a calculated free energy of binding of −28 ± 6 kcal/mol and a good fit to the reported SAXS data throughout the simulations. Multiple salt bridges, hydrogen bonds, and favorable van der Waals interactions are observed in the binding interface that promote complexation. The binding to FA4 is further supported by DLS competition experiments with the prototypical FA4 ligand, ibuprofen, showing displacement of detemir by ibuprofen. This study provides information on albumin–detemir binding on a molecular level, which could be utilized in a rational design of future lipidated albumin-binding peptides

Substitutions in SARS-CoV-2 Mpro Selected by Protease Inhibitor Boceprevir Confer Resistance to Nirmatrelvir

Nirmatrelvir, which targets the SARS-CoV-2 main protease (Mpro), is the first-in-line drug for prevention and treatment of severe COVID-19, and additional Mpro inhibitors are in development. However, the risk of resistance development threatens the future efficacy of such direct-acting antivirals. To gain knowledge on viral correlates of resistance to Mpro inhibitors, we selected resistant SARS-CoV-2 under treatment with the nirmatrelvir-related protease inhibitor boceprevir. SARS-CoV-2 selected during five escape experiments in VeroE6 cells showed cross-resistance to nirmatrelvir with up to 7.3-fold increased half-maximal effective concentration compared to original SARS-CoV-2, determined in concentration–response experiments. Sequence analysis revealed that escape viruses harbored Mpro substitutions L50F and A173V. For reverse genetic studies, these substitutions were introduced into a cell-culture-infectious SARS-CoV-2 clone. Infectivity titration and analysis of genetic stability of cell-culture-derived engineered SARS-CoV-2 mutants showed that L50F rescued the fitness cost conferred by A173V. In the concentration–response experiments, A173V was the main driver of resistance to boceprevir and nirmatrelvir. Structural analysis of Mpro suggested that A173V can cause resistance by making boceprevir and nirmatrelvir binding less favorable. This study contributes to a comprehensive overview of the resistance profile of the first-in-line COVID-19 treatment nirmatrelvir and can thus inform population monitoring and contribute to pandemic preparedness

Nirmatrelvir-resistant SARS-CoV-2 variants with high fitness in an infectious cell culture system

The oral protease inhibitor nirmatrelvir is of key importance for prevention of severe coronavirus disease 2019 (COVID-19). To facilitate resistance monitoring, we studied severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) escape from nirmatrelvir in cell culture. Resistant variants harbored combinations of substitutions in the SARS-CoV-2 main protease (Mpro). Reverse genetics revealed that E166V and L50F + E166V conferred high resistance in infectious culture, replicon, and Mpro systems. While L50F, E166V, and L50F + E166V decreased replication and Mpro activity, L50F and L50F + E166V variants had high fitness in the infectious system. Naturally occurring L50F compensated for fitness cost of E166V and promoted viral escape. Molecular dynamics simulations revealed that E166V and L50F + E166V weakened nirmatrelvir-Mpro binding. Polymerase inhibitor remdesivir and monoclonal antibody bebtelovimab retained activity against nirmatrelvir-resistant variants, and combination with nirmatrelvir enhanced treatment efficacy compared to individual compounds. These findings have implications for monitoring and ensuring treatments with efficacy against SARS-CoV-2 and emerging sarbecoviruses