Complex intramolecular mechanics of G-actin--an elastic network study

Computational Ways to Enhance Protein Inhibitor Design

Two new computational approaches are described to aid in the design of new peptide-based drugs by evaluating ensembles of protein structures from their dynamics and through the assessing of structures using empirical contact potential. These approaches build on the concept that conformational variability can aid in the binding process and, for disordered proteins, can even facilitate the binding of more diverse ligands. This latter consideration indicates that such a design process should be less restrictive so that multiple inhibitors might be effective. The example chosen here focuses on proteins/peptides that bind to hemagglutinin (HA) to block the large-scale conformational change for activation. Variability in the conformations is considered from sets of experimental structures, or as an alternative, from their simple computed dynamics; the set of designe peptides/small proteins from the David Baker lab designed to bind to hemagglutinin, is the large set considered and is assessed with the new empirical contact potentials.

Detection of Single Molecules Using Stochastic Resonance of Bistable Oligomers

Ultra-sensitive elements for nanoscale devices capable of detecting single molecules are in demand for many important applications. It is generally accepted that the inevitable stochastic disturbance of a sensing element by its surroundings will limit detection at the molecular level. However, a phenomenon exists (stochastic resonance) in which the environmental noise acts abnormally: it amplifies, rather than distorts, a weak signal. Stochastic resonance is inherent in non-linear bistable systems with criticality at which the bistability emerges. Our computer simulations have shown that the large-scale conformational dynamics of a short oligomeric fragment of thermosrespective polymer, poly-N-isopropylmethacrylamid, resemble the mechanical movement of nonlinear bistable systems. The oligomers we have studied demonstrate spontaneous vibrations and stochastic resonance activated by conventional thermal noise. We have observed reasonable shifts of the spontaneous vibrations and stochastic resonance modes when attaching an analyte molecule to the oligomer. Our simulations have shown that spontaneous vibrations and stochastic resonance of the bistable thermoresponsive oligomers are sensitive to both the analyte molecular mass and the binding affinity. All these effects indicate that the oligomers with mechanic-like bistability may be utilized as ultrasensitive operational units capable of detecting single molecules.

Simple mechanics of protein machines

While belonging to the nanoscale, protein machines are so complex that tracing even a small fraction of their cycle requires weeks of calculations on supercomputers. Surprisingly, many aspects of their operation can be however already reproduced by using very simple mechanical models of elastic networks. The analysis suggests that, similar to other self-organized complex systems, functional collective dynamics in such proteins is effectively reduced to a low-dimensional attractive manifold.

Coarse-Grained Protein Dynamics Studies Using Elastic Network Models

Elastic networks have been used as simple models of proteins to study their slow structural dynamics. They consist of point-like particles connected by linear Hookean springs and hence are convenient for linear normal mode analysis around a given reference structure. Furthermore, dynamic simulations using these models can provide new insights. As the computational cost associated with these models is considerably lower compared to that of all-atom models, they are also convenient for comparative studies between multiple protein structures. In this review, we introduce examples of coarse-grained molecular dynamics studies using elastic network models and their derivatives, focusing on the nonlinear phenomena, and discuss their applicability to large-scale macromolecular assemblies.

Fractal globules: a new approach to artificial molecular machines

The over-damped relaxation of elastic networks constructed by contact maps of hierarchically folded fractal (crumpled) polymer globules was investigated in detail. It was found that the relaxation dynamics of an anisotropic fractal globule is very similar to the behavior of biological molecular machines like motor proteins. When it is perturbed, the system quickly relaxes to a low-dimensional manifold, M, with a large basin of attraction and then slowly approaches equilibrium, not escaping M. Taking these properties into account, it is suggested that fractal globules, even those made by synthetic polymers, are artificial molecular machines that can transform perturbations into directed quasimechanical motion along a defined path.

Structure-Encoded Global Motions and Their Role in Mediating Protein-Substrate Interactions

Recent structure-based computational studies suggest that, in contrast to the classical description of equilibrium fluctuations as wigglings and jigglings, proteins have access to well-defined spectra of collective motions, called intrinsic dynamics, encoded by their structure under native state conditions. In particular, the global modes of motions (at the low frequency end of the spectrum) are shown by multiple studies to be highly robust to minor differences in the structure or to detailed interactions at the atomic level. These modes, encoded by the overall fold, usually define the mechanisms of interactions with substrates. They can be estimated by low-resolution models such as the elastic network models (ENMs) exclusively based on interresidue contact topology. The ability of ENMs to efficiently assess the global motions intrinsically favored by the overall fold as well as the relevance of these predictions to the dominant changes in structure experimentally observed for a given protein in the presence of different substrates suggest that the intrinsic dynamics plays a role in mediating protein-substrate interactions. These observations underscore the functional significance of structure-encoded dynamics, or the importance of the predisposition to favor functional global modes in the evolutionary selection of structures.

Multiscale impact of nucleotides and cations on the conformational equilibrium, elasticity and rheology of actin filaments and crosslinked networks

Cells are able to respond to mechanical forces and deformations. The actin cytoskeleton, a highly dynamic scaffolding structure, plays an important role in cell mechano-sensing. Thus, understanding rheological behaviors of the actin cytoskeleton is critical for delineating mechanical behaviors of cells. The actin cytoskeleton consists of interconnected actin filaments (F-actin) that form via self-assembly of actin monomers. It has been shown that molecular changes of the monomer subunits impact the rigidity of F-actin. However, it remains inconclusive whether or not the molecular changes can propagate to the network level and thus alter the rheological properties of actin networks. Here, we focus on how cation binding and nucleotide state tune the molecular conformation and rigidity of F-actin and a representative rheological behavior of actin networks, strain-stiffening. We employ a multiscale approach by combining established computational techniques: molecular dynamics, normal mode analysis and Brownian dynamics. Our findings indicate that different combinations of nucleotide (ATP, ADP or ADP-Pi) and cation [Formula: see text] or [Formula: see text] at one or multiple sites) binding change the molecular conformation of F-actin by varying inter- and intra-strand interactions which bridge adjacent subunits between and within F-actin helical strands. This is reflected in the rigidity of actin filaments against bending and stretching. We found that differences in extension and bending rigidity of F-actin induced by cation binding to the low-, intermediate- and high-affinity sites vary the strain-stiffening response of actin networks crosslinked by rigid crosslinkers, such as scruin, whereas they minimally impact the strain-stiffening response when compliant crosslinkers, such as filamin A or [Formula: see text]-actinin, are used.

Coarse-grain simulations of active molecular machines in lipid bilayers

A coarse-grain method for simulations of the dynamics of active protein inclusions in lipid bilayers is described. It combines the previously proposed hybrid simulations of bilayers [M.-J. Huang, R. Kapral, A. S. Mikhailov, and H.-Y. Chen, J. Chem. Phys. 137, 055101 (2012)], based on molecular dynamics for the lipids and multi-particle collision dynamics for the solvent, with an elastic-network description of active proteins. The method is implemented for a model molecular machine which performs active conformational motions induced by ligand binding and its release after reaction. The situation characteristic for peripheral membrane proteins is considered. Statistical investigations of the effects of single active or passive inclusions on the shape of the membrane are carried out. The results show that the peripheral machine produces asymmetric perturbations of the thickness of two leaflets of the membrane. It also produces a local saddle in the midplane height of the bilayer. Analysis of the power spectrum of the fluctuations of the membrane midplane shows that the conformational motion of the machine perturbs these membrane fluctuations. The hydrodynamic lipid flows induced by cyclic conformational changes in the machine are analyzed. It is shown that such flows are long-ranged and should provide an additional important mechanism for interactions between active inclusions in biological membranes.