Reaction Dynamics of ATP Hydrolysis in Actin Determined by ab Initio Molecular Dynamics Simulations

Mechanism of phosphate release from actin filaments

After ATP-actin monomers assemble filaments, the ATP's [Formula: see text]-phosphate is hydrolyzedwithin seconds and dissociates over minutes. We used all-atom molecular dynamics simulations to sample the release of phosphate from filaments and study residues that gate release. Dissociation of phosphate from Mg2+ is rate limiting and associated with an energy barrier of 20 kcal/mol, consistent with experimental rates of phosphate release. Phosphate then diffuses within an internal cavity toward a gate formed by R177, as suggested in prior computational studies and cryo-EM structures. The gate is closed when R177 hydrogen bonds with N111 and is open when R177 forms a salt bridge with D179. Most of the time, interactions of R177 with other residues occlude the phosphate release pathway. Machine learning analysis reveals that the occluding interactions fluctuate rapidly, underscoring the secondary role of backdoor gate opening in Pi release, in contrast with the previous hypothesis that gate opening is the primary event.

Study of the Myosin Relay Helix Peptide by Molecular Dynamics Simulations, Pump-Probe and 2D Infrared Spectroscopy

The Davydov model was conjectured to describe how an amide I excitation created during ATP hydrolysis in myosin might be significant in providing energy to drive myosin's chemomechanical cycle. The free energy surfaces of the myosin relay helix peptide dissolved in 2,2,2-trifluoroethanol (TFE), determined by metadynamics simulations, demonstrate local minima differing in free energy by only ~2 kT, corresponding to broken and stabilized hydrogen bonds, respectively. Experimental pump-probe and 2D infrared spectroscopy were performed on the peptide dissolved in TFE. The relative heights of two peaks seen in the pump-probe data and the corresponding relative volumes of diagonal peaks seen in the 2D-IR spectra at time delays between 0.5 ps and 1 ps differ noticeably from what is seen at earlier or later time delays or in the linear spectrum, indicating that a vibrational excitation may influence the conformational state of this helix. Thus, it is possible that the presence of an amide I excitation may be a direct factor in the conformational state taken on by the myosin relay helix following ATP hydrolysis in myosin.

Mechanism of Phosphate Release from Actin Filaments

After ATP-actin monomers assemble filaments, the ATP's γ-phosphate is hydrolyzed within seconds and dissociates over minutes. We used all-atom molecular dynamics simulations to sample the release of phosphate from filaments and study residues that gate release. Dissociation of phosphate from Mg2+ is rate limiting and associated with an energy barrier of 20 kcal/mol, consistent with experimental rates of phosphate release. Phosphate then diffuses in an internal cavity toward a gate formed by R177 suggested in prior computational studies and cryo-EM structures. The gate is closed when R177 hydrogen bonds with N111 and is open when R177 forms a salt bridge with D179. Most of the time interactions of R177 with other residues occludes the phosphate release pathway. Machine learning analysis reveals that the occluding interactions fluctuate rapidly, underscoring the secondary role of backdoor gate opening in Pi release, in contrast with the previous hypothesis that gate opening is the primary event.

Structural Significance of His73 in F-Actin Dynamics: Insights from Ab Initio Study

F-actin dynamics (polymerization and depolymerization) are associated with nucleotide exchange, providing the driving forces for dynamic cellular activities. As an important residue in the nucleotide state-sensing region in actin, His73 is often found to be methylated in natural actin and directly participates in F-actin dynamics by regulating nucleotide exchange. The interaction between His73 and its neighboring residue, Gly158, has significance for F-actin dynamics. However, this weak chemical interaction is difficult to characterize using classic molecular modeling methods. In this study, ab initio modeling was employed to explore the binding energy between His73 and Gly158. The results confirm that the methyl group on the His73 side chain contributes to the structural stability of atomistic networks in the nucleotide state-sensing region of actin monomers and confines the material exchange (Pi release) pathway within F-actin dynamics. Further binding energy analyses of actin structures under different nucleotide states showed that the potential model of His73/Gly158 hydrogen bond breaking in the material exchange mechanism is not obligatory within F-actin dynamics.

Towards a structural understanding of the remodeling of the actin cytoskeleton

