Potential Pro-Inflammatory Effect of Vitamin E Analogs through Mitigation of Tetrahydrocannabinol (THC) Binding to the Cannabinoid 2 Receptor.

Vitamin E acetate, which is used as a diluent of tetrahydrocannabinol (THC), has been reported as the primary causative agent of e-cigarette, or vaping, product use-associated lung injury (EVALI). Here, we employ in vitro assays, docking, and molecular dynamics (MD) computer simulations to investigate the interaction of vitamin E with the membrane-bound cannabinoid 2 receptor (CB2R), and its role in modulating the binding affinity of THC to CB2R. From the MD simulations, we determined that vitamin E interacts with both CB2R and membrane phospholipids. Notably, the synchronized effect of these interactions likely facilitates vitamin E acting as a lipid modulator for the cannabinoid system. Furthermore, MD simulation and trajectory analysis show that when THC binds to CB2R in the presence of vitamin E, the binding cavity widens, facilitating the entry of water molecules into it, leading to a reduced interaction of THC with CB2R. Additionally, the interaction between THC and vitamin E in solution is stabilized by several H bonds, which can directly limit the interaction of free THCs with CB2R. Overall, both the MD simulations and the in vitro dissociation assay results indicate that THC binding to CB2R is reduced in the presence of vitamin E. Our study discusses the role of vitamin E in limiting the effect of THCs and its implications on the reported pathology of EVALI.

Targeting SARS-CoV-2 M3CLpro by HCV NS3/4a Inhibitors: In Silico Modeling and In Vitro Screening.

Currently the entire human population is in the midst of a global pandemic caused by SARS-CoV-2 (Severe Acute Respiratory Syndrome CoronaVirus 2). This highly pathogenic virus has to date caused >71 million infections and >1.6 million deaths in >180 countries. Several vaccines and drugs are being studied as possible treatments or prophylactics of this viral infection. M3CLpro (coronavirus main cysteine protease) is a promising drug target as it has a significant role in viral replication. Here we use the X-ray crystal structure of M3CLpro in complex with boceprevir to study the dynamic changes of the protease upon ligand binding. The binding free energy was calculated for water molecules at different locations of the binding site, and molecular dynamics (MD) simulations were carried out for the M3CLpro/boceprevir complex, to thoroughly understand the chemical environment of the binding site. Several HCV NS3/4a protease inhibitors were tested in vitro against M3CLpro. Specifically, asunaprevir, narlaprevir, paritaprevir, simeprevir, and telaprevir all showed inhibitory effects on M3CLpro. Molecular docking and MD simulations were then performed to investigate the effects of these ligands on M3CLpro and to provide insights into the chemical environment of the ligand binding site. Our findings and observations are offered to help guide the design of possible potent protease inhibitors and aid in coping with the COVID-19 pandemic.

Discovery of Novel Small-Molecule Inhibitors of SARS-CoV-2 Main Protease as Potential Leads for COVID-19 Treatment.

The main protease of SARS-CoV-2 virus, Mpro, is an essential element for viral replication, and inhibitors targeting Mpro are currently being investigated in many drug development programs as a possible treatment for COVID-19. An in vitro pilot screen of a highly focused collection of compounds was initiated to identify new lead scaffolds for Mpro. These efforts identified a number of hits. The most effective of these was compound SIMR-2418 having an inhibitory IC50 value of 20.7 µM. Molecular modeling studies were performed to understand the binding characteristics of the identified compounds. The presence of a cyclohexenone warhead group facilitated covalent binding with the Cys145 residue of Mpro. Our results highlight the challenges of targeting Mpro protease and pave the way toward the discovery of potent lead molecules.

Molecular simulations of peptide amphiphiles.

This review describes recent progress in the area of molecular simulations of peptide assemblies, including peptide-amphiphiles and drug-amphiphiles. The ability to predict the structure and stability of peptide self-assemblies from the molecular level up is vital to the field of nanobiotechnology. Computational methods such as molecular dynamics offer the opportunity to characterize intermolecular forces between peptide-amphiphiles that are critical to the self-assembly process. Furthermore, these computational methods provide the ability to computationally probe the structure of these supramolecular assemblies at the molecular level, which is a challenge experimentally. Herein, we briefly highlight progress in the areas of all-atomistic and coarse-grained simulation studies investigating the self-assembly process of short peptides and peptide amphiphiles. We also discuss recent all-atomistic and coarse-grained simulations of the self-assembly of a drug-amphiphile into elongated filaments. Next, we discuss how these computational methods can provide further insight into the pathway of cylindrical nanofiber formation and predict their biocompatibility by studying the interaction of these peptide-amphiphile nanostructures with model cell membranes.

Sequence Dependence in Nucleosome Dynamics.

The basic packaging unit of eukaryotic chromatin is the nucleosome that contains 145-147 base pair duplex DNA wrapped around an octameric histone protein. While the DNA sequence plays a crucial role in controlling the positioning of the nucleosome, the molecular details behind the interplay between DNA sequence and nucleosome dynamics remain relatively unexplored. This study analyzes this interplay in detail by performing all-atom molecular dynamics simulations of nucleosomes, comparing the human α-satellite palindromic (ASP) and the strong positioning "Widom-601" DNA sequence at time scales of 12 µs. The simulations are performed at salt concentrations 10-20 times higher than physiological salt concentrations to screen the electrostatic interactions and promote unwrapping. These microsecond-long simulations give insight into the molecular-level sequence-dependent events that dictate the pathway of DNA unwrapping. We find that the "ASP" sequence forms a loop around SHL ± 5 for three sets of simulations. Coincident with loop formation is a cooperative increase in contacts with the neighboring N-terminal H2B tail and C-terminal H2A tail and the release of neighboring counterions. We find that the Widom-601 sequence exhibits a strong breathing motion of the nucleic acid ends. Coincident with the breathing motion is the collapse of the full N-terminal H3 tail and formation of an α-helix that interacts with the H3 histone core. We postulate that the dynamics of these histone tails and their modification with post-translational modifications (PTMs) may play a key role in governing this dynamics.

Redetermination and absolute configuration of (+)-7-epiclusianone.

The absolute configuration of 3-benzoyl-4-hy-droxy-6,6-dimethyl-1,5,7-tris-(3-methyl-but-2-en-yl)bicyclo-[3.3.1]non-3-ene-2,9-dione, C(33)H(42)O(4), isolated from Hypericum hypericoides, has been determined. The previous study [Xiao et al. (2007 ▶). J. Nat. Prod.70, 1779-1782] gave only the established relative configuration. The three stereogenic centers are now established as 1R, 5R and 7S on the basis of the refinement of the Flack absolute structure parameter against Cu Kα data and correspond to a specific rotation of [α](D) (20) = +66°. The enol-hy-droxy group forms an intra-molecular O-H⋯O hydrogen bond to close an S(6) ring.

Understanding the Dynamics of the Structural States of Cannabinoid Receptors and the Role of Different Modulators.

The cannabinoid receptors CB1R and CB2R are members of the G protein-coupled receptor (GPCR) family. These receptors have recently come to light as possible therapeutic targets for conditions affecting the central nervous system. However, because CB1R is known to have psychoactive side effects, its potential as a drug target is constrained. Therefore, targeting CB2R has become the primary focus of recent research. Using various molecular modeling studies, we analyzed the active, inactive, and intermediate states of both CBRs in this study. We conducted in-depth research on the binding properties of various groups of cannabinoid modulators, including agonists, antagonists, and inverse agonists, with all of the different conformational states of the CBRs. The binding effects of these modulators were studied on various CB structural features, including the movement of the transmembrane helices, the volume of the binding cavity, the internal fluids, and the important GPCR properties. Then, using in vitro experiments and computational modeling, we investigated how vitamin E functions as a lipid modulator to influence THC binding. This comparative examination of modulator binding to CBRs provides significant insight into the mechanisms of structural alterations and ligand affinity, which can directly help in the rational design of selective modulators that target either CB1R or CB2R.

