Biophysics-Guided Lead Discovery of HBV Capsid Assembly Modifiers.

Hepatitis B virus (HBV) is the leading cause of chronic liver pathologies worldwide. HBV nucleocapsid, a key structural component, is formed through the self-assembly of the capsid protein units. Therefore, interfering with the self-assembly process is a promising approach for the development of novel antiviral agents. Applied to HBV, this approach has led to several classes of capsid assembly modulators (CAMs). Here, we report structurally novel CAMs with moderate activity and low toxicity, discovered through a biophysics-guided approach combining docking, molecular dynamics simulations, and a series of assays with a particular emphasis on biophysical experiments. Several of the identified compounds induce the formation of aberrant capsids and inhibit HBV DNA replication in vitro, suggesting that they possess modest capsid assembly modulation effects. The synergistic computational and experimental approaches provided key insights that facilitated the identification of compounds with promising activities. The discovery of preclinical CAMs presents opportunities for subsequent optimization efforts, thereby opening new avenues for HBV inhibition.

Substrate recognition mechanism of the endoplasmic reticulum-associated ubiquitin ligase Doa10.

Doa10 (MARCHF6 in metazoans) is a large polytopic membrane-embedded E3 ubiquitin ligase in the endoplasmic reticulum (ER) that plays an important role in quality control of cytosolic and ER proteins. Although Doa10 is highly conserved across eukaryotes, it is not understood how Doa10 recognizes its substrates. Here, we define the substrate recognition mechanism of Doa10 by structural and functional analyses on Saccharomyces cerevisiae Doa10 and its model substrates. Cryo-EM analysis shows that Doa10 has unusual architecture with a large lipid-filled central cavity, and its conserved middle domain forms an additional water-filled lateral tunnel open to the cytosol. Our biochemical data and molecular dynamics simulations suggest that the entrance of the substrate's degron peptide into the lateral tunnel is required for efficient polyubiquitination. The N- and C-terminal membrane domains of Doa10 seem to form fence-like features to restrict polyubiquitination to those proteins that can access the central cavity and lateral tunnel. Our study reveals how extended hydrophobic sequences at the termini of substrate proteins are recognized by Doa10 as a signal for quality control.

Substrate recognition mechanism of the endoplasmic reticulum-associated ubiquitin ligase Doa10.

Doa10 (MARCH6 in metazoans) is a large polytopic membrane-embedded E3 ubiquitin ligase in the endoplasmic reticulum (ER) that plays an important role in quality control of cytosolic and ER proteins. Although Doa10 is highly conserved across eukaryotes, it is not understood how Doa10 recognizes its substrates. Here, we defined the substrate recognition mechanism of Doa10 by structural and functional analyses on Saccharomyces cerevisiae Doa10 and its well-defined degron Deg1. Cryo-EM analysis shows that Doa10 has unusual architecture with a large lipid-filled central cavity, and its conserved middle domain forms an additional water-filled lateral tunnel open to the cytosol. Our biochemical data and molecular dynamics simulations suggest that the entrance of the substrate's degron peptide into the lateral tunnel is required for efficient polyubiquitination. The N- and C-terminal membrane domains of Doa10 seem to form fence-like features to restrict polyubiquitination to those proteins that can access the central cavity and lateral tunnel.

Structural basis of mitochondrial protein import by the TIM23 complex.

Mitochondria import nearly all of their approximately 1,000-2,000 constituent proteins from the cytosol across their double-membrane envelope1-5. Genetic and biochemical studies have shown that the conserved protein translocase, termed the TIM23 complex, mediates import of presequence-containing proteins (preproteins) into the mitochondrial matrix and inner membrane. Among about ten different subunits of the TIM23 complex, the essential multipass membrane protein Tim23, together with the evolutionarily related protein Tim17, has long been postulated to form a protein-conducting channel6-11. However, the mechanism by which these subunits form a translocation path in the membrane and enable the import process remains unclear due to a lack of structural information. Here we determined the cryo-electron microscopy structure of the core TIM23 complex (heterotrimeric Tim17-Tim23-Tim44) from Saccharomyces cerevisiae. Contrary to the prevailing model, Tim23 and Tim17 themselves do not form a water-filled channel, but instead have separate, lipid-exposed concave cavities that face in opposite directions. Our structural and biochemical analyses show that the cavity of Tim17, but not Tim23, forms the protein translocation path, whereas Tim23 probably has a structural role. The results further suggest that, during translocation of substrate polypeptides, the nonessential subunit Mgr2 seals the lateral opening of the Tim17 cavity to facilitate the translocation process. We propose a new model for the TIM23-mediated protein import and sorting mechanism, a central pathway in mitochondrial biogenesis.