Actin filaments (F-actin) are a key component of eukaryotic cells. Whether serving as a scaffold for myosin or using their polymerization to push onto cellular components, their function is always related to force generation. To control and fine-tune force production, cells have a large array of actin-binding proteins (ABPs) dedicated to control every aspect of actin polymerization, filament localization, and their overall mechanical properties. Although great advances have been made in our biochemical understanding of the remodeling of the actin cytoskeleton, the structural basis of this process is still being deciphered. In this review, we summarize our current understanding of this process. We outline how ABPs control the nucleation and disassembly, and how these processes are affected by the nucleotide state of the filaments. In addition, we highlight recent advances in the understanding of actomyosin force generation, and describe recent advances brought forward by the developments of electron cryomicroscopy.

Multiscale simulation of actin filaments and actin-associated proteins

Actin is an important cytoskeletal protein that serves as a building block to form filament networks that span across the cell. These networks are orchestrated by a myriad of other cytoskeletal entities including the unbranched filament-forming protein formin and branched network-forming protein complex Arp2/3. Computational models have been able to provide insights into many important structural transitions that are involved in forming these networks, and into the nature of interactions essential for actin filament formation and for regulating the behavior of actin-associated proteins. In this review, we summarize a subset of such models that focus on the atomistic features and those that can integrate atomistic features into a larger picture in a multiscale fashion.

Structural evidence for the roles of divalent cations in actin polymerization and activation of ATP hydrolysis

Actin polymerization is a divalent cation-dependent process. Here we identify a cation binding site on the surface of actin in a 2.0-Å resolution X-ray structure of actin and find evidence of three additional sites in published high-resolution structures. These cations are stable in molecular dynamics (MD) simulations of the filament, suggesting a functional role in polymerization or filament rigidity. Polymerization activates the ATPase activity of the incorporating actin protomers. Careful analysis of water molecules that approach the ATP in the MD simulations revealed Gln137-activated water to be in a suitable position in F-actin, to initiate attack for ATP hydrolysis, and its occupancy was dependent on bound cations.

The structure of the actin filament is known at a resolution that has allowed the architecture of protein components to be unambiguously assigned. However, fully understanding the chemistry of the system requires higher resolution to identify the ions and water molecules involved in polymerization and ATP hydrolysis. Here, we find experimental evidence for the association of cations with the surfaces of G-actin in a 2.0-Å resolution X-ray structure of actin bound to a Cordon-Bleu WH2 motif and in previously determined high-resolution X-ray structures. Three of four reoccurring divalent cation sites were stable during molecular dynamics (MD) simulations of the filament, suggesting that these sites may play a functional role in stabilizing the filament. We modeled the water coordination at the ATP-bound Mg²⁺, which also proved to be stable during the MD simulations. Using this model of the filament with a hydrated ATP-bound Mg²⁺, we compared the cumulative probability of an activated hydrolytic water molecule approaching the γ-phosphorous of ATP, in comparison with G-actin, in the MD simulations. The cumulative probability increased in F-actin in line with the activation of actin’s ATPase activity on polymerization. However, inclusion of the cations in the filament lowered cumulative probability, suggesting the rate of hydrolysis may be linked to filament flexibility. Together, these data extend the possible roles of Mg²⁺ in polymerization and the mechanism of polymerization-induced activation of actin’s ATPase activity.

Allostery in the dengue virus NS3 helicase: Insights into the NTPase cycle from molecular simulations

The C-terminus domain of non-structural 3 (NS3) protein of the Flaviviridae viruses (e.g. HCV, dengue, West Nile, Zika) is a nucleotide triphosphatase (NTPase) -dependent superfamily 2 (SF2) helicase that unwinds double-stranded RNA while translocating along the nucleic polymer. Due to these functions, NS3 is an important target for antiviral development yet the biophysics of this enzyme are poorly understood. Microsecond-long molecular dynamic simulations of the dengue NS3 helicase domain are reported from which allosteric effects of RNA and NTPase substrates are observed. The presence of a bound single-stranded RNA catalytically enhances the phosphate hydrolysis reaction by affecting the dynamics and positioning of waters within the hydrolysis active site. Coupled with results from the simulations, electronic structure calculations of the reaction are used to quantify this enhancement to be a 150-fold increase, in qualitative agreement with the experimental enhancement factor of 10–100. Additionally, protein-RNA interactions exhibit NTPase substrate-induced allostery, where the presence of a nucleotide (e.g. ATP or ADP) structurally perturbs residues in direct contact with the phosphodiester backbone of the RNA. Residue-residue network analyses highlight pathways of short ranged interactions that connect the two active sites. These analyses identify motif V as a highly connected region of protein structure through which energy released from either active site is hypothesized to move, thereby inducing the observed allosteric effects. These results lay the foundation for the design of novel allosteric inhibitors of NS3.

