Gating mechanisms during actin filament elongation by formins

  1. Fikret Aydin
  2. Naomi Courtemanche
  3. Thomas D Pollard
  4. Gregory A Voth  Is a corresponding author
  1. The University of Chicago, United States
  2. University of Minnesota, United States
  3. Yale University, United States
9 figures, 2 tables and 5 additional files

Figures

Processive association of formin with the barbed end of a growing filament.

(A) End-on and side views of formin FH1FH2 domains interacting with the barbed end of an actin filament (gray) in closed (FH2 domain in red color) and open (FH2 domain in green color) states. The black curved lines represent FH1 domains. (B) The addition of an actin subunit (orange) from the bulk solution to the filament end, and stepping of an FH2 domain (green) onto the newly added actin subunit.

https://doi.org/10.7554/eLife.37342.002
Homology models of FH2 domains on the barbed end of actin filaments.

(A) Overall quality of the homology models of Cdc12, Bni1 and mDia1 FH2 domains after 200 ns of MD simulations. The graph compares the z-scores of the three homology models with the z-scores of all experimentally determined native proteins from Protein Data Bank. (B, C) Ribbon diagrams of a dimer of Cdc12 FH2 domains interacting with the barbed end of a space-filling model of an actin filament seven-mer in the state before the formin steps onto the newly added actin subunit A1 on the barbed end. The actin subunits are numbered from A1 to A7, starting from the barbed end. (B) Side view. (C) View from the barbed end without subunit A1 and with labels on the regions of FH2 domains. (D) The root mean square deviations (RMSD) of C-alpha atoms over time during the all-atom MD simulations of dimers of FH2 domains of Cdc12, Bni1 and mDia1 on an actin seven-mer. The first 160 ns of the trajectory for Bni1 is from Baker et al. (2015). The trajectories of the whole complex, FH2 domains and linkers are displayed in different colors. (E, F) Structural alignment of ribbon diagrams taken at the end of 200 ns all-atom simulations of the FH2 domains of Bni1, Cdc12 and mDia1 associated with the actin filament seven-mer. Structural features are labeled. (E) Views from the barbed end with actin subunits A2 and A3 in the upper panel. The FHL and FHT domains are shown separately by aligning the coiled-coil regions in the lower panel. Superimposed features are shown in gray and features that differ between the three formins are color-coded. The actin filaments are not aligned and only shown to guide the location of formins. (F) The superimposed ribbon diagrams of the separate parts of the FH2 domains of the three formins.

https://doi.org/10.7554/eLife.37342.003
Figure 2—source data 1

Cdc12 FH2 interacting with a seven-mer filament.

PDB file of the structure of Cdc12 FH2 domain interacting with an actin filament composed of seven actin subunits, obtained at 500 ns of the all-atom MD simulation.

https://doi.org/10.7554/eLife.37342.004
Figure 2—source data 2

Bni1 FH2 interacting with a seven-mer filament.

PDB file of the structure of Bni1 FH2 domain interacting with an actin filament composed of seven actin subunits, obtained at 500 ns of the all-atom MD simulation.

https://doi.org/10.7554/eLife.37342.005
Figure 2—source data 3

mDia1 FH2 interacting with a seven-mer filament.

PDB file of the structure of mDia1 FH2 domain interacting with an actin filament composed of seven actin subunits, obtained at 500 ns of the all-atom MD simulation.

https://doi.org/10.7554/eLife.37342.006
Figure 2—source data 4

Cdc12 FH2 interacting with a five-mer filament.

PDB file of the structure of Cdc12 FH2 domain interacting with an actin filament composed of five actin subunits, obtained at 350 ns of the all-atom MD simulation.

https://doi.org/10.7554/eLife.37342.007
Figure 2—source data 5

Bni1 FH2 interacting with a five-mer filament.

PDB file of the structure of Bni1 FH2 domain interacting with an actin filament composed of five actin subunits, obtained at 350 ns of the all-atom MD simulation.

