Search for: All records

We analyzed the spike protein S1/S2 cleavage of selected strains of a prototype coronavirus, mouse hepatitis virus (MHV) by the cellular protease furin, in order to understand the structural requirements underlying the sequence selectivity of the scissile segment. The probability of cleavage of selected MHV strains was first evaluated from furin cleavage scores predicted by the ProP computer software, and then cleavage was measured experimentally with a fluorogenic peptide cleavage assay consisting of S1/S2 peptide mimics and purified furin. We found that in vitro cleavability varied across MHV strains in line with predicted results—but with the notable exception of MHV-A59, which was not cleaved despite a high score predicted for its sequence. Using the known X-Ray structure of furin in complex with a substrate-like inhibitor as an initial structural reference, we carried out molecular dynamics (MD) simulations to learn the modes of binding of the peptides in the furin active site, and the suitability of the complex for initiation of the enzymatic cleavage. We identified the 3D structural requirements of the furin active site configuration that enable bound peptides to undergo cleavage, and the way in which the various strains tested experimentally are fulfilling these requirements. We find that despite some flexibility in the organization of the peptide bound to the active site of the enzyme, the presence of a histidine at P2 of MHV-A59 fails to properly orient the sidechain of His194 of the furin catalytic triad and therefore produces a distortion that renders the peptide/complex structural configuration in the active site incompatible with requirements for cleavage initiation. The Ser/Thr in P1 of MHV-2 and MHV-S has a similar effect of distorting the conformation of the furin active site residues produced by the elimination of the canonical salt-bridge formed by arginine in P1 position. This work informs a study of coronavirus infection and pathogenesis with respect to the function of the viral spike protein, and suggests an important process of viral adaptation and evolution within the spike S1/S2 structural loop. 

G protein-coupled receptors (GPCRs) play a key role in many cellular signaling mechanisms, and must select among multiple coupling possibilities in a ligand-specific manner in order to carry out a myriad of functions in diverse cellular contexts. Much has been learned about the molecular mechanisms of ligand-GPCR complexes from Molecular Dynamics (MD) simulations. However, to explore ligand-specific differences in the response of a GPCR to diverse ligands, as is required to understand ligand bias and functional selectivity, necessitates creating very large amounts of data from the needed large-scale simulations. This becomes a Big Data problem for the high dimensionality analysis of the accumulated trajectories. Here we describe a new machine learning (ML) approach to the problem that is based on transforming the analysis of GPCR function-related, ligand-specific differences encoded in the MD simulation trajectories into a representation recognizable by state-of-the-art deep learning object recognition technology. We illustrate this method by applying it to recognize the pharmacological classification of ligands bound to the 5-HT2A and D2 subtypes of class-A GPCRs from the serotonin and dopamine families. The ML-based approach is shown to perform the classification task with high accuracy, and we identify the molecular determinants of the classifications in the context of GPCR structure and function. This study builds a framework for the efficient computational analysis of MD Big Data collected for the purpose of understanding ligand-specific GPCR activity.