Modeling and Simulation of the NMDA Receptor at Coarse-Grained and Atomistic Levels.

N-Methyl-D-aspartate (NMDA) receptors are glutamate-gated excitatory channels that play essential roles in brain functions. While high-resolution structures were solved for an allosterically inhibited form of functional NMDA receptor, other key functional states (particularly the active open-channel state) have not yet been resolved at atomic resolutions. To decrypt the molecular mechanism of the NMDA receptor activation, structural modeling and simulation are instrumental in providing detailed information about the dynamics and energetics of the receptor in various functional states. In this chapter, we describe coarse-grained modeling of the NMDA receptor using an elastic network model and related modeling/analysis tools (e.g., normal mode analysis, flexibility and hotspot analysis, cryo-EM flexible fitting, and transition pathway modeling) based on available structures. Additionally, we show how to build an atomistic model of the active-state receptor with targeted molecular dynamics (MD) simulation and explore its energetics and dynamics with conventional MD simulation. Taken together, these modeling and simulation can offer rich structural and dynamic information which will guide experimental studies of the activation of this key receptor.

Interfacial interactions between spider silk protein and cellulose studied by molecular dynamics simulation.

CONTEXT: Due to their excellent biocompatibility and degradability, cellulose/spider silk protein composites hold a significant value in biomedical applications such as tissue engineering, drug delivery, and medical dressings. The interfacial interactions between cellulose and spider silk protein affect the properties of the composite. Therefore, it is important to understand the interfacial interactions between spider silk protein and cellulose to guide the design and optimization of composites. The study of the adsorption of protein on specific surfaces of cellulose crystal can be very complex using experimental methods. Molecular dynamics simulations allow the exploration of various physical and chemical changes at the atomic level of the material and enable an atomic description of the interactions between cellulose crystal planes and spider silk protein. In this study, molecular dynamics simulations were employed to investigate the interfacial interactions between spider silk protein (NTD) and cellulose surfaces. Findings of RMSD, RMSF, and secondary structure showed that the structure of NTD proteins remained unchanged during the adsorption process. Cellulose contact numbers and hydrogen bonding trends on different crystalline surfaces suggest that van der Waals forces and hydrogen bonding interactions drive the binding of proteins to cellulose. These findings reveal the interaction between cellulose and protein at the molecular level and provide theoretical guidance for the design and synthesis of cellulose/spider silk protein composites. METHODS: MD simulations were all performed using the GROMACS-5.1 software package and run with CHARMM36 carbohydrate force field. Molecular dynamics simulations were performed for 500 ns for the simulated system.

Unveiling the Molecular Interactions Between Human Transferrin and Limonene: Natural Compounds in Alzheimer's Disease Therapeutics.

Background: Neurodegeneration is a term describing an irreversible process of neuronal damage. In recent decades, research efforts have been directed towards deepening our knowledge of numerous neurodegenerative disorders, with a particular focus on conditions such as Alzheimer's disease (AD). Human transferrin (htf) is a key player in maintaining iron homeostasis within brain cells. Any disturbance in this equilibrium gives rise to the emergence of neurodegenerative diseases and associated pathologies, particularly AD. Limonene, a natural compound found in citrus fruits and various plants, has shown potential neuroprotective properties. Objective: In this study, our goal was to unravel the binding of limonene with htf, with the intention of comprehending the interaction mechanism of limonene with htf. Methods: Binding was scrutinized using fluorescence quenching and UV-Vis spectroscopic analyses. The binding mechanism of limonene was further investigated at the atomic level through molecular docking and extensive 200 ns molecular dynamic simulation (MD) studies. Results: Molecular docking uncovered that limonene interacted extensively with the deep cavity located within the htf binding pocket. MD results indicated that binding of limonene to htf did not induce substantial structural alterations, ultimately forming stable complex. The findings from fluorescence binding indicated a pronounced interaction between limonene and htf, limonene binds to htf with a binding constant (K) of 0.1×105 M-1. UV spectroscopy also advocated stable htf-limonene complex formation. Conclusions: The study deciphered the binding mechanism of limonene with htf, providing a platform to use limonene in AD therapeutics in context of iron homeostasis.

Treatment of flexibility of protein backbone in simulations of protein-ligand interactions using steered molecular dynamics.