https://doi.org/10.7554/eLife.37342.008
Figure 2—source data 6

mDia1 FH2 interacting with a five-mer filament.

PDB file of the structure of mDia1 FH2 domain interacting with an actin filament composed of five actin subunits, obtained at 350 ns of the all-atom MD simulation.

https://doi.org/10.7554/eLife.37342.009
Steric clashes of FH2 dimers with incoming actin subunits during all-atom simulations of seven-mer and five-mer filaments.

(A–C) FH2 dimers on seven-mer filaments including actin A1. (A) Time course of volume fractions of actin subunit (+A1) occupied by FHL during 500 ns simulations. Measurements start at 0 ns for Cdc12 and mDia1, and at 160 ns for Bni1. (B, C) Ribbon diagrams of Cdc12, Bni1 and mDia1 FH2 domains interacting with the barbed end of a space-filling model of an actin filament at the beginning (left) and at the 200 ns (middle) MD simulations. In the models on the right, red actin subunit +A1 was added at the end of the simulations. (B) Side views. (C) Views from the barbed end. (D, E) FH2 dimers on five-mer filaments without actin subunits A1 and A7. Time courses of volume fractions of incoming actin subunit (A1) occupied by FHT during (D) 350 ns simulations for the first replicas and (E) 200 ns simulations for the second replicas. The initial structures were generated by removing subunits A1 and A7 at the end of 200 ns all-atom simulations of seven-mer filaments from (A–C). FH2 dimers on a five-mer filament without actins A1 and A7 is shown in the inset.

https://doi.org/10.7554/eLife.37342.010
Figure 4 with 1 supplement
Effect of FH2 domains on the barbed end configurations of actin seven-mer filaments.

The panels show helical twist angles between subunits at the barbed ends of actin filaments associated with three different formin FH2 dimers during all-atom MD simulations. (A–C) The angles between the actin subunits (A1–A2 and A2–A3) as a function of time. The measurements start after initial equilibrations (Figure 2) at 100 ns for Cdc12 and mDia1 or at 160 ns for Bni1. (D–G) Comparison of the distributions of angles between actin subunits A1-A2 and A2-A3 during two different time intervals. (D, F) t = 100–150 ns for Cdc12 and mDia1, and t = 160–210 ns for Bni1. (E, G) t = 450–500 ns for all systems.

https://doi.org/10.7554/eLife.37342.011
Figure 4—source data 1

Twist angles between actin subunits as a function of time.

Helical twist angles between subunits of seven-mer filaments (A1–A2, A2–A3, A3–A4, A4–A5 and A5–A6) associated with three different formin FH2 dimers during all-atom MD simulations. The angles are outputted at every 40 ps during the time intervals given in Figure 4.

https://doi.org/10.7554/eLife.37342.013
Figure 4—figure supplement 1
Effect of FH2 domains on the configurations of actin seven-mer filaments.

Comparison of the distributions of angles between the actin pair of subunits during three different time intervals. (A,D,G) t = 100–120 ns for Cdc12 and mDia1, and t = 160–180 ns for Bni1. (B,E,H) t = 280–300 ns for all systems. (C,F,I) t = 480–500 ns for all systems.

https://doi.org/10.7554/eLife.37342.012
Figure 5 with 2 supplements
Effect of FH2 domains on the barbed end configurations of actin five-mer filaments before the addition of actin subunit A1.

(A) The helical twist angles between actin subunits A2 and A3 as a function of time during 350 ns of all-atom MD simulations (replica 1) of the FH2 domains of Bni1, Cdc12 and mDia1 associated with the actin filament five-mer consisting of subunits A2, A3, A4, A5 and A6. The initial structures in these simulations were generated by removing subunits A1 and A7 at the end of 200 ns all-atom simulations of the FH2 domains of Bni1, Cdc12 and mDia1 associated with the actin filament seven-mer. (B, C) Comparison of the distributions of angles between actin subunits A2 and A3 during the last time intervals of two independent simulations (replica 1 and replica 2). (B) The last 50 ns of the 350 ns simulations. (C) The last 20 ns of the 200 ns simulations. (D–F) The number of contacts between the lasso, knob and post regions of the FH2 domains and actin subunits A2 and A3 as a function of time during 350 ns of the all-atom MD simulations (replica 1).

