scholarly journals How Communication Pathways Bridge Local and Global Conformations in an IgG4 Antibody: a Molecular Dynamics Study

2021 ◽  
Author(s):  
Thomas Tarenzi ◽  
Marta Rigoli ◽  
Raffaello Potestio
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Md Simulations ◽  
Specific Sequence ◽  
Multi Scale ◽  
Antibody Structure ◽  
Complementarity Determining Regions ◽  
Global And Local ◽  
Cell Death Protein ◽  
Antibody Function

The affinity of an antibody for its antigen is primarily determined by the specific sequence and structural arrangement of the complementarity-determining regions (CDRs). Recently, however, evidence has accumulated that points toward a nontrivial relation between the CDR and distal sites on the antibody structure: variations in the binding strengths have been observed upon mutating amino acids separated from the paratope by several nanometers, thus suggesting the existence of a communication network within antibodies whose extension and relevance might be deeper than insofar expected. In this work, we test this hypothesis by means of molecular dynamics (MD) simulations of the IgG4 monoclonal antibody pembrolizumab, an approved drug that targets the programmed cell death protein 1 (PD-1). The molecule is simulated in both the apo and holo states, totalling 4 μs of MD trajectory. The analysis of these simulations shows that the bound antibody explores a restricted range of conformations with respect to the apo one, and that the global conformation of the molecule correlates with that of the CDR; a pivotal role in this relationship is played by the relatively short hinge, which mechanically couples Fab and Fc domains. These results support the hypothesis that pembrolizumab behaves as a complex machinery, with a multi-scale hierarchy of global and local conformational changes that communicate with one another. The analysis pipeline developed in this work is general, and it can help shed further light on the mechanistic aspects of antibody function.

scholarly journals Communication pathways bridge local and global conformations in an IgG4 antibody

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Thomas Tarenzi ◽  
Marta Rigoli ◽  
Raffaello Potestio
Keyword(s):  
Conformational Changes ◽  
Md Simulations ◽  
Specific Sequence ◽  
Multi Scale ◽  
Complementarity Determining Regions ◽  
Global And Local ◽  
Cell Death Protein ◽  
Antibody Function ◽  
Scale Hierarchy ◽  
Igg4 Antibody

AbstractThe affinity of an antibody for its antigen is primarily determined by the specific sequence and structural arrangement of the complementarity-determining regions (CDRs). Recent evidence, however, points toward a nontrivial relation between the CDR and distal sites: variations in the binding strengths have been observed upon mutating residues separated from the paratope by several nanometers, thus suggesting the existence of a communication network within antibodies, whose extension and relevance might be deeper than insofar expected. In this work, we test this hypothesis by means of molecular dynamics (MD) simulations of the IgG4 monoclonal antibody pembrolizumab, an approved drug that targets the programmed cell death protein 1 (PD-1). The molecule is simulated in both the apo and holo states, totalling 4 μs of MD trajectory. The analysis of these simulations shows that the bound antibody explores a restricted range of conformations with respect to the apo one, and that the global conformation of the molecule correlates with that of the CDR. These results support the hypothesis that pembrolizumab featues a multi-scale hierarchy of intertwined global and local conformational changes. The analysis pipeline developed in this work is general, and it can help shed further light on the mechanistic aspects of antibody function.

scholarly journals A Virtual Reality Ensemble Molecular Dynamics Workflow to Study Complex Conformational Changes in Proteins

2020 ◽  
Author(s):  
Jordi Juárez-Jiménez ◽  
Philip Tew ◽  
Michael o'connor ◽  
Salome Llabres ◽  
Rebecca Sage ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Virtual Reality ◽  
Conformational Changes ◽  
Large Scale ◽  
A Priori ◽  
Computational Cost ◽  
Md Simulations ◽  
Collective Variables ◽  
Sampling Protocols ◽  
State Models

<p>Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics and their biological function. Currently it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome ‘collective variable’ enhanced sampling protocols. Here we describe a framework that combines ensemble MD simulations and virtual-reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables <i>a priori</i>. We further show that eMD-VR generated pathways can be combined with Markov State Models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts. </p>

scholarly journals Combining Virtual Reality Visualization with Ensemble Molecular Dynamics to Study Complex Protein Conformational Changes

2020 ◽  
Author(s):  
Jordi Juárez-Jiménez ◽  
Philip Tew ◽  
Michael o'connor ◽  
Salome Llabres ◽  
Rebecca Sage ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Virtual Reality ◽  
Conformational Changes ◽  
Large Scale ◽  
Computational Cost ◽  
Md Simulations ◽  
Collective Variables ◽  
Protein Conformational Changes ◽  
Sampling Protocols ◽  
State Models