To ensure that an external force can break the interaction between a protein and a ligand, the steered molecular dynamics simulation requires a harmonic restrained potential applied to the protein backbone. A usual practice is that all or a certain number of protein's heavy atoms or Cα atoms are fixed, being restrained by a small force. This present study reveals that while fixing both either all heavy atoms and or all Cα atoms is not a good approach, while fixing a too small number of few atoms sometimes cannot prevent the protein from rotating under the influence of the bulk water layer, and the pulled molecule may smack into the wall of the active site. We found that restraining the Cα atoms under certain conditions is more relevant. Thus, we would propose an alternative solution in which only the Cα atoms of the protein at a distance larger than 1.2 nm from the ligand are restrained. A more flexible, but not too flexible, protein will be expected to lead to a more natural release of the ligand.

Conformations of a highly expressed Z19 α-zein studied with AlphaFold2 and MD simulations.

α-zeins are amphiphilic maize seed storage proteins with material properties suitable for a multitude of applications e.g., in renewable plastics, foods, therapeutics and additive manufacturing (3D-printing). To exploit their full potential, molecular-level insights are essential. The difficulties in experimental atomic-resolution characterization of α-zeins have resulted in a diversity of published molecular models. However, deep-learning α-zein models are largely unexplored. Therefore, this work studies an AlphaFold2 (AF2) model of a highly expressed α-zein using molecular dynamics (MD) simulations. The sequence of the α-zein cZ19C2 gave a loosely packed AF2 model with 7 α-helical segments connected by turns/loops. Compact tertiary structure was limited to a C-terminal bundle of three α-helices, each showing notable agreement with a published consensus sequence. Aiming to chart possible α-zein conformations in practically relevant solvents, rather than the native solid-state, the AF2 model was subjected to MD simulations in water/ethanol mixtures with varying ethanol concentrations. Despite giving structurally diverse endpoints, the simulations showed several patterns: In water and low ethanol concentrations, the model rapidly formed compact globular structures, largely preserving the C-terminal bundle. At ≥ 50 mol% ethanol, extended conformations prevailed, consistent with previous SAXS studies. Tertiary structure was partially stabilized in water and low ethanol concentrations, but was disrupted in ≥ 50 mol% ethanol. Aggregated results indicated minor increases in helicity with ethanol concentration. ß-sheet content was consistently low (∼1%) across all conditions. Beyond structural dynamics, the rapid formation of branched α-zein aggregates in aqueous environments was highlighted. Furthermore, aqueous simulations revealed favorable interactions between the protein and the crosslinking agent glycidyl methacrylate (GMA). The proximity of GMA epoxide carbons and side chain hydroxyl oxygens simultaneously suggested accessible reactive sites in compact α-zein conformations and pre-reaction geometries for methacrylation. The findings may assist in expanding the applications of these technologically significant proteins, e.g., by guiding chemical modifications.

Elucidating Daptomycin's Antibacterial Efficacy: Insights into the Tripartite Complex with Lipid II and Phospholipids in Bacterial Septum Membrane.

This study elucidated the mechanism of formation of a tripartite complex containing daptomycin (Dap), lipid II, and phospholipid phosphatidylglycerol in the bacterial septum membrane, which was previously reported as the cause of the antibacterial action of Dap against gram-positive bacteria via molecular dynamics and enhanced sampling methods. Others have suggested that this transient complex ushers in the inhibition of cell wall synthesis by obstructing the downstream polymerization and cross-linking processes involving lipid II, which is absent in the presence of cardiolipin lipid in the membrane. In this work, we observed that the complex was stabilized by Ca2+-mediated electrostatic interactions between Dap and lipid head groups, hydrophobic interaction, hydrogen bonds, and salt bridges between the lipopeptide and lipids and was associated with Dap concentration-dependent membrane depolarization, thinning of the bilayer, and increased lipid tail disorder. Residues Orn6 and Kyn13, along with the DXDG motif, made simultaneous contact with constituent lipids, hence playing a crucial role in the formation of the complex. Incorporating cardiolipin into the membrane model led to its competitively displacing lipid II away from the Dap, reducing the lifetime of the complex and the nonexistence of lipid tail disorder and membrane depolarization. No evidence of water permeation inside the membrane hydrophobic interior was noted in all of the systems studied. Additionally, it was shown that using hydrophobic contacts between Dap and lipids as collective variables for enhanced sampling gave rise to a free energy barrier for the translocation of the lipopeptide. A better understanding of Dap's antibacterial mechanism, as studied through this work, will help develop lipopeptide-based antibiotics for rising Dap-resistant bacteria.