https://doi.org/10.7554/eLife.37342.014
Figure 5—source data 1

Twist angles between actin subunits as a function of time.

Helical twist angles between subunits of five-mer filaments (A2–A3, A3–A4 and A4–A5) associated with three different formin FH2 dimers during all-atom MD of two independent simulations (replica 1 and replica 2). The angles are outputted at every 20 ps during the time intervals given in Figure 5.

https://doi.org/10.7554/eLife.37342.017
Figure 5—figure supplement 1
The interactions of FH2 domains (post regions) with an actin filament and the barbed end configuration.

The correlations for the number of contacts between the post regions of the FH2 domains and actin subunits (A2 and A3) and the distributions of twist angles of A2-A3 during 350 ns all-atom MD simulations of five-mer filaments.

https://doi.org/10.7554/eLife.37342.015
Figure 5—figure supplement 2
The interactions of FH2 domains (lasso and knob regions) with an actin filament and the barbed end configuration.

The correlations for the number of contacts between the lasso and knob regions of the FH2 domains and actin subunits (A2 and A3) and the distributions of twist angles of A2-A3 during 350 ns all-atom MD simulations of five-mer filaments.

https://doi.org/10.7554/eLife.37342.016
Free-energy sampling of the conformational mobility of FHT domains interacting with five-mer filaments.

(A) Two independent metabasin metadynamics (MBMetaD) simulations for each formin were carried out (run for 80 ns) to understand the conformational mobility of the FHT domain. The collective variables (CVs) defined in (B) were selected to describe the mobility of the FHT domain in the region of the incoming actin subunit A1 (shown by dashed black lines). The first CV (CV1) is the distance between the center of mass (COM) of FHT knob and the COM of actin subunit A2. CV2 is the distance between the COM of the FHT post and the COM of actin subunit A2. Initial distances from the knob and post to the COM of A2 are marked with ‘X’. If a larger area was favored to be explored, that would mean that the incoming actin monomer could be more easily accommodated into the barbed end as the rearrangement of the FHT domain requires a lower energy barrier. (B) Space-filling models of the barbed end of the filament (subunit A2 is yellow and subunit A3 is brown) and ribbon diagrams of the FH2 domains of Cdc12 (green with the knob and post of FHT in purple). (C) The distributions of angles between actin subunits A2 and A3 during two independent MBMetaD simulations.

https://doi.org/10.7554/eLife.37342.018
Figure 7 with 1 supplement
Interactions of FH2 domains with actin filament seven-mers.

The time-averaged number of contacts over the last 20 ns of the 200 ns simulations between the three actin subunits at the barbed ends of the filaments and the FH2 domains of mDia1 (red), Bni1 (blue) and Cdc12 (green). The actin subunits are numbered from 1 to 3, starting with the newly added A1 subunit. (A, B) A1 actin subunits; (C, D) A2 actin subunits; and (E, F) A3 actin subunits. (A, C, E) show the contacts of the FHL domains and (B, D, F) show the contacts of the FHT domains. A pair of residues was considered to be in contact if the distance between their C-alpha atoms was ≤12 Å. (G) View from the pointed end of the structural alignment of FH2 domains taken from the end of 200 ns all-atom simulations of Bni1, Cdc12 and mDia1 FH2 domains interacting with an actin filament seven-mer. The orange lines point to the knA and knB helices of mDia1, Cdc12 and Bni1 FH2 domains. (H, I) Total nonbonded interaction energy (sum of van der Waals and electrostatic interactions) between (H) FHL (leading) and (I) FHT (trailing) knob helices of mDia1, Bni1, Cdc12 formins and the barbed end of actin subunits A2 and A3.

