A wrench in the works of human acetylcholinesterase: soman induced conformational changes revealed by molecular dynamics simulations.

Irreversible inactivation of human acetylcholinesterase (hAChE) by organophosphorous pesticides (OPs) and chemical weapon agents (CWA) has severe morbidity and mortality consequences. We present data from quantum mechanics/molecular mechanics (QM/MM) and 80 classical molecular dynamics (MD) simulations of the apo and soman-adducted forms of hAChE to investigate the effects on the dynamics and protein structure when the catalytic Serine 203 is phosphonylated. We find that the soman phosphonylation of the active site Ser203 follows a water assisted addition-elimination mechanism with the elimination of the fluoride ion being the highest energy barrier at 6.

Computational Analysis of a Zn-Bound Tris(imidazolyl) Calix[6]arene Aqua Complex: Toward Incorporating Second-Coordination Sphere Effects into Carbonic Anhydrase Biomimetics.

Molecular dynamics simulations and quantum-mechanical calculations were performed to characterize a supramolecular tris(imidazolyl) calix[6]arene Zn(2+) aqua complex, as a biomimetic model for the catalyzed hydration of carbon dioxide to bicarbonate, H2O + CO2 → H(+) + HCO3(-). On the basis of potential-of-mean-force (PMF) calculations, stable conformations had distorted 3-fold symmetry and supported either one or zero encapsulated water molecules. The conformation with an encapsulated water molecule is calculated to be lower in free energy than the conformation with an empty cavity (ΔG = 1.

Dynamic linear response theory for conformational relaxation of proteins.

Dynamic linear response theory is adapted to the problem of computing the time evolution of the atomic coordinates of a protein in response to the unbinding of a ligand molecule from a binding pocket within the protein. When the ligand dissociates from the molecule, the protein molecule finds itself out of equilibrium and its configuration begins to change, ultimately coming to a new stable configuration corresponding to equilibrium in a force field that lacks the ligand-protein interaction terms. Dynamic linear response theory (LRT) relates the nonequilibrium motion of the protein atoms that ensues after the ligand molecule dissociates to equilibrium dynamics in the force field, or equivalently, on the potential energy surface (PES) relevant to the unliganded protein.

Langevin dynamics of molecules with internal rigid fragments in the harmonic regime.

An approximation scheme is developed to compute Brownian motion according to the Langevin equation for a molecular system moving in a harmonic force field (corresponding to a quadratic potential energy surface) and characterized by one or more rigid internal fragments. This scheme, which relies on elements of the rotation translation block (RTB) method for computing vibrational normal modes of large molecules developed by Sanejouand and co-workers [Biopolymers 34, 759 (1994); Proteins: Struct., Funct.