Protein Interface Regulating the Inserting Process of Imidazole Ionic Liquids into the Cell Membrane.

Ionic liquids (ILs) have shown promising potential in membrane protein extraction; however, the underlying mechanism remains unclear. Herein, we employed GPU-accelerated molecular dynamics (MD) simulations to investigate the dynamic insertion process of ILs into cell membranes containing membrane proteins. Our findings reveal that ILs spontaneously insert into the membrane, and the presence of membrane proteins significantly decelerates the rate of IL insertion into the membrane. Specifically, the relationship between the insertion rate and inserting free energy exhibits non-monotonic changes, which can be attributed to interfacial effects. The protein-water interface acts as trap for free ions and ionic clusters, while free ions preferentially insert into the membrane from the protein-lipid interface, which limits the insertion rate due to its narrowness. Thus, the insertion rate is governed by a combination of the free energy and interfacial effects. These findings provide valuable insights into the interfacial effects of protein-lipid bilayers and have implications for various biochemical-related applications.

Exploring the mechanism of Jinlida granules against type 2 diabetes mellitus by an integrative pharmacology strategy.

Jinlida granule (JLD) is a Traditional Chinese Medicine (TCM) formula used for the treatment of type 2 diabetes mellitus (T2DM). However, the mechanism of JLD treatment for T2DM is not fully revealed. In this study, we explored the mechanism of JLD against T2DM by an integrative pharmacology strategy. Active components and corresponding targets were retrieved from Traditional Chinese Medicine System Pharmacology (TCMSP), SwissADME and Bioinformatics Analysis Tool for Molecular Mechanisms of Traditional Chinese Medicine Database (BATMAN-TCM) database. T2DM-related targets were obtained from Drugbank and Genecards databases. The protein-protein interaction (PPI) network was constructed and analyzed with STRING (Search Toll for the Retrieval of Interacting Genes/proteins) and Cytoscape to get the key targets. Then, Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genomes (KEGG) enrichment analyses were performed with the Database for Annotation, Visualization and Integrated Discovery (DAVID). Lastly, the binding capacities and reliability between potential active components and the targets were verified with molecular docking and molecular dynamics simulation. In total, 185 active components and 337 targets of JLD were obtained. 317 targets overlapped with T2DM-related targets. RAC-alpha serine/threonine-protein kinase (AKT1), tumor necrosis factor (TNF), interleukin-6 (IL-6), cellular tumor antigen p53 (TP53), prostaglandin G/H synthase 2 (PTGS2), Caspase-3 (CASP3) and signal transducer and activator of transcription 3 (STAT3) were identified as seven key targets by the topological analysis of the PPI network. GO and KEGG enrichment analyses showed that the effects were primarily associated with gene expression, signal transduction, apoptosis and inflammation. The pathways were mainly enriched in PI3K-AKT signaling pathway and AGE-RAGE signaling pathway in diabetic complications. Molecular docking and molecular dynamics simulation verified the good binding affinity between the key components and targets. The predicted results may provide a theoretical basis for drug screening of JLD and a new insight for the therapeutic effect of JLD on T2DM.

Immunoinformatics design of a structural proteins driven multi-epitope candidate vaccine against different SARS-CoV-2 variants based on fynomer.

The ideal vaccines for combating diseases that may emerge in the future require more than simply inactivating a few pathogenic strains. This study aims to provide a peptide-based multi-epitope vaccine effective against various severe acute respiratory syndrome coronavirus 2 strains. To design the vaccine, a library of peptides from the spike, nucleocapsid, membrane, and envelope structural proteins of various strains was prepared. Then, the final vaccine structure was optimized using the fully protected epitopes and the fynomer scaffold. Using bioinformatics tools, the antigenicity, allergenicity, toxicity, physicochemical properties, population coverage, and secondary and three-dimensional structures of the vaccine candidate were evaluated. The bioinformatic analyses confirmed the high quality of the vaccine. According to further investigations, this structure is similar to native protein and there is a stable and strong interaction between vaccine and receptors. Based on molecular dynamics simulation, structural compactness and stability in binding were also observed. In addition, the immune simulation showed that the vaccine can stimulate immune responses similar to real conditions. Finally, codon optimization and in silico cloning confirmed efficient expression in Escherichia coli. In conclusion, the fynomer-based vaccine can be considered as a new style in designing and updating vaccines to protect against coronavirus disease.