<p>Molecular dynamics (MD) simulations are increasingly used to elucidate relationships between protein structure, dynamics and their biological function. Currently it is extremely challenging to perform MD simulations of large-scale structural rearrangements in proteins that occur on millisecond timescales or beyond, as this requires very significant computational resources, or the use of cumbersome ‘collective variable’ enhanced sampling protocols. Here we describe a framework that combines ensemble MD simulations and virtual-reality visualization (eMD-VR) to enable users to interactively generate realistic descriptions of large amplitude, millisecond timescale protein conformational changes in proteins. Detailed tests demonstrate that eMD-VR substantially decreases the computational cost of folding simulations of a WW domain, without the need to define collective variables <i>a priori</i>. We further show that eMD-VR generated pathways can be combined with Markov State Models to describe the thermodynamics and kinetics of large-scale loop motions in the enzyme cyclophilin A. Our results suggest eMD-VR is a powerful tool for exploring protein energy landscapes in bioengineering efforts. </p>

scholarly journals Molecular dynamics of the immune checkpoint programmed cell death protein I, PD-1: conformational changes of the BC-loop upon binding of the ligand PD-L1 and the monoclonal antibody nivolumab

2020 ◽  
Vol 21 (S17) ◽  
Author(s):  
Bernhard Roither ◽  
Chris Oostenbrink ◽  
Wolfgang Schreiner
Keyword(s):  
Molecular Dynamics ◽  
Signal Transduction ◽  
T Cells ◽  
Cell Death ◽  
Programmed Cell Death ◽  
Conformational Changes ◽  
Immune Checkpoint ◽  
Active Conformation ◽  
Dynamics Simulations ◽  
Cell Death Protein

Abstract Background The immune checkpoint receptor programmed cell death protein I (PD-1) has been identified as a key target in immunotherapy. PD-1 reduces the risk of autoimmunity by inducing apoptosis in antigen-specific T cells upon interaction with programmed cell death protein ligand I (PD-L1). Various cancer types overexpress PD-L1 to evade the immune system by inducing apoptosis in tumor-specific CD8+ T cells. The clinically used blocking antibody nivolumab binds to PD-1 and inhibits the immunosuppressive interaction with PD-L1. Even though PD-1 is already used as a drug target, the exact mechanism of the receptor is still a matter of debate. For instance, it is hypothesized that the signal transduction is based on an active conformation of PD-1. Results Here we present the results of the first molecular dynamics simulations of PD-1 with a complete extracellular domain with a focus on the role of the BC-loop of PD-1 upon binding PD-L1 or nivolumab. We could demonstrate that the BC-loop can form three conformations. Nivolumab binds to the BC-loop according to the conformational selection model whereas PD-L1 induces allosterically a conformational change of the BC-loop. Conclusion Due to the structural differences of the BC-loop, a signal transduction based on active conformation cannot be ruled out. These findings will have an impact on drug design and will help to refine immunotherapy blocking antibodies.

Molecular Basis Explanation of Imatinib Resistance of Bcr-Abl Due to T315I and P-Loop Mutations from Molecular Dynamics Simulations.

Blood ◽  
2007 ◽  
Vol 110 (11) ◽  
pp. 2917-2917
Author(s):  
Tai-Sung Lee ◽  
Steven Potts ◽  
Hagop Kantarjian ◽  
Jorge Cortes ◽  
Francis Giles ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Electrostatic Interactions ◽  
Large Scale ◽  
Resistance Mechanism ◽  
Md Simulations ◽  
Resistance Mechanisms ◽  
Static Structure ◽  
Imatinib Resistance ◽  
Dynamics Simulations

Abstract Molecular dynamics (MD) simulations on the complex of imatinib with the wild-type, T315I, and other 10 P-loop mutants of the tyrosine kinase Bcr-Abl have been performed to study the imatinib resistance mechanism at the atomic level. MD simulations show that large scale computational simulations could offer insight information that a static structure or simple homology modeling methods cannot provide for studying the Bcr-Abl imatinib resistance problem, especially in the case of conformational changes due to remote mutations. By utilizing the Molecular Mechanics/Poisson-Boltzmann surface area (MM-PBSA) techniques and analyzing the interactions between imatinib and individual residues, imatinib resistance mechanisms not previously thought have been revealed. Non-directly contacted P-loop mutations either unfavorably change the direct electrostatic interactions with imatinib, or cause the conformational changes influencing the contact energies between imatinib and other non-P-loop residues. We demonstrate that imatinib resistance of T315I mainly comes from the breakdown of the interactions between imatinib and E286 and M290, contradictory to previously suggested that the missing hydrogen bonding is the main contribution. We also demonstrate that except for the mutations of the direct contact residues, such as L248 and Y253, the unfavorable electrostatic interaction between P-loop and imatinib is the main reason for resistance for the P-loop mutations. Furthermore, in Y255H, protonation of the histidin is essential for rendering this mutation resistant to Gleevec. Our results demonstrate that MD is a powerful way to verify and predict clinical response or resistance to imatinib and other potential drugs.

