RESUMEN
Biology spans a continuum of length and time scales. Individual experimental methods only glimpse discrete pieces of this spectrum but can be combined to construct a more holistic view. In this Review, we detail the latest advancements in volume electron microscopy (vEM) and cryo-electron tomography (cryo-ET), which together can visualize biological complexity across scales from the organization of cells in large tissues to the molecular details inside native cellular environments. In addition, we discuss emerging methodologies for integrating three-dimensional electron microscopy (3DEM) imaging with multimodal data, including fluorescence microscopy, mass spectrometry, single-particle analysis, and AI-based structure prediction. This multifaceted approach fills gaps in the biological continuum, providing functional context, spatial organization, molecular identity, and native interactions. We conclude with a perspective on incorporating diverse data into computational simulations that further bridge and extend length scales while integrating the dimension of time.
Asunto(s)
Biología , Microscopía Electrónica , Microscopía por Crioelectrón/métodos , Tomografía con Microscopio Electrónico/métodos , Microscopía Fluorescente , Tiempo , Simulación por ComputadorRESUMEN
Bacteria encode myriad defences that target the genomes of infecting bacteriophage, including restriction-modification and CRISPR-Cas systems1. In response, one family of large bacteriophages uses a nucleus-like compartment to protect its replicating genomes by excluding host defence factors2-4. However, the principal composition and structure of this compartment remain unknown. Here we find that the bacteriophage nuclear shell assembles primarily from one protein, which we name chimallin (ChmA). Combining cryo-electron tomography of nuclear shells in bacteriophage-infected cells and cryo-electron microscopy of a minimal chimallin compartment in vitro, we show that chimallin self-assembles as a flexible sheet into closed micrometre-scale compartments. The architecture and assembly dynamics of the chimallin shell suggest mechanisms for its nucleation and growth, and its role as a scaffold for phage-encoded factors mediating macromolecular transport, cytoskeletal interactions, and viral maturation.
Asunto(s)
Bacterias , Bacteriófagos , Compartimento Celular , Proteínas Virales , Ensamble de Virus , Bacterias/citología , Bacterias/inmunología , Bacterias/metabolismo , Bacterias/virología , Bacteriófagos/química , Bacteriófagos/inmunología , Bacteriófagos/fisiología , Bacteriófagos/ultraestructura , Microscopía por Crioelectrón , Proteínas Virales/química , Proteínas Virales/metabolismo , Proteínas Virales/ultraestructuraRESUMEN
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.
Asunto(s)
COVID-19 , Humanos , SARS-CoV-2/metabolismo , Enzima Convertidora de Angiotensina 2/metabolismo , Unión Proteica , Peptidil-Dipeptidasa A/metabolismoRESUMEN
High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in nonstructural protein 14 (nsp14), which excises nucleotides including antiviral drugs misincorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here, we determined a 1.6-Å resolution crystal structure of severe acute respiratory syndrome CoV 2 (SARS-CoV-2) ExoN in complex with its essential cofactor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. We also show that the ExoN activity can rescue a stalled RNA primer poisoned with sofosbuvir and allow RdRp to continue its extension in the presence of the chain-terminating drug, biochemically recapitulating proofreading in SARS-CoV-2 replication. Molecular dynamics simulations further show remarkable flexibility of multidomain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA binding to support its exonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anticoronaviral drugs or strategies to attenuate the viral virulence.
Asunto(s)
Exorribonucleasas/química , Simulación de Dinámica Molecular , Conformación de Ácido Nucleico , Dominios Proteicos , ARN Viral/química , SARS-CoV-2/enzimología , Proteínas no Estructurales Virales/química , Sitios de Unión/genética , COVID-19/virología , Dominio Catalítico , Cristalografía por Rayos X , Exorribonucleasas/genética , Exorribonucleasas/metabolismo , Humanos , Lisina/química , Lisina/genética , Lisina/metabolismo , Mutación Missense , Unión Proteica , ARN Viral/genética , ARN Viral/metabolismo , SARS-CoV-2/fisiología , Proteínas no Estructurales Virales/genética , Proteínas no Estructurales Virales/metabolismoRESUMEN
Schwann cells (SCs) undergo phenotypic transformation and then orchestrate nerve repair following a peripheral nervous system injury. The low-density lipoprotein receptor-related protein-1 (LRP1) is significantly upregulated in SCs in response to acute injury, activating cJun and promoting SC survival. Matrix-metalloproteinase-9 (MMP-9) is an LRP1 ligand that binds LRP1 through its hemopexin domain (PEX) and activates SC survival signaling and migration. To identify novel peptide mimetics within the hemopexin domain of MMP-9, we examined the crystal structure of PEX, synthesized four peptides, and examined their potential to bind and activate LRP1. We demonstrate that a 22 amino acid peptide, peptide 2, was the only peptide that activated Akt and ERK1/2 signaling in SCs, similar to a glutathione s-transferase (GST)-fused holoprotein, GST-PEX. Intraneural injection of peptide 2, but not vehicle, into crush-injured sciatic nerves activated cJun greater than 2.5-fold in wild-type mice, supporting that peptide 2 can activate the SC repair signaling in vivo. Peptide 2 also bound to Fc-fusion proteins containing the ligand-binding motifs of LRP1, clusters of complement-like repeats (CCRII and CCRIV). Pulldown and computational studies of alanine mutants of peptide 2 showed that positively charged lysine and arginine amino acids within the peptide are critical for stability and binding to CCRII. Collectively, these studies demonstrate that a novel peptide derived from PEX can serve as an LRP1 agonist and possesses qualities previously associated with LRP1 binding and SC signaling in vitro and in vivo.
Asunto(s)
Hemopexina , Metaloproteinasa 9 de la Matriz , Ratones , Animales , Hemopexina/metabolismo , Metaloproteinasa 9 de la Matriz/metabolismo , Ligandos , Transducción de Señal/fisiología , Péptidos/farmacología , Péptidos/metabolismo , Células de Schwann/metabolismo , Proteína 1 Relacionada con Receptor de Lipoproteína de Baja Densidad/metabolismoRESUMEN
Determining the viability of a new drug molecule is a time- and resource-intensive task that makes computer-aided assessments a vital approach to rapid drug discovery. Here we develop a machine learning algorithm, iMiner, that generates novel inhibitor molecules for target proteins by combining deep reinforcement learning with real-time 3D molecular docking using AutoDock Vina, thereby simultaneously creating chemical novelty while constraining molecules for shape and molecular compatibility with target active sites. Moreover, through the use of various types of reward functions, we have introduced novelty in generative tasks for new molecules such as chemical similarity to a target ligand, molecules grown from known protein bound fragments, and creation of molecules that enforce interactions with target residues in the protein active site. The iMiner algorithm is embedded in a composite workflow that filters out Pan-assay interference compounds, Lipinski rule violations, uncommon structures in medicinal chemistry, and poor synthetic accessibility with options for cross-validation against other docking scoring functions and automation of a molecular dynamics simulation to measure pose stability. We also allow users to define a set of rules for the structures they would like to exclude during the training process and postfiltering steps. Because our approach relies only on the structure of the target protein, iMiner can be easily adapted for the future development of other inhibitors or small molecule therapeutics of any target protein.
RESUMEN
Surfactant monolayers at sea spray aerosol (SSA) surfaces regulate various atmospheric processes including gas transfer, cloud interactions, and radiative properties. Most experimental studies of SSA employ a simplified surfactant mixture of long-chain fatty acids (LCFAs) as a proxy for the sea surface microlayer or SSA surface. However, medium-chain fatty acids (MCFAs) make up nearly 30% of the FA fraction in nascent SSA. Given that LCFA monolayers are easily disrupted upon the introduction of chemical heterogeneity (such as mixed chain lengths), simple FA proxies are unlikely to represent realistic SSA interfaces. Integrating experimental and computational techniques, we characterize the impact that partially soluble MCFAs have on the properties of atmospherically relevant LCFA mixtures. We explore the extent to which the MCFA lauric acid (LA) is surface stabilized by varying acidity, salinity, and monolayer composition. We also discuss the impacts of pH on LCFA-assisted LA retention, where the presence of LCFAs may shift the surface-adsorption equilibria of laurateâthe conjugate baseâtoward higher surface activities. Molecular dynamic simulations suggest a mechanism for the enhanced surface retention of laurate. We conclude that increased FA heterogeneity at SSA surfaces promotes surface activity of soluble FA species, altering monolayer phase behavior and impacting climate-relevant atmospheric processes.
RESUMEN
To advance the mission of in silico cell biology, modeling the interactions of large and complex biological systems becomes increasingly relevant. The combination of molecular dynamics (MD) simulations and Markov state models (MSMs) has enabled the construction of simplified models of molecular kinetics on long timescales. Despite its success, this approach is inherently limited by the size of the molecular system. With increasing size of macromolecular complexes, the number of independent or weakly coupled subsystems increases, and the number of global system states increases exponentially, making the sampling of all distinct global states unfeasible. In this work, we present a technique called independent Markov decomposition (IMD) that leverages weak coupling between subsystems to compute a global kinetic model without requiring the sampling of all combinatorial states of subsystems. We give a theoretical basis for IMD and propose an approach for finding and validating such a decomposition. Using empirical few-state MSMs of ion channel models that are well established in electrophysiology, we demonstrate that IMD models can reproduce experimental conductance measurements with a major reduction in sampling compared with a standard MSM approach. We further show how to find the optimal partition of all-atom protein simulations into weakly coupled subunits.
Asunto(s)
Cadenas de Markov , Proteínas/metabolismo , Simulación por Computador , Cinética , Modelos Moleculares , Simulación de Dinámica Molecular , Conformación Proteica , Proteínas/químicaRESUMEN
To investigate the evolution of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in the immune population, we coincupi bated the authentic virus with a highly neutralizing plasma from a COVID-19 convalescent patient. The plasma fully neutralized the virus for seven passages, but, after 45 d, the deletion of F140 in the spike N-terminal domain (NTD) N3 loop led to partial breakthrough. At day 73, an E484K substitution in the receptor-binding domain (RBD) occurred, followed, at day 80, by an insertion in the NTD N5 loop containing a new glycan sequon, which generated a variant completely resistant to plasma neutralization. Computational modeling predicts that the deletion and insertion in loops N3 and N5 prevent binding of neutralizing antibodies. The recent emergence in the United Kingdom, South Africa, Brazil, and Japan of natural variants with similar changes suggests that SARS-CoV-2 has the potential to escape an effective immune response and that vaccines and antibodies able to control emerging variants should be developed.
Asunto(s)
Sustitución de Aminoácidos , Enzima Convertidora de Angiotensina 2/inmunología , Anticuerpos Neutralizantes/inmunología , Anticuerpos Antivirales/inmunología , COVID-19/inmunología , SARS-CoV-2/genética , Glicoproteína de la Espiga del Coronavirus/inmunología , Enzima Convertidora de Angiotensina 2/química , Enzima Convertidora de Angiotensina 2/genética , Animales , Anticuerpos Neutralizantes/química , Anticuerpos Neutralizantes/genética , Anticuerpos Neutralizantes/farmacología , Anticuerpos Antivirales/química , Anticuerpos Antivirales/genética , Anticuerpos Antivirales/farmacología , Sitios de Unión , COVID-19/genética , COVID-19/virología , Chlorocebus aethiops , Convalecencia , Expresión Génica , Humanos , Evasión Inmune , Sueros Inmunes/química , Modelos Moleculares , Mutación , Pruebas de Neutralización , Unión Proteica , Conformación Proteica en Hélice alfa , Conformación Proteica en Lámina beta , Dominios y Motivos de Interacción de Proteínas , SARS-CoV-2/efectos de los fármacos , SARS-CoV-2/inmunología , SARS-CoV-2/patogenicidad , Glicoproteína de la Espiga del Coronavirus/química , Glicoproteína de la Espiga del Coronavirus/genética , Células VeroRESUMEN
Janus kinases (JAK), a group of proteins in the nonreceptor tyrosine kinase (NRTKs) family, play a crucial role in growth, survival, and angiogenesis. They are activated by cytokines through the Janus kinase-signal transducer and activator of a transcription (JAK-STAT) signaling pathway. JAK-STAT signaling pathways have significant roles in the regulation of cell division, apoptosis, and immunity. Identification of the V617F mutation in the Janus homology 2 (JH2) domain of JAK2 leading to myeloproliferative disorders has stimulated great interest in the drug discovery community to develop JAK2-specific inhibitors. However, such inhibitors should be selective toward JAK2 over other JAKs and display an extended residence time. Recently, novel JAK2/STAT5 axis inhibitors (N-(1H-pyrazol-3-yl)pyrimidin-2-amino derivatives) have displayed extended residence times (hours or longer) on target and adequate selectivity excluding JAK3. To facilitate a deeper understanding of the kinase-inhibitor interactions and advance the development of such inhibitors, we utilize a multiscale Markovian milestoning with Voronoi tessellations (MMVT) approach within the Simulation-Enabled Estimation of Kinetic Rates v.2 (SEEKR2) program to rank order these inhibitors based on their kinetic properties and further explain the selectivity of JAK2 inhibitors over JAK3. Our approach investigates the kinetic and thermodynamic properties of JAK-inhibitor complexes in a user-friendly, fast, efficient, and accurate manner compared to other brute force and hybrid-enhanced sampling approaches.
Asunto(s)
Antineoplásicos , Inhibidores de las Cinasas Janus , Citocinas , Janus Quinasa 2 , Inhibidores de Proteínas Quinasas/farmacología , Inhibidores de Proteínas Quinasas/química , Transducción de Señal , Factores de Transcripción STAT/antagonistas & inhibidoresRESUMEN
We present SEEKR2 (simulation-enabled estimation of kinetic rates version 2)âthe latest iteration in the family of SEEKR programs for using multiscale simulation methods to computationally estimate the kinetics and thermodynamics of molecular processes, in particular, ligand-receptor binding. SEEKR2 generates equivalent, or improved, results compared to the earlier versions of SEEKR but with significant increases in speed and capabilities. SEEKR2 has also been built with greater ease of usability and with extensible features to enable future expansions of the method. Now, in addition to supporting simulations using NAMD, calculations may be run with the fast and extensible OpenMM simulation engine. The Brownian dynamics portion of the calculation has also been upgraded to Browndye 2. Furthermore, this version of SEEKR supports hydrogen mass repartitioning, which significantly reduces computational cost, while showing little, if any, loss of accuracy in the predicted kinetics.
Asunto(s)
Simulación de Dinámica Molecular , Cinética , Ligandos , Unión Proteica , TermodinámicaRESUMEN
The discovery of new drugs is a time consuming and expensive process. Methods such as virtual screening, which can filter out ineffective compounds from drug libraries prior to expensive experimental study, have become popular research topics. As the computational drug discovery community has grown, in order to benchmark the various advances in methodology, organizations such as the Drug Design Data Resource have begun hosting blinded grand challenges seeking to identify the best methods for ligand pose-prediction, ligand affinity ranking, and free energy calculations. Such open challenges offer a unique opportunity for researchers to partner with junior students (e.g., high school and undergraduate) to validate basic yet fundamental hypotheses considered to be uninteresting to domain experts. Here, we, a group of high school-aged students and their mentors, present the results of our participation in Grand Challenge 4 where we predicted ligand affinity rankings for the Cathepsin S protease, an important protein target for autoimmune diseases. To investigate the effect of incorporating receptor dynamics on ligand affinity rankings, we employed the Relaxed Complex Scheme, a molecular docking method paired with molecular dynamics-generated receptor conformations. We found that Cathepsin S is a difficult target for molecular docking and we explore some advanced methods such as distance-restrained docking to try to improve the correlation with experiments. This project has exemplified the capabilities of high school students when supported with a rigorous curriculum, and demonstrates the value of community-driven competitions for beginners in computational drug discovery.
Asunto(s)
Benchmarking , Diseño de Fármacos , Sitios de Unión , Niño , Humanos , Ligandos , Simulación del Acoplamiento Molecular , Unión Proteica , Conformación Proteica , TermodinámicaRESUMEN
CC chemokine receptor 2 (CCR2) is a part of the chemokine receptor family, an important class of therapeutic targets. These class A G-protein coupled receptors (GPCRs) are involved in mammalian signaling pathways and control cell migration toward endogenous CC chemokine ligands, named for the adjacent cysteine motif on their N terminus. Chemokine receptors and their associated ligands are involved in a wide range of diseases and thus have become important drug targets. CCR2, in particular, promotes the metastasis of cancer cells and is also implicated in autoimmunity-driven type-1 diabetes, diabetic nephropathy, multiple sclerosis, asthma, atherosclerosis, neuropathic pain, and rheumatoid arthritis. Although promising, CCR2 antagonists have been largely unsuccessful to date. Here, we investigate the effect of an orthosteric and an allosteric antagonist on CCR2 dynamics by coupling long-timescale molecular dynamics simulations with Markov-state model theory. We find that the antagonists shift CCR2 into several stable inactive conformations that are distinct from the crystal structure conformation and disrupt a continuous internal water and sodium ion pathway, preventing transitions to an active-like state. Several metastable conformations present a cryptic drug-binding pocket near the allosteric site that may be amenable to targeting with small molecules. Without antagonists, the apo dynamics reveal intermediate conformations along the activation pathway that provide insight into the basal dynamics of CCR2 and may also be useful for future drug design.
Asunto(s)
Receptores CCR2/química , Sitio Alostérico , Secuencias de Aminoácidos , Sitios de Unión , Humanos , Ligandos , Simulación de Dinámica Molecular , Unión Proteica , Conformación Proteica , Receptores CCR2/antagonistas & inhibidores , Receptores CCR2/genética , Receptores CCR2/metabolismoRESUMEN
COVID-19 has resulted in huge numbers of infections and deaths worldwide and brought the most severe disruptions to societies and economies since the Great Depression. Massive experimental and computational research effort to understand and characterize the disease and rapidly develop diagnostics, vaccines, and drugs has emerged in response to this devastating pandemic and more than 130 000 COVID-19-related research papers have been published in peer-reviewed journals or deposited in preprint servers. Much of the research effort has focused on the discovery of novel drug candidates or repurposing of existing drugs against COVID-19, and many such projects have been either exclusively computational or computer-aided experimental studies. Herein, we provide an expert overview of the key computational methods and their applications for the discovery of COVID-19 small-molecule therapeutics that have been reported in the research literature. We further outline that, after the first year the COVID-19 pandemic, it appears that drug repurposing has not produced rapid and global solutions. However, several known drugs have been used in the clinic to cure COVID-19 patients, and a few repurposed drugs continue to be considered in clinical trials, along with several novel clinical candidates. We posit that truly impactful computational tools must deliver actionable, experimentally testable hypotheses enabling the discovery of novel drugs and drug combinations, and that open science and rapid sharing of research results are critical to accelerate the development of novel, much needed therapeutics for COVID-19.
Asunto(s)
Tratamiento Farmacológico de COVID-19 , Simulación por Computador , Diseño de Fármacos , Descubrimiento de Drogas/métodos , Reposicionamiento de Medicamentos , Antivirales/uso terapéutico , COVID-19/virología , Ensayos Clínicos como Asunto , Humanos , Pandemias , SARS-CoV-2/efectos de los fármacosRESUMEN
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.
Asunto(s)
Simulación de Dinámica Molecular , SARS-CoV-2/química , Virión/química , COVID-19 , Análisis de Componente Principal , Proteínas Virales/químicaRESUMEN
The coronavirus disease 2019 (COVID-19) pandemic has swept over the world in the past months, causing significant loss of life and consequences to human health. Although numerous drug and vaccine development efforts are underway, there are many outstanding questions on the mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral association to angiotensin-converting enzyme 2 (ACE2), its main host receptor, and host cell entry. Structural and biophysical studies indicate some degree of flexibility in the viral extracellular spike glycoprotein and at the receptor-binding domain (RBD)-receptor interface, suggesting a role in infection. Here, we perform explicitly solvated, all-atom, molecular dynamics simulations of the glycosylated, full-length, membrane-bound ACE2 receptor in both an apo and spike RBD-bound state to probe the intrinsic dynamics of the ACE2 receptor in the context of the cell surface. A large degree of fluctuation in the full-length structure is observed, indicating hinge bending motions at the linker region connecting the head to the transmembrane helix while still not disrupting the ACE2 homodimer or ACE2-RBD interfaces. This flexibility translates into an ensemble of ACE2 homodimer conformations that could sterically accommodate binding of the spike trimer to more than one ACE2 homodimer and suggests a mechanical contribution of the host receptor toward the large spike conformational changes required for cell fusion. This work presents further structural and functional insights into the role of ACE2 in viral infection that can potentially be exploited for the rational design of effective SARS-CoV-2 therapeutics.
Asunto(s)
Enzima Convertidora de Angiotensina 2/metabolismo , COVID-19/enzimología , COVID-19/virología , SARS-CoV-2/fisiología , Enzima Convertidora de Angiotensina 2/química , Humanos , Simulación de Dinámica Molecular , Multimerización de ProteínaRESUMEN
Changeux et al. (Changeux et al. C. R. Biol. 343:33-39.) recently suggested that the SARS-CoV-2 spike protein may interact with nicotinic acetylcholine receptors (nAChRs) and that such interactions may be involved in pathology and infectivity. This hypothesis is based on the fact that the SARS-CoV-2 spike protein contains a sequence motif similar to known nAChR antagonists. Here, we use molecular simulations of validated atomically detailed structures of nAChRs and of the spike to investigate the possible binding of the Y674-R685 region of the spike to nAChRs. We examine the binding of the Y674-R685 loop to three nAChRs, namely the human α4ß2 and α7 subtypes and the muscle-like αßγδ receptor from Tetronarce californica. Our results predict that Y674-R685 has affinity for nAChRs. The region of the spike responsible for binding contains a PRRA motif, a four-residue insertion not found in other SARS-like coronaviruses. The conformational behavior of the bound Y674-R685 is highly dependent on the receptor subtype; it adopts extended conformations in the α4ß2 and α7 complexes but is more compact when bound to the muscle-like receptor. In the α4ß2 and αßγδ complexes, the interaction of Y674-R685 with the receptors forces the loop C region to adopt an open conformation, similar to other known nAChR antagonists. In contrast, in the α7 complex, Y674-R685 penetrates deeply into the binding pocket in which it forms interactions with the residues lining the aromatic box, namely with TrpB, TyrC1, and TyrC2. Estimates of binding energy suggest that Y674-R685 forms stable complexes with all three nAChR subtypes. Analyses of simulations of the glycosylated spike show that the Y674-R685 region is accessible for binding. We suggest a potential binding orientation of the spike protein with nAChRs, in which they are in a nonparallel arrangement to one another.
Asunto(s)
Receptores Nicotínicos/metabolismo , Glicoproteína de la Espiga del Coronavirus/metabolismo , Glicosilación , Humanos , Simulación de Dinámica Molecular , Péptidos/química , Péptidos/metabolismo , Unión Proteica , Receptores Nicotínicos/química , Glicoproteína de la Espiga del Coronavirus/química , TermodinámicaRESUMEN
The 55-residue OCRE domains of the splicing factors RBM5 and RBM10 contain 15 tyrosines in compact, globular folds. At 25 °C, all 15 tyrosines show symmetric 1 H NMR spectra, with averaged signals for the pairs of δ- and ϵ-ring hydrogens. At 4 °C, two tyrosines were identified as showing 1 H NMR line-broadening due to lowered frequency of the ring-flipping. For the other 13 tyrosine rings, it was not evident, from the 1 H NMR data alone, whether they were either all flipping at high frequencies, or whether slowed flipping went undetected due to small chemical-shift differences between pairs of exchanging ring hydrogen atoms. Here, we integrate 1 H NMR spectroscopy and molecular dynamics (MD) simulations to determine the tyrosine ring-flip frequencies. In the RBM10-OCRE domain, we found that, for 11 of the 15 tyrosines, these frequencies are in the range 2.0×106 to 1.3×108 â s-1 , and we established an upper limit of <1.0×106 â s-1 for the remaining four residues. The experimental data and the MD simulation are mutually supportive, and their combined use extends the analysis of aromatic ring-flip events beyond the limitations of routine 1 H NMR line-shape analysis into the nanosecond frequency range.
Asunto(s)
Proteínas de Ciclo Celular/química , Proteínas de Unión al ADN/química , Simulación de Dinámica Molecular , Factores de Empalme de ARN/química , Proteínas de Unión al ARN/química , Proteínas Supresoras de Tumor/química , Tirosina/química , Secuencias de Aminoácidos , Humanos , Espectroscopía de Protones por Resonancia MagnéticaRESUMEN
Human immunodeficiency virus type 1 (HIV-1) Vif recruits a cellular ubiquitin ligase complex to degrade antiviral APOBEC3 enzymes (APOBEC3C-H) and PP2A phosphatase regulators (PPP2R5A to PPP2R5E). While APOBEC3 antagonism is the canonical function of HIV-1 Vif, this viral accessory protein is also known to trigger G2/M cell cycle arrest. Vif initiates G2/M arrest by degrading multiple PPP2R5 family members, an activity prevalent among diverse HIV-1 and simian immunodeficiency virus (SIV) isolates. Here, computational protein-protein docking was used to delineate a Vif/CBF-ß/PPP2R5 complex in which Vif is predicted to bind the same PPP2R5 surface as physiologic phosphatase targets. This model was tested using targeted mutagenesis of amino acid residues within or adjacent to the putative interface to show loss or retention, respectively, of Vif-induced PPP2R5 degradation activity. Additionally, expression of a peptide that mimics cellular targets of PPP2R5s robustly inhibited Vif-mediated degradation of PPP2R5A but not APOBEC3G. Moreover, live-cell imaging studies examining Vif-mediated degradation of PPP2R5A and APOBEC3G within the same cell revealed that PPP2R5A degradation kinetics are comparable to those of APOBEC3G with a half-life of roughly 6 h postinfection, demonstrating that Vif can concurrently mediate the degradation of distinct cellular substrates. Finally, experiments with a panel of patient-derived Vif isolates indicated that PPP2R5A degradation activity is common in patient-derived isolates. Taken together, these results support a model in which PPP2R5 degradation and global changes in the cellular phosphoproteome are likely to be advantageous for viral pathogenesis.IMPORTANCE A critical function of HIV-1 Vif is to counteract the family of APOBEC3 innate immune proteins. It is also widely accepted that Vif induces G2/M cell cycle arrest in several different cell types. Recently, it has been shown that Vif degrades multiple PPP2R5 phosphoregulators to induce the G2/M arrest phenotype. Here, computational approaches are used to test a structural model of the Vif/PPP2R5 complex. In addition, imaging studies are used to show that Vif degrades these PPP2R5 substrates in roughly the same time frame as APOBEC3 degradation and that this activity is prevalent in patient-derived Vif isolates. These studies are important by further defining PPP2R5 proteins as a bona fide substrate of HIV-1 Vif.
Asunto(s)
Desaminasa APOBEC-3G/química , VIH-1/genética , Proteína Fosfatasa 2/química , Productos del Gen vif del Virus de la Inmunodeficiencia Humana/química , Desaminasa APOBEC-3G/genética , Desaminasa APOBEC-3G/metabolismo , Secuencia de Aminoácidos , Sustitución de Aminoácidos , Sitios de Unión , Expresión Génica , Regulación de la Expresión Génica , Vectores Genéticos/química , Vectores Genéticos/metabolismo , Células HEK293 , Infecciones por VIH/virología , VIH-1/aislamiento & purificación , VIH-1/metabolismo , Células HeLa , Interacciones Huésped-Patógeno/genética , Humanos , Cinética , Modelos Moleculares , Mutación , Unión Proteica , Dominios y Motivos de Interacción de Proteínas , Proteína Fosfatasa 2/genética , Proteína Fosfatasa 2/metabolismo , Estructura Secundaria de Proteína , Proteolisis , Proteínas Recombinantes/química , Proteínas Recombinantes/genética , Proteínas Recombinantes/metabolismo , Transducción de Señal , Especificidad por Sustrato , Productos del Gen vif del Virus de la Inmunodeficiencia Humana/genética , Productos del Gen vif del Virus de la Inmunodeficiencia Humana/metabolismoRESUMEN
Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the finite element method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.