https://doi.org/10.7554/eLife.37342.019
Figure 7—figure supplement 1
Interactions of each region in FH2 domains with an actin filament. The number of contacts between the lasso, linker, knob and post regions of the FH2 domains and actin subunits (A1, A2 and A3) as a function of time during last 50 ns of 500 ns all-atom MD simulations.
https://doi.org/10.7554/eLife.37342.020
Construction and simulations of coarse-grained (CG) models of FH2 domains.

(A) Assignment of coarse-grained (CG) sites of the Cdc12 FH2 domain dimer by the EDCG method (FHL is green and FHT is orange). (B) Connection of CG sites by harmonic bonds defined by heteroENM and matched to the fluctuations observed in all-atom simulations. (C) CG model of the Cdc12 FH2 domains on a seven-mer actin filament. FH2 domains are shown in green, actin monomers are colored as follow: A1 (gray), A2 (yellow), A3 (orange), A4 (cyan), A5 (pink), A6 (lime) and A7 (purple). (D) Plot comparing root mean square fluctuations (RMSFs) of the CG sites of the Cdc12 FH2 dimer obtained from AA simulations (AA-CG mapped) with RMSFs of these sites in CG simulations using a range of cutoff distances for the harmonic springs that connect every CG site. The cutoff distance of 50 Å (blue) gives the best match with the RMSFs of AA simulations (yellow). (E, F) Time courses of the root mean squared displacements of the FH2 dimers from three independent CG simulations spanning two microseconds at low and high dielectric constant. (E) Dielectric constant is 1. (F) Dielectric constant is 3.

https://doi.org/10.7554/eLife.37342.023
Conformational changes of FH2 domains on the barbed ends of actin filament seven-mers during coarse-grained (CG) simulations.

(A) Side and (B) bottom views of structures of Cdc12, Bni1, and mDia1 FH2 domains interacting with the barbed ends of actin filaments at the end (t = 2000 ns) of the CG simulations (dielectric constant of 3, replicate 2). (C) Time courses of the center-of-mass displacements of FHL domains from their initial positions (shown as percentage change) in the direction perpendicular to the long axis of the actin filaments. A positive percentage indicates a shift towards the actin filament (favors the ‘closed’ state), whereas a negative percentage indicates a shift away from the actin filament (favors the ‘open’ state). The shifting of a representative FHL domain in two different directions is schematically shown at the bottom.

https://doi.org/10.7554/eLife.37342.024

Tables

Table 1
Stability of salt bridges between formin FH2 domains and actin subunits.

Stability was measured as the percent of the time that a salt bridge formed between knob helices of the three FH2 domains and actin barbed end grooves of actin subunits A2 or A3 during (A) the last 20 ns of AA simulations spanning 200 ns of seven-mer filaments, (B) the last 50 ns of the simulations spanning 500 ns of seven-mer filaments and (C) the last 50 ns of the simulations spanning 350 ns of five-mer filaments.

https://doi.org/10.7554/eLife.37342.021
Cdc12 knob helices/Actin barbed end groove
FH2 residueFH2 domainActin residueActin subunit(A) Percent formed(B) Percent formed(C) Percent formed
K1068FHLD25A278.276.986.5
E1090FHLR147A2-63.113.1
E1093FHLR147A299.4--
K1099FHLE167A2--41.2
K1105FHLE167A2-52.677.1
K1107FHLE167A241.5--
K1068FHTD25A395.270.782.4
K1072FHTD25A364.371.485.3
K1072FHTD24A3-15.612.8
E1093FHTR147A394.897.399.8
K1099FHTE167A3-47.218.1
Bni1 knob helices/Actin barbed end groove
E1463FHLR147A283.689.588.6
R1423FHLE167A20.4010.541-
E1463FHTR147A390.088.593.8
K1467FHTE167A397.495.7-
mDia1 knob helices/Actin barbed end groove
R851FHLD25A297.893.5-
E871FHLR147A28.4218.8-
K879FHLE167A225.33.69-
K838FHTE167A398.698.625.7
Table 2
Salt bridges between formin FH2 domains and actin subunits.

