Exploration micromechanism of VP35 IID interaction and recognition dsRNA: A molecular dynamics simulation

ABSTRACT Multifunctional viral protein (VP35) encoded by the highly pathogenic Ebola viruses (EBOVs) can antagonize host double‐stranded RNA (dsRNA) sensors and immune response because of the simultaneous recognition of dsRNA backbone and blunt ends. Mutation of select hydrophobic conserved basic residues within the VP35 inhibitory domain (IID) abrogates its dsRNA‐binding activity, and impairs VP35‐mediated interferon (IFN) antagonism. Herein the detailed binding mechanism between dsRNA and WT, single mutant, and double mutant were investigated by all‐atom molecular dynamics (MD) simulation and binding energy calculation. R312A/R322A double mutations results in a completely different binding site and orientation upon the structure analyses. The calculated binding free energy results reveal that R312A, R322A, and K339A single mutations decrease the binding free energies by 17.82, 13.18, and 13.68 kcal mol −1 , respectively. The binding energy decomposition indicates that the strong binding affinity of the key residues is mainly due to the contributions of electrostatic interactions in the gas phase, where come from the positively charged side chain and the negatively charged dsRNA backbone. R312A, R322A, and K339A single mutations have no significant effect on VP35 IID conformation, but the mutations influence the contributions of electrostatic interactions in the gas phase. The calculated results reveal that end‐cap residues which mainly contribute VDW interactions can recognize and capture dsRNA blunt ends, and the central basic residues (R312, R322, and K339) which mainly contribute favorable electrostatic interactions with dsRNA backbone can fix dsRNA binding site and orientation. Proteins 2017; 85:1008–1023. © 2017 Wiley Periodicals, Inc.