Non-structural protein 3 (NS3) is a Flaviviridae (e.g. Hepatitis C, dengue, and Zika viruses) helicase that unwinds double stranded RNA while translocating along the nucleic polymer during viral genome replication. As a member of superfamily 2 (SF2) helicases, NS3 utilizes the free energy of nucleotide triphosphate (NTP) binding, hydrolysis, and product unbinding to perform its functions. While much is known about SF2 helicases, the pathways and mechanisms through which free energy is transduced between the NTP hydrolysis active site and RNA binding cleft remains elusive. Here we present a multiscale computational study to characterize the allosteric effects induced by the RNA and NTPase substrates (ATP, ADP, and P_i) as well as the pathways of short-range, residue-residue interactions that connect the two active sites. Results from this body of molecular dynamics simulations and electronic structure calculations are highlighted in context to the NTPase enzymatic cycle, allowing for development of testable hypotheses for validation of these simulations. Our insights, therefore, provide novel details about the biophysics of NS3 and guide the next generation of experimental studies.

Theoretical studies of the function switch and mechanism of AceK as a highly active ATPase

As a multi-function enzyme, AceK integrates kinase, phosphatase and ATPase activities in a single active site and these functions are delicately regulated..

Computational characterization of the chemical step in the GTP hydrolysis by Ras-GAP for the wild-type and G13V mutated Ras

The free energy profiles for the chemical reaction of the guanosine triphosphate hydrolysis GTP + H2O → GDP + Pi by Ras-GAP for the wild-type and G13V mutated Ras were computed by using molecular dynamics protocols with the QM(ab initio)/MM potentials. The results are consistent with the recent measurements of reaction kinetics in Ras-GAP showing about two-order reduction of the rate constant upon G13V mutation in Ras: the computed activation barrier on the free energy profile is increased by 3 kcal/mol upon the G13V replacement. The major reason for a higher energy barrier is a shift of the "arginine finger" (R789 from GAP) from the favorable position in the active site. The results of simulations provide support for the mechanism of the reference reaction according to which the Q61 side chain directly participates in chemical transformations at the proton transfer stage.

Unraveling the mystery of ATP hydrolysis in actin filaments

Actin performs its myriad cellular functions by the growth and disassembly of its filamentous form. The hydrolysis of ATP in the actin filament has been shown to modulate properties of the filament, thus making it a pivotal regulator of the actin life cycle. Actin has evolved to selectively hydrolyze ATP in the filamentous form, F-actin, with an experimentally observed rate increase over the monomeric form, G-actin, of 4.3 × 10⁴. The cause of this dramatic increase in rate is investigated in this paper using extensive QM/MM simulations of both G- and F-actin. To compute the free energy of hydrolysis in both systems, metadynamics is employed along two collective variables chosen to describe the reaction coordinates of hydrolysis. F-actin is modeled as a monomer with restraints applied to coarse-grained variables enforced to keep it in a filament-like conformation. The simulations reveal a barrier height reduction for ATP hydrolysis in F-actin as compared to G-actin of 8 ± 1 kcal/mol, in good agreement with the experimentally measured barrier height reduction of 7 ± 1 kcal/mol. The barrier height reduction is influenced by an enhanced rotational diffusion of water in F-actin as compared to G-actin and shorter water wires between Asp154 and the nucleophilic water in F-actin, leading to more rapid proton transport.

Mechanism of deep-sea fish α-actin pressure tolerance investigated by molecular dynamics simulations

The pressure tolerance of monomeric α-actin proteins from the deep-sea fish Coryphaenoides armatus and C. yaquinae was compared to that of non-deep-sea fish C. acrolepis, carp, and rabbit/human/chicken actins using molecular dynamics simulations at 0.1 and 60 MPa. The amino acid sequences of actins are highly conserved across a variety of species. The actins from C. armatus and C. yaquinae have the specific substitutions Q137K/V54A and Q137K/L67P, respectively, relative to C. acrolepis, and are pressure tolerant to depths of at least 6000 m. At high pressure, we observed significant changes in the salt bridge patterns in deep-sea fish actins, and these changes are expected to stabilize ATP binding and subdomain arrangement. Salt bridges between ATP and K137, formed in deep-sea fish actins, are expected to stabilize ATP binding even at high pressure. At high pressure, deep-sea fish actins also formed a greater total number of salt bridges than non-deep-sea fish actins owing to the formation of inter-helix/strand and inter-subdomain salt bridges. Free energy analysis suggests that deep-sea fish actins are stabilized to a greater degree by the conformational energy decrease associated with pressure effect.