Glycoproteomic landscape and structural dynamics of TIM family immune checkpoints enabled by mucinase SmE

Mucin-domain glycoproteins are densely O-glycosylated and play critical roles in a host of biological functions. In particular, the T cell immunoglobulin and mucin-domain containing family of proteins (TIM-1, -3, -4) decorate immune cells and act as key regulators in cellular immunity. However, their dense O-glycosylation remains enigmatic, primarily due to the challenges associated with studying mucin domains. Here, we demonstrate that the mucinase SmE has a unique ability to cleave at residues bearing very complex glycans. SmE enables improved mass spectrometric analysis of several mucins, including the entire TIM family. With this information in-hand, we perform molecular dynamics (MD) simulations of TIM-3 and -4 to understand how glycosylation affects structural features of these proteins. Finally, we use these models to investigate the functional relevance of glycosylation for TIM-3 function and ligand binding. Overall, we present a powerful workflow to better understand the detailed molecular structures and functions of the mucinome.

SARS-CoV-2 evolved variants optimize binding to cellular glycocalyx

Viral variants of concern continue to arise for SARS-CoV-2, potentially impacting both methods for detection and mechanisms of action. Here, we investigate the effect of an evolving spike positive charge in SARS-CoV-2 variants and subsequent interactions with heparan sulfate and the angiotensin converting enzyme 2 (ACE2) in the glycocalyx. We show that the positively charged Omicron variant evolved enhanced binding rates to the negatively charged glycocalyx. Moreover, we discover that while the Omicron spike-ACE2 affinity is comparable to that of the Delta variant, the Omicron spike interactions with heparan sulfate are significantly enhanced, giving rise to a ternary complex of spike-heparan sulfate-ACE2 with a large proportion of double-bound and triple-bound ACE2. Our findings suggest that SARS-CoV-2 variants evolve to be more dependent on heparan sulfate in viral attachment and infection. This discovery enables us to engineer a second-generation lateral-flow test strip that harnesses both heparin and ACE2 to reliably detect all variants of concern, including Omicron.

Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors

Vaccines and drugs have helped reduce disease severity and blunt the spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, ongoing virus transmission, continuous evolution, and increasing selective pressures have the potential to yield viral variants capable of resisting these interventions. Here, we investigate the susceptibility of natural variants of the main protease [Mpro; 3C-like protease (3CLpro)] of SARS-CoV-2 to protease inhibitors. Multiple single amino acid changes in Mpro confer resistance to nirmatrelvir (the active component of Paxlovid). An additional clinical-stage inhibitor, ensitrelvir (Xocova), shows a different resistance mutation profile. Importantly, phylogenetic analyses indicate that several of these resistant variants have pre-existed the introduction of these drugs into the human population and are capable of spreading. These results encourage the monitoring of resistance variants and the development of additional protease inhibitors and other antiviral drugs with different mechanisms of action and resistance profiles for combinatorial therapy.

Targeted protein S-nitrosylation of ACE2 inhibits SARS-CoV-2 infection

Prevention of infection and propagation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a high priority in the Coronavirus Disease 2019 (COVID-19) pandemic. Here we describe S-nitrosylation of multiple proteins involved in SARS-CoV-2 infection, including angiotensin-converting enzyme 2 (ACE2), the receptor for viral entry. This reaction prevents binding of ACE2 to the SARS-CoV-2 spike protein, thereby inhibiting viral entry, infectivity and cytotoxicity. Aminoadamantane compounds also inhibit coronavirus ion channels formed by envelope (E) protein. Accordingly, we developed dual-mechanism aminoadamantane nitrate compounds that inhibit viral entry and, thus, the spread of infection by S-nitrosylating ACE2 via targeted delivery of the drug after E protein channel blockade. These non-toxic compounds are active in vitro and in vivo in the Syrian hamster COVID-19 model and, thus, provide a novel avenue to pursue therapy.

Glycoproteomic landscape and structural dynamics of TIM family immune checkpoints enabled by mucinase SmE