ESAT-6 undergoes self-association at phagosomal pH and an ESAT-6-specific nanobody restricts M. tuberculosis growth in macrophages.

Mycobacterium tuberculosis (Mtb) is known to survive within macrophages by compromising the integrity of the phagosomal compartment in which it resides. This activity primarily relies on the ESX-1 secretion system, predominantly involving the protein duo ESAT-6 and CFP-10. CFP-10 likely acts as a chaperone, while ESAT-6 likely disrupts phagosomal membrane stability via a largely unknown mechanism. we employ a series of biochemical analyses, protein modeling techniques, and a novel ESAT-6-specific nanobody to gain insight into the ESAT-6's mode of action. First, we measure the binding kinetics of the tight 1:1 complex formed by ESAT-6 and CFP-10 at neutral pH. Subsequently, we demonstrate a rapid self-association of ESAT-6 into large complexes under acidic conditions, leading to the identification of a stable tetrameric ESAT-6 species. Using molecular dynamics simulations, we pinpoint the most probable interaction interface. Furthermore, we show that cytoplasmic expression of an anti-ESAT-6 nanobody blocks Mtb replication, thereby underlining the pivotal role of ESAT-6 in intracellular survival. Together, these data suggest that ESAT-6 acts by a pH-dependent mechanism to establish two-way communication between the cytoplasm and the Mtb-containing phagosome.

Transient interactions modulate the affinity of NF-κB transcription factors for DNA.

The dimeric nuclear factor kappa B (NF-κB) transcription factors (TFs) regulate gene expression by binding to a variety of κB DNA elements with conserved G:C-rich flanking sequences enclosing a degenerate central region. Toward defining mechanistic principles of affinity regulated by degeneracy, we observed an unusual dependence of the affinity of RelA on the identity of the central base pair, which appears to be noncontacted in the complex crystal structures. The affinity of κB sites with A or T at the central position is ~10-fold higher than with G or C. The crystal structures of neither the complexes nor the free κB DNAs could explain the differences in affinity. Interestingly, differential dynamics of several residues were revealed in molecular dynamics simulation studies, where simulation replicates totaling 148 µs were performed on NF-κB:DNA complexes and free κB DNAs. Notably, Arg187 and Arg124 exhibited selectivity in transient interactions that orchestrated a complex interplay among several DNA-interacting residues in the central region. Binding and simulation studies with mutants supported these observations of transient interactions dictating specificity. In combination with published reports, this work provides insights into the nuanced mechanisms governing the discriminatory binding of NF-κB family TFs to κB DNA elements and sheds light on cancer pathogenesis of cRel, a close homolog of RelA.

Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations.

Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.

Identification of a functional missense variant in the matrix metallopeptidase 10 (MMP10) gene in two families with premature myocardial infarction.

A positive family history is a major independent risk factor for atherosclerosis, and genetic variation is an important aspect of cardiovascular disease research. We identified a heterozygous missense variant p.L245P in the MMP10 gene in two families with premature myocardial infarction using whole-exome sequencing. The aim of this study was to investigate the consequences of this variant using in-silico and functional in-vitro assays. Molecular dynamics simulations were used to analyze protein interactions, calculate free binding energy, and measure the volume of the substrate-binding cleft of MMP10-TIMP1 models. The p.L245P variant showed an altered protein surface, different intra- and intermolecular interactions of MMP10-TIMP1, a lower total free binding energy between MMP10-TIMP1, and a volume-minimized substrate-binding cleft of MMP10 compared to the wild-type. For the functional assays, human THP-1 cells were transfected with plasmids containing MMP10 cDNA carrying the p.L245P and wild-type variant and differentiated into macrophages. Macrophage adhesion and migration assays were then conducted, and pro-inflammatory chemokine levels were evaluated. The p.L245P variant led to macrophages that were more adherent, less migratory, and secreted higher levels of the pro-inflammatory chemokines CXCL1 and CXCL8 than wild-type macrophages. Thus, the p.L245P variant in MMP10 may influence the pathogenesis of atherosclerosis in families with premature myocardial infarction by altering protein - protein interactions, macrophage adhesion and migration, and expression of pro-inflammatory chemokines, which may increase plaque rupture. These results could contribute to the development of selective MMP10 inhibitors and reduce the risk of atherosclerosis in families with a history of premature myocardial infarction.

