skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Sign Gradient Descent Algorithms for Kinetostatic Protein Folding
This paper proposes a sign gradient descent (SGD) algorithm for predicting the three-dimensional folded protein molecule structures under the kinetostatic compliance method (KCM). In the KCM framework, which can be used to simulate the range of motion of peptide-based nanorobots/nanomachines, protein molecules are modeled as a large number of rigid nano-linkages that form a kinematic mechanism under motion constraints imposed by chemical bonds while folding under the kinetostatic effect of nonlinear interatomic force fields. In a departure from the conventional successive kinetostatic fold compliance framework, the proposed SGD-based iterative algorithm in this paper results in convergence to the local minima of the free energy of protein molecules corresponding to their final folded conformations in a faster and more robust manner. KCM-based folding dynamics simulations of the backbone chains of protein molecules demonstrate the effectiveness of the proposed algorithm.  more » « less
Award ID(s):
2153744
PAR ID:
10523223
Author(s) / Creator(s):
;
Publisher / Repository:
IEEE
Date Published:
ISBN:
979-8-3503-3039-7
Page Range / eLocation ID:
1 to 6
Format(s):
Medium: X
Location:
2023 International Conference on Manipulation, Automation and Robotics at Small Scales (MARSS 2023)
Sponsoring Org:
National Science Foundation
More Like this
  1. Numerical simulations of protein folding enable the design of protein-based nanomachines and nanorobots by predicting folded three-dimensional protein structures with high accuracy and revealing the protein conformation transitions during folding and unfolding. In the kinetostatic compliance method (KCM) for folding simulations, protein molecules are represented as ensembles of rigid nano-linkages connected by chemical bonds, and the folding process is driven by the kinetostatic influence of nonlinear interatomic force fields until the system converges to a free-energy minimum of the protein. Despite its strengths, the conventional KCM framework demands an excessive number of iterations to reach folded protein conformations, with each iteration requiring costly computations of interatomic force fields. To address these limitations, this work introduces a family of sign gradient descent (SGD) algorithms for predicting folded protein structures. Unlike the heuristic-based iterations of the conventional KCM framework, the proposed SGD algorithms rely on the sign of the free-energy gradient to guide the kinetostatic folding process. Owing to their faster and more robust convergence, the proposed SGD-based algorithms reduce the computational burden of interatomic force field evaluations required to reach folded conformations. Their effectiveness is demonstrated through numerical simulations of KCM-based folding of protein backbone chains. 
    more » « less
  2. The kinetostatic compliance method (KCM) models protein molecules as nanomechanisms consisting of numerous rigid peptide plane linkages. These linkages articulate with respect to each other through changes in the molecule dihedral angles, resulting in a kinematic mechanism with hyper degrees of freedom. Within the KCM framework, nonlinear interatomic forces drive protein folding by guiding the molecule’s dihedral angle vector towards its lowest energy state in a kinetostatic manner. This paper proposes a numerical integrator that is well suited to KCM-based protein folding and overcomes the limitations of traditional explicit Euler methods with fixed step size. Our proposed integration scheme is based on pseudo-transient continuation with an adaptive step size updating rule that can efficiently compute protein folding pathways, namely, the transient three-dimensional configurations of protein molecules during folding. Numerical simulations utilizing the KCM approach on protein backbones confirm the effectiveness of the proposed integrator. 
    more » « less
  3. Abstract This paper addresses the problem of algorithmic prediction of protein folding pathways, namely, the transient three-dimensional conformations of protein molecules during folding, under constrained rates of entropy change. We formulate the physics-based prediction of folding pathways as a control synthesis problem, where the control inputs guide the protein folding simulations. These folding control inputs are obtained from largescale trust-region subproblems (TRS) utilizing a computationally efficient algorithm with no need for outer iterations. The proposed control synthesis approach, which leverages the solutions obtained from a special generalized eigenvalue problem, avoids potentially cumbersome and unpredictable iterative computations at each protein conformation. Moreover, the TRS-based control inputs align the closed-loop dynamics closely with the kinetostatic compliance method (KCM) reference vector field while satisfying ellipsoidal constraints on the folding control inputs. Finally, we provide conditions for existence and uniqueness of the resulting closed-loop solutions, which are the protein folding pathways under constraints on the rate of entropy change. Numerical simulations utilizing the KCM approach on protein backbones confirm the effectiveness of the proposed framework. 
    more » « less
  4. Understanding the process of protein unfolding plays a crucial role in various applications such as design of folding-based protein engines. Using the well-established kinetostatic compliance (KCM)-based method for modeling of protein conformation dynamics and a recent nonlinear control theoretic approach to KCM-based protein folding, this letter formulates protein unfolding as a destabilizing control analysis/synthesis problem. In light of this formulation, it is shown that the Chetaev instability framework can be used to investigate the KCM-based unfolding dynamics. In particular, a Chetaev function for analysis of unfolding dynamics under the effect of optical tweezers and a class of control Chetaev functions for synthesizing control inputs that elongate protein strands from their folded conformations are presented. Based on the presented control Chetaev function, an unfolding input is derived from the Artstein-Sontag universal formula and the results are compared against optical tweezer-based unfolding. 
    more » « less
  5. This paper investigates the problem of prediction of protein molecule folding pathways under entropy-loss constraints by formulating a control synthesis problem whose solutions are obtained by solving large-scale quadratic programming (QP) optimizations with nonlinear constraints. The utilized non-iterative and computationally efficient algorithm, which is based on solving generalized eigenvalue problems, prevents an unpredictable and potentially large number of iterations at each protein conformation for computing the folding control inputs. The synthesized control inputs remain close to the renowned kinetostatic compliance method (KCM) reference vector field while satisfying proper quadratic inequality constraints that limit the rate of molecule entropy-loss during folding. 
    more » « less