Abstract:We developed an efficient Monte Carlo Simulated Annealing (MCSA) program for modeling protein loops with high speed. The total conformational energy in each step of MCSA simulation consists of two parts: the nonbonded atomic interaction represented by a simple soft-sphere potential and the harmonic distance constraint to ensure the smooth connection of the loop segment to the rest of the protein structure. The soft-sphere potential was a simplified potential that has been successfully used by the authors in modeling the carbohydrate part of glycoprotein systems [H. Zhang, Y. Yang, L. Lai, and Y. Tang (1996), Carbohydrate Research, Vol. 284, pp. 25–34]. It only considers the purely repulsive steric interactions to avoid artificial attractive forces between atoms in the absence of solvent molecules. The N-terminal of the loop segment was connected to the bulk protein part, and two dummy main-chain atoms N and Cα immediately following the C-terminal of the loop segment were constrained to their real positions in the protein structure, which not only assures the correct geometry of loop-protein connection but also is more rigorous than the previous work. To improve the speed, two strategies, the local region method and grid-mapping method, were devised to accelerate the computation of environmental interaction that is responsible for the major part of the computing consumption. The grid-mapping method can reduce computational time dramatically. Conformations with rational steric packing and smooth connection to the rest of the protein structure were generated by the MCSA program, and then were refined by the empirical force field CHARMm [B. R. Brook, R. E. Braccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus (1983), Journal of Computational Chemistry, Vol. 4, pp. 187–217]. Bovine pancreatic trypsin inhibitor (BPTI) was used as an example to test the ability of loop modeling of the method, and five loops in BPTI were calculated. Conformations close to the crystal structure were generated for all of them. With the criteria of CHARMm energy, near-native conformations can be selected, for example, the backbone rms deviation 0.93 A from the crystal structure was gotten for the longest 9-residue loop. © 1997 John Wiley & Sons, Inc.

A Fast and Efficient Program for Modeling Protein Loops

Calculation of Protein Surface Loops Using Monte-Carlo Simulated Annealing Simulation

CGLFold: a contact-assisted <i>de novo</i> protein structure prediction using global exploration and loop perturbation sampling algorithm

Accurate Structure Prediction for Protein Loops Based on Molecular Dynamics Simulations with RSFF2C

Prediction of protein loop conformations using multiscale modeling methods with physical energy scoring functions

Monte Carlo Sampling of Protein Folding by Combining an All-Atom Physics-Based Model with a Native State Bias

Comprehensive assessment of protein loop modeling programs on large-scale datasets: prediction accuracy and efficiency

Efficiently Explore The Energy Landscape Of Proteins In Molecular Dynamics Simulations By Amplifying Collective Motions

Protein Conformational Transitions from All-Atom Adaptively Biased Path Optimization

A generative probabilistic model for Loop modeling

Highly accurate and efficient deep learning paradigm for full-atom protein loop modeling with KarmaLoop

Adaptive Simulated Annealing Based Protein Loop Modeling of Neurotoxins

Modeling Protein Loop Structure by Cyclic Coordinate Descent-based Approach.

LoopSage: An energy-based Monte Carlo approach for the loop extrusion modeling of chromatin

Computational Model for Protein Unfolding Simulation.

Clustering and Percolation in Protein Loop Structures

Discovering Protein Conformational Flexibility Through Artificial Intelligence Aided Molecular Dynamics.

Sequence-based Protein Structure Optimization Using Enhanced Simulated Annealing Algorithm on a Coarse-Grained Model

LoopSage: An Energy-Based Monte Carlo approach for the Loop Extrusion Modelling of Chromatin

Multiple Simulated Annealing-Molecular Dynamics (MSA-MD) for Conformational Space Search of Peptide and Miniprotein

Simulating Protein Conformations through Global Optimization