CHARMM at 45: Enhancements in Accessibility, Functionality, and Speed.

Since its inception nearly a half century ago, CHARMM has been playing a central role in computational biochemistry and biophysics. Commensurate with the developments in experimental research and advances in computer hardware, the range of methods and applicability of CHARMM have also grown. This review summarizes major developments that occurred after 2009 when the last review of CHARMM was published. They include the following: new faster simulation engines, accessible user interfaces for convenient workflows, and a vast array of simulation and analysis methods that encompass quantum mechanical, atomistic, and coarse-grained levels, as well as extensive coverage of force fields. In addition to providing the current snapshot of the CHARMM development, this review may serve as a starting point for exploring relevant theories and computational methods for tackling contemporary and emerging problems in biomolecular systems. CHARMM is freely available for academic and nonprofit research at https://academiccharmm.org/program.

In Silico Screening, Molecular Dynamics Simulation and Binding Free Energy Identify Single-Point Mutations That Destabilize p53 and Reduce Binding to DNA.

Considering p53's pivotal role as a tumor suppressor protein, proactive identification and characterization of potentially harmful p53 mutations are crucial before they appear in the population. To address this, four computational prediction tools-SIFT, Polyphen-2, PhD-SNP, and MutPred2-utilizing sequence-based and machine-learning algorithms, were employed to identify potentially deleterious p53 nsSNPs (nonsynonymous single nucleotide polymorphisms) that may impact p53 structure, dynamics, and binding with DNA. These computational methods identified three variants, namely, C141Y, C238S, and L265P, as detrimental to p53 stability. Furthermore, molecular dynamics (MD) simulations revealed that all three variants exhibited heightened structural flexibility compared to the native protein, especially the C141Y and L265P mutations. Consequently, due to the altered structure of mutant p53, the DNA-binding affinity of all three variants decreased by approximately 1.8 to 9.7 times compared to wild-type p53 binding with DNA (14 µM). Notably, the L265P mutation exhibited an approximately ten-fold greater reduction in binding affinity. Consequently, the presence of the L265P mutation in p53 could pose a substantial risk to humans. Given that p53 regulates abnormal tumor growth, this research carries significant implications for surveillance efforts and the development of anticancer therapies.

Exploring Cannabinoids as Potential Inhibitors of SARS-CoV-2 Papain-like Protease: Insights from Computational Analysis and Molecular Dynamics Simulations.

The emergence of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has triggered a global COVID-19 pandemic, challenging healthcare systems worldwide. Effective therapeutic strategies against this novel coronavirus remain limited, underscoring the urgent need for innovative approaches. The present research investigates the potential of cannabis compounds as therapeutic agents against SARS-CoV-2 through their interaction with the virus's papain-like protease (PLpro) protein, a crucial element in viral replication and immune evasion. Computational methods, including molecular docking and molecular dynamics (MD) simulations, were employed to screen cannabis compounds against PLpro and analyze their binding mechanisms and interaction patterns. The results showed cannabinoids with binding affinities ranging from -6.1 kcal/mol to -4.6 kcal/mol, forming interactions with PLpro. Notably, Cannabigerolic and Cannabidiolic acids exhibited strong binding contacts with critical residues in PLpro's active region, indicating their potential as viral replication inhibitors. MD simulations revealed the dynamic behavior of cannabinoid-PLpro complexes, highlighting stable binding conformations and conformational changes over time. These findings shed light on the mechanisms underlying cannabis interaction with SARS-CoV-2 PLpro, aiding in the rational design of antiviral therapies. Future research will focus on experimental validation, optimizing binding affinity and selectivity, and preclinical assessments to develop effective treatments against COVID-19.

Discovery of Nirmatrelvir Resistance Mutations in SARS-CoV-2 3CLpro: A Computational-Experimental Approach.

