The Hybrid Idea of (Energy Minimization) Optimization Methods Applied to Study PrionProtein Structures Focusing on the beta2-alpha2 Loop

Abstract

In molecular mechanics, current generation potential energy functions provide a reasonably good compromise between accuracy and effectiveness. This paper firstly reviewed several most commonly used classical potential energy functions and their optimization methods used for energy minimization. To minimize a potential energy function, about 95% efforts are spent on the Lennard-Jones potential of van der Waals interactions; we also give a detailed review on some effective computational optimization methods in the Cambridge Cluster Database to solve the problem of Lennard-Jones clusters. From the reviews, we found the hybrid idea of optimization methods is effective, necessary and efficient for solving the potential energy minimization problem and the Lennard-Jones clusters problem. An application to prion protein structures is then done by the hybrid idea. We focus on the β2-α2 loop of prion protein structures, and we found (i) the species that has the clearly and highly ordered β2-α2 loop usually owns a 3 10 -helix in this loop, (ii) a "π-circle" Y128-F175-Y218-Y163-F175-Y169-R164-Y128(-Y162) is around the β2-α2 loop.

Role of Oxidized Gly25, Gly29, and Gly33 Residues on the Interactions of Aβ1–42 with Lipid Membranes

“…The systems were energy minimized using the steepest descent algorithm to remove all atomic clashes. , This was followed by an equilibration under NVT conditions for 1 ns, where the reference temperature of 310 K was regulated with the velocity-rescale thermostat and the time constant was set to 0.1 ps. During minimization and initial equilibration stages, the heavy atoms of protein and lipids were subjected to position restraints with a force constant of 1000 kJ/mol· nm 2 .…”

A Hairpin Motif in the Amyloid-β Peptide Is Important for Formation of Disease-Related Oligomers

“…The initial structures were placed in a cubic box; water molecules and 150 mM NaCl were added, and the systems were neutralized by inserting extra Na + ions. After that, the systems underwent energy minimization using the steepest descent algorithm to remove clashes between atoms, followed by two equilibration steps, each for 1 ns, under canonical ( NVT ) and isobaric–isothermic ( NPT ) ensemble conditions. During the equilibration, the pressure was kept at 1.0 bar using the Parrinello–Rahman pressure coupling method, , and the temperature was kept at 298 K using the velocity-rescale thermostat method .…”

Comparative molecular dynamics simulations of pathogenic and non-pathogenic huntingtin protein monomers and dimers

“…The initial (already collapsed) monomer structures were solvated in a ∼2,000 nm 3 cubic box, including 150 mM Na + and Cl − ions and overall system neutralization, leading to about 200,000 atoms. The energy of each system was initially minimized using the steepest descent algorithm ( Müller and Brown, 1979 ; Zhang, 2015 ), followed by 1.0 ns equilibration and then the production runs under NPT ensemble conditions. The pressure and temperature were maintained at 1.0 bar and 298 K using the Parrinello-Rahman pressure coupling method ( Parrinello and Rahman, 1981 ; Parrinello and Rahman, 1982 ) and a velocity-rescale thermostat method ( Bussi et al, 2007 ), respectively.…”