scholarly journals Insights into the substrate binding mechanism of SULT1A1 through molecular dynamics with excited normal modes simulations

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Balint Dudas ◽  
Daniel Toth ◽  
David Perahia ◽  
Arnaud B. Nicot ◽  
Erika Balog ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Conformational Changes ◽  
Molecular Mechanisms ◽  
Normal Modes ◽  
Md Simulations ◽  
Conformational Space ◽  
Metabolizing Enzymes ◽  
Protein Ligand Interactions ◽  
Ligand Interactions ◽  
Excited Normal Modes

AbstractSulfotransferases (SULTs) are phase II drug-metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3′-phosphoadenosine 5′-phosphosulfate (PAPS) to a substrate. It has been previously suggested that a considerable shift of SULT structure caused by PAPS binding could control the capability of SULT to bind large substrates. We employed molecular dynamics (MD) simulations and the recently developed approach of MD with excited normal modes (MDeNM) to elucidate molecular mechanisms guiding the recognition of diverse substrates and inhibitors by SULT1A1. MDeNM allowed exploring an extended conformational space of PAPS-bound SULT1A1, which has not been achieved up to now by using classical MD. The generated ensembles combined with docking of 132 SULT1A1 ligands shed new light on substrate and inhibitor binding mechanisms. Unexpectedly, our simulations and analyses on binding of the substrates estradiol and fulvestrant demonstrated that large conformational changes of the PAPS-bound SULT1A1 could occur independently of the co-factor movements that could be sufficient to accommodate large substrates as fulvestrant. Such structural displacements detected by the MDeNM simulations in the presence of the co-factor suggest that a wider range of drugs could be recognized by PAPS-bound SULT1A1 and highlight the utility of including MDeNM in protein–ligand interactions studies where major rearrangements are expected.

scholarly journals NAPS update: network analysis of molecular dynamics data and protein–nucleic acid complexes

2019 ◽  
Vol 47 (W1) ◽  
pp. W462-W470 ◽  
Author(s):  
Broto Chakrabarty ◽  
Varun Naganathan ◽  
Kanak Garg ◽  
Yash Agarwal ◽  
Nita Parekh
Keyword(s):  
Molecular Dynamics ◽  
Nucleic Acid ◽  
Network Analysis ◽  
Conformational Changes ◽  
Protein Complexes ◽  
Md Simulations ◽  
Network Construction ◽  
Dynamic Correlation ◽  
Residue Interaction ◽  
Protein Nucleic Acid

Abstract Network theory is now a method of choice to gain insights in understanding protein structure, folding and function. In combination with molecular dynamics (MD) simulations, it is an invaluable tool with widespread applications such as analyzing subtle conformational changes and flexibility regions in proteins, dynamic correlation analysis across distant regions for allosteric communications, in drug design to reveal alternative binding pockets for drugs, etc. Updated version of NAPS now facilitates network analysis of the complete repertoire of these biomolecules, i.e., proteins, protein–protein/nucleic acid complexes, MD trajectories, and RNA. Various options provided for analysis of MD trajectories include individual network construction and analysis of intermediate time-steps, comparative analysis of these networks, construction and analysis of average network of the ensemble of trajectories and dynamic cross-correlations. For protein–nucleic acid complexes, networks of the whole complex as well as that of the interface can be constructed and analyzed. For analysis of proteins, protein–protein complexes and MD trajectories, network construction based on inter-residue interaction energies with realistic edge-weights obtained from standard force fields is provided to capture the atomistic details. Updated version of NAPS also provides improved visualization features, interactive plots and bulk execution. URL: http://bioinf.iiit.ac.in/NAPS/

scholarly journals Effects of number of parallel runs and frequency of bias-strength replacement in generalized ensemble molecular dynamics simulations

2019 ◽  
Vol 1 ◽  
pp. e4 ◽  
Author(s):  
Takuya Shimato ◽  
Kota Kasahara ◽  
Junichi Higo ◽  
Takuya Takahashi
Keyword(s):  
Molecular Dynamics ◽  
Potential Energy ◽  
Molecular Dynamics Simulations ◽  
Conformational Changes ◽  
High Frequency ◽  
Md Simulations ◽  
Energy Space ◽  
Bias Potential ◽  
Ensemble Approach ◽  
Dynamics Simulations

Background The generalized ensemble approach with the molecular dynamics (MD) method has been widely utilized. This approach usually has two features. (i) A bias potential, whose strength is replaced during a simulation, is applied. (ii) Sampling can be performed by many parallel runs of simulations. Although the frequency of the bias-strength replacement and the number of parallel runs can be adjusted, the effects of these settings on the resultant ensemble remain unclear. Method In this study, we performed multicanonical MD simulations for a foldable mini-protein (Trp-cage) and two unstructured peptides (8- and 20-residue poly-glutamic acids) with various settings. Results As a result, running many short simulations yielded robust results for the Trp-cage model. Regarding the frequency of the bias-potential replacement, although using a high frequency enhanced the traversals in the potential energy space, it did not promote conformational changes in all the systems.

Sign in / Sign up

Export Citation Format

Share Document