The COVID-19 pandemic has emphasized the urgency for effective antiviral therapies against SARS-CoV-2. Targeting the main protease (3CLpro) of the virus has emerged as a promising approach, and nirmatrelvir (PF-07321332), the active component of Pfizer's oral drug Paxlovid, has demonstrated remarkable clinical efficacy. However, the emergence of resistance mutations poses a challenge to its continued success. In this study, we employed alchemical free energy perturbation (FEP) alanine scanning to identify nirmatrelvir-resistance mutations within SARS-CoV-2 3CLpro. FEP identified several mutations, which were validated through in vitro IC50 experiments and found to result in 8- and 72-fold increases in nirmatrelvir IC50 values. Additionally, we constructed SARS-CoV-2 omicron replicons containing these mutations, and one of the mutants (S144A/E166A) displayed a 20-fold increase in EC50, confirming the role of FEP in identifying drug-resistance mutations. Our findings suggest that FEP can be a valuable tool in proactively monitoring the emergence of resistant strains and guiding the design of future inhibitors with reduced susceptibility to drug resistance. As nirmatrelvir is currently widely used for treating COVID-19, this research has important implications for surveillance efforts and antiviral development.

Synthesis and Characterization of the Donor-Acceptor Conjugated Polymer PBDB-T Implementing Group IV Element Germanium.

Here, we synthesized and characterized a novel two-dimensional (2D) conjugated electron donor-acceptor (D-A) copolymer (PBDB-T-Ge), wherein the substituent of triethyl germanium was added to the electron donor unit of the polymer. The Turbo-Grignard reaction was used to implement the group IV element into the polymer, resulting in a yield of 86%. This corresponding polymer, PBDB-T-Ge, exhibited a down-shift in the highest occupied molecular orbital (HOMO) level to -5.45 eV while the lowest unoccupied molecular orbital (LUMO) level was -3.64 eV. The peaks in UV-Vis absorption and the PL emission of PBDB-T-Ge were observed at 484 nm and 615 nm, respectively.

A computationally designed ACE2 decoy has broad efficacy against SARS-CoV-2 omicron variants and related viruses in vitro and in vivo.

SARS-CoV-2, especially B.1.1.529/omicron and its sublineages, continues to mutate to evade monoclonal antibodies and antibodies elicited by vaccination. Affinity-enhanced soluble ACE2 (sACE2) is an alternative strategy that works by binding the SARS-CoV-2 S protein, acting as a 'decoy' to block the interaction between the S and human ACE2. Using a computational design strategy, we designed an affinity-enhanced ACE2 decoy, FLIF, that exhibited tight binding to SARS-CoV-2 delta and omicron variants. Our computationally calculated absolute binding free energies (ABFE) between sACE2:SARS-CoV-2 S proteins and their variants showed excellent agreement to binding experiments. FLIF displayed robust therapeutic utility against a broad range of SARS-CoV-2 variants and sarbecoviruses, and neutralized omicron BA.5 in vitro and in vivo. Furthermore, we directly compared the in vivo therapeutic efficacy of wild-type ACE2 (non-affinity enhanced ACE2) against FLIF. A few wild-type sACE2 decoys have shown to be effective against early circulating variants such as Wuhan in vivo. Our data suggest that moving forward, affinity-enhanced ACE2 decoys like FLIF may be required to combat evolving SARS-CoV-2 variants. The approach described herein emphasizes how computational methods have become sufficiently accurate for the design of therapeutics against viral protein targets. Affinity-enhanced ACE2 decoys remain highly effective at neutralizing omicron subvariants.

Photophysical Insights of Halogenated Dipyrrolonaphthyridine-Diones as Potential Photodynamic Therapy Agents†.

We report the synthesis and photophysical characterization of novel halogenated dipyrrolonaphthyridine-diones (X2 -DPNDs, X = Cl, Br, and I), as candidates for photodynamic therapy (PDT) application. Apart from the heavy atom-induced spin-orbit coupling (SOC) dynamics in the investigated X2 -DPNDs, it was found that the position of the halogen atom (relative to the nitrogen of the pyrrole ring) also influenced the triplet excited state behavior. Interestingly, the faster/efficiency sensitization of 3 O2 to 1 O2 using X2 -DPND correlates with the rate of triplet population, kISC >1.6 × 108 s-1 for I2 -DPND vs kISC >2.9 × 109 s-1 for Cl2 -DPND and Br2 -DPND (where τISC  = 343 ± 3 ps for I2 -DPND and τISC  = 5-6 ns for Cl2 -DPND and Br2 -DPND are the lowest time constants/values for ISC). Furthermore, the heavy atom-induced SOC in Cl2 -DPND and Br2 -DPND did not lead to a reduction of the corresponding fluorescence (ca 75% vs 67% for the parent DPND). The attractive photophysical characteristics of Cl2 /Br2 -DPND put them on the landscape as not only promising PDT agents but also as fluorescence probes. The present study is a stepping stone in the development of novel organic photosystems for synergistic photomedicinal applications.

