1
|
Bosch A, Guzman HV, Pérez R. Adsorption-Driven Deformation and Footprints of the RBD Proteins in SARS-CoV-2 Variants on Biological and Inanimate Surfaces. J Chem Inf Model 2024; 64:5977-5990. [PMID: 39083670 PMCID: PMC11323246 DOI: 10.1021/acs.jcim.4c00460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Revised: 06/03/2024] [Accepted: 06/04/2024] [Indexed: 08/02/2024]
Abstract
Respiratory viruses, carried through airborne microdroplets, frequently adhere to surfaces, including plastics and metals. However, our understanding of the interactions between viruses and materials remains limited, particularly in scenarios involving polarizable surfaces. Here, we investigate the role of the receptor-binding domain (RBD) of the spike protein mutations on the adsorption of SARS-CoV-2 to hydrophobic and hydrophilic surfaces employing molecular simulations. To contextualize our findings, we contrast the interactions on inanimate surfaces with those on native biological interfaces, specifically the angiotensin-converting enzyme 2. Notably, we identify a 2-fold increase in structural deformations for the protein's receptor binding motif (RBM) onto inanimate surfaces, indicative of enhanced shock-absorbing mechanisms. Furthermore, the distribution of adsorbed amino acids (landing footprints) on the inanimate surface reveals a distinct regional asymmetry relative to the biological interface, with roughly half of the adsorbed amino acids arranged in opposite sites. In spite of the H-bonds formed at the hydrophilic substrate, the simulations consistently show a higher number of contacts and interfacial area with the hydrophobic surface, where the wild-type RBD adsorbs more strongly than the Delta or Omicron RBDs. In contrast, the adsorption of Delta and Omicron to hydrophilic surfaces was characterized by a distinctive hopping-pattern. The novel shock-absorbing mechanisms identified in the virus adsorption on inanimate surfaces show the embedded high-deformation capacity of the RBD without losing its secondary structure, which could lead to current experimental strategies in the design of virucidal surfaces.
Collapse
Affiliation(s)
- Antonio
M. Bosch
- Departamento
de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
| | - Horacio V. Guzman
- Departamento
de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
- Department
of Theoretical Physics, Jožef Stefan
Institute, SI-1000 Ljubljana, Slovenia
| | - Rubén Pérez
- Departamento
de Física Teórica de la Materia Condensada, Universidad Autónoma de Madrid, E-28049 Madrid, Spain
- Condensed
Matter Physics Center (IFIMAC), Universidad
Autónoma de Madrid, E-28049 Madrid, Spain
| |
Collapse
|
2
|
Brown CM, Marrink SJ. Modeling membranes in situ. Curr Opin Struct Biol 2024; 87:102837. [PMID: 38744147 DOI: 10.1016/j.sbi.2024.102837] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Revised: 03/26/2024] [Accepted: 04/23/2024] [Indexed: 05/16/2024]
Abstract
Molecular dynamics simulations of cellular membranes have come a long way-from simple model lipid bilayers to multicomponent systems capturing the crowded and complex nature of real cell membranes. In this opinionated minireview, we discuss the current challenge to simulate the dynamics of membranes in their native environment, in situ, with the prospect of reaching the level of whole cells and cell organelles using an integrative modeling framework.
Collapse
Affiliation(s)
- Chelsea M Brown
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, the Netherlands. https://twitter.com/chelseabrowncg
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Groningen, the Netherlands. s.j.marrinkrug.nl
| |
Collapse
|
3
|
Biriukov D, Vácha R. Pathways to a Shiny Future: Building the Foundation for Computational Physical Chemistry and Biophysics in 2050. ACS PHYSICAL CHEMISTRY AU 2024; 4:302-313. [PMID: 39069976 PMCID: PMC11274290 DOI: 10.1021/acsphyschemau.4c00003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2024] [Revised: 03/15/2024] [Accepted: 03/18/2024] [Indexed: 07/30/2024]
Abstract
In the last quarter-century, the field of molecular dynamics (MD) has undergone a remarkable transformation, propelled by substantial enhancements in software, hardware, and underlying methodologies. In this Perspective, we contemplate the future trajectory of MD simulations and their possible look at the year 2050. We spotlight the pivotal role of artificial intelligence (AI) in shaping the future of MD and the broader field of computational physical chemistry. We outline critical strategies and initiatives that are essential for the seamless integration of such technologies. Our discussion delves into topics like multiscale modeling, adept management of ever-increasing data deluge, the establishment of centralized simulation databases, and the autonomous refinement, cross-validation, and self-expansion of these repositories. The successful implementation of these advancements requires scientific transparency, a cautiously optimistic approach to interpreting AI-driven simulations and their analysis, and a mindset that prioritizes knowledge-motivated research alongside AI-enhanced big data exploration. While history reminds us that the trajectory of technological progress can be unpredictable, this Perspective offers guidance on preparedness and proactive measures, aiming to steer future advancements in the most beneficial and successful direction.
Collapse
Affiliation(s)
- Denys Biriukov
- CEITEC
− Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- National
Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
| | - Robert Vácha
- CEITEC
− Central European Institute of Technology, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- National
Centre for Biomolecular Research, Faculty of Science, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
- Department
of Condensed Matter Physics, Faculty of Science, Masaryk University, Kotlářská 267/2, 611 37 Brno, Czech
Republic
| |
Collapse
|
4
|
Cornet J, Coulonges N, Pezeshkian W, Penissat-Mahaut M, Desgrez-Dautet H, Marrink SJ, Destainville N, Chavent M, Manghi M. There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning. SOFT MATTER 2024; 20:4998-5013. [PMID: 38884641 DOI: 10.1039/d4sm00089g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/18/2024]
Abstract
We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.
Collapse
Affiliation(s)
- Julie Cornet
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, France.
| | - Nelly Coulonges
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, France.
- Institut de Pharmacologie et Biologie Structurale, Université de Toulouse, CNRS, Université Toulouse III - Paul Sabatier, 31400, Toulouse, France.
| | - Weria Pezeshkian
- Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen, Denmark
| | - Maël Penissat-Mahaut
- Institut de Pharmacologie et Biologie Structurale, Université de Toulouse, CNRS, Université Toulouse III - Paul Sabatier, 31400, Toulouse, France.
| | - Hermes Desgrez-Dautet
- Laboratoire de Microbiologie et Génétique Moléculaires (LMGM), Centre de Biologie Intégrative (CBI), Université de Toulouse, CNRS, UPS, France
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | | | - Matthieu Chavent
- Institut de Pharmacologie et Biologie Structurale, Université de Toulouse, CNRS, Université Toulouse III - Paul Sabatier, 31400, Toulouse, France.
- Laboratoire de Microbiologie et Génétique Moléculaires (LMGM), Centre de Biologie Intégrative (CBI), Université de Toulouse, CNRS, UPS, France
| | - Manoel Manghi
- Laboratoire de Physique Théorique, Université de Toulouse, CNRS, UPS, France.
| |
Collapse
|
5
|
Balabdaoui F, Wierzbicki T, Bao E. Reconstruction of the real 3D shape of the SARS-CoV-2 virus. Biophys J 2024; 123:1297-1310. [PMID: 38715359 PMCID: PMC11140469 DOI: 10.1016/j.bpj.2024.04.019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2023] [Revised: 02/23/2024] [Accepted: 04/19/2024] [Indexed: 05/19/2024] Open
Abstract
The photographs of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus taken by electron transmission microscopy and cryoelectron microscopy provide only a 2D silhouette. The viruses appear to look like distorted circles. The present paper questions the real shape of the SARS-CoV-2 virus and makes an attempt to give an answer. Is this a general ellipsoid, a spheroid with rotational symmetry, a sphere, or something else? The answer requires the application of tools from three different disciplines: structural mechanics, microbiology, and statistics. A total of 590 virus photographs taken from 22 recently published papers were examined. From this experimental data pool, the histogram of diameter ratios was built from the 283 measurements where the virus images could be approximated as ellipses. The curve peaks at the diameter ratio of 1.22. The transformation equation for the spatial shape to the planar shade was derived for a fixed light source of the microscope. This equation involves an unknown orientation of the viruses with respect to the microscope. Two sets of models were developed, one with a uniform distribution of the virus orientation and the other with the orientation defined by the normalized beta distribution. In both sets of models, the unknown diameter ratio of the spheroidal virus was regarded as a random realization from translated gamma distributions. The parameters of the distribution of the kernel functions were determined by minimizing the mean square difference between the predicted and measured 2D histograms. The information included in the measured histograms was found to be insufficient to find an unknown distribution of the virus's orientation. Simply too many unknown parameters render the solution physically unrealistic. The minimization procedure with a uniform probability of virus orientation predicted the peak of the aspect ratio of the 3D spheroid at 1.32. Based on this result, models of the virus will be developed in the continuation of this research for a full dynamic analysis.
Collapse
Affiliation(s)
| | | | - Emma Bao
- Duke University, Durham, North Carolina
| |
Collapse
|
6
|
SureshKumar H, Appadurai R, Srivastava A. Glycans modulate lipid binding in Lili-Mip lipocalin protein: insights from molecular simulations and protein network analyses. Glycobiology 2024; 34:cwad094. [PMID: 38015986 DOI: 10.1093/glycob/cwad094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Revised: 11/23/2023] [Accepted: 11/24/2023] [Indexed: 11/30/2023] Open
Abstract
The unique viviparous Pacific Beetle cockroaches provide nutrition to their embryo by secreting milk proteins Lili-Mip, a lipid-binding glycoprotein that crystallises in-vivo. The resolved in-vivo crystal structure of variably glycosylated Lili-Mip shows a classical Lipocalin fold with an eight-stranded antiparallel beta-barrel enclosing a fatty acid. The availability of physiologically unaltered glycoprotein structure makes Lili-Mip a very attractive model system to investigate the role of glycans on protein structure, dynamics, and function. Towards that end, we have employed all-atom molecular dynamics simulations on various glycosylated stages of a bound and free Lili-Mip protein and characterised the impact of glycans and the bound lipid on the dynamics of this glycoconjugate. Our work provides important molecular-level mechanistic insights into the role of glycans in the nutrient storage function of the Lili-Mip protein. Our analyses show that the glycans stabilise spatially proximal residues and regulate the low amplitude opening motions of the residues at the entrance of the binding pocket. Glycans also preserve the native orientation and conformational flexibility of the ligand. However, we find that either deglycosylation or glycosylation with high-mannose and paucimannose on the core glycans, which better mimic the natural insect glycosylation state, significantly affects the conformation and dynamics. A simple but effective distance- and correlation-based network analysis of the protein also reveals the key residues regulating the barrel's architecture and ligand binding characteristics in response to glycosylation.
Collapse
Affiliation(s)
- Harini SureshKumar
- Molecular Biophysics Unit, Indian Institute of Science, C. V. Raman Road, Bangalore, KA 560012, India
| | - Rajeswari Appadurai
- Molecular Biophysics Unit, Indian Institute of Science, C. V. Raman Road, Bangalore, KA 560012, India
| | - Anand Srivastava
- Molecular Biophysics Unit, Indian Institute of Science, C. V. Raman Road, Bangalore, KA 560012, India
| |
Collapse
|
7
|
Methorst J, van Hilten N, Hoti A, Stroh KS, Risselada HJ. When Data Are Lacking: Physics-Based Inverse Design of Biopolymers Interacting with Complex, Fluid Phases. J Chem Theory Comput 2024; 20:1763-1776. [PMID: 38413010 PMCID: PMC10938504 DOI: 10.1021/acs.jctc.3c00874] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2023] [Revised: 01/03/2024] [Accepted: 01/03/2024] [Indexed: 02/29/2024]
Abstract
Biomolecular research traditionally revolves around comprehending the mechanisms through which peptides or proteins facilitate specific functions, often driven by their relevance to clinical ailments. This conventional approach assumes that unraveling mechanisms is a prerequisite for wielding control over functionality, which stands as the ultimate research goal. However, an alternative perspective emerges from physics-based inverse design, shifting the focus from mechanisms to the direct acquisition of functional control strategies. By embracing this methodology, we can uncover solutions that might not have direct parallels in natural systems, yet yield crucial insights into the isolated molecular elements dictating functionality. This provides a distinctive comprehension of the underlying mechanisms.In this context, we elucidate how physics-based inverse design, facilitated by evolutionary algorithms and coarse-grained molecular simulations, charts a promising course for innovating the reverse engineering of biopolymers interacting with intricate fluid phases such as lipid membranes and liquid protein phases. We introduce evolutionary molecular dynamics (Evo-MD) simulations, an approach that merges evolutionary algorithms with the Martini coarse-grained force field. This method directs the evolutionary process from random amino acid sequences toward peptides interacting with complex fluid phases such as biological lipid membranes, offering significant promises in the development of peptide-based sensors and drugs. This approach can be tailored to recognize or selectively target specific attributes such as membrane curvature, lipid composition, membrane phase (e.g., lipid rafts), and protein fluid phases. Although the resulting optimal solutions may not perfectly align with biological norms, physics-based inverse design excels at isolating relevant physicochemical principles and thermodynamic driving forces governing optimal biopolymer interaction within complex fluidic environments. In addition, we expound upon how physics-based evolution using the Evo-MD approach can be harnessed to extract the evolutionary optimization fingerprints of protein-lipid interactions from native proteins. Finally, we outline how such an approach is uniquely able to generate strategic training data for predictive neural network models that cover the whole relevant physicochemical domain. Exploring challenges, we address key considerations such as choosing a fitting fitness function to delineate the desired functionality. Additionally, we scrutinize assumptions tied to system setup, the targeted protein structure, and limitations posed by the utilized (coarse-grained) force fields and explore potential strategies for guiding evolution with limited experimental data. This discourse encapsulates the potential and remaining obstacles of physics-based inverse design, paving the way for an exciting frontier in biomolecular research.
Collapse
Affiliation(s)
- Jeroen Methorst
- Leiden
Institute of Chemistry, Leiden University, 2333 CC Leiden, The Netherlands
- Department
of Physics, Technische Universität
Dortmund, 44227 Dortmund, Germany
| | - Niek van Hilten
- Leiden
Institute of Chemistry, Leiden University, 2333 CC Leiden, The Netherlands
| | - Art Hoti
- Leiden
Institute of Chemistry, Leiden University, 2333 CC Leiden, The Netherlands
| | - Kai Steffen Stroh
- Department
of Physics, Technische Universität
Dortmund, 44227 Dortmund, Germany
| | - Herre Jelger Risselada
- Leiden
Institute of Chemistry, Leiden University, 2333 CC Leiden, The Netherlands
- Department
of Physics, Technische Universität
Dortmund, 44227 Dortmund, Germany
| |
Collapse
|
8
|
Alvares CMS, Semino R. Force matching and iterative Boltzmann inversion coarse grained force fields for ZIF-8. J Chem Phys 2024; 160:094115. [PMID: 38445731 DOI: 10.1063/5.0190807] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Accepted: 02/06/2024] [Indexed: 03/07/2024] Open
Abstract
Despite the intense activity at electronic and atomistic resolutions, coarse grained (CG) modeling of metal-organic frameworks remains largely unexplored. One of the main reasons for this is the lack of adequate CG force fields. In this work, we present iterative Boltzmann inversion and force matching (FM) force fields for modeling ZIF-8 at three different coarse grained resolutions. Their ability to reproduce structure, elastic tensor, and thermal expansion is evaluated and compared with that of MARTINI force fields considered in previous work [Alvares et al., J. Chem. Phys. 158, 194107 (2023)]. Moreover, MARTINI and FM are evaluated for their ability to depict the swing effect, a subtle phase transition ZIF-8 undergoes when loaded with guest molecules. Overall, we found that all our force fields reproduce structure reasonably well. Elastic constants and volume expansion results are analyzed, and the technical and conceptual challenges of reproducing them are explained. Force matching exhibits promising results for capturing the swing effect. This is the first time these CG methods, widely applied in polymer and biomolecule communities, are deployed to model porous solids. We highlight the challenges of fitting CG force fields for these materials.
Collapse
Affiliation(s)
| | - Rocio Semino
- Sorbonne Université, CNRS, Physico-chimie des Electrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
| |
Collapse
|
9
|
Bekdash R, Yoshida K, Nair MS, Qiu L, Ahdout J, Tsai HY, Uryu K, Soni RK, Huang Y, Ho DD, Yazawa M. Developing inhibitory peptides against SARS-CoV-2 envelope protein. PLoS Biol 2024; 22:e3002522. [PMID: 38483887 PMCID: PMC10939250 DOI: 10.1371/journal.pbio.3002522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2023] [Accepted: 01/25/2024] [Indexed: 03/17/2024] Open
Abstract
Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has affected approximately 800 million people since the start of the Coronavirus Disease 2019 (COVID-19) pandemic. Because of the high rate of mutagenesis in SARS-CoV-2, it is difficult to develop a sustainable approach for prevention and treatment. The Envelope (E) protein is highly conserved among human coronaviruses. Previous studies reported that SARS-CoV-1 E deficiency reduced viral propagation, suggesting that E inhibition might be an effective therapeutic strategy for SARS-CoV-2. Here, we report inhibitory peptides against SARS-CoV-2 E protein named iPep-SARS2-E. Leveraging E-induced alterations in proton homeostasis and NFAT/AP-1 pathway in mammalian cells, we developed screening platforms to design and optimize the peptides that bind and inhibit E protein. Using Vero-E6 cells, human-induced pluripotent stem cell-derived branching lung organoid and mouse models with SARS-CoV-2, we found that iPep-SARS2-E significantly inhibits virus egress and reduces viral cytotoxicity and propagation in vitro and in vivo. Furthermore, the peptide can be customizable for E protein of other human coronaviruses such as Middle East Respiratory Syndrome Coronavirus (MERS-CoV). The results indicate that E protein can be a potential therapeutic target for human coronaviruses.
Collapse
Affiliation(s)
- Ramsey Bekdash
- Department of Rehabilitation and Regenerative Medicine, Columbia University, New York, New York, United States of America
- Columbia Stem Cell Initiative, Columbia University, New York, New York, United States of America
- Department of Pharmacology, Columbia University, New York, New York, United States of America
| | - Kazushige Yoshida
- Department of Rehabilitation and Regenerative Medicine, Columbia University, New York, New York, United States of America
- Columbia Stem Cell Initiative, Columbia University, New York, New York, United States of America
| | - Manoj S. Nair
- Aaron Diamond AIDS Research Center, Columbia University, New York, New York, United States of America
| | - Lauren Qiu
- Department of Rehabilitation and Regenerative Medicine, Columbia University, New York, New York, United States of America
- Columbia Stem Cell Initiative, Columbia University, New York, New York, United States of America
- Department of Biological Science, Columbia University, New York, New York, United States of America
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
| | - Johnathan Ahdout
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
| | - Hsiang-Yi Tsai
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
| | - Kunihiro Uryu
- EMSCOPIC, New York, New York, United States of America
| | - Rajesh K. Soni
- Proteomics and Macromolecular Crystallography Shared Resource, Columbia University, New York, New York, United States of America
| | - Yaoxing Huang
- Aaron Diamond AIDS Research Center, Columbia University, New York, New York, United States of America
| | - David D. Ho
- Aaron Diamond AIDS Research Center, Columbia University, New York, New York, United States of America
- Department of Microbiology and Immunology, Columbia University, New York, New York, United States of America
- Division of Infectious Diseases, Department of Medicine, Columbia University, New York, New York, United States of America
| | - Masayuki Yazawa
- Department of Rehabilitation and Regenerative Medicine, Columbia University, New York, New York, United States of America
- Columbia Stem Cell Initiative, Columbia University, New York, New York, United States of America
- Department of Pharmacology, Columbia University, New York, New York, United States of America
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
| |
Collapse
|
10
|
Samsudin F, Zuzic L, Marzinek JK, Bond PJ. Mechanisms of allostery at the viral surface through the eyes of molecular simulation. Curr Opin Struct Biol 2024; 84:102761. [PMID: 38142635 DOI: 10.1016/j.sbi.2023.102761] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 11/13/2023] [Accepted: 12/04/2023] [Indexed: 12/26/2023]
Abstract
The outermost surface layer of any virus is formed by either a capsid shell or envelope. Such layers have traditionally been thought of as immovable structures, but it is becoming apparent that they cannot be viewed exclusively as static architectures protecting the viral genome. A limited number of proteins on the virion surface must perform a multitude of functions in order to orchestrate the viral life cycle, and allostery can regulate their structures at multiple levels of organization, spanning individual molecules, protomers, large oligomeric assemblies, or entire viral surfaces. Here, we review recent contributions from the molecular simulation field to viral surface allostery, with a particular focus on the trimeric spike glycoprotein emerging from the coronavirus surface, and the icosahedral flaviviral envelope complex. As emerging viral pathogens continue to pose a global threat, an improved understanding of viral dynamics and allosteric regulation will prove crucial in developing novel therapeutic strategies.
Collapse
Affiliation(s)
- Firdaus Samsudin
- Bioinformatics Institute (A∗STAR), 30 Biopolis Street, 07-01 Matrix, 138671, Singapore
| | - Lorena Zuzic
- Department of Chemistry, Langelandsgade 140, Aarhus University, Aarhus 8000, Denmark
| | - Jan K Marzinek
- Bioinformatics Institute (A∗STAR), 30 Biopolis Street, 07-01 Matrix, 138671, Singapore
| | - Peter J Bond
- Bioinformatics Institute (A∗STAR), 30 Biopolis Street, 07-01 Matrix, 138671, Singapore; Department of Biological Sciences, 16 Science Drive 4, National University of Singapore, 117558, Singapore.
| |
Collapse
|
11
|
Pezeshkian W, Ipsen JH. Mesoscale simulation of biomembranes with FreeDTS. Nat Commun 2024; 15:548. [PMID: 38228588 PMCID: PMC10792169 DOI: 10.1038/s41467-024-44819-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2023] [Accepted: 01/05/2024] [Indexed: 01/18/2024] Open
Abstract
We present FreeDTS software for performing computational research on biomembranes at the mesoscale. In this software, a membrane is represented by a dynamically triangulated surface equipped with vertex-based inclusions to integrate the effects of integral and peripheral membrane proteins. Several algorithms are included in the software to simulate complex membranes at different conditions such as framed membranes with constant tension, vesicles and high-genus membranes with various fixed volumes or constant pressure differences and applying external forces to membrane regions. Furthermore, the software allows the user to turn off the shape evolution of the membrane and focus solely on the organization of proteins. As a result, we can take realistic membrane shapes obtained from, for example, cryo-electron tomography and backmap them into a finer simulation model. In addition to many biomembrane applications, this software brings us a step closer to simulating realistic biomembranes with molecular resolution. Here we provide several interesting showcases of the power of the software but leave a wide range of potential applications for interested users.
Collapse
Affiliation(s)
- Weria Pezeshkian
- Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100, Copenhagen, Denmark.
| | - John H Ipsen
- MEMPHYS/PhyLife, Department of Physics, Chemistry and Pharmacy (FKF), University of Southern Denmark, Campusvej 55, 5230, Odense M, Denmark
| |
Collapse
|
12
|
Nonn A, Kiss B, Pezeshkian W, Tancogne-Dejean T, Cerrone A, Kellermayer M, Bai Y, Li W, Wierzbicki T. Inferring mechanical properties of the SARS-CoV-2 virus particle with nano-indentation tests and numerical simulations. J Mech Behav Biomed Mater 2023; 148:106153. [PMID: 37865016 DOI: 10.1016/j.jmbbm.2023.106153] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 09/14/2023] [Accepted: 09/27/2023] [Indexed: 10/23/2023]
Abstract
The pandemic caused by the SARS-CoV-2 virus has claimed more than 6.5 million lives worldwide. This global challenge has led to accelerated development of highly effective vaccines tied to their ability to elicit a sustained immune response. While numerous studies have focused primarily on the spike (S) protein, less is known about the interior of the virus. Here we propose a methodology that combines several experimental and simulation techniques to elucidate the internal structure and mechanical properties of the SARS-CoV-2 virus. The mechanical response of the virus was analyzed by nanoindentation tests using a novel flat indenter and evaluated in comparison to a conventional sharp tip indentation. The elastic properties of the viral membrane were estimated by analytical solutions, molecular dynamics (MD) simulations on a membrane patch and by a 3D Finite Element (FE)-beam model of the virion's spike protein and membrane molecular structure. The FE-based inverse engineering approach provided a reasonable reproduction of the mechanical response of the virus from the sharp tip indentation and was successfully verified against the flat tip indentation results. The elastic modulus of the viral membrane was estimated in the range of 7-20 MPa. MD simulations showed that the presence of proteins significantly reduces the fracture strength of the membrane patch. However, FE simulations revealed an overall high fracture strength of the virus, with a mechanical behavior similar to the highly ductile behavior of engineering metallic materials. The failure mechanics of the membrane during sharp tip indentation includes progressive damage combined with localized collapse of the membrane due to severe bending. Furthermore, the results support the hypothesis of a close association of the long membrane proteins (M) with membrane-bound hexagonally packed ribonucleoproteins (RNPs). Beyond improved understanding of coronavirus structure, the present findings offer a knowledge base for the development of novel prevention and treatment methods that are independent of the immune system.
Collapse
Affiliation(s)
- Aida Nonn
- CMM Lab, Faculty of Mechanical Engineering, OTH Regensburg, 93053, Regensburg, Germany.
| | - Bálint Kiss
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, H-1094, Hungary; ELKH-SE Biophysical Virology Research Group, Budapest, H-1094, Hungary
| | - Weria Pezeshkian
- Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100, Copenhagen, Denmark
| | | | - Albert Cerrone
- Computational Hydraulics Laboratory, University of Notre Dame, Notre Dame, IN, 46556, USA
| | - Miklos Kellermayer
- Department of Biophysics and Radiation Biology, Semmelweis University, Budapest, H-1094, Hungary; ELKH-SE Biophysical Virology Research Group, Budapest, H-1094, Hungary
| | - Yuanli Bai
- Department of Mechanical and Aerospace of Engineering, University of Central Florida, 4000 Central Florida Blvd., Orlando, FL, 32816, USA
| | - Wei Li
- Impact and Crashworthiness Lab, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Tomasz Wierzbicki
- Impact and Crashworthiness Lab, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| |
Collapse
|
13
|
Kovalenko I, Kholina E, Fedorov V, Khruschev S, Vasyuchenko E, Meerovich G, Strakhovskaya M. Interaction of Methylene Blue with Severe Acute Respiratory Syndrome Coronavirus 2 Envelope Revealed by Molecular Modeling. Int J Mol Sci 2023; 24:15909. [PMID: 37958892 PMCID: PMC10650479 DOI: 10.3390/ijms242115909] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2023] [Revised: 10/24/2023] [Accepted: 10/31/2023] [Indexed: 11/15/2023] Open
Abstract
Methylene blue has multiple antiviral properties against Severe Acute Respiratory Syndrome-related Coronavirus 2 (SARS-CoV-2). The ability of methylene blue to inhibit different stages of the virus life cycle, both in light-independent and photodynamic processes, is used in clinical practice. At the same time, the molecular aspects of the interactions of methylene blue with molecular components of coronaviruses are not fully understood. Here, we use Brownian dynamics to identify methylene blue binding sites on the SARS-CoV-2 envelope. The local lipid and protein composition of the coronavirus envelope plays a crucial role in the binding of this cationic dye. Viral structures targeted by methylene blue include the S and E proteins and negatively charged lipids. We compare the obtained results with known experimental data on the antiviral effects of methylene blue to elucidate the molecular basis of its activity against coronaviruses.
Collapse
Affiliation(s)
- Ilya Kovalenko
- Faculty of Biology, Lomonosov Moscow State University, Moscow 119234, Russia; (I.K.); (E.K.); (V.F.); (S.K.); (E.V.)
- Scientific and Educational Mathematical Center «Sofia Kovalevskaya Northwestern Center for Mathematical Research», Pskov State University, Pskov 180000, Russia
| | - Ekaterina Kholina
- Faculty of Biology, Lomonosov Moscow State University, Moscow 119234, Russia; (I.K.); (E.K.); (V.F.); (S.K.); (E.V.)
| | - Vladimir Fedorov
- Faculty of Biology, Lomonosov Moscow State University, Moscow 119234, Russia; (I.K.); (E.K.); (V.F.); (S.K.); (E.V.)
| | - Sergei Khruschev
- Faculty of Biology, Lomonosov Moscow State University, Moscow 119234, Russia; (I.K.); (E.K.); (V.F.); (S.K.); (E.V.)
| | - Ekaterina Vasyuchenko
- Faculty of Biology, Lomonosov Moscow State University, Moscow 119234, Russia; (I.K.); (E.K.); (V.F.); (S.K.); (E.V.)
| | - Gennady Meerovich
- Prokhorov General Physics Institute of the Russian Academy of Sciences, Moscow 119991, Russia
- Institute for Physics and Engineering in Biomedicine, National Research Nuclear University “MEPHI”, Moscow 115409, Russia
| | - Marina Strakhovskaya
- Faculty of Biology, Lomonosov Moscow State University, Moscow 119234, Russia; (I.K.); (E.K.); (V.F.); (S.K.); (E.V.)
| |
Collapse
|
14
|
Borges-Araújo L, Patmanidis I, Singh AP, Santos LHS, Sieradzan AK, Vanni S, Czaplewski C, Pantano S, Shinoda W, Monticelli L, Liwo A, Marrink SJ, Souza PCT. Pragmatic Coarse-Graining of Proteins: Models and Applications. J Chem Theory Comput 2023; 19:7112-7135. [PMID: 37788237 DOI: 10.1021/acs.jctc.3c00733] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.
Collapse
Affiliation(s)
- Luís Borges-Araújo
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Ilias Patmanidis
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Akhil P Singh
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
| | - Lucianna H S Santos
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Adam K Sieradzan
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Stefano Vanni
- Department of Biology, University of Fribourg, Chemin du Musée 10, Fribourg CH-1700, Switzerland
- Institut de Pharmacologie Moléculaire et Cellulaire, Université Côte d'Azur, Inserm, CNRS, 06560 Valbonne, France
| | - Cezary Czaplewski
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Sergio Pantano
- Biomolecular Simulations Group, Institut Pasteur de Montevideo, Montevideo 11400, Uruguay
| | - Wataru Shinoda
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-naka, Kita, Okayama 700-8530, Japan
| | - Luca Monticelli
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, Wita Stwosza 63, 80-308 Gdansk, Poland
| | - Siewert J Marrink
- Groningen Biomolecular Sciences and Biotechnology Institute and Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands
| | - Paulo C T Souza
- Molecular Microbiology and Structural Biochemistry (MMSB, UMR 5086), CNRS, University of Lyon, 7 Passage du Vercors, 69007 Lyon, France
| |
Collapse
|
15
|
Ren M, Ma Z, Zhao L, Wang Y, An H, Sun F. Self-Association of ACE-2 with Different RBD Amounts: A Dynamic Simulation Perspective on SARS-CoV-2 Infection. J Chem Inf Model 2023; 63:4423-4432. [PMID: 37382878 DOI: 10.1021/acs.jcim.3c00041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/30/2023]
Abstract
Transmissibility of SARS-CoV-2 initially relies on its trimeric Spike-RBDs to tether the ACE-2 on host cells, and enhanced self-association of ACE-2 engaged with Spike facilitates the viral infection. Two primary packing modes of Spike-ACE2 heteroproteins exist potentially due to discrepant amounts of RBDs loading on ACE-2, but the resultant self-association difference is inherently unclear. We used extensive coarse-grained dynamic simulations to characterize the self-association efficiency, the conformation relevance, and the molecular mechanism of ACE-2 with different RBD amounts. It was revealed that the ACE-2 hanging two/full RBDs (Mode-A) rapidly dimerized into the heteroprotein complex in a compact "linear" conformation, while the bare ACE-2 showed weakened self-association and a protein complex. The RBD-tethered ectodomains of ACE-2 presented a more upright conformation relative to the membrane, and the intermolecular ectodomains were predominantly packed by the neck domains, which was obligated to the rapid protein self-association in a compact pattern. Noted is the fact that the ACE-2 tethered by a single RBD (Mode-B) retained considerable self-association efficiency and clustering capability, which unravels the interrelation of ACE-2 colocalization and protein cross-linkage. The molecular perspectives in this study expound the self-association potency of ACE-2 with different RBD amounts and the viral activity implications, which can greatly enhance our comprehension of SARS-CoV-2 infection details.
Collapse
Affiliation(s)
- Meina Ren
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin 300401, China
| | - Ziyi Ma
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin 300401, China
| | - Lina Zhao
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin 300401, China
| | - Yanjiao Wang
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin 300401, China
| | - Hailong An
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin 300401, China
| | - Fude Sun
- Key Laboratory of Molecular Biophysics, Hebei Province, Institute of Biophysics, School of Health Science & Biomedical Engineering, Hebei University of Technology, Tianjin 300401, China
| |
Collapse
|
16
|
Lynch DL, Pavlova A, Fan Z, Gumbart JC. Understanding Virus Structure and Dynamics through Molecular Simulations. J Chem Theory Comput 2023. [PMID: 37192279 DOI: 10.1021/acs.jctc.3c00116] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
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.
Collapse
Affiliation(s)
- Diane L Lynch
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Anna Pavlova
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Zixing Fan
- Interdisciplinary Bioengineering Graduate Program, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| |
Collapse
|