Temperature-dependent iron motion in extremophile rubredoxins - no need for 'corresponding states'.

Extremophile organisms are known that can metabolize at temperatures down to - 25 °C (psychrophiles) and up to 122 °C (hyperthermophiles). Understanding viability under extreme conditions is relevant for human health, biotechnological applications, and our search for life elsewhere in the universe. Information about the stability and dynamics of proteins under environmental extremes is an important factor in this regard. Here we compare the dynamics of small Fe-S proteins - rubredoxins - from psychrophilic and hyperthermophilic microorganisms, using three different nuclear techniques as well as molecular dynamics calculations to quantify motion at the Fe site. The theory of 'corresponding states' posits that homologous proteins from different extremophiles have comparable flexibilities at the optimum growth temperatures of their respective organisms. Although 'corresponding states' would predict greater flexibility for rubredoxins that operate at low temperatures, we find that from 4 to 300 K, the dynamics of the Fe sites in these homologous proteins are essentially equivalent.

Transcriptional and structural analysis of non-specific lipid transfer proteins modulated by fungal endophytes in Antarctic plants under drought.

Lipid transfer proteins (LTPs) play crucial roles in various biological processes in plants, such as pollen tube adhesion, phospholipid transfer, cuticle synthesis, and response to abiotic stress. While a few members of the non-specific LTPs (nsLTPs) have been identified, their structural characteristics remain largely unexplored. Given the observed improvement in the performance of Antarctic plants facing water deficit when associated with fungal endophytes, this study aimed to assess the role of these symbiotic organisms in the transcriptional modulation of putative nsLTPs. The study focused on identifying and characterizing two nsLTP in the Antarctic plant Colobanthus quitensis that exhibit responsiveness to drought stress. Furthermore, we investigated the influence of Antarctic endophytic fungi on the expression profiles of these nsLTPs, as these fungi have been known to enhance plant physiological and biochemical performance under water deficit conditions. Through 3D modeling, docking, and molecular dynamics simulations with different substrates, the conducted structural and ligand-protein interaction analyses showed that differentially expressed nsLTPs displayed the ability to interact with various ligands, with a higher affinity towards palmitoyl-CoA. Overall, our findings suggest a regulatory mechanism for the expression of these two nsLTPs in Colobanthus quitensis under drought stress, further modulated by the presence of endophytic fungi.

Screening approach by a combination of computational and in vitro experiments: identification of fluvastatin sodium as a potential SIRT2 inhibitor.

BACKGROUND: Sirtuins (SIRTs) are NAD+-dependent deacetylases that play various roles in numerous pathophysiological processes, holding promise as therapeutic targets worthy of further investigation. Among them, the SIRT2 subtype is closely associated with tumorigenesis and malignancies. Dysregulation of SIRT2 activation can regulate the expression levels of related genes in cancer cells, leading to tumor occurrence and metastasis. METHODS: In this study, we used computer simulations to screen for novel SIRT2 inhibitors from the FDA database, based on which 10 compounds with high docking scores and good interactions were selected for in vitro anti-pancreatic cancer metastasis testing and enzyme binding inhibition experiments. The results showed that fluvastatin sodium may possess inhibitory activity against SIRT2. Subsequently, fluvastatin sodium was subjected to molecular docking experiments with various SIRT isoforms, and the combined results from Western blotting experiments indicated its potential as a SIRT2 inhibitor. Next, molecular docking, molecular dynamics (MD) simulations, and binding free energy calculations were performed, revealing the binding mode of fluvastatin sodium at the SIRT2 active site, further validating the stability and interaction of the ligand-protein complex under physiological conditions. RESULTS: Overall, this study provides a systematic virtual screening workflow for the discovery of SIRT2 activity inhibitors, identifies the potential inhibitory effect of fluvastatin sodium as a lead compound on SIRT2, and opens up a new direction for developing highly active and selectively targeted SIRT2 inhibitors.

Virtual Screening of Small Molecules Targeting BCL2 with Machine Learning, Molecular Docking, and MD Simulation.

This study aimed to identify potential BCL-2 small molecule inhibitors using deep neural networks (DNN) and random forest (RF), algorithms as well as molecular docking and molecular dynamics (MD) simulations to screen a library of small molecules. The RF model classified 61% (2355/3867) of molecules as 'Active'. Further analysis through molecular docking with Vina identified CHEMBL3940231, CHEMBL3938023, and CHEMBL3947358 as top-scored small molecules with docking scores of -11, -10.9, and 10.8 kcal/mol, respectively. MD simulations validated these compounds' stability and binding affinity to the BCL2 protein.