Understanding Virus Structure and Dynamics through Molecular Simulations.

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

SARS-CoV-2 spike opening dynamics and energetics reveal the individual roles of glycans and their collective impact.

The trimeric spike (S) glycoprotein, which protrudes from the SARS-CoV-2 viral envelope, binds to human ACE2, initiated by at least one protomer's receptor binding domain (RBD) switching from a "down" (closed) to an "up" (open) state. Here, we used large-scale molecular dynamics simulations and two-dimensional replica exchange umbrella sampling calculations with more than a thousand windows and an aggregate total of 160 µs of simulation to investigate this transition with and without glycans. We find that the glycosylated spike has a higher barrier to opening and also energetically favors the down state over the up state. Analysis of the S-protein opening pathway reveals that glycans at N165 and N122 interfere with hydrogen bonds between the RBD and the N-terminal domain in the up state, while glycans at N165 and N343 can stabilize both the down and up states. Finally, we estimate how epitope exposure for several known antibodies changes along the opening path. We find that the BD-368-2 antibody's epitope is continuously exposed, explaining its high efficacy.

The Mechanism of Action of Hepatitis B Virus Capsid Assembly Modulators Can Be Predicted from Binding to Early Assembly Intermediates.

Interfering with the self-assembly of virus nucleocapsids is a promising approach for the development of novel antiviral agents. Applied to hepatitis B virus (HBV), this approach has led to several classes of capsid assembly modulators (CAMs) that target the virus by either accelerating nucleocapsid assembly or misdirecting it into noncapsid-like particles, thereby inhibiting the HBV replication cycle. Here, we have assessed the structures of early nucleocapsid assembly intermediates, bound with and without CAMs, using molecular dynamics simulations. We find that distinct conformations of the intermediates are induced depending on whether the bound CAM accelerates or misdirects assembly. Specifically, the assembly intermediates with bound misdirecting CAMs appear to be flattened relative to those with bound accelerators. Finally, the potency of CAMs within the same class was studied. We find that an increased number of contacts with the capsid protein and favorable binding energies inferred from free energy perturbation calculations are indicative of increased potency.

TRPV1 Activation by Anandamide via a Unique Lipid Pathway.

The capsaicin receptor, transient receptor potential vanilloid type 1 (TRPV1), is a polymodal channel that has been implicated in the perception of pain and can be modulated by a variety of cannabinoid ligands. Here we report TRPV1 channel activation by the endocannabinoid, anandamide (AEA), in a unique, peripheral binding site via extended MD simulations. These results aim to expand the understanding of TRPV1 and assist in the development of new TRPV1 modulators.

Machine Learning Reveals the Critical Interactions for SARS-CoV-2 Spike Protein Binding to ACE2.

SARS-CoV and SARS-CoV-2 bind to the human ACE2 receptor in practically identical conformations, although several residues of the receptor-binding domain (RBD) differ between them. Herein, we have used molecular dynamics (MD) simulations, machine learning (ML), and free-energy perturbation (FEP) calculations to elucidate the differences in binding by the two viruses. Although only subtle differences were observed from the initial MD simulations of the two RBD-ACE2 complexes, ML identified the individual residues with the most distinctive ACE2 interactions, many of which have been highlighted in previous experimental studies. FEP calculations quantified the corresponding differences in binding free energies to ACE2, and examination of MD trajectories provided structural explanations for these differences. Lastly, the energetics of emerging SARS-CoV-2 mutations were studied, showing that the affinity of the RBD for ACE2 is increased by N501Y and E484K mutations but is slightly decreased by K417N.

ACE2 glycans preferentially interact with SARS-CoV-2 over SARS-CoV.

We report a distinct difference in the interactions of the glycans of the host-cell receptor, ACE2, with SARS-CoV-2 and SARS-CoV S-protein receptor-binding domains (RBDs). Our analysis demonstrates that the ACE2 glycan at N322 enhances interactions with the SARS-CoV-2 RBD while the ACE2 glycan at N90 may offer protection against infections of both coronaviruses depending on its composition. The interactions of the ACE2 glycan at N322 with SARS-CoV RBD are blocked by the presence of the RBD glycan at N357 of the SARS-CoV RBD. The absence of this glycosylation site on SARS-CoV-2 RBD may enhance its binding with ACE2.

