1
|
Pirojsirikul T, Lee VS, Nimmanpipug P. Unraveling Bacterial Single-Stranded Sequence Specificities: Insights from Molecular Dynamics and MMPBSA Analysis of Oligonucleotide Probes. Mol Biotechnol 2024; 66:582-591. [PMID: 38374320 DOI: 10.1007/s12033-024-01082-0] [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: 10/03/2022] [Accepted: 01/10/2024] [Indexed: 02/21/2024]
Abstract
We utilized molecular dynamics (MD) simulations and Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) free energy calculations to investigate the specificity of two oligonucleotide probes, namely probe B and probe D, in detecting single-stranded DNA (ssDNA) within three bacteria families: Enterobacteriaceae, Pasteurellaceae, and Vibrionaceae. Due to the limited understanding of molecular mechanisms in the previous research, we have extended the discussion to focus specifically on investigating the binding process of bacteria-probe DNA duplexes, with an emphasis on analyzing the binding free energy. The role of electrostatic contributions in the specificity between the oligonucleotide probes and the bacterial ssDNAs was investigated and found to be crucial. Our calculations yielded results that were highly consistent with the experimental data. Through our study, we have successfully exhibited the benefits of utilizing in-silico approaches as a powerful virtual-screening tool, particularly in research areas that demand a thorough comprehension of molecular interactions.
Collapse
Affiliation(s)
- Teerapong Pirojsirikul
- Division of Physical Science, Faculty of Science, Prince of Songkla University, Songkhla, 90110, Thailand.
| | - Vannajan Sanghiran Lee
- Department of Chemistry, Center of Theoretical and Computational Physics, Faculty of Science, University of Malaya, 50603, Kuala Lumpur, Malaysia
| | - Piyarat Nimmanpipug
- Department of Chemistry, Faculty of Science, Chiang Mai University, Chiang Mai, 50200, Thailand
| |
Collapse
|
2
|
Jonniya NA, Poddar S, Mahapatra S, Kar P. Computer-aided Affinity Enhancement of a Cross-reactive Antibody against Dengue Virus Envelope Domain III. Cell Biochem Biophys 2023; 81:737-755. [PMID: 37735329 DOI: 10.1007/s12013-023-01175-8] [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] [Accepted: 09/11/2023] [Indexed: 09/23/2023]
Abstract
The dengue virus (DENV), composed of four distinct but serologically related Flaviviruses, causes the most important emerging viral disease, with nearly 400 million infections yearly. Currently, there are no approved therapies. Although DENV infection induces lifelong immunity against the same serotype, the antibodies raised contribute to severe disease in heterotypic infections. Therefore, understanding the mechanism of DENV neutralization by antibodies is crucial in the design of vaccines against all serotypes. This study reports a comparative structural and energetic analysis of the monoclonal antibody (mAb) 4E11 in complex with its target domain III of the envelope protein for all four DENV serotypes. We use extensive replica molecular dynamics simulations in conjunction with the binding free energy calculations. Further single point and double mutations were designed through computational site-directed mutagenesis and observed that the re-engineered antibody exhibits high affinity to binding and broadly neutralizing activity against serotypes. Our results showed improved binding affinity by the gain of enthalpy, which could be attributed to the stabilization of salt-bridge and hydrogen bond interactions at the antigen-antibody interface. The findings provide valuable results in understanding the structural dynamics and energetic contributions that will be helpful to the design of high-affinity antibodies against dengue infections.
Collapse
Affiliation(s)
- Nisha Amarnath Jonniya
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore, 453552, Madhya Pradesh, India
- Department of Pharmacology, University of California San Diego, 9500 Gilman Drive, La Jolla, CA, 92093, USA
| | - Sayan Poddar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore, 453552, Madhya Pradesh, India
| | - Subhasmita Mahapatra
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore, 453552, Madhya Pradesh, India
| | - Parimal Kar
- Department of Biosciences and Biomedical Engineering, Indian Institute of Technology Indore, Khandwa Road, Indore, 453552, Madhya Pradesh, India.
| |
Collapse
|
3
|
Zhang X, Wu R, Liu Q, Zhou L. Interfacial water molecules contribute to antibody binding to the receptor-binding domain of SARS-CoV-2 spike protein. J Biomol Struct Dyn 2023; 41:14929-14938. [PMID: 37042961 DOI: 10.1080/07391102.2023.2199336] [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: 11/23/2022] [Accepted: 02/18/2023] [Indexed: 04/13/2023]
Abstract
Antibodies that recognize the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), especially the neutralizing antibodies, carry great hope in the treatment and final elimination of COVID-19. Driven by a synchronized global effort, thousands of antibodies against the spike protein have been identified during the past two years, with the structural information available at atomistic detail for hundreds of these antibodies. We developed an improved molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method including explicitly treated interfacial water to calculate the binding free energy between representative antibodies and the receptor binding domain (RBD) domain of SARS-COV-2 spike proteins. We discovered that explicit treatment of water molecules located at the interface between RBD and antibody effectively improves the results for the WT and variants of concern (VOC) systems. Interfacial water molecules, together with surface and internal water molecules, behave drastically from bulk water and exert peculiar impacts on protein dynamics and energy, and thus warrant explicit treatment to complement implicit solvent models. Our results illustrate the importance of including interfacial water molecules to approach efficient and reliable prediction of binding free energy.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Xin Zhang
- Institute of Molecular Physiology, Shenzhen Bay Laboratory, Shenzhen, Guangdong, P.R. China
| | - Ruiping Wu
- Institute of Molecular Physiology, Shenzhen Bay Laboratory, Shenzhen, Guangdong, P.R. China
- Department of Chemistry, Southern University of Science and Technology, Shenzhen, Guangdong, P. R. China
| | - Qinglian Liu
- Department of Physiology and Biophysics, School of Medicine, Virginia Commonwealth University, Richmond, Virginia, USA
| | - Lei Zhou
- Institute of Molecular Physiology, Shenzhen Bay Laboratory, Shenzhen, Guangdong, P.R. China
| |
Collapse
|
4
|
Fernández-Quintero ML, Ljungars A, Waibl F, Greiff V, Andersen JT, Gjølberg TT, Jenkins TP, Voldborg BG, Grav LM, Kumar S, Georges G, Kettenberger H, Liedl KR, Tessier PM, McCafferty J, Laustsen AH. Assessing developability early in the discovery process for novel biologics. MAbs 2023; 15:2171248. [PMID: 36823021 PMCID: PMC9980699 DOI: 10.1080/19420862.2023.2171248] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 01/18/2023] [Indexed: 02/25/2023] Open
Abstract
Beyond potency, a good developability profile is a key attribute of a biological drug. Selecting and screening for such attributes early in the drug development process can save resources and avoid costly late-stage failures. Here, we review some of the most important developability properties that can be assessed early on for biologics. These include the influence of the source of the biologic, its biophysical and pharmacokinetic properties, and how well it can be expressed recombinantly. We furthermore present in silico, in vitro, and in vivo methods and techniques that can be exploited at different stages of the discovery process to identify molecules with liabilities and thereby facilitate the selection of the most optimal drug leads. Finally, we reflect on the most relevant developability parameters for injectable versus orally delivered biologics and provide an outlook toward what general trends are expected to rise in the development of biologics.
Collapse
Affiliation(s)
- Monica L. Fernández-Quintero
- Center for Molecular Biosciences Innsbruck (CMBI), Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Anne Ljungars
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Franz Waibl
- Center for Molecular Biosciences Innsbruck (CMBI), Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Victor Greiff
- Department of Immunology, University of Oslo, Oslo, Norway
| | - Jan Terje Andersen
- Department of Immunology, University of Oslo, Oslo University Hospital Rikshospitalet, Oslo, Norway
- Institute of Clinical Medicine and Department of Pharmacology, University of Oslo, Oslo, Norway
| | | | - Timothy P. Jenkins
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Bjørn Gunnar Voldborg
- National Biologics Facility, Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Lise Marie Grav
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| | - Sandeep Kumar
- Biotherapeutics Discovery, Boehringer Ingelheim Pharmaceuticals Inc, Ridgefield, CT, USA
| | - Guy Georges
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Hubert Kettenberger
- Roche Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, Germany
| | - Klaus R. Liedl
- Center for Molecular Biosciences Innsbruck (CMBI), Department of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innsbruck, Austria
| | - Peter M. Tessier
- Department of Chemical Engineering, Pharmaceutical Sciences and Biomedical Engineering, Biointerfaces Institute, University of Michigan, Ann Arbor, Michigan, USA
| | - John McCafferty
- Department of Medicine, Cambridge Institute of Therapeutic Immunology and Infectious Disease, University of Cambridge, Cambridge, UK
- Maxion Therapeutics, Babraham Research Campus, Cambridge, UK
| | - Andreas H. Laustsen
- Department of Biotechnology and Biomedicine, Technical University of Denmark, Kongens Lyngby, Denmark
| |
Collapse
|
5
|
Ahmad B, Choi S. Unraveling the Tomaralimab Epitope on the Toll-like Receptor 2 via Molecular Dynamics and Deep Learning. ACS OMEGA 2022; 7:28226-28237. [PMID: 35990491 PMCID: PMC9386714 DOI: 10.1021/acsomega.2c02559] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Accepted: 07/22/2022] [Indexed: 06/15/2023]
Abstract
Tomaralimab (OPN-305) is the first humanized immunoglobulin G4 monoclonal antibody against TLR2 and is designed to prevent inflammation that is driven by inappropriate or excessive activation of innate immune pathways. Here, we constructed a homology model of Tomaralimab and its complex with TLR2 at different mapped epitopes and unraveled their behavior at the atomistic level. Furthermore, we predicted a novel epitope (leucine-rich region 9-12) near the lipopeptide-binding site that can be targeted and studied for the utility of therapeutic antibodies. A geometric deep learning algorithm was used to envisage Tomaralimab binding affinity changes upon mutation. There was a significant difference in binding affinity for Tomaralimab following epitope-mutated alanine substitutions of Val266, Pro294, Arg295, Asn319, Pro326, and His372. Using deep learning-based ΔΔG prediction, we computationally contrasted human TLR2-TLR2, TLR2-TLR1, and TLR2-TLR6 dimerization. These results reveal the mechanism that underlies Tomaralimab binding to TLR2 and should help to design structure-based mimics or bispecific antibodies that can be used to inhibit both lipopeptide-binding and TLR2 dimerization.
Collapse
Affiliation(s)
- Bilal Ahmad
- Department
of Molecular Science and Technology, Ajou
University, Suwon 16499, Korea
- S&K
Therapeutics, Ajou University
Campus Plaza 418, 199 Worldcup-ro, Yeongtong-gu, Suwon 16502, Korea
| | - Sangdun Choi
- Department
of Molecular Science and Technology, Ajou
University, Suwon 16499, Korea
- S&K
Therapeutics, Ajou University
Campus Plaza 418, 199 Worldcup-ro, Yeongtong-gu, Suwon 16502, Korea
| |
Collapse
|
6
|
Wong MTY, Kelm S, Liu X, Taylor RD, Baker T, Essex JW. Higher Affinity Antibodies Bind With Lower Hydration and Flexibility in Large Scale Simulations. Front Immunol 2022; 13:884110. [PMID: 35707541 PMCID: PMC9190259 DOI: 10.3389/fimmu.2022.884110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Accepted: 04/22/2022] [Indexed: 11/25/2022] Open
Abstract
We have carried out a long-timescale simulation study on crystal structures of nine antibody-antigen pairs, in antigen-bound and antibody-only forms, using molecular dynamics with enhanced sampling and an explicit water model to explore interface conformation and hydration. By combining atomic level simulation and replica exchange to enable full protein flexibility, we find significant numbers of bridging water molecules at the antibody-antigen interface. Additionally, a higher proportion of interactions excluding bulk waters and a lower degree of antigen bound CDR conformational sampling are correlated with higher antibody affinity. The CDR sampling supports enthalpically driven antibody binding, as opposed to entropically driven, in that the difference between antigen bound and unbound conformations do not correlate with affinity. We thus propose that interactions with waters and CDR sampling are aspects of the interface that may moderate antibody-antigen binding, and that explicit hydration and CDR flexibility should be considered to improve antibody affinity prediction and computational design workflows.
Collapse
Affiliation(s)
- Mabel T. Y. Wong
- School of Chemistry, University of Southampton, Southampton, United Kingdom
| | | | | | | | | | - Jonathan W. Essex
- School of Chemistry, University of Southampton, Southampton, United Kingdom
- *Correspondence: Jonathan W. Essex,
| |
Collapse
|
7
|
Greene D, Barton M, Luchko T, Shiferaw Y. Computational Analysis of Binding Interactions between the Ryanodine Receptor Type 2 and Calmodulin. J Phys Chem B 2021; 125:10720-10735. [PMID: 34533024 DOI: 10.1021/acs.jpcb.1c03896] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to a variety of cardiac arrhythmias, such as catecholaminergic polymorphic ventricular tachycardia (CPVT). RyR2 is regulated by calmodulin (CaM), and mutations that disrupt their interaction can cause aberrant calcium release, leading to an arrhythmia. It was recently shown that increasing the RyR2-CaM binding affinity could rescue a defective CPVT-related RyR2 channel to near wild-type behavior. However, the interactions that determine the binding affinity at the RyR2-CaM binding interface are not well understood. In this study, we identify the key domains and interactions, including several new interactions, involved in the binding of CaM to RyR2. Also, our comparison between the wild-type and V3599K mutant suggests how the RyR2-CaM binding affinity can be increased via a change in the central and N-terminal lobe binding contacts for CaM. This computational approach provides new insights into the effect of a mutation at the RyR2-CaM binding interface, and it may find utility in drug design for the future treatment of cardiac arrhythmias.
Collapse
Affiliation(s)
- D'Artagnan Greene
- Department of Physics, California State University, Northridge, California 91330, United States
| | - Michael Barton
- Department of Physics, California State University, Northridge, California 91330, United States
| | - Tyler Luchko
- Department of Physics, California State University, Northridge, California 91330, United States
| | - Yohannes Shiferaw
- Department of Physics, California State University, Northridge, California 91330, United States
| |
Collapse
|
8
|
Protein-Protein Interactions: Insight from Molecular Dynamics Simulations and Nanoparticle Tracking Analysis. Molecules 2021; 26:molecules26185696. [PMID: 34577167 PMCID: PMC8472368 DOI: 10.3390/molecules26185696] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2021] [Revised: 09/13/2021] [Accepted: 09/16/2021] [Indexed: 11/23/2022] Open
Abstract
Protein-protein interaction plays an essential role in almost all cellular processes and biological functions. Coupling molecular dynamics (MD) simulations and nanoparticle tracking analysis (NTA) assay offered a simple, rapid, and direct approach in monitoring the protein-protein binding process and predicting the binding affinity. Our case study of designed ankyrin repeats proteins (DARPins)—AnkGAG1D4 and the single point mutated AnkGAG1D4-Y56A for HIV-1 capsid protein (CA) were investigated. As reported, AnkGAG1D4 bound with CA for inhibitory activity; however, it lost its inhibitory strength when tyrosine at residue 56 AnkGAG1D4, the most key residue was replaced by alanine (AnkGAG1D4-Y56A). Through NTA, the binding of DARPins and CA was measured by monitoring the increment of the hydrodynamic radius of the AnkGAG1D4-gold conjugated nanoparticles (AnkGAG1D4-GNP) and AnkGAG1D4-Y56A-GNP upon interaction with CA in buffer solution. The size of the AnkGAG1D4-GNP increased when it interacted with CA but not AnkGAG1D4-Y56A-GNP. In addition, a much higher binding free energy (∆GB) of AnkGAG1D4-Y56A (−31 kcal/mol) obtained from MD further suggested affinity for CA completely reduced compared to AnkGAG1D4 (−60 kcal/mol). The possible mechanism of the protein-protein binding was explored in detail by decomposing the binding free energy for crucial residues identification and hydrogen bond analysis.
Collapse
|
9
|
Roy U. Structural and molecular analyses of functional epitopes and escape mutants in Japanese encephalitis virus envelope protein domain III. Immunol Res 2021; 68:81-89. [PMID: 32445181 PMCID: PMC7243247 DOI: 10.1007/s12026-020-09130-y] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
The Japanese encephalitis virus (JEV) is one of the vector borne causes of encephalitis found in southeastern Asia. This positive single-stranded RNA virus is a member of the Flaviviridae family, which notably includes dengue, tick-borne, West Nile, Zika as well as yellow fever, and transmits to humans by infected mosquitos. The main site of interactions for antibodies against this virus is the envelope protein domain III (ED3). The present report investigates the time-dependent structural and conformational changes of JEV ED3 functional epitopes and escape mutants by computer simulations. The results indicate the presence of significant structural differences between the functional epitopes and the escape mutants. Mutation-induced structural/conformational instabilities of this type can decrease the antibody neutralization activity. Among the different escape mutants studied here, Ser40Lys/Asp41Arg appear to be most unstable, while Ser40Glu/Asp41Leu exhibit the lowest structural variations. The highest level of escape mutation observed in Ser40Lys is linked to the relatively higher values of root mean square deviation/fluctuation found in the molecular dynamics simulation of this protein. Secondary-structure deviations and depletion of H bonding are other contributing factors to the protein’s increased instability. Overall, the proteins with residue 41 mutations are found to be structurally more ordered than those with residue 40 mutations. The detailed time-based structural assessment of the mutant epitopes described here may contribute to the development of novel vaccines and antiviral drugs necessary to defend against future outbreaks of JEV escape mutants.
Collapse
Affiliation(s)
- Urmi Roy
- Department of Chemistry & Biomolecular Science, Clarkson University, 8 Clarkson Avenue, Potsdam, NY, 13699-5820, USA.
| |
Collapse
|
10
|
Yamashita T. Toward rational antibody design: recent advancements in molecular dynamics simulations. Int Immunol 2019; 30:133-140. [PMID: 29346652 DOI: 10.1093/intimm/dxx077] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2017] [Accepted: 01/11/2018] [Indexed: 01/02/2023] Open
Abstract
Because antibodies have become an important therapeutic tool, rational antibody design is a challenging issue involving various science and technology fields. From the computational aspect, many types of design-assist methods have been developed, but their accuracy is not fully satisfactory. Because of recent advancements in computational power, molecular dynamics (MD) simulation has become a helpful tool to trace the motion of proteins and to characterize their properties. Thus, MD simulation has been applied to various systems involving antigen-antibody complexes and has been shown to provide accurate insight into antigen-antibody interactions and dynamics at an atomic resolution. Therefore, it is highly possible that MD simulation will play several roles complementing the conventional antibody design. In this review, we address several important features of MD simulation in the context of rational antibody design.
Collapse
Affiliation(s)
- Takefumi Yamashita
- Laboratory for Systems Biology and Medicine, Research Center for Advanced Science and Technology, The University of Tokyo, Japan
| |
Collapse
|
11
|
Soong JX, Chan SK, Lim TS, Choong YS. Optimisation of human V H domain antibodies specific to Mycobacterium tuberculosis heat shock protein (HSP16.3). J Comput Aided Mol Des 2019; 33:375-385. [PMID: 30689080 DOI: 10.1007/s10822-019-00186-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Accepted: 01/22/2019] [Indexed: 11/29/2022]
Abstract
Mycobacterium tuberculosis (Mtb) 16.3 kDa heat shock protein 16.3 (HSP16.3) is a latency-associated antigen that can be targeted for latent tuberculosis (TB) diagnostic and therapeutic development. We have previously developed human VH domain antibodies (dAbs; clone E3 and F1) specific against HSP16.3. In this work, we applied computational methods to optimise and design the antibodies in order to improve the binding affinity with HSP16.3. The VH domain antibodies were first docked to the dimer form of HSP16.3 and further sampled using molecular dynamics simulation. The calculated binding free energy of the HSP16.3-dAb complexes showed non-polar interactions were responsible for the antigen-antibody association. Per-residue free energy decomposition and computational alanine scanning have identified one hotspot residue for E3 (Y391) and 4 hotspot residues for F1 (M394, Y396, R397 and M398). These hotspot residues were then mutated and evaluated by binding free energy calculations. Phage ELISA assay was carried out on the potential mutants (E3Y391W, F1M394E, F1R397N and F1M398Y). The experimental assay showed improved binding affinities of E3Y391W and F1M394E against HSP16.3 compared with the wild type E3 and F1. This case study has thus showed in silico methods are able to assist in optimisation or improvement of antibody-antigen binding.
Collapse
Affiliation(s)
- Jia Xin Soong
- Institute for Research in Molecular Medicine (INFORMM), Universiti Sains Malaysia, Minden, 11800, Penang, Malaysia
| | - Soo Khim Chan
- Institute for Research in Molecular Medicine (INFORMM), Universiti Sains Malaysia, Minden, 11800, Penang, Malaysia
| | - Theam Soon Lim
- Institute for Research in Molecular Medicine (INFORMM), Universiti Sains Malaysia, Minden, 11800, Penang, Malaysia.,Analytical Biochemistry Research Centre, Universiti Sains Malaysia, Minden, 11800, Penang, Malaysia
| | - Yee Siew Choong
- Institute for Research in Molecular Medicine (INFORMM), Universiti Sains Malaysia, Minden, 11800, Penang, Malaysia.
| |
Collapse
|
12
|
Kaur H, Sain N, Mohanty D, Salunke DM. Deciphering evolution of immune recognition in antibodies. BMC STRUCTURAL BIOLOGY 2018; 18:19. [PMID: 30563492 PMCID: PMC6299584 DOI: 10.1186/s12900-018-0096-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/17/2018] [Accepted: 11/14/2018] [Indexed: 11/29/2022]
Abstract
Background Antibody, the primary effector molecule of the immune system, evolves after initial encounter with the antigen from a precursor form to a mature one to effectively deal with the antigen. Antibodies of a lineage diverge through antigen-directed isolated pathways of maturation to exhibit distinct recognition potential. In the context of evolution in immune recognition, diversity of antigen cannot be ignored. While there are reports on antibody lineage, structural perspective with respect to diverse recognition potential in a lineage has never been studied. Hence, it is crucial to evaluate how maturation leads to topological tailoring within a lineage enabling them to interact with significantly distinct antigens. Results A data-driven approach was undertaken for the study. Global experimental mouse and human antibody-antigen complex structures from PDB were compiled into a coherent database of germline-linked antibodies bound with distinct antigens. Structural analysis of all lineages showed variations in CDRs of both H and L chains. Observations of conformational adaptation made from analysis of static structures were further evaluated by characterizing dynamics of interaction in two lineages, mouse VH1–84 and human VH5–51. Sequence and structure analysis of the lineages explained that somatic mutations altered the geometries of individual antibodies with common structural constraints in some CDRs. Additionally, conformational landscape obtained from molecular dynamics simulations revealed that incoming pathogen led to further conformational divergence in the paratope (as observed across datasets) even while maintaining similar overall backbone topology. MM-GB/SA analysis showed binding energies to be in physiological range. Results of the study are coherent with experimental observations. Conclusions The findings of this study highlight basic structural principles shaping the molecular evolution of a lineage for significantly diverse antigens. Antibodies of a lineage follow different developmental pathways while preserving the imprint of the germline. From the study, it can be generalized that structural diversification of the paratope is an outcome of natural selection of a conformation from an available ensemble, which is further optimized for antigen interaction. The study establishes that starting from a common lineage, antibodies can mature to recognize a wide range of antigens. This hypothesis can be further tested and validated experimentally. Electronic supplementary material The online version of this article (10.1186/s12900-018-0096-1) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Harmeet Kaur
- Regional Centre for Biotechnology, Biotech Science Cluster, Faridabad, Haryana, 121001, India.,Manipal Academy of Higher Education, Manipal, Karnataka, 576104, India
| | - Neetu Sain
- National Institute of Immunology, New Delhi, Delhi, 110067, India
| | - Debasisa Mohanty
- National Institute of Immunology, New Delhi, Delhi, 110067, India
| | - Dinakar M Salunke
- Regional Centre for Biotechnology, Biotech Science Cluster, Faridabad, Haryana, 121001, India. .,International Centre for Genetic Engineering and Biotechnology, New Delhi, Delhi, 110067, India.
| |
Collapse
|
13
|
Hadden JA, Perilla JR. All-atom virus simulations. Curr Opin Virol 2018; 31:82-91. [PMID: 30181049 PMCID: PMC6456034 DOI: 10.1016/j.coviro.2018.08.007] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2018] [Revised: 08/04/2018] [Accepted: 08/13/2018] [Indexed: 12/11/2022]
Abstract
The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.
Collapse
Affiliation(s)
- Jodi A Hadden
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States.
| | - Juan R Perilla
- Department of Chemistry and Biochemistry, University of Delaware, Newark, DE 19716, United States
| |
Collapse
|
14
|
Wang C, Greene D, Xiao L, Qi R, Luo R. Recent Developments and Applications of the MMPBSA Method. Front Mol Biosci 2018; 4:87. [PMID: 29367919 PMCID: PMC5768160 DOI: 10.3389/fmolb.2017.00087] [Citation(s) in RCA: 332] [Impact Index Per Article: 55.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/18/2017] [Accepted: 11/30/2017] [Indexed: 12/23/2022] Open
Abstract
The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.
Collapse
Affiliation(s)
- Changhao Wang
- Chemical and Materials Physics Graduate Program, University of California, Irvine, Irvine, CA, United States
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Physics and Astronomy, University of California, Irvine, Irvine, CA, United States
| | - D'Artagnan Greene
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
| | - Li Xiao
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, United States
| | - Ruxi Qi
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
| | - Ray Luo
- Chemical and Materials Physics Graduate Program, University of California, Irvine, Irvine, CA, United States
- Department of Molecular Biology and Biochemistry, University of California, Irvine, Irvine, CA, United States
- Department of Biomedical Engineering, University of California, Irvine, Irvine, CA, United States
- Department of Chemical Engineering and Materials Science, University of California, Irvine, Irvine, CA, United States
| |
Collapse
|
15
|
Abstract
With the drive toward high throughput molecular dynamics (MD) simulations involving ever-greater numbers of simulation replicates run for longer, biologically relevant timescales (microseconds), the need for improved computational methods that facilitate fully automated MD workflows gains more importance. Here we report the development of an automated workflow tool to perform AMBER GPU MD simulations. Our workflow tool capitalizes on the capabilities of the Kepler platform to deliver a flexible, intuitive, and user-friendly environment and the AMBER GPU code for a robust and high-performance simulation engine. Additionally, the workflow tool reduces user input time by automating repetitive processes and facilitates access to GPU clusters, whose high-performance processing power makes simulations of large numerical scale possible. The presented workflow tool facilitates the management and deployment of large sets of MD simulations on heterogeneous computing resources. The workflow tool also performs systematic analysis on the simulation outputs and enhances simulation reproducibility, execution scalability, and MD method development including benchmarking and validation.
Collapse
|
16
|
Lees WD, Stejskal L, Moss DS, Shepherd AJ. Investigating Substitutions in Antibody-Antigen Complexes Using Molecular Dynamics: A Case Study with Broad-spectrum, Influenza A Antibodies. Front Immunol 2017; 8:143. [PMID: 28261207 PMCID: PMC5309259 DOI: 10.3389/fimmu.2017.00143] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2016] [Accepted: 01/30/2017] [Indexed: 11/20/2022] Open
Abstract
In studying the binding of host antibodies to the surface antigens of pathogens, the structural and functional characterization of antibody–antigen complexes by X-ray crystallography and binding assay is important. However, the characterization requires experiments that are typically time consuming and expensive: thus, many antibody–antigen complexes are under-characterized. For vaccine development and disease surveillance, it is often vital to assess the impact of amino acid substitutions on antibody binding. For example, are there antibody substitutions capable of improving binding without a loss of breadth, or antigen substitutions that lead to antigenic escape? The questions cannot be answered reliably from sequence variation alone, exhaustive substitution assays are usually impractical, and alanine scans provide at best an incomplete identification of the critical residue–residue interactions. Here, we show that, given an initial structure of an antibody bound to an antigen, molecular dynamics simulations using the energy method molecular mechanics with Generalized Born surface area (MM/GBSA) can model the impact of single amino acid substitutions on antibody–antigen binding energy. We apply the technique to three broad-spectrum antibodies to influenza A hemagglutinin and examine both previously characterized and novel variant strains observed in the human population that may give rise to antigenic escape. We find that in some cases the impact of a substitution is local, while in others it causes a reorientation of the antibody with wide-ranging impact on residue–residue interactions: this explains, in part, why the change in chemical properties of a residue can be, on its own, a poor predictor of overall change in binding energy. Our estimates are in good agreement with experimental results—indeed, they approximate the degree of agreement between different experimental techniques. Simulations were performed on commodity computer hardware; hence, this approach has the potential to be widely adopted by those undertaking infectious disease research. Novel aspects of this research include the application of MM/GBSA to investigate binding between broadly binding antibodies and a viral glycoprotein; the development of an approach for visualizing substrate–ligand interactions; and the use of experimental assay data to rescale our predictions, allowing us to make inferences about absolute, as well as relative, changes in binding energy.
Collapse
Affiliation(s)
- William D Lees
- Institute of Structural and Molecular Biology, Birkbeck College , London , UK
| | - Lenka Stejskal
- Institute of Structural and Molecular Biology, Birkbeck College , London , UK
| | - David S Moss
- Institute of Structural and Molecular Biology, Birkbeck College , London , UK
| | - Adrian J Shepherd
- Institute of Structural and Molecular Biology, Birkbeck College , London , UK
| |
Collapse
|