Percentage of the time that salt bridges formed between the lasso and linker regions of Cdc12, Bni1 and mDia1 FH2 domains and actin subunits during the last 20 ns of AA simulations spanning 200 ns of seven-mer filaments.

https://doi.org/10.7554/eLife.37342.022
FH2 residueFH2 domainActin residueActin subunit% Formed
Cdc12/actin
R990FHL-lassoD363A2100.
K992FHT-lassoE125A199.6
K1038FHL-linkerE99A271.5
K1041FHL-linkerE100A252.5
K1045FHL-linkerE2A272.3
K1045FHL-linkerD3A293.2
K1046FHL-linkerE4A250.3
K1038FHT-linkerD363A120.6
Bni1/actin
K1357FHL-lassoD363A285.2
K1357FHT-lassoD363A123.0
K1359FHT-lassoE125A199.2
R1402FHL-linkerE99A299.4
E1403FHL-linkerK359A295.0
K1410FHL-linkerE4A224.0
K1410FHL-linkerE100A281.8
K1412FHL-linkerD3A257.5
R1402FHT-linkerD363A198.4
mDia1/actin
R764FHL-lassoD363A223.4
K807FHL-linkerE125A249.9
K813FHL-linkerE4A237.9
K813FHL-linkerE99A220.4
K826FHL-linkerD3A270.9
K826FHL-linkerE2A243.3
K827FHL-linkerE2A213.2
K828FHL-linkerE2A211.4
K807FHT-linkerE125A123.0
K810FHT-linkerE83A143.7
K810FHT-linkerE125A146.7
K813FHT-linkerD51A199.0
E816FHT-linkerR37A196.8
K828FHT-linkerE2A325.1

Additional files

Supplementary file 1

Definitions of features of Cdc12, Bni1 and mDia1 FH2 domains.

The columns show feature name and residues comprising corresponding feature.

https://doi.org/10.7554/eLife.37342.025
Supplementary file 2

Template modeling scores (TM-scores) of the FH2 domain homology models based on different templates.

TM-score is used to determine the similarity of protein structures. Based on statistics, TM-score between 0.0 and 0.17 indicates random structural similarity and TM-score between 0.5 and 1.00 indicates having about the same fold [43]. TM-scores were obtained through this website: https://zhanglab.ccmb.med.umich.edu/TM-align/

https://doi.org/10.7554/eLife.37342.026
Supplementary file 3

Correlations between interactions of FH2 domain with an actin filament and the barbed end configuration.

Pearson correlation coefficients for the number of contacts between the lasso, knob and post regions of the FH2 domains and actin subunits (A2 and A3) and the distributions of twist angles of A2-A3 (given in Figure 5—figure supplement 2 and Figure 5—figure supplement 1) during 350 ns all-atom MD simulations of five-mer filaments.

https://doi.org/10.7554/eLife.37342.027
Supplementary file 4

Average interactions of different FH2 regions with an actin filament.

The average number of contacts (with standard deviations and t-statistics) between the lasso, knob and post regions of the FH2 domains and actin subunits (A2 and A3) from 350 ns of the all-atom MD simulations of five-mer filaments.

https://doi.org/10.7554/eLife.37342.028
Transparent reporting form
https://doi.org/10.7554/eLife.37342.029

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Fikret Aydin
  2. Naomi Courtemanche
  3. Thomas D Pollard
  4. Gregory A Voth
(2018)
Gating mechanisms during actin filament elongation by formins
eLife 7:e37342.
https://doi.org/10.7554/eLife.37342