SNP-Associated Substitutions of Amino Acid Residues in the dNTP Selection Subdomain Decrease Polß Polymerase Activity.

In the cell, DNA polymerase ß (Polß) is involved in many processes aimed at maintaining genome stability and is considered the main repair DNA polymerase participating in base excision repair (BER). Polß can fill DNA gaps formed by other DNA repair enzymes. Single-nucleotide polymorphisms (SNPs) in the POLB gene can affect the enzymatic properties of the resulting protein, owing to possible amino acid substitutions. For many SNP-associated Polß variants, an association with cancer, owing to changes in polymerase activity and fidelity, has been shown. In this work, kinetic analyses and molecular dynamics simulations were used to examine the activity of naturally occurring polymorphic variants G274R, G290C, and R333W. Previously, the amino acid substitutions at these positions have been found in various types of tumors, implying a specific role of Gly-274, Gly-290, and Arg-333 in Polß functioning. All three polymorphic variants had reduced polymerase activity. Two substitutions-G274R and R333W-led to the almost complete disappearance of gap-filling and primer elongation activities, a decrease in the deoxynucleotide triphosphate-binding ability, and a lower polymerization constant, due to alterations of local contacts near the replaced amino acid residues. Thus, variants G274R, G290C, and R333W may be implicated in an elevated level of unrepaired DNA damage.

Computational Estimation of Residues Involving Resistance to the SARS-CoV-2 Main Protease Inhibitor Ensitrelvir Based on Virtual Alanine Scan of the Active Site.

Ensitrelvir is a noncovalent inhibitor of the main protease (Mpro) of severe acute respiratory syndrome coronavirus 2. Acquisition of drug resistance in virus-derived proteins is a serious therapeutic concern, and drug resistance occurs due to amino acid mutations. In this study, we computationally constructed 24 mutants, in which one residue around the active site was replaced with alanine and performed molecular dynamics simulations to the complex of Mpro and ensitrelvir to predict the residues involved in drug resistance. We evaluated the changes in the entire protein structure and ligand configuration in each of these mutants and estimated which residues were involved in ensitrelvir recognition. This method is called a virtual alanine scan. In nine mutants (S1A, T26A, H41A, M49A, L141A, H163A, E166A, V186A, and R188A), although the entire protein structure and catalytic dyad (cysteine (Cys)145 and histidine (His)41) were not significantly moved, the ensitrelvir configuration changed. Thus, it is considered that these mutants did not recognize ensitrelvir while maintaining Mpro enzymatic activities, and Ser1, Thr26, His41, Met49, Leu141, His163, Glu166, Val186, and Arg188 may be related to ensitrelvir resistance. The ligand shift noted in M49A was similar to that observed in M49I, which has been shown to be experimentally ensitrelvir resistant. These findings suggest that our research approach can predict mutations that incite drug resistance.

New tetrahydroisoquinoline-4-carbonitrile derivatives as potent agents against cyclin-dependent kinases, crystal structures, and computational studies.

The synthesis of two new hexahydroisoquinoline-4-carbonitrile derivatives (3a and 3b) is reported along with spectroscopic data and their crystal structures. In compound 3a, the intramolecular O-H···O hydrogen bond constraints the acetyl and hydroxyl groups to be syn. In the crystal, inversion dimers are generated by C-H···O hydrogen bonds and are connected into layers parallel to (10-1) by additional C-H···O hydrogen bonds. The layers are stacked with Cl···S contacts 0.17 Å less than the sum of the respective van der Waals radii. The conformation of the compound 3b is partially determined by the intramolecular O-H···O hydrogen bond. A puckering analysis of the tetrahydroisoquinoline unit was performed. In the crystal, O-H···O and C-H···O hydrogen bonds together with C-H···π(ring) interactions form layers parallel to (01-1) which pack with normal van der Waals interactions. To understand the binding efficiency and stability of the title molecules, molecular docking, and 100 ns dynamic simulation analyses were performed with CDK5A1. To rationalize their structure-activity relationship(s), a DFT study at the B3LYP/6-311++G** theoretical level was also done. The 3D Hirshfled surfaces were also taken to investigate the crystal packings of both compounds. In addition, their ADMET properties were explored.Communicated by Ramaswamy H. Sarma.