Dataset Information

A stepwise docking molecular dynamics approach for simulating antibody recognition with substantial conformational changes.

ABSTRACT: Conformational changes or rearrangements are common events during inter-biomolecular recognition. Tracking these changes are essential for exploring the allosteric mechanism and it is usually achieved by molecular dynamics simulation in silico. We previously identified a broad-neutralizing antibody against H5 influenza virus, 13D4, and solved the crystal structures of the free 13D4 Fab and its complex with hemagglutinin (HA). Structural comparison of the unbound and bound 13D4 Fabs showed that the heavy chain complementarity-determining region 3 (HCDR3) undergoes a substantial conformational rearrangement when it recognizes the receptor-binding site (RBS). Here, we used molecular dynamics (MD) to simulate the conformational changes that occur during antibody recognition. We showed that neither conventional MD nor steered MD could recapitulate the loop fitting of the RBS structure contour. Consequently, to simulate these challenging conformational changes, we engaged a stepwise docking MD method that allowed for the gradual docking of the ligand to receptor. This new method recapitulates the bound shape of the HCDR3 and provides the best approximation of the shape rendered by the co-crystal structure, with an RMSD of 0.926 Å. This strategy affords a flexible MD approach for simulating complicated conformational changes that occur during molecular recognition, and helps to provide an understanding of the involved allosteric mechanism.

SUBMITTER: Huang Y

PROVIDER: S-EPMC8816672 | biostudies-literature | 2022

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

A stepwise docking molecular dynamics approach for simulating antibody recognition with substantial conformational changes.

Huang Yang Y Li Zizhen Z Hong Qiyang Q Zhou Lizhi L Ma Yue Y Hu Yisha Y Xin Jiabao J Li Tingting T Kong Zhibo Z Zheng Qingbing Q Chen Yixin Y Zhao Qinjian Q Gu Ying Y Zhang Jun J Wang Yingbin Y Yu Hai H Li Shaowei S Xia Ningshao N

Computational and structural biotechnology journal 20220118

Conformational changes or rearrangements are common events during inter-biomolecular recognition. Tracking these changes are essential for exploring the allosteric mechanism and it is usually achieved by molecular dynamics simulation <i>in silico</i>. We previously identified a broad-neutralizing antibody against H5 influenza virus, 13D4, and solved the crystal structures of the free 13D4 Fab and its complex with hemagglutinin (HA). Structural comparison of the unbound and bound 13D4 Fabs showed ...[more]

PMID: 35198128

Similar Datasets

Project description:Conformational sampling of complex biomolecules is an emerging frontier in drug discovery. Advances in lab-based structural biology and related computational approaches like AlphaFold have made great strides in obtaining static protein structures for biologically relevant targets. However, biology is in constant motion, and many important biological processes rely on conformationally driven events. Conventional molecular dynamics (MD) simulations run on standard hardware are impractical for many drug design projects, where conformationally driven biological events can take microseconds to milliseconds or longer. An alternative approach is to focus the search on a limited region of conformational space defined by a putative reaction coordinate (i.e., path collective variable). The search space is typically limited by applying restraints, which can be guided by insights about the underlying biological process of interest. The challenge is striking a balance between the degree to which the system is constrained and still allowing for natural motions along the path. A plethora of restraints exist to limit the size of conformational search space, although each has drawbacks when simulating complex biological motions. In this work, we present a three-stage procedure to construct realistic path collective variables (PCVs) and introduce a new kind of barrier restraint that is particularly well suited for complex conformationally driven biological events, such as allosteric modulations and conformational signaling. The PCV presented here is all-atom (as opposed to C-alpha or backbone only) and is derived from all-atom MD trajectory frames. The new restraint relies on a barrier function (specifically, the scaled reciprocal function), which we show is particularly beneficial in the context of molecular dynamics, where near-hard-wall restraints are needed with zero tolerance to restraint violation. We have implemented our PCV and barrier restraint within a hybrid sampling framework that combines well-tempered metadynamics and extended-Lagrangian adaptive biasing force (meta-eABF). We use three particular examples of high pharmaceutical interest to demonstrate the value of this approach: (1) sampling the distance from ubiquitin to a protein of interest within the supramolecular cullin-RING ligase complex, (2) stabilizing the wild-type conformation of the oncogenic mutant JAK2-V617F pseudokinase domain, and (3) inducing an activated state of the stimulator of interferon genes (STING) protein observed upon ligand binding. For examples 2 and 3, we present statistical analysis of meta-eABF free energy estimates and, for each case, code for reproducing this work.

Dataset Information

A stepwise docking molecular dynamics approach for simulating antibody recognition with substantial conformational changes.

Publications

A stepwise docking molecular dynamics approach for simulating antibody recognition with substantial conformational changes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets