Abstract:An efficient methodology, further referred to as ICM, for versatile modeling operations and global energy optimization on arbitrarily fixed multimolecular systems is described. It is aimed at protein structure prediction, homology modeling, molecular docking, nuclear magnetic resonance (NMR) structure determination, and protein design. The method uses and further develops a previously introduced approach to model biomolecular structures in which bond lengths, bond angles, and torsion angles are considered as independent variables, any subset of them being fixed. Here we simplify and generalize the basic description of the system, introduce the variable dihedral phase angle, and allow arbitrary connections of the molecules and conventional definition of the torsion angles. Algorithms for calculation of energy derivatives with respect to internal variables in the topological tree of the system and for rapid evaluation of accessible surface are presented. Multidimensional variable restraints are proposed to represent the statistical information about the torsion angle distributions in proteins. To incorporate complex energy terms as solvation energy and electrostatics into a structure prediction procedure, a “double‐energy” Monte Carlo minimization procedure in which these terms are omitted during the minimization stage of the random step and included for the comparison with the previous conformation in a Markov chain is proposed and justified. The ICM method is applied successfully to a molecular docking problem. The procedure finds the correct parallel arrangement of two rigid helixes from a leucine zipper domain as the lowest‐energy conformation (0.5 Å root mean square, rms, deviation from the native structure) starting from completely random configuration. Structures with antiparallel helixes or helixes staggered by one helix turn had energies higher by about 7 or 9 kcal/mol, respectively. Soft docking was also attempted. A docking procedure allowing side‐chain flexibility also converged to the parallel configuration starting from the helixes optimized individually. To justdy an internal coordinate approach to the structure prediction as opposed to a Cartesian one, energy hypersurfaces around the native structure of the squash seeds trypsin inhibitor were studied. Torsion angle minimization from the optimal conformation randomly distorted up to the rms deviation of 2.2 Å or angular rms deviation of l0° restored the native conformation in most cases. In contrast, Cartesian coordinate minimization did not reach the minimum from deviations as small as 0.3 Å or 2°. We conclude that the most promising detailed approach to the protein‐folding problem would consist of some coarse global sampling strategy combined with the local energy minimization in the torsion coordinate space. © 1994 by John Wiley & Sons, Inc.

TRAJELIX: A Computational Tool for the Geometric Characterization of Protein Helices During Molecular Dynamics Simulations

Conformational Properties and Elastic Behavior of Protein-Like Lattice Polymers

Estimating conformational heterogeneity of tryptophan synthase with a template‐based Alphafold2 approach

PROTERAN: animated terrain evolution for visual analysis of patterns in protein folding trajectory.

Bendix: intuitive helix geometry analysis and abstraction

TriMem: A parallelized hybrid Monte Carlo software for efficient simulations of lipid membranes

Geo-Measures: A PyMOL plugin for protein structure ensembles analysis

Molecular Dynamics and Ion Mobility Spectrometry Study of Model Β-Hairpin Peptide, Trpzip1.

HTMoL: full-stack solution for remote access, visualization, and analysis of Molecular Dynamics trajectory data

Protalign: a 3-dimensional protein alignment assessment tool

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

ICM—A new method for protein modeling and design: Applications to docking and structure prediction from the distorted native conformation

Disentangling folding from energetic traps in simulations of disordered proteins

Folding pathways of a helix-turn-helix model protein

A Structural-informatics Approach for Tracing β-Sheets: Building Pseudo-Cα Traces for β-Strands in Intermediate-resolution Density Maps

Assessing the helical stability of polyXYs at the boundaries of Intrinsically Disordered Regions with MD simulations

Exploring the trimerization process of a transmembrane helix with an ionizable residue by molecular dynamics simulations: a case study of transmembrane domain 5 of LMP-1

Structural Refinement of Proteins by Restrained Molecular Dynamics Simulations with Non-Interacting Molecular Fragments

Riemannian Geometry and Molecular Surfaces I: Spectrum of the Laplacian

MoDEL (Molecular Dynamics Extended Library): A Database of Atomistic Molecular Dynamics Trajectories

Efficient and Reproducible Folding Simulations of the Trp-cage Protein with Multiscale Molecular Dynamics