Tuning Proton Transfer Thermodynamics in SARS-CoV-2 Main Protease: Implications for Catalysis and Inhibitor Design.

The catalytic reaction in SARS-CoV-2 main protease is activated by a proton transfer (PT) from Cys145 to His41. The same PT is likely also required for the covalent binding of some inhibitors. Here we use a multiscale computational approach to investigate the PT thermodynamics in the apo enzyme and in complex with two potent inhibitors, N3 and the α-ketoamide 13b. We show that with the inhibitors the free energy cost to reach the charge-separated state of the active-site dyad is lower, with N3 inducing the most significant reduction. We also show that a few key sites (including specific water molecules) significantly enhance or reduce the thermodynamic feasibility of the PT reaction, with selective desolvation of the active site playing a crucial role. The approach presented is a cost-effective procedure to identify the enzyme regions that control the activation of the catalytic reaction and is thus also useful to guide the design of inhibitors.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease.

The main protease (Mpro) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of Mpro, a cysteine protease, have been determined, facilitating structure-based drug design. Mpro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, Mpro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nucleophile Cys145 have been debated in previous studies of SARS-CoV Mpro, but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 Mpro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of Mpro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α-ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N δ (HD) and N ϵ (HE) protonation of His41 and His164, respectively, the α-ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 Mpro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

Tuning Proton Transfer Thermodynamics in SARS-Cov-2 Main Protease: Implications for Catalysis and Inhibitor Design.

In this comutational work a hybrid quantum mechanics/molecular mechanics approach, the MD-PMM approach, is used to investigate the proton transfer reaction the activates the catalytic activity of SARS-CoV-2 main protease. The proton transfer thermodynamics is investigated for the apo ensyme (i.e., without any bound substrate or inhibitor) and in the presence of a inhibitor, N3, which was previously shown to covalently bind SARS-CoV-2 main protease.

Inhibitor binding influences the protonation states of histidines in SARS-CoV-2 main protease.

The main protease (M pro ) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an attractive target for antiviral therapeutics. Recently, many high-resolution apo and inhibitor-bound structures of M pro , a cysteine protease, have been determined, facilitating structure-based drug design. M pro plays a central role in the viral life cycle by catalyzing the cleavage of SARS-CoV-2 polyproteins. In addition to the catalytic dyad His41-Cys145, M pro contains multiple histidines including His163, His164, and His172. The protonation states of these histidines and the catalytic nu-cleophile Cys145 have been debated in previous studies of SARS-CoV M pro , but have yet to be investigated for SARS-CoV-2. In this work we have used molecular dynamics simulations to determine the structural stability of SARS-CoV-2 M pro as a function of the protonation assignments for these residues. We simulated both the apo and inhibitor-bound enzyme and found that the conformational stability of the binding site, bound inhibitors, and the hydrogen bond networks of M pro are highly sensitive to these assignments. Additionally, the two inhibitors studied, the peptidomimetic N3 and an α -ketoamide, display distinct His41/His164 protonation-state-dependent stabilities. While the apo and the N3-bound systems favored N δ (HD) and N ϵ (HE) protonation of His41 and His164, respectively, the α -ketoamide was not stably bound in this state. Our results illustrate the importance of using appropriate histidine protonation states to accurately model the structure and dynamics of SARS-CoV-2 M pro in both the apo and inhibitor-bound states, a necessary prerequisite for drug-design efforts.

A Closer Look at Anandamide Interaction With TRPV1.

The transient receptor potential subfamily vanilloid type 1 ion channel (TRPV1), located in the peripheral nervous system has been implicated in the perception of pain and possesses the ability to be modulated by various cannabinoid ligands. Because of its location, TRPV1 is an ideal target for the development of novel pain therapeutics. Literature precedent suggests a wide range of cannabinoid ligands can activate TRPV1, but the location and mode of entry is not well understood. Understanding the modes in which cannabinoids can enter and bind to TRPV1 can aid in rational drug design. The first endogenous ligand identified for TRPV1 was the endocannabinoid, anandamide (AEA). The Molecular Dynamics (MD) studies discussed here investigate the entry mode of AEA into TRPV1. During the course of the 10+ microsecond MD simulations, two distinct binding modes were observed: AEA binding in the tunnel formed by the S1-S4 region, and AEA binding in the vanilloid binding pocket, with preference for the former. Unbiased MD simulations have revealed multiple spontaneous binding events into the S1-S4 region, with only one event of AEA binding the vanilloid binding pocket. These results suggest that AEA enters TRPV1 via a novel location between helices S1-S4 via the lipid bilayer.