The Interaction of Supramolecular Anticancer Drug Amphiphiles with Phospholipid Membranes.

The shape of drug delivery vehicles impacts both the circulation time and the effectiveness of the vehicle. Peptide-based drug amphiphiles (DAs) are promising new candidates as drug delivery vehicles that can self-assemble into shapes such as nanofilament and nanotube (diameter ~ 6-10 nm). The number of conjugated drugs affects the IC50 of these DAs, which is correlated to the effective cellular uptake. Characterizing and optimizing the interaction of these DAs and their assemblies with the cellular membrane is experimentally challenging. Long-time molecular dynamics can determine if the DA molecular structure affects the translocation across and interaction with the cellular membrane. Here, we report long-time atomistic simulation on Anton 2 (up to 25 µs) of these DAs with model cellular membranes. Results indicate that the interaction of these DAs with model cellular membranes is dependent on the number of conjugated drugs. We find that, with increased drug loading, the hydrophobic drug (camptothecin) builds up in the outer hydrophobic core of the membrane, pulling in positively charged peptide groups. Next, we computationally probe the interaction of differing shapes of these model drug delivery vehicles-nanofilament and nanotube-with the same model membranes, finding that the interaction of these nanostructures with the membrane is strongly repulsive. Results suggest that the hydrogen bond density between the nanostructure and the membrane may play a key role in modulating the interaction between the nanostructure and the membrane. Taken together, these results offer important insights for the rational design of peptide-based drug delivery vehicles.

Rational Coarse-Grained Molecular Dynamics Simulations of Supramolecular Anticancer Nanotubes.

Peptide self-assembly has been used to design an array of nanostructures that possess functional biomedical applications. Experimental studies have reported nanofilament and nanotube formation from peptide-based drug amphiphiles (DAs). These DAs have shown to possess an inherently high drug loading with a tunable release mechanism. Herein, we report rational coarse-grained molecular dynamics simulations of the self-assembly process and the structure and stability of preassembled nanotubes at longer timescales (µs). We find that aggregation between these DAs at the submicrosecond timescale is driven by directional aromatic interactions between the drugs. The drugs form a large and high-density nucleus that is stable throughout microsecond timescales. Simulations of nanotubes characterize the drug-drug stacking and find correlations at nanometer length scales. These simulations can inform the rational molecular design of drug amphiphiles.

Effect of Nucleotide State on the Protofilament Conformation of Tubulin Octamers.

At the molecular level, the dynamic instability (random growth and shrinkage) of the microtubule (MT) is driven by the nucleotide state (GTP vs GDP) in the ß subunit of the tubulin dimers at the MT cap. Here, we use large-scale molecular dynamics (MD) simulations and normal-mode analysis (NMA) to characterize the effect of a single GTP cap layer on tubulin octamers composed of two neighboring protofilaments (PFs). We utilize recently reported high-resolution structures of dynamic MTs to simulate a GDP octamer both with and without a single GTP cap layer. We perform multiple replicas of long-time atomistic MD simulations (3 replicas, 0.3 µs for each replica, 0.9 µs for each octamer system, and 1.8 µs total) of both octamers. We observe that a single GTP cap layer induces structural differences in neighboring PFs, finding that one PF possesses a gradual curvature, compared to the second PF which possesses a kinked conformation. This results in either curling or splaying between these PFs. We suggest that this is due to asymmetric strengths of longitudinal contacts between the two PFs. Furthermore, using NMA, we calculate mechanical properties of these octamer systems and find that octamer system with a single GTP cap layer possesses a lower flexural rigidity.