Mucin-domain glycoproteins are densely O-glycosylated and play critical roles in a host of biological functions. In particular, the T cell immunoglobulin and mucin-domain containing family of proteins (TIM-1, -3, -4) decorate immune cells and act as key checkpoint inhibitors in cancer. However, their dense O-glycosylation remains enigmatic both in terms of glycoproteomic landscape and structural dynamics, primarily due to the challenges associated with studying mucin domains. Here, we present a mucinase (SmE) and demonstrate its ability to selectively cleave along the mucin glycoprotein backbone, similar to others of its kind. Unlike other mucinases, though, SmE harbors the unique ability to cleave at residues bearing extremely complex glycans which enabled improved mass spectrometric analysis of several mucins, including the entire TIM family. With this information in-hand, we performed molecular dynamics (MD) simulations of TIM-3 and -4 to demonstrate how glycosylation affects structural features of these proteins. Overall, we present a powerful workflow to better understand the detailed molecular structures of the mucinome.

Transmissible SARS-CoV-2 variants with resistance to clinical protease inhibitors

Vaccines and drugs have helped reduce disease severity and blunt the spread of SARS-CoV-2. However, ongoing virus transmission, continuous evolution, and increasing selective pressures have the potential to yield viral variants capable of resisting these interventions. Here, we investigate the susceptibility of natural variants of the main protease (Mpro/3CLpro) of SARS-CoV-2 to protease inhibitors. Multiple single amino acid changes in Mpro confer resistance to nirmatrelvir (the active component of Paxlovid). An additional clinical-stage inhibitor, ensitrelvir (Xocova), shows a different resistance mutation profile. Importantly, phylogenetic analyses indicate that several of these resistant variants have pre-existed the introduction of these drugs into the human population and are capable of spreading. These results encourage the monitoring of resistance variants and the development of additional protease inhibitors and other antiviral drugs with different mechanisms of action and resistance profiles for combinatorial therapy.

Spike-heparan sulfate interactions in SARS-CoV-2 infection

Recent biochemical, biophysical, and genetic studies have shown that heparan sulfate, a major component of the cellular glycocalyx, participates in infection of SARS-CoV-2 by facilitating the so-called open conformation of the spike protein, which is required for binding to ACE2. This review highlights the involvement of heparan sulfate in the SARS-CoV-2 infection cycle and argues that there is a high degree of coordination between host cell heparan sulfate and asparagine-linked glycans on the spike in enabling ACE2 binding and subsequent infection. The discovery that spike protein binding and infection depends on both viral and host glycans provides insights into the evolution, spread and potential therapies for SARS-CoV-2 and its variants.

Tuaimenal A, a Meroterpene from the Irish Deep-Sea Soft Coral Duva florida, Displays Inhibition of the SARS-CoV-2 3CLpro Enzyme

Cold water benthic environments are a prolific source of structurally diverse molecules with a range of bioactivities against human disease. Specimens of a previously chemically unexplored soft coral, Duva florida, were collected during a deep-sea cruise that sampled marine invertebrates along the Irish continental margin in 2018. Tuaimenal A (1), a cyclized merosesquiterpenoid representing a new carbon scaffold with a highly substituted chromene core, was discovered through exploration of the soft coral secondary metabolome via NMR-guided fractionation. The absolute configuration was determined through vibrational circular dichroism. Functional biochemical assays and in silico docking experiments found tuaimenal A selectively inhibits the viral main protease (3CLpro) of SARS-CoV-2.

Targeted protein S-nitrosylation of ACE2 as potential treatment to prevent spread of SARS-CoV-2 infection

Prevention of infection and propagation of SARS-CoV-2 is of high priority in the COVID-19 pandemic. Here, we describe S-nitrosylation of multiple proteins involved in SARS-CoV-2 infection, including angiotensin converting enzyme 2 (ACE2), the receptor for viral entry. This reaction prevents binding of ACE2 to the SARS-CoV-2 Spike protein, thereby inhibiting viral entry, infectivity, and cytotoxicity. Aminoadamantane compounds also inhibit coronavirus ion channels formed by envelope (E) protein. Accordingly, we developed dual-mechanism aminoadamantane nitrate compounds that inhibit viral entry and thus spread of infection by S-nitrosylating ACE2 via targeted delivery of the drug after E-protein channel blockade. These non-toxic compounds are active in vitro and in vivo in the Syrian hamster COVID-19 model, and thus provide a novel avenue for therapy.

GlycoGrip: Cell Surface-Inspired Universal Sensor for Betacoronaviruses

