1
|
Carou-Senra P, Rodríguez-Pombo L, Awad A, Basit AW, Alvarez-Lorenzo C, Goyanes A. Inkjet Printing of Pharmaceuticals. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2024; 36:e2309164. [PMID: 37946604 DOI: 10.1002/adma.202309164] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Revised: 10/23/2023] [Indexed: 11/12/2023]
Abstract
Inkjet printing (IJP) is an additive manufacturing process that selectively deposits ink materials, layer-by-layer, to create 3D objects or 2D patterns with precise control over their structure and composition. This technology has emerged as an attractive and versatile approach to address the ever-evolving demands of personalized medicine in the healthcare industry. Although originally developed for nonhealthcare applications, IJP harnesses the potential of pharma-inks, which are meticulously formulated inks containing drugs and pharmaceutical excipients. Delving into the formulation and components of pharma-inks, the key to precise and adaptable material deposition enabled by IJP is unraveled. The review extends its focus to substrate materials, including paper, films, foams, lenses, and 3D-printed materials, showcasing their diverse advantages, while exploring a wide spectrum of therapeutic applications. Additionally, the potential benefits of hardware and software improvements, along with artificial intelligence integration, are discussed to enhance IJP's precision and efficiency. Embracing these advancements, IJP holds immense potential to reshape traditional medicine manufacturing processes, ushering in an era of medical precision. However, further exploration and optimization are needed to fully utilize IJP's healthcare capabilities. As researchers push the boundaries of IJP, the vision of patient-specific treatment is on the horizon of becoming a tangible reality.
Collapse
Affiliation(s)
- Paola Carou-Senra
- Departamento de Farmacología, Farmacia y Tecnología Farmacéutica, I+D Farma Group (GI-1645), Facultad de Farmacia, Instituto de Materiales (iMATUS) and Health Research Institute of Santiago de Compostela (IDIS), Universidade de Santiago de Compostela, Santiago de Compostela, 15782, Spain
| | - Lucía Rodríguez-Pombo
- Departamento de Farmacología, Farmacia y Tecnología Farmacéutica, I+D Farma Group (GI-1645), Facultad de Farmacia, Instituto de Materiales (iMATUS) and Health Research Institute of Santiago de Compostela (IDIS), Universidade de Santiago de Compostela, Santiago de Compostela, 15782, Spain
| | - Atheer Awad
- Department of Clinical, Pharmaceutical and Biological Sciences, University of Hertfordshire, College Lane, Hatfield, AL10 9AB, UK
| | - Abdul W Basit
- Department of Pharmaceutics, UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London, WC1N 1AX, UK
- FABRX Ltd., Henwood House, Henwood, Ashford, Kent, TN24 8DH, UK
- FABRX Artificial Intelligence, Carretera de Escairón 14, Currelos (O Saviñao), CP 27543, Spain
| | - Carmen Alvarez-Lorenzo
- Departamento de Farmacología, Farmacia y Tecnología Farmacéutica, I+D Farma Group (GI-1645), Facultad de Farmacia, Instituto de Materiales (iMATUS) and Health Research Institute of Santiago de Compostela (IDIS), Universidade de Santiago de Compostela, Santiago de Compostela, 15782, Spain
| | - Alvaro Goyanes
- Departamento de Farmacología, Farmacia y Tecnología Farmacéutica, I+D Farma Group (GI-1645), Facultad de Farmacia, Instituto de Materiales (iMATUS) and Health Research Institute of Santiago de Compostela (IDIS), Universidade de Santiago de Compostela, Santiago de Compostela, 15782, Spain
- Department of Pharmaceutics, UCL School of Pharmacy, University College London, 29-39 Brunswick Square, London, WC1N 1AX, UK
- FABRX Ltd., Henwood House, Henwood, Ashford, Kent, TN24 8DH, UK
- FABRX Artificial Intelligence, Carretera de Escairón 14, Currelos (O Saviñao), CP 27543, Spain
| |
Collapse
|
2
|
Giron CC, Laaksonen A, Barroso da Silva FL. Differences between Omicron SARS-CoV-2 RBD and other variants in their ability to interact with cell receptors and monoclonal antibodies. J Biomol Struct Dyn 2023; 41:5707-5727. [PMID: 35815535 DOI: 10.1080/07391102.2022.2095305] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2022] [Accepted: 06/23/2022] [Indexed: 12/23/2022]
Abstract
SARS-CoV-2 remains a health threat with the continuous emergence of new variants. This work aims to expand the knowledge about the SARS-CoV-2 receptor-binding domain (RBD) interactions with cell receptors and monoclonal antibodies (mAbs). By using constant-pH Monte Carlo simulations, the free energy of interactions between the RBD from different variants and several partners (Angiotensin-Converting Enzyme-2 (ACE2) polymorphisms and various mAbs) were predicted. Computed RBD-ACE2-binding affinities were higher for two ACE2 polymorphisms (rs142984500 and rs4646116) typically found in Europeans which indicates a genetic susceptibility. This is amplified for Omicron (BA.1) and its sublineages BA.2 and BA.3. The antibody landscape was computationally investigated with the largest set of mAbs so far in the literature. From the 32 studied binders, groups of mAbs were identified from weak to strong binding affinities (e.g. S2K146). These mAbs with strong binding capacity and especially their combination are amenable to experimentation and clinical trials because of their high predicted binding affinities and possible neutralization potential for current known virus mutations and a universal coronavirus.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Carolina Corrêa Giron
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, SP, Brazil
- Universidade Federal do Triângulo Mineiro, Hospital de Clínicas, Uberaba, MG, Brazil
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, Stockholm, Sweden
- State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing, PR China
- Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Iasi, Romania
- Department of Engineering Sciences and Mathematics, Division of Energy Science, Luleå University of Technology, Luleå, Sweden
- Department of Chemical and Geological Sciences, University of Cagliari, Monserrato, Italy
| | - Fernando Luís Barroso da Silva
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, SP, Brazil
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
3
|
Blanco PM, Narambuena CF, Madurga S, Mas F, Garcés JL. Unusual Aspects of Charge Regulation in Flexible Weak Polyelectrolytes. Polymers (Basel) 2023; 15:2680. [PMID: 37376324 DOI: 10.3390/polym15122680] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2023] [Revised: 06/07/2023] [Accepted: 06/10/2023] [Indexed: 06/29/2023] Open
Abstract
This article reviews the state of the art of the studies on charge regulation (CR) effects in flexible weak polyelectrolytes (FWPE). The characteristic of FWPE is the strong coupling of ionization and conformational degrees of freedom. After introducing the necessary fundamental concepts, some unconventional aspects of the the physical chemistry of FWPE are discussed. These aspects are: (i) the extension of statistical mechanics techniques to include ionization equilibria and, in particular, the use of the recently proposed Site Binding-Rotational Isomeric State (SBRIS) model, which allows the calculation of ionization and conformational properties on the same foot; (ii) the recent progresses in the inclusion of proton equilibria in computer simulations; (iii) the possibility of mechanically induced CR in the stretching of FWPE; (iv) the non-trivial adsorption of FWPE on ionized surfaces with the same charge sign as the PE (the so-called "wrong side" of the isoelectric point); (v) the influence of macromolecular crowding on CR.
Collapse
Affiliation(s)
- Pablo M Blanco
- Physical Chemistry Unit, Materials Science and Physical Chemistry Department & Research Institute of Theoretical and Computational Chemistry (IQTCUB), Barcelona University (UB), 08028 Barcelona, Catalonia, Spain
| | - Claudio F Narambuena
- Grupo de Bionanotecnologia y Sistemas Complejos, Infap-CONICET & Facultad Regional San Rafael, Universidad Tecnológica Nacional, San Rafael 5600, Argentina
| | - Sergio Madurga
- Physical Chemistry Unit, Materials Science and Physical Chemistry Department & Research Institute of Theoretical and Computational Chemistry (IQTCUB), Barcelona University (UB), 08028 Barcelona, Catalonia, Spain
| | - Francesc Mas
- Physical Chemistry Unit, Materials Science and Physical Chemistry Department & Research Institute of Theoretical and Computational Chemistry (IQTCUB), Barcelona University (UB), 08028 Barcelona, Catalonia, Spain
| | - Josep L Garcés
- Chemistry Department, Technical School of Agricultural Engineering & AGROTECNIO, Lleida University (UdL), 25003 Lleida, Catalonia, Spain
| |
Collapse
|
4
|
Neamtu A, Mocci F, Laaksonen A, Barroso da Silva FL. Towards an optimal monoclonal antibody with higher binding affinity to the receptor-binding domain of SARS-CoV-2 spike proteins from different variants. Colloids Surf B Biointerfaces 2023; 221:112986. [PMID: 36375294 PMCID: PMC9617679 DOI: 10.1016/j.colsurfb.2022.112986] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Revised: 09/13/2022] [Accepted: 10/27/2022] [Indexed: 11/13/2022]
Abstract
A highly efficient and robust multiple scales in silico protocol, consisting of atomistic Molecular Dynamics (MD), coarse-grain (CG) MD, and constant-pH CG Monte Carlo (MC), has been developed and used to study the binding affinities of selected antigen-binding fragments of the monoclonal antibody (mAbs) CR3022 and several of its here optimized versions against 11 SARS-CoV-2 variants including the wild type. Totally 235,000 mAbs structures were initially generated using the RosettaAntibodyDesign software, resulting in top 10 scored CR3022-like-RBD complexes with critical mutations and compared to the native one, all having the potential to block virus-host cell interaction. Of these 10 finalists, two candidates were further identified in the CG simulations to be the best against all SARS-CoV-2 variants. Surprisingly, all 10 candidates and the native CR3022 exhibited a higher affinity for the Omicron variant despite its highest number of mutations. The multiscale protocol gives us a powerful rational tool to design efficient mAbs. The electrostatic interactions play a crucial role and appear to be controlling the affinity and complex building. Studied mAbs carrying a more negative total net charge show a higher affinity. Structural determinants could be identified in atomistic simulations and their roles are discussed in detail to further hint at a strategy for designing the best RBD binder. Although the SARS-CoV-2 was specifically targeted in this work, our approach is generally suitable for many diseases and viral and bacterial pathogens, leukemia, cancer, multiple sclerosis, rheumatoid, arthritis, lupus, and more.
Collapse
Affiliation(s)
- Andrei Neamtu
- Department of Physiology, "Grigore T. Popa" University of Medicine and Pharmacy of Iasi, Str. Universitatii nr. 16, 700051 Iasi, România; TRANSCEND Centre - Regional Institute of Oncology (IRO) Iasi, Str. General Henri Mathias Berthelot, Nr. 2-4 Iași, România
| | - Francesca Mocci
- University of Cagliari, Department of Chemical and Geological Sciences, Campus Monserrato, SS 554 bivio per Sestu, 09042 Monserrato, Italy
| | - Aatto Laaksonen
- Centre of Advanced Research in Bionanoconjugates and Biopolymers, PetruPoni Institute of Macromolecular Chemistry Aleea Grigore Ghica-Voda, 41 A, 700487 Iasi, Romania; University of Cagliari, Department of Chemical and Geological Sciences, Campus Monserrato, SS 554 bivio per Sestu, 09042 Monserrato, Italy; Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden; State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing 210009, PR China; Department of Engineering Sciences and Mathematics, Division of Energy Science, Luleå University of Technology, SE-97187 Luleå, Sweden
| | - Fernando L Barroso da Silva
- Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. café, s/no - campus da USP, BR-14040-903 Ribeirão Preto, SP, Brazil; Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA.
| |
Collapse
|
5
|
Barroso da Silva FL, Giron CC, Laaksonen A. Electrostatic Features for the Receptor Binding Domain of SARS-COV-2 Wildtype and Its Variants. Compass to the Severity of the Future Variants with the Charge-Rule. J Phys Chem B 2022; 126:6835-6852. [PMID: 36066414 DOI: 10.1021/acs.jpcb.2c04225] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Electrostatic intermolecular interactions are important in many aspects of biology. We have studied the main electrostatic features involved in the interaction of the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein with the human receptor Angiotensin-converting enzyme 2 (ACE2). As the principal computational tool, we have used the FORTE approach, capable to model proton fluctuations and computing free energies for a very large number of protein-protein systems under different physical-chemical conditions, here focusing on the RBD-ACE2 interactions. Both the wild-type and all critical variants are included in this study. From our large ensemble of extensive simulations, we obtain, as a function of pH, the binding affinities, charges of the proteins, their charge regulation capacities, and their dipole moments. In addition, we have calculated the pKas for all ionizable residues and mapped the electrostatic coupling between them. We are able to present a simple predictor for the RBD-ACE2 binding based on the data obtained for Alpha, Beta, Gamma, Delta, and Omicron variants, as a linear correlation between the total charge of the RBD and the corresponding binding affinity. This "RBD charge rule" should work as a quick test of the degree of severity of the coming SARS-CoV-2 variants in the future.
Collapse
Affiliation(s)
- Fernando L Barroso da Silva
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. café, s/no-campus da USP, BR-14040-903 Ribeirão Preto, SP, Brazil.,Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Carolina Corrêa Giron
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. café, s/no-campus da USP, BR-14040-903 Ribeirão Preto, SP, Brazil.,Hospital de Clínicas, Universidade Federal do Triângulo Mineiro, Av. Getúlio Guaritá, 38025-440 Uberaba, MG, Brazil
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden.,State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing, 210009, P. R. China.,Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Aleea Grigore Ghica-Voda, 41A, 700487 Iasi, Romania.,Department of Engineering Sciences and Mathematics, Division of Energy Science, Luleå University of Technology, SE-97187 Luleå, Sweden.,Department of Chemical and Geological Sciences, Campus Monserrato, University of Cagliari, SS 554 bivio per Sestu, 09042 Monserrato, Italy
| |
Collapse
|
6
|
Poveda-Cuevas SA, Etchebest C, da Silva FLB. Self-association features of NS1 proteins from different flaviviruses. Virus Res 2022; 318:198838. [PMID: 35662566 DOI: 10.1016/j.virusres.2022.198838] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Revised: 05/26/2022] [Accepted: 05/27/2022] [Indexed: 12/22/2022]
Abstract
Flaviviruses comprise a large group of arboviral species that are distributed in several countries of the tropics, neotropics, and some temperate zones. Since they can produce neurological pathologies or vascular damage, there has been intense research seeking better diagnosis and treatments for their infections in the last decades. The flavivirus NS1 protein is a relevant clinical target because it is involved in viral replication, immune evasion, and virulence. Being a key factor in endothelial and tissue-specific modulation, NS1 has been largely studied to understand the molecular mechanisms exploited by the virus to reprogram host cells. A central part of the viral maturation processes is the NS1 oligomerization because many stages rely on these protein-protein assemblies. In the present study, the self-associations of NS1 proteins from Zika, Dengue, and West Nile viruses are examined through constant-pH coarse-grained biophysical simulations. Free energies of interactions were estimated for different oligomeric states and pH conditions. Our results show that these proteins can form both dimers and tetramers under conditions near physiological pH even without the presence of lipids. Moreover, pH plays an important role mainly controlling the regimes where van der Waals interactions govern their association. Finally, despite the similarity at the sequence level, we found that each flavivirus has a well-characteristic protein-protein interaction profile. These specific features can provide new hints for the development of binders both for better diagnostic tools and the formulation of new therapeutic drugs.
Collapse
Affiliation(s)
- Sergio A Poveda-Cuevas
- Universidade de São Paulo, Programa Interunidades em Bioinformática, Rua do Matão, 1010, BR-05508-090 São Paulo, São Paulo, Brazil; Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. do Café, s/no-Campus da USP, BR-14040-903 Ribeirão Preto, São Paulo, Brazil; University of São Paulo and Université de Paris International Laboratory in Structural Bioinformatics, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. do Café, s/no-Campus da USP, Bloco B, BR-14040-903 Ribeirão Preto, São Paulo, Brazil.; Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, United States
| | - Catherine Etchebest
- Université Paris Cité, Biologie Intégrée du Globule Rouge, Equipe 2, INSERM, F-75015 Paris, France; University of São Paulo and Université de Paris International Laboratory in Structural Bioinformatics, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. do Café, s/no-Campus da USP, Bloco B, BR-14040-903 Ribeirão Preto, São Paulo, Brazil
| | - Fernando L Barroso da Silva
- Universidade de São Paulo, Programa Interunidades em Bioinformática, Rua do Matão, 1010, BR-05508-090 São Paulo, São Paulo, Brazil; Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. do Café, s/no-Campus da USP, BR-14040-903 Ribeirão Preto, São Paulo, Brazil; University of São Paulo and Université de Paris International Laboratory in Structural Bioinformatics, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. do Café, s/no-Campus da USP, Bloco B, BR-14040-903 Ribeirão Preto, São Paulo, Brazil..
| |
Collapse
|
7
|
Lunkad R, Barroso da Silva FL, Košovan P. Both Charge-Regulation and Charge-Patch Distribution Can Drive Adsorption on the Wrong Side of the Isoelectric Point. J Am Chem Soc 2022; 144:1813-1825. [DOI: 10.1021/jacs.1c11676] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Raju Lunkad
- Department of Physical and Macromolecular Chemistry, Charles University, Hlavova 8, 128 43 Prague, Czech Republic
| | - Fernando L. Barroso da Silva
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
- Department of Biomolecular Sciences, School of Pharmaceutical Sciences at Ribeirão Preto, University of São Paulo, Ribeirão Preto, São Paulo 14040-900, Brazil
| | - Peter Košovan
- Department of Physical and Macromolecular Chemistry, Charles University, Hlavova 8, 128 43 Prague, Czech Republic
| |
Collapse
|
8
|
Reilley DJ, Wang J, Dokholyan NV, Alexandrova AN. Titr-DMD-A Rapid, Coarse-Grained Quasi-All-Atom Constant pH Molecular Dynamics Framework. J Chem Theory Comput 2021; 17:4538-4549. [PMID: 34165292 PMCID: PMC10662685 DOI: 10.1021/acs.jctc.1c00338] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The pH-dependence of enzyme fold stability and catalytic activity is a fundamentally dynamic, structural property which is difficult to study. The challenges and expense of investigating dynamic, atomic scale behavior experimentally means that computational methods, particularly constant pH molecular dynamics (CpHMD), are well situated tools for this. However, these methods often struggle with affordable sampling of sufficiently long time scales while also obtaining accurate pKa prediction and verifying the structures they generate. We introduce Titr-DMD, an affordable CpHMD method that combines the quasi-all-atom coarse-grained discrete molecular dynamics (DMD) method for conformational sampling with Propka for pKa prediction, to circumvent these issues. The combination enables rapid sampling on limited computational resources, while simulations are still performed on the atomic scale. We benchmark the method on a set of proteins with experimentally attested pKa and on the pH triggered conformational change in a staphylococcal nuclease mutant, a rare experimental study of such behavior. Our results show Titr-DMD to be an effective and inexpensive method to study pH-coupled protein dynamics.
Collapse
Affiliation(s)
- David J Reilley
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Jian Wang
- Department of Pharmacology, Department of Biochemistry and Molecular Biology, Penn State University College of Medicine, Hershey, Pennsylvania 17033, United States
| | - Nikolay V Dokholyan
- Department of Pharmacology, Department of Biochemistry and Molecular Biology, Penn State University College of Medicine, Hershey, Pennsylvania 17033, United States
- Departments of Chemistry and Biomedical Engineering, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
- California NanoSystems Institute, Los Angeles, California 90095-1569, United States
| |
Collapse
|
9
|
Giron CC, Laaksonen A, Barroso da Silva FL. Up State of the SARS-COV-2 Spike Homotrimer Favors an Increased Virulence for New Variants. FRONTIERS IN MEDICAL TECHNOLOGY 2021; 3:694347. [PMID: 35047936 PMCID: PMC8757851 DOI: 10.3389/fmedt.2021.694347] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Accepted: 05/31/2021] [Indexed: 12/23/2022] Open
Abstract
The COVID-19 pandemic has spread worldwide. However, as soon as the first vaccines-the only scientifically verified and efficient therapeutic option thus far-were released, mutations combined into variants of SARS-CoV-2 that are more transmissible and virulent emerged, raising doubts about their efficiency. This study aims to explain possible molecular mechanisms responsible for the increased transmissibility and the increased rate of hospitalizations related to the new variants. A combination of theoretical methods was employed. Constant-pH Monte Carlo simulations were carried out to quantify the stability of several spike trimeric structures at different conformational states and the free energy of interactions between the receptor-binding domain (RBD) and angiotensin-converting enzyme II (ACE2) for the most worrying variants. Electrostatic epitopes were mapped using the PROCEEDpKa method. These analyses showed that the increased virulence is more likely to be due to the improved stability to the S trimer in the opened state, in which the virus can interact with the cellular receptor, ACE2, rather than due to alterations in the complexation RBD-ACE2, since the difference observed in the free energy values was small (although more attractive in general). Conversely, the South African/Beta variant (B.1.351), compared with the SARS-CoV-2 wild type (wt), is much more stable in the opened state with one or two RBDs in the up position than in the closed state with three RBDs in the down position favoring the infection. Such results contribute to understanding the natural history of disease and indicate possible strategies for developing new therapeutic molecules and adjusting the vaccine doses for higher B-cell antibody production.
Collapse
Affiliation(s)
- Carolina Corrêa Giron
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, Brazil
- Hospital de Clínicas, Universidade Federal do Triângulo Mineiro, Uberaba, Brazil
| | - Aatto Laaksonen
- Arrhenius Laboratory, Department of Materials and Environmental Chemistry, Stockholm University, Stockholm, Sweden
- State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing, China
- Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Iasi, Romania
- Division of Energy Science, Department of Engineering Sciences and Mathematics, Luleå University of Technology, Luleå, Sweden
| | - Fernando Luís Barroso da Silva
- Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Ribeirão Preto, Brazil
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, United States
| |
Collapse
|
10
|
Poveda-Cuevas SA, Barroso da Silva FL, Etchebest C. How the Strain Origin of Zika Virus NS1 Protein Impacts Its Dynamics and Implications to Their Differential Virulence. J Chem Inf Model 2021; 61:1516-1530. [PMID: 33651942 DOI: 10.1021/acs.jcim.0c01377] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Viruses can impact and affect human populations in a severe way. The appropriate differentiation among several species or strains of viruses is one of the biggest challenges for virology and infectiology studies. The detection of measurables-quantified discrepancies allows for more accurate clinical diagnoses and treatments for viral diseases. In the present study, we have used a computational approach to explore the dynamical properties of the nonstructural protein 1 from two strains of Zika virus. Our results show that despite a high sequence similarity, the two viral proteins from different origins can exhibit significant dissimilar structural dynamics, which complement their reported differential virulence. The present study opens up new ways in the understanding of the infectivity for these biological entities.
Collapse
Affiliation(s)
- Sergio A Poveda-Cuevas
- Programa Interunidades em Bioinformática, Universidade de São Paulo, Rua do Matão, 1010, BR, 05508-090 São Paulo, São Paulo, Brazil.,Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. do Café, s/no-Campus da USP, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil.,University of São Paulo and Université de Paris International Laboratory in Structural Bioinformatics, Av. do Café, s/no-Campus da USP, Bloco B, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
| | - Fernando Luís Barroso da Silva
- Programa Interunidades em Bioinformática, Universidade de São Paulo, Rua do Matão, 1010, BR, 05508-090 São Paulo, São Paulo, Brazil.,Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. do Café, s/no-Campus da USP, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil.,Department of Chemical and Biomolecular Engineering, North Carolina State University, 27695 Raleigh, North Carolina, United States.,University of São Paulo and Université de Paris International Laboratory in Structural Bioinformatics, Av. do Café, s/no-Campus da USP, Bloco B, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
| | - Catherine Etchebest
- Institut National de la Transfusion Sanguine, 6 Rue Alexandre Cabanel, 75015 Paris, France.,Institut National de la Santé et de la Recherche Médicale, UMR_S 1134, Biologie Intégrée du Globule Rouge, Equipe 2, INSERM, Dynamique des Structures et des Interactions Moléculaires, F-75015 Paris, France.,Université de Paris, 5 Rue Thomas Mann, 75013 Paris, France.,University of São Paulo and Université de Paris International Laboratory in Structural Bioinformatics, Av. do Café, s/no-Campus da USP, Bloco B, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
| |
Collapse
|
11
|
Corrêa Giron C, Laaksonen A, Barroso da Silva FL. On the interactions of the receptor-binding domain of SARS-CoV-1 and SARS-CoV-2 spike proteins with monoclonal antibodies and the receptor ACE2. Virus Res 2020; 285:198021. [PMID: 32416259 PMCID: PMC7228703 DOI: 10.1016/j.virusres.2020.198021] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2020] [Revised: 05/01/2020] [Accepted: 05/11/2020] [Indexed: 01/12/2023]
Abstract
A new betacoronavirus named SARS-CoV-2 has emerged as a new threat to global health and economy. A promising target for both diagnosis and therapeutics treatments of the new disease named COVID-19 is the coronavirus (CoV) spike (S) glycoprotein. By constant-pH Monte Carlo simulations and the PROCEEDpKa method, we have mapped the electrostatic epitopes for four monoclonal antibodies and the angiotensin-converting enzyme 2 (ACE2) on both SARS-CoV-1 and the new SARS-CoV-2 S receptor binding domain (RBD) proteins. We also calculated free energy of interactions and shown that the S RBD proteins from both SARS viruses binds to ACE2 with similar affinities. However, the affinity between the S RBD protein from the new SARS-CoV-2 and ACE2 is higher than for any studied antibody previously found complexed with SARS-CoV-1. Based on physical chemical analysis and free energies estimates, we can shed some light on the involved molecular recognition processes, their clinical aspects, the implications for drug developments, and suggest structural modifications on the CR3022 antibody that would improve its binding affinities for SARS-CoV-2 and contribute to address the ongoing international health crisis.
Collapse
MESH Headings
- Angiotensin-Converting Enzyme 2
- Antibodies, Monoclonal/immunology
- Antibodies, Monoclonal/metabolism
- Antibodies, Viral/immunology
- Antibodies, Viral/metabolism
- Betacoronavirus/chemistry
- Betacoronavirus/immunology
- Computer Simulation
- Epitope Mapping
- Humans
- Models, Molecular
- Monte Carlo Method
- Peptidyl-Dipeptidase A/chemistry
- Peptidyl-Dipeptidase A/metabolism
- Protein Binding
- Protein Conformation
- Protein Interaction Domains and Motifs
- Protein Interaction Mapping
- Receptors, Virus/chemistry
- Receptors, Virus/metabolism
- Severe acute respiratory syndrome-related coronavirus/chemistry
- Severe acute respiratory syndrome-related coronavirus/immunology
- SARS-CoV-2
- Spike Glycoprotein, Coronavirus/chemistry
- Spike Glycoprotein, Coronavirus/immunology
- Spike Glycoprotein, Coronavirus/metabolism
- Thermodynamics
Collapse
Affiliation(s)
- Carolina Corrêa Giron
- Universidade Federal do Triângulo Mineiro, Departamento de Saúde Coletiva, Rua Vigário Carlos, 38025-350 Uberaba, MG, Brazil; Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. café, s/no - campus da USP, BR-14040-903 Ribeirão Preto SP, Brazil
| | - Aatto Laaksonen
- Department of Materials and Environmental Chemistry, Arrhenius Laboratory, Stockholm University, SE-106 91 Stockholm, Sweden; State Key Laboratory of Materials-Oriented and Chemical Engineering, Nanjing Tech University, Nanjing, 210009, PR China; Centre of Advanced Research in Bionanoconjugates and Biopolymers, Petru Poni Institute of Macromolecular Chemistry, Aleea Grigore Ghica-Voda, 41A, 700487 Iasi, Romania; Department of Engineering Sciences and Mathematics, Division of Energy Science, Luleå University of Technology, SE-97187 Luleå, Sweden
| | - Fernando L Barroso da Silva
- Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. café, s/no - campus da USP, BR-14040-903 Ribeirão Preto SP, Brazil; Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, United States.
| |
Collapse
|
12
|
Zanetti-Polzi L, Daidone I, Amadei A. Fully Atomistic Multiscale Approach for p Ka Prediction. J Phys Chem B 2020; 124:4712-4722. [PMID: 32427481 DOI: 10.1021/acs.jpcb.0c01752] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The ionization state of titratable amino acids strongly affects proteins structure and functioning in a large number of biological processes. It is therefore essential to be able to characterize the pKa of ionizable groups inside proteins and to understand its microscopic determinants in order to gain insights into many functional properties of proteins. A big effort has been devoted to the development of theoretical approaches for the prediction of deprotonation free energies, yet the accurate theoretical/computational calculation of pKa values is recognized as a current challenge. A methodology based on a hybrid quantum/classical approach is here proposed for the computation of deprotonation free energies. The method is applied to calculate the pKa of formic acid, methylammonium, and methanethiol, providing results in good agreement with the corresponding experimental estimates. The pKa is also calculated for aspartic acid and lysine as single residues in solution and for three aspartic/glutamic acids inside a well-characterized protein: hen egg white lysozyme. While for small molecules the method is able to deal with multiple protonation states of all titratable groups, this becomes computationally very expensive for proteins. The calculated pKa values for the single amino acids (except for the zwitterionic aspartic acid) and inside the protein display a systematic shift with respect to the experimental values that suggests that the fine balance between hydrophobic and polar interactions might be not accurately reproduced by the usual classical force-fields, thus affecting the computation of deprotonation free energies. The calculated pKa shifts inside the protein are in good agreement with the corresponding experimental ones (within 1 pKa unit), well reproducing the pKa changes due to the protein environment even in the case of large pKa shifts.
Collapse
Affiliation(s)
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, Via Vetoio, I-67010 L'Aquila, Italy
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome "Tor Vergata", Via della Ricerca Scientifica, I-00185 Rome, Italy
| |
Collapse
|
13
|
Frigori RB, Barroso da Silva FL, Carvalho PPD, Alves NA. Occurrence of Biased Conformations as Precursors of Assembly States in Fibril Elongation of Amyloid-β Fibril Variants: An In Silico Study. J Phys Chem B 2020; 124:2798-2805. [DOI: 10.1021/acs.jpcb.0c01360] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Rafael B. Frigori
- Universidade Tecnológica Federal do Paraná, Rua Cristo Rei 19, Toledo 85902-490, Paraná, Brazil
| | - Fernando L. Barroso da Silva
- Departamento de Ciências Biomoleculares, FCFRP, Universidade de São Paulo, Avenida do Café, s/no, Ribeirão Preto 14040-903, São Paulo, Brazil
| | - Patrícia P. D. Carvalho
- Departamento de Fı́sica, FFCLRP, Universidade de São Paulo, Avenida Bandeirantes, 3900, Ribeirão Preto 14040-901, São Paulo, Brazil
| | - Nelson A. Alves
- Departamento de Fı́sica, FFCLRP, Universidade de São Paulo, Avenida Bandeirantes, 3900, Ribeirão Preto 14040-901, São Paulo, Brazil
| |
Collapse
|
14
|
Barroso da Silva FL, Carloni P, Cheung D, Cottone G, Donnini S, Foegeding EA, Gulzar M, Jacquier JC, Lobaskin V, MacKernan D, Mohammad Hosseini Naveh Z, Radhakrishnan R, Santiso EE. Understanding and Controlling Food Protein Structure and Function in Foods: Perspectives from Experiments and Computer Simulations. Annu Rev Food Sci Technol 2020; 11:365-387. [PMID: 31951485 DOI: 10.1146/annurev-food-032519-051640] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The structure and interactions of proteins play a critical role in determining the quality attributes of many foods, beverages, and pharmaceutical products. Incorporating a multiscale understanding of the structure-function relationships of proteins can provide greater insight into, and control of, the relevant processes at play. Combining data from experimental measurements, human sensory panels, and computer simulations through machine learning allows the construction of statistical models relating nanoscale properties of proteins to the physicochemical properties, physiological outcomes, and tastes of foods. This review highlights several examples of advanced computer simulations at molecular, mesoscale, and multiscale levels that shed light on the mechanisms at play in foods, thereby facilitating their control. It includes a practical simulation toolbox for those new to in silico modeling.
Collapse
Affiliation(s)
- Fernando Luís Barroso da Silva
- School of Pharmaceutical Sciences at Ribeirão Preto, University of São Paulo, BR-14040-903, Ribeirão Preto, São Paulo, Brazil
| | - Paolo Carloni
- Institute for Computational Biomedicine (IAS-5/INM-9), Forschungszentrum Jülich, 52425 Jülich, Germany.,Department of Physics, RWTH Aachen University, 52062 Aachen, Germany
| | - David Cheung
- School of Chemistry, National University of Ireland Galway, Galway, Ireland
| | - Grazia Cottone
- Department of Physics and Chemistry, University of Palermo, 90128 Palermo, Italy
| | - Serena Donnini
- Department of Biological and Environmental Science, University of Jyväskylä, Jyväskylä 40014, Finland
| | - E Allen Foegeding
- Department of Food, Bioprocessing, & Nutrition Sciences, North Carolina State University, Raleigh, North Carolina 27695, USA
| | - Muhammad Gulzar
- UCD School of Agriculture and Food Science, University College Dublin, Dublin 4, Ireland
| | | | | | - Donal MacKernan
- UCD School of Physics, University College Dublin, Dublin 4, Ireland
| | | | - Ravi Radhakrishnan
- Department of Bioengineering, University of Pennsylvania, Philadelphia, Pennsylvania, 19104, USA
| | - Erik E Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
15
|
Poveda-Cuevas SA, Etchebest C, Barroso da Silva FL. Identification of Electrostatic Epitopes in Flavivirus by Computer Simulations: The PROCEEDpKa Method. J Chem Inf Model 2019; 60:944-963. [DOI: 10.1021/acs.jcim.9b00895] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Sergio A. Poveda-Cuevas
- Universidade de São Paulo, Programa Interunidades em Bioinformática, Rua do Matão, 1010, BR, 05508-090 São Paulo, São Paulo, Brazil
- Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. Café, s/no−Campus da USP, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
- University of São Paulo-Université Sorbonne Paris Cité International Laboratory in Structural Bioinformatics, Av. do Café, s/no−FCFRP, Bloco B, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
| | - Catherine Etchebest
- Université de Paris, Biologie Intégrée du Globule Rouge, UMR_S1134, BIGR, INSERM, F-75015 Paris, France
- Equipe 2, Dynamique des Structures et des Interactions Moléculaires, Université Paris Diderot−Paris 7, INTS, 6 Rue Alexandre Cabanel, 75015 Paris, France
- Laboratoire d’Excellence GR-Ex, Paris, France
- University of São Paulo-Université Sorbonne Paris Cité International Laboratory in Structural Bioinformatics, Av. do Café, s/no−FCFRP, Bloco B, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
| | - Fernando L. Barroso da Silva
- Universidade de São Paulo, Programa Interunidades em Bioinformática, Rua do Matão, 1010, BR, 05508-090 São Paulo, São Paulo, Brazil
- Universidade de São Paulo, Departamento de Ciências Biomoleculares, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. Café, s/no−Campus da USP, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
- University of São Paulo-Université Sorbonne Paris Cité International Laboratory in Structural Bioinformatics, Av. do Café, s/no−FCFRP, Bloco B, BR, 14040-903 Ribeirão Preto, São Paulo, Brazil
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27695, United States
| |
Collapse
|
16
|
Mendonça DC, Macedo JN, Guimarães SL, Barroso da Silva FL, Cassago A, Garratt RC, Portugal RV, Araujo APU. A revised order of subunits in mammalian septin complexes. Cytoskeleton (Hoboken) 2019; 76:457-466. [PMID: 31608568 DOI: 10.1002/cm.21569] [Citation(s) in RCA: 54] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2019] [Revised: 09/18/2019] [Accepted: 09/19/2019] [Indexed: 12/15/2022]
Abstract
Septins are GTP binding proteins considered to be novel components of the cytoskeleton. They polymerize into filaments based on hexameric or octameric core particles in which two copies of either three or four different septins, respectively, assemble into a specific sequence. Viable combinations of the 13 human septins are believed to obey substitution rules in which the different septins involved must come from distinct subgroups. The hexameric assembly, for example, has been reported to be SEPT7-SEPT6-SEPT2-SEPT2-SEPT6-SEPT7. Here, we have replaced SEPT2 by SEPT5 according to the substitution rules and used transmission electron microscopy to demonstrate that the resulting recombinant complex assembles into hexameric particles which are inverted with respect that predicted previously. MBP-SEPT5 constructs and immunostaining show that SEPT5 occupies the terminal positions of the hexamer. We further show that this is also true for the assembly including SEPT2, in direct contradiction with that reported previously. Consequently, both complexes expose an NC interface, as reported for yeast, which we show to be more susceptible to high salt concentrations. The correct assembly for the canonical combination of septins 2-6-7 is therefore established to be SEPT2-SEPT6-SEPT7-SEPT7-SEPT6-SEPT2, implying the need for revision of the mechanisms involved in filament assembly.
Collapse
Affiliation(s)
| | - Joci N Macedo
- São Carlos Institute of Physics, USP, São Carlos, SP, Brazil
- Federal Institute of Education, Science and Technology of Rondonia
| | | | - Fernando L Barroso da Silva
- Faculty of Pharmaceutical Sciences, USP, Ribeirão Preto, SP, Brazil
- UMR_S 1134, Université Paris Diderot, Paris, France
| | - Alexandre Cassago
- Brazilian Nanotechnology National Laboratory, CNPEM, Campinas, SP, Brazil
| | | | - Rodrigo V Portugal
- Brazilian Nanotechnology National Laboratory, CNPEM, Campinas, SP, Brazil
| | - Ana P U Araujo
- São Carlos Institute of Physics, USP, São Carlos, SP, Brazil
| |
Collapse
|
17
|
Pasquali S, Frezza E, Barroso da Silva FL. Coarse-grained dynamic RNA titration simulations. Interface Focus 2019; 9:20180066. [PMID: 31065339 DOI: 10.1098/rsfs.2018.0066] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/13/2019] [Indexed: 12/11/2022] Open
Abstract
Electrostatic interactions play a pivotal role in many biomolecular processes. The molecular organization and function in biological systems are largely determined by these interactions. Owing to the highly negative charge of RNA, the effect is expected to be more pronounced in this system. Moreover, RNA base pairing is dependent on the charge of the base, giving rise to alternative secondary and tertiary structures. The equilibrium between uncharged and charged bases is regulated by the solution pH, which is therefore a key environmental condition influencing the molecule's structure and behaviour. By means of constant-pH Monte Carlo simulations based on a fast proton titration scheme, coupled with the coarse-grained model HiRE-RNA, molecular dynamic simulations of RNA molecules at constant pH enable us to explore the RNA conformational plasticity at different pH values as well as to compute electrostatic properties as local pK a values for each nucleotide.
Collapse
Affiliation(s)
- S Pasquali
- Laboratoire de Cristallographie et RMN Biologiques, CNRS UMR 8015, Université Paris Descartes, Paris 75006, France
| | - E Frezza
- Laboratoire de Cristallographie et RMN Biologiques, CNRS UMR 8015, Université Paris Descartes, Paris 75006, France
| | - F L Barroso da Silva
- Departamento de Física e Química, Faculdade de Ciência s Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. do café, s/no, Ribeirão Preto, SP BR-14040-903, Brazil.,Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
18
|
Barroso da Silva FL, Sterpone F, Derreumaux P. OPEP6: A New Constant-pH Molecular Dynamics Simulation Scheme with OPEP Coarse-Grained Force Field. J Chem Theory Comput 2019; 15:3875-3888. [DOI: 10.1021/acs.jctc.9b00202] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Fernando Luís Barroso da Silva
- Departamento de Física e Química, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Universidade de São Paulo, Av. do café, s/no, Ribeirão Preto, São Paulo BR-14040-903, Brazil
- Laboratoire de Biochimie Theórique, UPR 9080 CNRS, Institut de Biologie Physico Chimique, Université Paris Diderot − Paris 7 et Université Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
| | - Fabio Sterpone
- Laboratoire de Biochimie Theórique, UPR 9080 CNRS, Institut de Biologie Physico Chimique, Université Paris Diderot − Paris 7 et Université Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| |
Collapse
|
19
|
Poveda-Cuevas S, Etchebest C, Barroso da Silva FL. Insights into the ZIKV NS1 Virology from Different Strains through a Fine Analysis of Physicochemical Properties. ACS OMEGA 2018; 3:16212-16229. [PMID: 31458257 PMCID: PMC6643396 DOI: 10.1021/acsomega.8b02081] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/17/2018] [Accepted: 11/14/2018] [Indexed: 05/02/2023]
Abstract
The flavivirus genus has several organisms responsible for generating various diseases in humans. Recently, especially in tropical regions, Zika virus (ZIKV) has raised great health concerns due to the high number of cases affecting the area during the last years that has been accompanied by a rise in the cases of the Guillain-Barré syndrome and fetal and neonatal microcephaly. Diagnosis is still difficult since the clinical symptoms between ZIKV and other flaviviruses (e.g., dengue and yellow fever) are highly similar. The understanding of their common physicochemical properties that are pH-dependent and biomolecular interaction features and their differences sheds light on the relation strain-virulence and might suggest alternative strategies toward differential serological diagnostics and therapeutic intervention. Due to their immunogenicity, the primary focus of this study was on the ZIKV nonstructural proteins 1 (NS1). By means of computational studies and semiquantitative theoretical analyses, we calculated the main physicochemical properties of this protein from different strains that are directly responsible for the biomolecular interactions and, therefore, can be related to the differential infectivity of the strains. We also mapped the electrostatic differences at both the sequence and structural levels for the strains from Uganda to Brazil, which could suggest possible molecular mechanisms for the increase of the virulence of ZIKV in Brazil. Exploring the interfaces used by NS1 to self-associate in some different oligomeric states and interact with membranes and the antibody, we could map the strategy used by the ZIKV during its evolutionary process. This indicates possible molecular mechanisms that can be correlated with the different immunological responses. By comparing with the known antibody structure available for the West Nile virus, we demonstrated that this antibody would have difficulties to neutralize the NS1 from the Brazilian strain. The present study also opens up perspectives to computationally design high-specificity antibodies.
Collapse
Affiliation(s)
- Sergio
A. Poveda-Cuevas
- Programa
Interunidades em Bioinformática, Universidade de São Paulo, São Paulo 05508-090, Brazil
- Departamento
de Física e Química, Faculdade de Ciências Farmacêuticas
de Ribeirão Preto, Universidade de
São Paulo, Ribeirão Preto, São Paulo 14040-903, Brazil
- University
of São Paulo and Université Sorbonne Paris Cité
Joint International Laboratory in Structural Bioinformatics
| | - Catherine Etchebest
- Institut
National de la Transfusion Sanguine, Paris 75015, France
- Biologie
Intégrée du Globule Rouge, Equipe 2, Dynamique des Structures
et des Interactions Moléculaires, Institut National de la Santé et de la Recherche Médicale,
UMR_S 1134, Paris 75015, France
- Université
Sorbonne Paris Cité and Université Paris Diderot, 75013 Paris, France
- University
of São Paulo and Université Sorbonne Paris Cité
Joint International Laboratory in Structural Bioinformatics
| | - Fernando L. Barroso da Silva
- Programa
Interunidades em Bioinformática, Universidade de São Paulo, São Paulo 05508-090, Brazil
- Departamento
de Física e Química, Faculdade de Ciências Farmacêuticas
de Ribeirão Preto, Universidade de
São Paulo, Ribeirão Preto, São Paulo 14040-903, Brazil
- University
of São Paulo and Université Sorbonne Paris Cité
Joint International Laboratory in Structural Bioinformatics
- Department
of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, United States
- E-mail: and
| |
Collapse
|
20
|
Srivastava D, Santiso E, Gubbins K, Barroso da Silva FL. Computationally Mapping pK a Shifts Due to the Presence of a Polyelectrolyte Chain around Whey Proteins. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:11417-11428. [PMID: 28859478 DOI: 10.1021/acs.langmuir.7b02271] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
Experimental studies have shown the formation of soluble complexes in the pure repulsive Coulombic regime even when the net charges of the protein and the polyelectrolyte have the same sign ( De Kruif et al. Curr. Opin. Colloid Interface Sci. 2004 , 9 , 340 ; De Vries et al. J. Chem. Phys. 2003 , 118 , 4649 ; Grymonpre et al. Biomacromolecules 2001 , 2 , 422 ; Hattori et al. Langmuir 2000 , 16 , 9738 ). This attractive phenomenon has often been described as "complexation on the wrong side of pI". While one theory assumes the existence of "charged patches" on the protein surface from ion-dipole interactions, thus allowing a polyelectrolyte to bind to an oppositely heterogeneous charged protein region, another theoretical view considers the induced-charge interactions to be the dominant factor in these complexations. This charge regulation mechanism can be described by proton fluctuations resulting from mutual rearrangements of the distributions of the charged groups, due to perturbations of the acid-base equilibrium. Using constant-pH Monte Carlo simulations and several quantitative and visual analysis tools, we investigate the significance of each of these interactions for two whey proteins, α-lactalbumin (α-LA) and lysozyme (LYZ). Through physical chemistry parameters, free energies of interactions, and the mapping of amino acid pKa shifts and polyelectrolyte trajectories, we show the charge regulation mechanism to be the most important contributor in protein-polyelectrolyte complexation regardless of pH, dipole moment, and protein capacitance in a low salt regime.
Collapse
Affiliation(s)
- Deepti Srivastava
- Department of Chemical and Biomolecular Engineering, North Carolina State University , Raleigh, North Carolina 27695, United States
| | - Erik Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University , Raleigh, North Carolina 27695, United States
| | - Keith Gubbins
- Department of Chemical and Biomolecular Engineering, North Carolina State University , Raleigh, North Carolina 27695, United States
| | - Fernando Luís Barroso da Silva
- Department of Chemical and Biomolecular Engineering, North Carolina State University , Raleigh, North Carolina 27695, United States
- Department of Physics and Chemistry, School of Pharmaceutical Sciences at Ribeirão Preto, University of São Paulo , 14040-903 Ribeirão Preto, SP, Brazil
| |
Collapse
|
21
|
Barroso daSilva FL, Dias LG. Development of constant-pH simulation methods in implicit solvent and applications in biomolecular systems. Biophys Rev 2017; 9:699-728. [PMID: 28921104 PMCID: PMC5662048 DOI: 10.1007/s12551-017-0311-5] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2017] [Accepted: 08/01/2017] [Indexed: 12/20/2022] Open
Abstract
pH is a critical parameter for biological and technological systems directly related with electrical charges. It can give rise to peculiar electrostatic phenomena, which also makes them more challenging. Due to the quantum nature of the process, involving the forming and breaking of chemical bonds, quantum methods should ideally by employed. Nevertheless, due to the very large number of ionizable sites, different macromolecular conformations, salt conditions, and all other charged species, the CPU time cost simply becomes prohibitive for computer simulations, making this a quite complex problem. Simplified methods based on Monte Carlo sampling have been devised and will be reviewed here, highlighting the updated state-of-the-art of this field, advantages, and limitations of different theoretical protocols for biomolecular systems (proteins and nucleic acids). Following a historical perspective, the discussion will be associated with the applications to protein interactions with other proteins, polyelectrolytes, and nanoparticles.
Collapse
Affiliation(s)
- Fernando Luís Barroso daSilva
- Departamento de Física e Química, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. do café, s/no. - Universidade de São Paulo, BR-14040-903, Ribeirão Preto, SP, Brazil.
- UCD School of Physics, UCD Institute for Discovery, University College Dublin, Belfield, Dublin 4, Ireland.
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC, USA.
| | - Luis Gustavo Dias
- Departamento de Química, Faculdade de Filosofia, Ciências e Letras de Ribeirão Preto, Av. Bandeirantes, 3900 - Universidade de São Paulo, BR-14040-901, Ribeirão Preto, SP, Brazil
| |
Collapse
|
22
|
Barroso da Silva FL, Derreumaux P, Pasquali S. Protein-RNA complexation driven by the charge regulation mechanism. Biochem Biophys Res Commun 2017; 498:264-273. [PMID: 28709871 DOI: 10.1016/j.bbrc.2017.07.027] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2017] [Revised: 06/27/2017] [Accepted: 07/05/2017] [Indexed: 11/27/2022]
Abstract
Electrostatic interactions play a pivotal role in many (bio)molecular association processes. The molecular organization and function in biological systems are largely determined by these interactions from pure Coulombic contributions to more peculiar mesoscopic forces due to ion-ion correlation and proton fluctuations. The latter is a general electrostatic mechanism that gives attraction particularly at low electrolyte concentrations. This charge regulation mechanism due to titrating amino acid and nucleotides residues is discussed here in a purely electrostatic framework. By means of constant-pH Monte Carlo simulations based on a fast coarse-grained titration proton scheme, a new computer molecular model was devised to study protein-RNA interactions. The complexation between the RNA silencing suppressor p19 viral protein and the 19-bp small interfering RNA was investigated at different solution pH and salt conditions. The outcomes illustrate the importance of the charge regulation mechanism that enhances the association between these macromolecules in a similar way as observed for other protein-polyelectrolyte systems typically found in colloidal science. Due to the highly negative charge of RNA, the effect is more pronounced in this system as predicted by the Kirkwood-Shumaker theory. Our results contribute to the general physico-chemical understanding of macromolecular complexation and shed light on the extensive role of RNA in the cell's life.
Collapse
Affiliation(s)
- Fernando Luís Barroso da Silva
- Departamento de Física e Química, Faculdade de Ciências Farmacêuticas de Ribeirão Preto, Av. do café, s/no. - Universidade de São Paulo, BR-14040-903 Ribeirão Preto, SP, Brazil; Laboratoire de Biochimie Theórique, UPR 9080 CNRS, Institut de Biologie Physico Chimique, Université Paris Diderot - Paris 7 et Université Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France.
| | - Philippe Derreumaux
- Laboratoire de Biochimie Theórique, UPR 9080 CNRS, Institut de Biologie Physico Chimique, Université Paris Diderot - Paris 7 et Université Sorbonne Paris Cité, 13 rue Pierre et Marie Curie, 75005 Paris, France
| | - Samuela Pasquali
- Laboratoire de Cristallographie et RMN Biologiques, UMR 8015 CNRS, Faculté des sciences pharmaceutiques et biologiques, Universtié Paris Descartes et Université Sorbonne Paris Cité, 4 Avenue de l'Observatoire, 75006 Paris, France
| |
Collapse
|