Abstract:Despite the progress in prediction of protein complexes over the last decade, recent blind protein complex structure prediction challenges revealed limited success rates (less than 20% models with DockQ score > 0.4) on targets that exhibit significant conformational change upon binding. To overcome limitations in capturing backbone motions, we developed a new, aggressive sampling method that incorporates temperature replica exchange Monte Carlo (T-REMC) and conformational sampling techniques within docking protocols in Rosetta. Our method, ReplicaDock 2.0, mimics induced-fit mechanism of protein binding to sample backbone motions across putative interface residues on-the-fly, thereby recapitulating binding-partner induced conformational changes. Furthermore, ReplicaDock 2.0 clocks in at 150-500 CPU hours per target (protein-size dependent); a runtime that is significantly faster than Molecular Dynamics based approaches. For a benchmark set of 88 proteins with moderate to high flexibility (unbound-to-bound iRMSD over 1.2 Å), ReplicaDock 2.0 successfully docks 61% of moderately flexible complexes and 35% of highly flexible complexes. Additionally, we demonstrate that by biasing backbone sampling particularly towards residues comprising flexible loops or hinge domains, highly flexible targets can be predicted to under 2 Å accuracy. This indicates that additional gains are possible when mobile protein segments are known. Significance Statement Proteins bind each other in a highly specific and regulated manner, and these associated dynamics of binding are intimately linked to their function. Conventional techniques of structure determination such as cryo-EM, X-ray crystallography and NMR are time-consuming and arduous. Using a temperature-replica exchange Monte Carlo approach that mimics the kinetic mechanism of “induced fit” binding, we improved prediction of protein complex structures, particularly for targets that exhibit considerable conformational changes upon binding (Interface root mean square deviation (unbound-bound) > 1.2 Å. Capturing these binding-induced conformational changes in proteins can aid us in better understanding biological mechanisms and suggest intervention strategies for disease mechanisms.

Replica Exchange with Solute Tempering for Protein Conformational Sampling

Replica exchange molecular dynamics method for protein folding simulation.

Simulated Solute Tempering 2

Replica Exchange Molecular Dynamics: A Practical Application Protocol with Solutions to Common Problems and a Peptide Aggregation and Self-Assembly Example.

Replica exchange with guided annealing for accelerated sampling of disordered protein conformations

Enhanced Conformational Sampling using Replica Exchange with Collective-Variable Tempering

Boosting the Conformational Sampling by Combining Replica Exchange with Solute Tempering and Well-Sliced Metadynamics

A Second Look at Canonical Sampling of Biomolecules using Replica Exchange Simulation

Mapping transiently formed and sparsely populated conformations on a complex energy landscape

Self-assembly of a model supramolecular polymer studied by replica exchange with solute tempering

A Selective Integrated Tempering Method.

Efficiently determining membrane-bound conformations of peripheral membrane proteins using replica exchange with hybrid tempering

Constant pH Replica Exchange Molecular Dynamics in Biomolecules Using a Discrete Protonation Model

Replica exchange of expanded ensembles: A generalized ensemble approach with enhanced flexibility and parallelizability

Induced fit with replica exchange improves protein complex structure prediction

Parallel Implementation of the Replica Exchange Molecular Dynamics Algorithm on Blue Gene/L

Prediction of protein–protein complexes using replica exchange with repulsive scaling

Combining Enhanced Sampling with Experiment Directed Simulation of the GYG peptide

Replica Exchange Simulations of Transient Encounter Complexes in Protein–protein Association

Binding Energy Distribution Analysis Method: Hamiltonian Replica Exchange with Torsional Flattening for Binding Mode Prediction and Binding Free Energy Estimation