Inspired by the role of cell-surface glycoproteins as coreceptors for pathogens, we report the development of GlycoGrip: a glycopolymer-based lateral flow assay for detecting SARS-CoV-2 and its variants. GlycoGrip utilizes glycopolymers for primary capture and antispike antibodies labeled with gold nanoparticles for signal-generating detection. A lock-step integration between experiment and computation has enabled efficient optimization of GlycoGrip test strips which can selectively, sensitively, and rapidly detect SARS-CoV-2 and its variants in biofluids. Employing the power of the glycocalyx in a diagnostic assay has distinct advantages over conventional immunoassays as glycopolymers can bind to antigens in a multivalent capacity and are highly adaptable for mutated strains. As new variants of SARS-CoV-2 are identified, GlycoGrip will serve as a highly reconfigurable biosensor for their detection. Additionally, via extensive ensemble-based docking simulations which incorporate protein and glycan motion, we have elucidated important clues as to how heparan sulfate and other glycocalyx components may bind the spike glycoprotein during SARS-CoV-2 host-cell infection. GlycoGrip is a promising and generalizable alternative to costly, labor-intensive RT-PCR, and we envision it will be broadly useful, including for rural or low-income populations that are historically undertested and under-reported in infection statistics.

CIFDock: A novel CHARMM-based flexible receptor-flexible ligand docking protocol

Docking studies play a critical role in the current workflow of drug discovery. However, limitations may often arise through factors including inadequate ligand sampling, a lack of protein flexibility, scoring function inadequacies (e.g., due to metals, co-factors, etc.), and difficulty in retaining explicit water molecules. Herein, we present a novel CHARMM-based induced fit docking (CIFDock) workflow that can circumvent these limitations by employing all-atom force fields coupled to enhanced sampling molecular dynamics procedures. Self-guided Langevin dynamics simulations are used to effectively sample relevant ligand conformations, side chain orientations, crystal water positions, and active site residue motion. Protein flexibility is further enhanced by dynamic sampling of side chain orientations using an expandable rotamer library. Steps in the procedure consisting of fixing individual components (e.g., the ligand) while sampling the other components (e.g., the residues in the active site of the protein) allow for the complex to adapt to conformational changes. Ultimately, all components of the complex-the protein, ligand, and waters-are sampled simultaneously and unrestrained with SGLD to capture any induced fit effects. This modular flexible docking procedure is automated using CHARMM scripting, interfaced with SLURM array processing, and parallelized to use the desired number of processors. We validated the CIFDock procedure by performing cross-docking studies using a data set comprised of 21 pharmaceutically relevant proteins. Five variants of the CHARMM-based SWISSDOCK scoring functions were created to quantify the results of the final generated poses. Results obtained were comparable to, or in some cases improved upon, commercial docking program data.

A glycan gate controls opening of the SARS-CoV-2 spike protein

SARS-CoV-2 infection is controlled by the opening of the spike protein receptor binding domain (RBD), which transitions from a glycan-shielded 'down' to an exposed 'up' state to bind the human angiotensin-converting enzyme 2 receptor and infect cells. While snapshots of the 'up' and 'down' states have been obtained by cryo-electron microscopy and cryo-electron tomagraphy, details of the RBD-opening transition evade experimental characterization. Here over 130 µs of weighted ensemble simulations of the fully glycosylated spike ectodomain allow us to characterize more than 300 continuous, kinetically unbiased RBD-opening pathways. Together with ManifoldEM analysis of cryo-electron microscopy data and biolayer interferometry experiments, we reveal a gating role for the N-glycan at position N343, which facilitates RBD opening. Residues D405, R408 and D427 also participate. The atomic-level characterization of the glycosylated spike activation mechanism provided herein represents a landmark study for ensemble pathway simulations and offers a foundation for understanding the fundamental mechanisms of SARS-CoV-2 viral entry and infection.

A glycan gate controls opening of the SARS-CoV-2 spike protein

SARS-CoV-2 infection is controlled by the opening of the spike protein receptor binding domain (RBD), which transitions from a glycan-shielded "down" to an exposed "up" state in order to bind the human ACE2 receptor and infect cells. While snapshots of the "up" and "down" states have been obtained by cryoEM and cryoET, details of the RBD opening transition evade experimental characterization. Here, over 130 μs of weighted ensemble (WE) simulations of the fully glycosylated spike ectodomain allow us to characterize more than 300 continuous, kinetically unbiased RBD opening pathways. Together with ManifoldEM analysis of cryo-EM data and biolayer interferometry experiments, we reveal a gating role for the N-glycan at position N343, which facilitates RBD opening. Residues D405, R408, and D427 also participate. The atomic-level characterization of the glycosylated spike activation mechanism provided herein achieves a new high-water mark for ensemble pathway simulations and offers a foundation for understanding the fundamental mechanisms of SARS-CoV-2 viral entry and infection.

Modeling Boronic Acid Based Fluorescent Saccharide Sensors: Computational Investigation of d-Fructose Binding to Dimethylaminomethylphenylboronic Acid

Designing organic saccharide sensors for use in aqueous solution is a nontrivial endeavor. Incorporation of hydrogen bonding groups on a sensor's receptor unit to target saccharides is an obvious strategy but not one that is likely to ensure analyte-receptor interactions over analyte-solvent or receptor-solvent interactions. Phenylboronic acids are known to reversibly and covalently bind saccharides (diols in general) with highly selective affinity in aqueous solution. Therefore, recent work has sought to design such sensors and understand their mechanism for allowing fluorescence with bound saccharides. In past work, binding orientations of several saccharides were determined to dimethylaminomethylphenylboronic acid (DMPBA) receptors with an anthracene fluorophore; however, the binding orientation of d-fructose to such a sensor could not be determined. In this work, we investigate the potential binding modes by generating 20 possible bidentate and six possible tridentate modes between fructose and DMPBA, a simplified receptor model. Gas phase and implicit solvent geometry optimizations, with a myriad functional/basis set pairs, were carried out to identify the lowest energy bidentate and tridentate binding modes of d-fructose to DMPBA. An interesting hydrogen transfer was observed during selected bidentate gas phase optimizations; this transfer suggests a strong sharing of the hydrogen atom between the boronate hydroxyl and amine nitrogen.

The Good, the Bad, and the Ugly: "HiPen", a New Dataset for Validating (S)QM/MM Free Energy Simulations

Indirect (S)QM/MM free energy simulations (FES) are vital to efficiently incorporating sufficient sampling and accurate (QM) energetic evaluations when estimating free energies of practical/experimental interest. Connecting between levels of theory, i.e., calculating Δ A l o w → h i g h , remains to be the most challenging step within an indirect FES protocol. To improve calculations of Δ A l o w → h i g h , we must: (1) compare the performance of all FES methods currently available; and (2) compile and maintain datasets of Δ A l o w → h i g h calculated for a wide-variety of molecules so that future practitioners may replicate or improve upon the current state-of-the-art. Towards these two aims, we introduce a new dataset, "HiPen", which tabulates Δ A g a s M M → 3 o b (the free energy associated with switching from an M M to an S C C - D F T B molecular description using the 3ob parameter set in gas phase), calculated for 22 drug-like small molecules. We compare the calculation of this value using free energy perturbation, Bennett's acceptance ratio, Jarzynski's equation, and Crooks' equation. We also predict the reliability of each calculated Δ A g a s M M → 3 o b by evaluating several convergence criteria including sample size hysteresis, overlap statistics, and bias metric ( Π ). Within the total dataset, three distinct categories of molecules emerge: the "good" molecules, for which we can obtain converged Δ A g a s M M → 3 o b using Jarzynski's equation; "bad" molecules which require Crooks' equation to obtain a converged Δ A g a s M M → 3 o b ; and "ugly" molecules for which we cannot obtain reliably converged Δ A g a s M M → 3 o b with either Jarzynski's or Crooks' equations. We discuss, in depth, results from several example molecules in each of these categories and describe how dihedral discrepancies between levels of theory cause convergence failures even for these gas phase free energy simulations.

Computing converged free energy differences between levels of theory via nonequilibrium work methods: Challenges and opportunities

We demonstrate that Jarzynski's equation can be used to reliably compute free energy differences between low and high level representations of systems. The need for such a calculation arises when employing the so-called "indirect" approach to free energy simulations with mixed quantum mechanical/molecular mechanical (QM/MM) Hamiltonians; a popular technique for circumventing extensive simulations involving quantum chemical computations. We have applied this methodology to several small and medium sized organic molecules, both in the gas phase and explicit solvent. Test cases include several systems for which the standard approach; that is, free energy perturbation between low and high level description, fails to converge. Finally, we identify three major areas in which the difference between low and high level representations make the calculation of ΔAlow→high difficult: bond stretching and angle bending, different preferred conformations, and the response of the MM region to the charge distribution of the QM region. © 2016 Wiley Periodicals, Inc.

Methods for Efficiently and Accurately Computing Quantum Mechanical Free Energies for Enzyme Catalysis