Computationally Designed ACE2 Decoy Receptor Binds SARS-CoV-2 Spike (S) Protein with Tight Nanomolar Affinity.

Even with the availability of vaccines, therapeutic options for COVID-19 still remain highly desirable, especially in hospitalized patients with moderate or severe disease. Soluble ACE2 (sACE2) is a promising therapeutic candidate that neutralizes SARS CoV-2 infection by acting as a decoy. Using computational mutagenesis, we designed a number of sACE2 derivatives carrying three to four mutations. The top-predicted sACE2 decoy based on the in silico mutagenesis scan was subjected to molecular dynamics and free-energy calculations for further validation. After illuminating the mechanism of increased binding for our designed sACE2 derivative, the design was verified experimentally by flow cytometry and BLI-binding experiments. The computationally designed sACE2 decoy (ACE2-FFWF) bound the receptor-binding domain of SARS-CoV-2 tightly with low nanomolar affinity and ninefold affinity enhancement over the wild type. Furthermore, cell surface expression was slightly greater than wild-type ACE2, suggesting that the design is well-folded and stable. Having an arsenal of high-affinity sACE2 derivatives will help to buffer against the emergence of SARS CoV-2 variants. Here, we show that computational methods have become sufficiently accurate for the design of therapeutics for current and future viral pandemics.

In-silico screening and microsecond molecular dynamics simulations to identify single point mutations that destabilize ß-hexosaminidase A causing Tay-Sachs disease.

ß-hexosaminidase A (HexA) protein is responsible for the degradation of GM2 gangliosides in the central and peripheral nervous systems. Tay-Sachs disease occurs when HexA within Hexosaminidase does not properly function and harmful GM2 gangliosides begin to build up within the neurons. In this study, in silico methods such as SIFT, PolyPhen-2, PhD-SNP, and MutPred were utilized to analyze the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) on HexA in order to identify possible pathogenetic and deleterious variants. Molecular dynamics (MD) simulations showed that two mutants, P25S and W485R, experienced an increase in structural flexibility compared to the native protein. Particularly, there was a decrease in the overall number and frequencies of hydrogen bonds for the mutants compared to the wildtype. MM/GBSA calculations were performed to help assess the change in binding affinity between the wildtype and mutant structures and a mechanism-based inhibitor, NGT, which is known to help increase the residual activity of HexA. Both of the mutants experienced a decrease in the binding affinity from -23.8 kcal/mol in wildtype to -20.9 and -18.7 kcal/mol for the P25S and W485R variants of HexA, respectively.

Nitrone and Alkyne Cascade Reactions for Regio- and Diastereoselective 1-Pyrroline Synthesis.

The synthesis of 1-pyrrolines from N-alkenylnitrones and alkynes has been explored as a retrosynthetic alternative to traditional approaches. These cascade reactions are formal [4+1] cycloadditions that proceed through a proposed dipolar cycloaddition and N-alkenylisoxazoline [3,3']-sigmatropic rearrangement. A variety of cyclic alkynes and terminal alkynes have been shown to undergo the transformation with N-alkenylnitrones under mild conditions to provide the corresponding spirocyclic and densely substituted 1-pyrrolines with high regio- and diastereoselectivity. Mechanistic studies provide insight into the balance of steric and electronic effects that promote the cascade process and control the diastereo- and regioisomeric preferences of the 1-pyrroline products. Diastereoselective derivatization of the 1-pyrrolines prepared by the cascade reaction demonstrate the divergent synthetic utility of the new method.