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.

Prediction of transmembrane helix configurations by replica-exchange simulations

Structure Prediction of Membrane Proteins

Replica exchange molecular dynamics method for protein folding simulation.

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

Validation of Molecular Dynamics Simulations for Prediction of Three-Dimensional Structures of Small Proteins

Induced fit with replica exchange improves protein complex structure prediction

Predicting the helix-helix interactions from correlated residue mutations.

Simulating Protein Fold Switching by Replica-Exchange-with-Tunneling.

Computational Prediction of Heteromeric Protein Complex Disassembly Order with Hybrid Monte Carlo/Molecular Dynamics Simulation

Free-energy landscapes of transmembrane homodimers by bias-exchange adaptively biased molecular dynamics

Transmembrane Helices Topology Prediction: Using a Simplified Transmembrane HMM

[A Novel Segment-Training Algorithm for Transmembrane Helices Prediction].

Protein Secondary Structure Prediction Based on Ramachandran Maps

Structure Prediction for the Helical Skeletons Detected from the Low Resolution Protein Density Map

Construction of an Intermediate-Resolution Lattice Model and Re-Examination of the Helix-Coil Transition: a Dynamic Monte Carlo Simulation.

Evaluation for prediction methods of α-helix transmembrane protein topology

Structure calculation strategies for helical membrane proteins; a comparison study

Prediction of atomic structure from sequence for double helical DNA oligomers

ANALYSIS AND PREDICTION OF PROTEIN TRANSMEMBRANE HELICAL REGIONS BY WAVELET TRANSFORMATION

TM Finder: A prediction program for transmembrane protein segments using a combination of hydrophobicity and nonpolar phase helicity scales

Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry