Nucleotide-dependence of G-actin conformation from multiple molecular dynamics simulations and observation of a putatively polymerization-competent superclosed state

Actin R256 Mono-methylation Is a Conserved Post-translational Modification Involved in Transcription

Nuclear actin has been elusive due to the lack of knowledge about molecular mechanisms. From actin-containing chromatin remodeling complexes, we discovered an arginine mono-methylation mark on an evolutionarily conserved R256 residue of actin (R256me1). Actin R256 mutations in yeast affect nuclear functions and cause diseases in human. Interestingly, we show that an antibody specific for actin R256me1 preferentially stains nuclear actin over cytoplasmic actin in yeast, mouse, and human cells. We also show that actin R256me1 is regulated by protein arginine methyl transferase-5 (PRMT5) in HEK293 cells. A genome-wide survey of actin R256me1 mark provides a landscape for nuclear actin correlated with transcription. Further, gene expression and protein interaction studies uncover extensive correlations between actin R256me1 and active transcription. The discovery of actin R256me1 mark suggests a fundamental mechanism to distinguish nuclear actin from cytoplasmic actin through post-translational modification (PTM) and potentially implicates an actin PTM mark in transcription and human diseases.

Nuclear actin and actin PTMs are poorly understood. Kumar et al. discover a system of actin PTMs similar to histone PTMs, including a conserved mark on nuclear actin (R256me1) with potential implications for transcription and human diseases.

Charge-dependent interactions of monomeric and filamentous actin with lipid bilayers

The cytoskeletal protein actin polymerizes into filaments that are essential for the mechanical stability of mammalian cells. In vitro experiments showed that direct interactions between actin filaments and lipid bilayers are possible and that the net charge of the bilayer as well as the presence of divalent ions in the buffer play an important role. In vivo, colocalization of actin filaments and divalent ions are suppressed, and cells rely on linker proteins to connect the plasma membrane to the actin network. Little is known, however, about why this is the case and what microscopic interactions are important. A deeper understanding is highly beneficial, first, to obtain understanding in the biological design of cells and, second, as a possible basis for the building of artificial cortices for the stabilization of synthetic cells. Here, we report the results of coarse-grained molecular dynamics simulations of monomeric and filamentous actin in the vicinity of differently charged lipid bilayers. We observe that charges on the lipid head groups strongly determine the ability of actin to adsorb to the bilayer. The inclusion of divalent ions leads to a reversal of the binding affinity. Our in silico results are validated experimentally by reconstitution assays with actin on lipid bilayer membranes and provide a molecular-level understanding of the actin-membrane interaction.

Mechanical differences between ATP and ADP actin states: A molecular dynamics study

This paper aims to give a comprehensive atomistic modeling of the nanomechanical behavior of actin monomer. Actin is a ubiquitous and essential component of cytoskeleton which forms many different cellular structures. Despite for several years great effort has been devoted to the investigation of mechanical properties of the actin filament, studies on the nanomechanical behavior of actin monomer are still lacking. These scales are, however, important for a complete understanding of the role of actin as an important component in the cytoskeleton structure. Based on the accuracy of atomistic modeling methods such as molecular dynamics simulations, steered molecular dynamics method is performed to assess tension of monomeric G-actin molecule under different types of mechanical loading including axial and lateral. As a result, stress-strain curves are obtained in aqueous solution, with either ATP or ADP bound in the nucleotide binding pocket. The obtained results yield evaluation of the tensile stiffness of a single actin monomer in lateral and normal direction. In order to compare the behavior of ATP and ADP G-actins, the number of hydrogen bonds and nonbonded interactions between the nucleotide and the protein are analyzed. Moreover, The effect of virtual spring of steered molecular dynamics on the mechanical behavior of actin monomer is investigated. The results reveal increasing the virtual spring constant leads to convergence of the stiffness. Moreover, in this paper, a generalized model is proposed to extend the obtained results for the monomeric G-actin scale to the actin filament. Our modeling estimated a persistence length of actin filament 15.41 µm, close to experimental measurements. Moreover, In this paper, the breaking force actin-actin bond is evaluated using steered molecular dynamics simulation. By applying a tensile force, actin-actin bond ruptured at 4197.5 pN.

ATP and ADP actin states

This minireview is dedicated to the memory of Henryk Eisenberg and honors his major contributions to many areas of biophysics and to the analysis of macromolecular states and interactions in particular. This work reviews the ATP and ADP states of a ubiquitous protein, actins, and considers the present evidence for and against unique, nucleotide-dependent conformations of this protein. The effects of ATP and ADP on specific structural elements of actins, its loops and clefts, as revealed by mutational, crosslinking, spectroscopic, and EPR methods are discussed. It is concluded that the existing evidence points to dynamic equilibria of these structural elements among various conformational states in both ATP- and ADP-actins, with the nucleotides impacting the equilibria distributions.

CHARMM Force-Fields with Modified Polyphosphate Parameters Allow Stable Simulation of the ATP-Bound Structure of Ca(2+)-ATPase

Adenosine triphosphate (ATP) is an indispensable energy source in cells. In a wide variety of biological phenomena like glycolysis, muscle contraction/relaxation, and active ion transport, chemical energy released from ATP hydrolysis is converted to mechanical forces to bring about large-scale conformational changes in proteins. Investigation of structure-function relationships in these proteins by molecular dynamics (MD) simulations requires modeling of ATP in solution and ATP bound to proteins with accurate force-field parameters. In this study, we derived new force-field parameters for the triphosphate moiety of ATP based on the high-precision quantum calculations of methyl triphosphate. We tested our new parameters on membrane-embedded sarcoplasmic reticulum Ca(2+)-ATPase and four soluble proteins. The ATP-bound structure of Ca(2+)-ATPase remains stable during MD simulations, contrary to the outcome in shorter simulations using original parameters. Similar results were obtained with the four ATP-bound soluble proteins. The new force-field parameters were also tested by investigating the range of conformations sampled during replica-exchange MD simulations of ATP in explicit water. Modified parameters allowed a much wider range of conformational sampling compared with the bias toward extended forms with original parameters. A diverse range of structures agrees with the broad distribution of ATP conformations in proteins deposited in the Protein Data Bank. These simulations suggest that the modified parameters will be useful in studies of ATP in solution and of the many ATP-utilizing proteins.

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.

Nucleotide regulation of the structure and dynamics of G-actin

Actin, a highly conserved cytoskeletal protein found in all eukaryotic cells, facilitates cell motility and membrane remodeling via a directional polymerization cycle referred to as treadmilling. The nucleotide bound at the core of each actin subunit regulates this process. Although the biochemical kinetics of treadmilling has been well characterized, the atomistic details of how the nucleotide affects polymerization remain to be definitively determined. There is increasing evidence that the nucleotide regulation (and other characteristics) of actin cannot be fully described from the minimum energy structure, but rather depends on a dynamic equilibrium between conformations. In this work we explore the conformational mobility of the actin monomer (G-actin) in a coarse-grained subspace using umbrella sampling to bias all-atom molecular-dynamics simulations along the variables of interest. The results reveal that ADP-bound actin subunits are more conformationally mobile than ATP-bound subunits. We used a multiscale analysis method involving coarse-grained and atomistic representations of these simulations to characterize how the nucleotide affects the low-energy states of these systems. The interface between subdomains SD2-SD4, which is important for polymerization, is stabilized in an actin filament-like (F-actin) conformation in ATP-bound G-actin. Additionally, the nucleotide modulates the conformation of the SD1-SD3 interface, a region involved in the binding of several actin-binding proteins.

Coarse-graining provides insights on the essential nature of heterogeneity in actin filaments

Experiments have shown that actin is structurally polymorphic, but knowledge of the details of molecular level heterogeneity in both the dynamics of a single subunit and the interactions between subunits is still lacking. Here, using atomistic molecular dynamics simulations of the actin filament, we identify domains of atoms that move in a correlated fashion, quantify interactions between these domains using coarse-grained (CG) analysis methods, and perform CG simulations to explore the importance of filament heterogeneity. The persistence length and torsional stiffness calculated from molecular dynamics simulation data agree with experimental values. We additionally observe that distinct actin conformations coexist in actin filaments. The filaments also exhibit random twist angles that are broadly distributed. CG analysis reveals that interactions between equivalent CG pairs vary from one subunit to another. To explore the importance of heterogeneity on filament dynamics, we perform CG simulations using different methods of parameterization to show that only by including heterogeneous interactions can we reproduce the twist angles and related properties. Free energy calculations further suggest that in general the actin filament is best represented as a set of subunits with differing CG sites and interactions, and the incorporating heterogeneity into the CG interactions is more important than including that in the CG sites. Our work therefore presents a systematic method to explore molecular level detail in this large and complex biopolymer.

Effects of actin-binding proteins on the thermal stability of monomeric actin

Differential scanning calorimetry (DSC) was applied to investigate the thermal unfolding of rabbit skeletal muscle G-actin in its complexes with actin-binding proteins, cofilin, twinfilin, and profilin. The results show that the effects of these proteins on the thermal stability of G-actin depend on the nucleotide, ATP or ADP, bound in the nucleotide-binding cleft between actin subdomains 2 and 4. Interestingly, cofilin binding stabilizes both ATP-G-actin and ADP-G-actin, whereas twinfilin increases the thermal stability of the ADP-G-actin but not that of the ATP-G-actin. By contrast, profilin strongly decreases the thermal stability of the ATP-G-actin but has no appreciable effect on the ADP-G-actin. Comparison of these DSC results with literature data reveals a relationship between the effects of actin-binding proteins on the thermal unfolding of G-actin, stabilization or destabilization, and their effects on the rate of nucleotide exchange in the nucleotide-binding cleft, decrease or increase. These results suggest that the thermal stability of G-actin depends, at least partially, on the conformation of the nucleotide-binding cleft: the actin molecule is more stable when the cleft is closed, while an opening of the cleft leads to significant destabilization of G-actin. Thus, DSC studies of the thermal unfolding of G-actin can provide new valuable information about the conformational changes induced by actin-binding proteins in the actin molecule.

Structural states and dynamics of the D-loop in actin

Conformational changes induced by ATP hydrolysis on actin are involved in the regulation of complex actin networks. Previous structural and biochemical data implicate the DNase I binding loop (D-loop) of actin in such nucleotide-dependent changes. Here, we investigated the structural and conformational states of the D-loop (in solution) using cysteine scanning mutagenesis and site-directed labeling. The reactivity of D-loop cysteine mutants toward acrylodan and the mobility of spin labels on these mutants do not show patterns of an α-helical structure in monomeric and filamentous actin, irrespective of the bound nucleotide. Upon transition from monomeric to filamentous actin, acrylodan emission spectra and electron paramagnetic resonance line shapes of labeled mutants are blue-shifted and more immobilized, respectively, with the central residues (residues 43-47) showing the most drastic changes. Moreover, complex electron paramagnetic resonance line shapes of spin-labeled mutants suggest several conformational states of the D-loop. Together with a new (to our knowledge) actin crystal structure that reveals the D-loop in a unique hairpin conformation, our data suggest that the D-loop equilibrates in F-actin among different conformational states irrespective of the nucleotide state of actin.

Inhibition of interdomain motion in g-actin by the natural product latrunculin: a molecular dynamics study

As part of the cytoskeleton, actin is essential for the morphology, motility, and division of eukaryotic cells. Recent X-ray fiber diffraction studies have shown that the conformation of monomeric actin is flattened upon incorporation into the filament by a relative rotation of its two major domains. The antiproliferative activity of latrunculin, a macrolide toxin produced by sponges, seems to be related to its binding to monomeric actin and inhibition of polymerization. Yet, the mechanism of inhibition is not known in detail. Here, multiple explicit water molecular dynamics simulations show that latrunculin binding hinders the conformational transition related to actin polymerization. In particular, the presence of latrunculin at the interface of the two major domains of monomeric actin reduces the correlated displacement of Domain 2 with respect to Domain 1. Moreover, higher rotational flexibility between the two major domains is observed in the absence of ATP as compared to ATP-bound actin, offering a possible explanation as to why actin polymerizes more favorably in the absence of nucleotides.

Key structural features of the actin filament Arp2/3 complex branch junction revealed by molecular simulation

We investigated the structure, properties and dynamics of the actin filament branch junction formed by actin-related protein (Arp) 2/3 complex using all-atom molecular dynamics (MD) simulations based on a model fit to a reconstruction from electron tomograms. Simulations of the entire structure consisting of 31 protein subunits together with solvent molecules containing ∼3 million atoms were performed for an aggregate time of 175 ns. One 75-ns simulation of the original reconstruction was compared to two 50-ns simulations of alternate structures, showing that the hypothesized branch junction structure is very stable. Our simulations revealed that the interface between Arp2/3 complex and the mother actin filament features a large number of salt bridges and hydrophobic contacts, many of which are dynamic and formed/broken on the timescale of the simulation. The simulations suggest that the DNase binding loops in Arp3, and possibly Arp2, form stabilizing contacts with the mother filament. Unbiased comparison of models sampled from the MD simulation trajectory with the primary experimental electron tomography data identified regions were snapshots from the simulation provide atomic details of the model structures and also pinpoints regions where the initial modeling based on the electron tomogram reconstruction may be suboptimal.

Structural modeling and molecular dynamics simulation of the actin filament

Actin is a major structural protein of the eukaryotic cytoskeleton and enables cell motility. Here, we present a model of the actin filament (F-actin) that not only incorporates the global structure of the recently published model by Oda et al. but also conserves internal stereochemistry. A comparison is made using molecular dynamics simulation of the model with other recent F-actin models. A number of structural determents such as the protomer propeller angle, the number of hydrogen bonds, and the structural variation among the protomers are analyzed. The MD comparison is found to reflect the evolution in quality of actin models over the last 6 years. In addition, simulations of the model are carried out in states with both ADP or ATP bound and local hydrogen-bonding differences characterized.

Multiple Conformations of F-actin

Actin works within eukaryotic cells to facilitate a variety of cellular processes, which are driven by the assembly of G-actin (monomeric form) into F-actin (fibrous form), and the disassembly of F-actin into G-actin. F-actin adopts multiple conformations, which are specified by interactions with various actin-binding proteins. Knowledge of the multiple conformations of actin is the key for understanding its cellular functions. Recently, we published a refined model for F-actin. In this review, based on this model, we discuss the origin, mechanism, and possible physiological significance of the multiple conformations of F-actin.

Specific cleavage of the DNase-I binding loop dramatically decreases the thermal stability of actin

Differential scanning calorimetry was used to investigate the thermal unfolding of actin specifically cleaved within the DNaseI-binding loop between residues Met47-Gly48 or Gly42-Val43 by two bacterial proteases, subtilisin or ECP32/grimelysin (ECP), respectively. The results obtained show that both cleavages strongly decreased the thermal stability of monomeric actin with either ATP or ADP as a bound nucleotide. An even more pronounced difference in the thermal stability between the cleaved and intact actin was observed when both actins were polymerized into filaments. Similar to intact F-actin, both cleaved F-actins were significantly stabilized by phalloidin and aluminum fluoride; however, in all cases, the thermal stability of the cleaved F-actins was much lower than that of intact F-actin, and the stability of ECP-cleaved F-actin was lower than that of subtilisin-cleaved F-actin. These results confirm that the DNaseI-binding loop is involved in the stabilization of the actin structure, both in monomers and in the filament subunits, and suggest that the thermal stability of actin depends, at least partially, on the conformation of the nucleotide-binding cleft. Moreover, an additional destabilization of the unstable cleaved actin upon ATP/ADP replacement provides experimental evidence for the highly dynamic actin structure that cannot be simply open or closed, but rather should be considered as being able to adopt multiple conformations.

A nucleotide state-sensing region on actin

The nucleotide state of actin (ATP, ADP-P(i), or ADP) is known to impact its interactions with other actin molecules upon polymerization as well as with multiple actin binding proteins both in the monomeric and filamentous states of actin. Recently, molecular dynamics simulations predicted that a sequence located at the interface of subdomains 1 and 3 (W-loop; residues 165-172) changes from an unstructured loop to a beta-turn conformation upon ATP hydrolysis (Zheng, X., Diraviyam, K., and Sept, D. (2007) Biophys. J. 93, 1277-1283). This region participates directly in the binding to other subunits in F-actin as well as to cofilin, profilin, and WH2 domain proteins and, therefore, could contribute to the nucleotide sensitivity of these interactions. The present study demonstrates a reciprocal communication between the W-loop region and the nucleotide binding cleft on actin. Point mutagenesis of residues 167, 169, and 170 and their site-specific labeling significantly affect the nucleotide release from the cleft region, whereas the ATP/ADP switch alters the fluorescence of probes located in the W-loop. In the ADP-P(i) state, the W-loop adopts a conformation similar to that in the ATP state but different from the ADP state. Binding of latrunculin A to the nucleotide cleft favors the ATP-like conformation of the W-loop, whereas ADP-ribosylation of Arg-177 forces the W-loop into a conformation distinct from those in the ADP and ATP-states. Overall, our experimental data suggest that the W-loop of actin is a nucleotide sensor, which may contribute to the nucleotide state-dependent changes in F-actin and nucleotide state-modulated interactions of both G- and F-actin with actin-binding proteins.

Optimal identification of semi-rigid domains in macromolecules from molecular dynamics simulation

Biological function relies on the fact that biomolecules can switch between different conformations and aggregation states. Such transitions involve a rearrangement of parts of the biomolecules involved that act as dynamic domains. The reliable identification of such domains is thus a key problem in biophysics. In this work we present a method to identify semi-rigid domains based on dynamical data that can be obtained from molecular dynamics simulations or experiments. To this end the average inter-atomic distance-deviations are computed. The resulting matrix is then clustered by a constrained quadratic optimization problem. The reliability and performance of the method are demonstrated for two artificial peptides. Furthermore we correlate the mechanical properties with biological malfunction in three variants of amyloidogenic transthyretin protein, where the method reveals that a pathological mutation destabilizes the natural dimer structure of the protein. Finally the method is used to identify functional domains of the GroEL-GroES chaperone, thus illustrating the efficiency of the method for large biomolecular machines.