Enzyme activity is inherently linked to free energies of transition states, ligand binding, protonation/deprotonation, etc.; these free energies, and thus enzyme function, can be affected by residue mutations, allosterically induced conformational changes, and much more. Therefore, being able to predict free energies associated with enzymatic processes is critical to understanding and predicting their function. Free energy simulation (FES) has historically been a computational challenge as it requires both the accurate description of inter- and intramolecular interactions and adequate sampling of all relevant conformational degrees of freedom. The hybrid quantum mechanical molecular mechanical (QM/MM) framework is the current tool of choice when accurate computations of macromolecular systems are essential. Unfortunately, robust and efficient approaches that employ the high levels of computational theory needed to accurately describe many reactive processes (ie, ab initio, DFT), while also including explicit solvation effects and accounting for extensive conformational sampling are essentially nonexistent. In this chapter, we will give a brief overview of two recently developed methods that mitigate several major challenges associated with QM/MM FES: the QM non-Boltzmann Bennett's acceptance ratio method and the QM nonequilibrium work method. We will also describe usage of these methods to calculate free energies associated with (1) relative properties and (2) along reaction paths, using simple test cases with relevance to enzymes examples.

The conserved active site histidine-glutamate pair of ferrochelatase coordinately catalyzes porphyrin metalation

Ferrochelatase catalyzes the insertion of ferrous iron into protoporphyrin IX to generate heme. Despite recent research on the reaction mechanism of ferrochelatase, the precise roles and localization of individual active site residues in catalysis, particularly those involved in the insertion of the ferrous iron into the protoporphyrin IX substrate, remain controversial. One outstanding question is from which side of the macrocycle of the bound porphyin substrate is the ferrous iron substrate inserted. Pre-steady state kinetic experiments done under single-turnover conditions conclusively demonstrate that metal ion insertion is pH-dependent, and that the conserved active site His-Glu pair coordinately catalyzes the metal ion insertion reaction. Further, p[Formula: see text] calculations and molecular dynamic simulations indicate that the active site His is deprotonated and the protonation state of the Glu relates to the conformational state of ferrochelatase. Specifically, the conserved Glu in the open conformation of ferrochelatase is deprotonated, while it remains protonated in the closed conformation. These findings support not only the role of the His-Glu pair in catalyzing metal ion insertion, as these residues need to be deprotonated to bind the incoming metal ion, but also the importance of the relationship between the protonation state of the Glu residue and the conformation of ferrochelatase. Finally, the results of this study are consistent with our previous proposal that the unwinding of the [Formula: see text]-helix, the major structural determinant of the closed to open conformational transition in ferrochelatase, is associated with the Glu residue binding the Fe[Formula: see text] substrate from a mitochondrial Fe[Formula: see text] donor.

Efficiently computing pathway free energies: New approaches based on chain-of-replica and Non-Boltzmann Bennett reweighting schemes

BACKGROUND

Accurately modeling condensed phase processes is one of computation's most difficult challenges. Include the possibility that conformational dynamics may be coupled to chemical reactions, where multiscale (i.e., QM/MM) methods are needed, and this task becomes even more daunting.

METHODS

Free energy simulations (i.e., molecular dynamics), multiscale modeling, and reweighting schemes.

RESULTS

Herein, we present two new approaches for mitigating the aforementioned challenges. The first is a new chain-of-replica method (off-path simulations, OPS) for computing potentials of mean force (PMFs) along an easily defined reaction coordinate. This development is coupled with a new distributed, highly-parallel replica framework (REPDstr) within the CHARMM package. Validation of these new schemes is carried out on two processes that undergo conformational changes. First is the simple torsional rotation of butane, while a much more challenging glycosidic rotation (in vacuo and solvated) is the second. Additionally, a new approach that greatly improves (i.e., possibly an order of magnitude) the efficiency of computing QM/MM PMFs is introduced and compared to standard schemes. Our efforts are grounded in the recently developed method for efficiently computing QM-based free energies (i.e., QM-Non-Boltzmann Bennett, QM-NBB). Again, we validate this new technique by computing the QM/MM PMF of butane's torsional rotation.

CONCLUSIONS

The OPS-REPDstr method is a promising new approach that overcomes many limitations of standard pathway simulations in CHARMM. The combination of QM-NBB with pathway techniques is very promising as it offers significant advantages over current procedures.

GENERAL SIGNIFICANCE

Efficiently computing potentials of mean force is a major, unresolved, area of interest. This article is part of a Special Issue entitled Recent developments of molecular dynamics.