GPR6 Structural Insights: Homology Model Construction and Docking Studies.

GPR6 is an orphan G protein-coupled receptor that has been associated with the cannabinoid family because of its recognition of a sub-set of cannabinoid ligands. The high abundance of GPR6 in the central nervous system, along with high constitutive activity and a link to several neurodegenerative diseases make GPR6 a promising biological target. In fact, diverse research groups have demonstrated that GPR6 represents a possible target for the treatment of neurodegenerative disorders such as Parkinson's disease, Alzheimer's disease, and Huntington's disease. Several patents have claimed the use of a wide range of pyrazine derivatives as GPR6 inverse agonists for the treatment of Parkinson's disease symptoms and other dyskinesia syndromes. However, the full pharmacological importance of GPR6 has not yet been fully explored due to the lack of high potency, readily available ligands targeting GPR6. The long-term goal of the present study is to develop such ligands. In this paper, we describe our initial steps towards this goal. A human GPR6 homology model was constructed using a suite of computational techniques. This model permitted the identification of unique GPR6 structural features and the exploration of the GPR6 binding crevice. A subset of patented pyrazine analogs were docked in the resultant GPR6 inactive state model to validate the model, rationalize the structure-activity relationships from the reported patents and identify the key residues in the binding crevice for ligand recognition. We will take this structural knowledge into the next phase of GPR6 project, in which scaffold hopping will be used to design new GPR6 ligands.

The Nucleotide-Free State of the Cannabinoid CB2/Gi Complex.

In this issue of Cell, two papers report agonist-bound cryo-EM structures of the cannabinoid receptor, CB2, in complex with Gi. Importantly, beyond providing information that could help distinguish CB2 ligand binding from CB1, these structures support the existence of a nucleotide-free state during G-protein signaling.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning.

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

Towards A Molecular Understanding of The Cannabinoid Related Orphan Receptor GPR18: A Focus on Its Constitutive Activity.

The orphan G-protein coupled receptor (GPCR), GPR18, has been recently proposed as a potential member of the cannabinoid family as it recognizes several endogenous, phytogenic, and synthetic cannabinoids. Potential therapeutic applications for GPR18 include intraocular pressure, metabolic disorders, and cancer. GPR18 has been reported to have high constitutive activity, i.e., activation/signaling occurs in the absence of an agonist. This activity can be reduced significantly by the A3.39N mutation. At the intracellular (IC) ends of (transmembrane helices) TMH3 and TMH6 in GPCRs, typically, a pair of oppositely charged amino acids form a salt bridge called the "ionic lock". Breaking of this salt bridge creates an IC opening for coupling with G protein. The GPR18 "ionic lock" residues (R3.50/S6.33) can form only a hydrogen bond. In this paper, we test the hypothesis that the high constitutive activity of GPR18 is due to the weakness of its "ionic lock" and that the A3.39N mutation strengthens this lock. To this end, we report molecular dynamics simulations of wild-type (WT) GPR18 and the A3.39N mutant in fully hydrated (POPC) phophatidylcholine lipid bilayers. Results suggest that in the A3.39N mutant, TMH6 rotates and brings R3.50 and S6.33 closer together, thus strengthening the GPR18 "ionic lock".

Molecular Dynamics Methodologies for Probing Cannabinoid Ligand/Receptor Interaction.

The cannabinoid type 1 and 2 G-protein-coupled receptors are currently important pharmacological targets with significant drug discovery potential. These receptors have been shown to display functional selectivity or biased agonism, a property currently thought to have substantial therapeutic potential. Although recent advances in crystallization techniques have provided a wealth of structural information about this important class of membrane-embedded proteins, these structures lack dynamical information. In order to fully understand the interplay of structure and function for this important class of proteins, complementary techniques that address the dynamical aspects of their function are required such as NMR as well as a variety of other spectroscopies. Complimentary to these experimental approaches is molecular dynamics, which has been effectively used to help unravel, at the atomic level, the dynamics of ligand binding and activation of these membrane-bound receptors. Here, we discuss and present several representative examples of the application of molecular dynamics simulations to the understanding of the signatures of ligand-binding and -biased signaling at the cannabinoid type 1 and 2 receptors.