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.

Dealing with induced fit, conformational selection and secondary poses in molecular dynamics simulations for reliable free energy predictions

Enhanced molecular dynamics sampling of drug target conformations

Multistate Method to Efficiently Account for Tautomerism and Protonation in Alchemical Free-Energy Calculations

Predicting the bioactive conformations of macrocycles: a molecular dynamics-based docking procedure with DynaDock

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

Sampling of conformational ensemble for virtual screening using molecular dynamics simulations and normal mode analysis

Enhanced Sampling in Molecular Dynamics Simulations: How Many MD Snapshots can be Needed to Reproduce the Biological Behavior?

Molecular dynamics-based approaches for enhanced sampling of long-time, large-scale conformational changes in biomolecules

Validating Small-Molecule Force Fields for Macrocyclic Compounds Using NMR Data in Different Solvents

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

Statistical accuracy of molecular dynamics-based methods for sampling conformational ensembles of disordered proteins

A Discard-and-Restart MD algorithm for the sampling of realistic protein transition states and enhance structure-based drug discovery

A Different Kind of Restraint Suitable for Molecular Dynamics Simulations

Enhanced sampling techniques in molecular dynamics simulations of biological systems

Induced fit with replica exchange improves protein complex structure prediction

Identifying ligand binding sites and poses using GPU-accelerated Hamiltonian replica exchange molecular dynamics

CDOCKER and $$\Lambda$$-Dynamics for Prospective Prediction in D3R Grand Challenge 2

Computational Methods for Studying Conformational Behaviors of Cyclic Peptides.

Impact of replicas and simulation length on in silico behaviors of a protein domain