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

Molecular dynamics simulations support a multistep pathway for activation of branched actin filament nucleation by Arp2/3 complex

Actin-related protein 2/3 complex (Arp2/3 complex) catalyzes the nucleation of branched actin filaments that push against membranes in processes like cellular motility and endocytosis. During activation by WASP proteins, the complex must bind WASP and engage the side of a pre-existing (mother) filament before a branched filament is nucleated. Recent high-resolution structures of activated Arp2/3 complex revealed two major sets of activating conformational changes. How these activating conformational changes are triggered by interactions of Arp2/3 complex with actin filaments and WASP remains unclear. Here we use a recent high-resolution structure of Arp2/3 complex at a branch junction to design all-atom molecular dynamics simulations that elucidate the pathway between the active and inactive states. We ran a total of ∼4.6 microseconds of both unbiased and steered all-atom molecular dynamics simulations starting from three different binding states, including Arp2/3 complex within a branch junction, bound only to a mother filament, and alone in solution. These simulations indicate that the contacts with the mother filament are mostly insensitive to the massive rigid body motion that moves Arp2 and Arp3 into a short pitch helical (filament-like) arrangement, suggesting actin filaments alone do not stimulate the short pitch conformational change. In contrast, contacts with the mother filament stabilize subunit flattening in Arp3, an intrasubunit change that converts Arp3 from a conformation that mimics an actin monomer to one that mimics a filamentous actin subunit. Our results support a multistep activation pathway that has important implications for understanding how WASP-mediated activation allows Arp2/3 complex to assemble force-producing actin networks.

Evolutionary diversification reveals distinct somatic versus germline cytoskeletal functions of the Arp2 branched actin nucleator protein

Branched actin networks are critical in many cellular processes, including cell motility and division. Arp2, a protein within the 7-membered Arp2/3 complex, is responsible for generating branched actin. Given its essential roles, Arp2 evolves under stringent sequence conservation throughout eukaryotic evolution. We unexpectedly discovered recurrent evolutionary diversification of Arp2 in Drosophila, yielding independently arising paralogs Arp2D in obscura species and Arp2D2 in montium species. Both paralogs are unusually testis-enriched in expression relative to Arp2. We investigated whether their sequence divergence from canonical Arp2 led to functional specialization by replacing Arp2 in D. melanogaster with either Arp2D or Arp2D2. Despite their divergence, we surprisingly found both complement Arp2's essential function in the soma, suggesting they have preserved the ability to polymerize branched actin even in a non-native species. However, we found that Arp2D-expressing males are subfertile and display many defects throughout sperm development. We pinpointed two highly diverged structural regions in Arp2D that contribute to these defects: subdomain 2 and the C-terminus. We expected that germline function would be rescued by replacing Arp2D's long and charged C-terminus with Arp2's short C-terminus, yet surprisingly, the essential somatic function of Arp2D was lost. Therefore, while Arp2D's structural divergence is incompatible with D. melanogaster sperm development, its unique C-terminus has evolved a critical role in actin polymerization. Our findings suggest canonical Arp2's function differs between somatic versus germline contexts, and Arp2 paralogs have recurrently evolved and specialized for actin branching in the testis.

Steps of actin filament branch formation by Arp2/3 complex investigated with coarse-grained molecular dynamics

The nucleation of actin filament branches by the Arp2/3 complex involves activation through nucleation promotion factors (NPFs), recruitment of actin monomers, and binding of the complex to the side of actin filaments. Because of the large system size and processes that involve flexible regions and diffuse components, simulations of branch formation using all-atom molecular dynamics are challenging. We applied a coarse-grained model that retains amino-acid level information and allows molecular dynamics simulations in implicit solvent, with globular domains represented as rigid bodies and flexible regions allowed to fluctuate. We used recent electron microscopy structures of the inactive Arp2/3 complex bound to NPF domains and to mother actin filament for the activated Arp2/3 complex. We studied interactions of Arp2/3 complex with the activating VCA domain of the NPF Wiskott-Aldrich syndrome protein, actin monomers, and actin filament. We found stable configurations with one or two actin monomers bound along the branch filament direction and with CA domain of VCA associated to the strong and weak binding sites of the Arp2/3 complex, supporting prior structural studies and validating our approach. We reproduced delivery of actin monomers and CA to the Arp2/3 complex under different conditions, providing insight into mechanisms proposed in previous studies. Simulations of active Arp2/3 complex bound to a mother actin filament indicate the contribution of each subunit to the binding. Addition of the C-terminal tail of Arp2/3 complex subunit ArpC2, which is missing in the cryo-EM structure, increased binding affinity, indicating a possible stabilizing role of this tail.

Deciphering the Molecular Terms of Arp2/3 Allosteric Regulation from All-Atom Simulations and Dynamical Network Theory

The Arp2/3 molecular machine stimulates the generation of branched actin networks at the cytosolic surface of cellular membranes. Arp2/3 is thus pivotal for cell motility and migration, and its aberrant function is implicated in cancer invasion and metastasis. Here, all-atom multi μs-long molecular dynamics simulations and dynamical NetWork Analysis (NWA) unprecedentedly disclose the molecular terms of Arp2/3 regulation (activation/inhibition) by positive/negative allosteric modulators. After identifying the crucial structural elements underlying Arp2/3's conformational transition toward its active actin-polymerization-competent state, we decrypt the activating signaling paths heading from the allosteric effector (ATP) binding sites to these pivotal regions, also elucidating how small-molecule inhibitors scramble this signal-exchange. As a result, while ATP-induced signaling triggers a harmonious conformational transition toward active Arp2/3, the inhibitors disturb these information channels, desynchronizing Arp2/3 functional movements, thus hindering its activation. Our outcomes supply a conceptual basis for devising small-molecule inhibitors to block infiltrative cancer migration.

Cryo-electron tomography structure of Arp2/3 complex in cells reveals new insights into the branch junction

The actin-related protein (Arp)2/3 complex nucleates branched actin filament networks pivotal for cell migration, endocytosis and pathogen infection. Its activation is tightly regulated and involves complex structural rearrangements and actin filament binding, which are yet to be understood. Here, we report a 9.0 Å resolution structure of the actin filament Arp2/3 complex branch junction in cells using cryo-electron tomography and subtomogram averaging. This allows us to generate an accurate model of the active Arp2/3 complex in the branch junction and its interaction with actin filaments. Notably, our model reveals a previously undescribed set of interactions of the Arp2/3 complex with the mother filament, significantly different to the previous branch junction model. Our structure also indicates a central role for the ArpC3 subunit in stabilizing the active conformation.

A Burst of Genetic Innovation in Drosophila Actin-Related Proteins for Testis-Specific Function

Many cytoskeletal proteins perform fundamental biological processes and are evolutionarily ancient. For example, the superfamily of actin-related proteins (Arps) specialized early in eukaryotic evolution for diverse cellular roles in the cytoplasm and the nucleus. Despite its strict conservation across eukaryotes, we find that the Arp superfamily has undergone dramatic lineage-specific diversification in Drosophila. Our phylogenomic analyses reveal four independent Arp gene duplications that occurred in the common ancestor of the obscura group of Drosophila and have been mostly preserved in this lineage. All four obscura-specific Arp paralogs are predominantly expressed in the male germline and have evolved under positive selection. We focus our analyses on the divergent Arp2D paralog, which arose via a retroduplication event from Arp2, a component of the Arp2/3 complex that polymerizes branched actin networks. Computational modeling analyses suggest that Arp2D can replace Arp2 in the Arp2/3 complex and bind actin monomers. Together with the signature of positive selection, our findings suggest that Arp2D may augment Arp2's functions in the male germline. Indeed, we find that Arp2D is expressed during and following male meiosis, where it localizes to distinct locations such as actin cones-specialized cytoskeletal structures that separate bundled spermatids into individual mature sperm. We hypothesize that this unprecedented burst of genetic innovation in cytoskeletal proteins may have been driven by the evolution of sperm heteromorphism in the obscura group of Drosophila.

Plastic Deformation and Fragmentation of Strained Actin Filaments

The assembly of actin filaments and filament networks generate forces that drive cell and vesicle movement. These structures and the comprising actin filaments must be mechanically stable to sustain these forces and maintain their structural integrity. Filaments in these dynamic structures must also be disassembled to recycle and replenish the pool of actin monomers available for polymerization. Actin-severing proteins such as cofilin and contractile myosin motor proteins fragment these nominally stable structures. We developed a mesoscopic-length-scale actin filament model to investigate force-induced filament fragmentation. We show that fragmentation in our model occurs at curvatures similar to previous measurements of fragmentation within (cofil)actin and actin-cofilactin boundaries. Boundaries between bare and cofilin-decorated segments are brittle and fragment at small bending and twisting deformations. Extending filaments disperses strain uniformly over subunit interfaces, and filaments fragment with no detectable partial rupture or plastic deformation. In contrast, bending or twisting filaments imposes nonuniform interface strain and leads to partial interface rupture, accelerating filament fragmentation. As a result, the rupture force under compressive loads is an order of magnitude lower than under tensile loads. Partial interface rupture may be a primary mechanism of accelerating actin filament fragmentation by other actin-destabilizing proteins.

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.

Conformational changes in Arp2/3 complex induced by ATP, WASp-VCA, and actin filaments

We used fluorescence spectroscopy and EM to determine how binding of ATP, nucleation-promoting factors, actin monomers, and actin filaments changes the conformation of Arp2/3 complex during the process that nucleates an actin filament branch. We mutated subunits of Schizosaccharomyces pombe Arp2/3 complex for labeling with fluorescent dyes at either the C termini of Arp2 and Arp3 or ArpC1 and ArpC3. We measured Förster resonance energy transfer (FRET) efficiency (ET_eff) between the dyes in the presence of the various ligands. We also computed class averages from electron micrographs of negatively stained specimens. ATP binding made small conformational changes of the nucleotide-binding cleft of the Arp2 subunit. WASp-VCA, WASp-CA, and WASp-actin-VCA changed the ET_eff between the dyes on the Arp2 and Arp3 subunits much more than between dyes on ArpC1 and ArpC3. Ensemble FRET detected an additional structural change that brought ArpC1 and ArpC3 closer together when Arp2/3 complex bound actin filaments. VCA binding to Arp2/3 complex causes a conformational change that favors binding to the side of an actin filament, which allows further changes required to nucleate a daughter filament.

MCMap-A Computational Tool for Mapping Energy Landscapes of Transient Protein-Protein Interactions

MCMap is a tool particularly well-suited for analyzing energy landscapes of transient macromolecular complexes. The program applies a Monte Carlo strategy, where the ligand moves randomly in the electrostatic field of the receptor. By applying importance sampling, the major interaction sites are mapped, resulting in a global distribution of ligand-receptor complexes. This approach displays the dynamic character of transiently interacting protein complexes where not a single complex but an ensemble of complexes better describes the protein interactions. The software provides a broad range of analysis options which allow for relating the simulations to experimental data and for interpreting them on a structural level. The application of MCMap is exemplified by the electron-transfer complex of cytochrome c peroxidase and cytochrome c from baker's yeast. The functionality of MCMap and the visualization of simulation data are in particular demonstrated by studying the dependence of the association on ionic strength and on the oxidation state of the binding partner. Furthermore, microscopically, a repulsion of a second ligand can be seen in the ternary complex upon the change of the oxidation state of the bound cytochrome c. The software is made available as open source software together with the example and can be downloaded free of charge from http://www.bisb.uni-bayreuth.de/index.php?page=downloads.

Looking "Under the Hood" of Cellular Mechanotransduction with Computational Tools: A Systems Biomechanics Approach across Multiple Scales

Signal modulation has been developed in living cells throughout evolution to promote utilizing the same machinery for multiple cellular functions. Chemical and mechanical modules of signal transmission and transduction are interconnected and necessary for organ development and growth. However, due to the high complexity of the intercommunication of physical intracellular connections with biochemical pathways, there are many missing details in our overall understanding of mechanotransduction processes, i.e., the process by which mechanical signals are converted to biochemical cascades. Cell-matrix adhesions are mechanically coupled to the nucleus through the cytoskeleton. This modulated and tightly integrated network mediates the transmission of mechanochemical signals from the extracellular matrix to the nucleus. Various experimental and computational techniques have been utilized to understand the basic mechanisms of mechanotransduction, yet many aspects have remained elusive. Recently, in silico experiments have made important contributions to the field of mechanobiology. Herein, computational modeling efforts devoted to understanding integrin-mediated mechanotransduction pathways are reviewed, and an outlook is presented for future directions toward using suitable computational approaches and developing novel techniques for addressing important questions in the field of mechanotransduction.

Two Deafness-Causing Actin Mutations (DFNA20/26) Have Allosteric Effects on the Actin Structure

Point mutations in γ-cytoplasmic actin have been shown to result in autosomal-dominant, nonsyndromic, early-onset deafness. Two mutations at the same site, K118M and K118N, provide a unique opportunity to compare the effects of two dissimilar amino acid substitutions that produce a similar phenotype in humans. K118 resides in a helix that runs from K113 to T126, and mutations that alter the position, dynamics, and/or biochemistry of this helix can result in a wide range of pathologies. Using a combination of computational and experimental studies, both employing yeast actin, we find that these mutations at K118 result in changes in the structure and dynamics of the DNase-I loop, alterations in the structure of the H73 loop as well as the side-chain orientations of W79 and W86, changes in nucleotide exchange rates, and significant shifts in the twist of the actin monomer. Interestingly, in the case of K118N, the twist of the monomer is nearly identical to that of the F-actin protomer, and in vitro polymerization assays show that this mutation results in faster polymerization. Taken together, these results indicate that mutations at this site give rise to a series of small changes that can be tolerated in vivo but result in misregulation of actin assembly and dynamics.

Crystals of the Arp2/3 complex in two new space groups with structural information about actin-related protein 2 and potential WASP binding sites

Co-crystals of the bovine Arp2/3 complex with the CA motif from N-WASP in two new space groups were analyzed by X-ray diffraction. The crystals in the orthorhombic space group P212121 contained one complex per asymmetric unit, with unit-cell parameters a = 105.48, b = 156.71, c = 177.84 Å, and diffracted to 3.9 Å resolution. The crystals in the tetragonal space group P41 contained two complexes per asymmetric unit, with unit-cell parameters a = b = 149.93, c = 265.91 Å, and diffracted to 5.0 Å resolution. The electron-density maps of both new crystal forms had densities for small segments of subdomains 1 and 2 of Arp2. Both maps had density at the binding site on Arp3 for the C-terminal EWE tripeptide from N-WASP and a binding site proposed for the C motif of N-WASP in the barbed-end groove of Arp2. The map from the tetragonal crystal form had density near the barbed end of Arp3 that may correspond to the C helix of N-WASP. The noise levels and the low resolution of the maps made the assignment of specific molecular structures for any of these CA peptides impossible.

Expression of cortactin in human gliomas and its effect on migration and invasion of glioma cells

The aim of the present study was to investigate the role of cortactin in the infiltrative behavior of glioma cells and the potential mechanism of cortactin in promoting the migration and invasion of glioma cells. The expression of cortactin was detected by immunohistochemistry in 40 human glioma specimens and 8 non-tumor brain specimens. U251, LN229 and SNB19 glioma cells were employed for the in vitro study and assigned into the siRNA-cortactin (transfected with siRNA specific to cortactin), siRNA-NC (transfected with negative control RNA sequence) and siRNA-N (transfected with empty vector) groups. The expression of cortactin in different treated glioma cell groups was detected using western blot analysis and RT-qPCR. The migration and invasion of glioma cells under different treatments were assessed using a wound-healing assay and Transwell-chamber invasion assay, respectively. The lamellipodia of glioma cells following treatment were observed by immunofluorescence (IF) and changes of lamellipodia over time were imaged using an inverted microscope. The distribution of cortactin and the actin-related protein 2/3 (Arp2/3) complex in glioma cells were observed after IF detection. The expression of cortactin in the glioma specimens was significantly higher than that in non-tumor brain tissue (P<0.05) and positively correlated with the malignancy of glioma specimens (r=0.912, P=0.00). The cortactin expression in glioma cells was markedly inhibited (P<0.05) and their migration and invasion ability was also impaired significantly following treatment with siRNA (P<0.05) compared with the other two groups. The size and persistence time of lamellipodia were reduced after cortactin expression was inhibited in glioma cells. Cortactin and the Arp2/3 complex were co-localized in the front of glioma cells, where actin was polymerized and lamellipodia formed. Thus, the results revealed that, cortactin is crucial in invasion and migration of glioma cells, which may promote the migration and invasion of glioma cells by regulating lamellipodia formation, a process requiring the combination of cortactin and the Arp2/3 complex.

Molecular dynamics simulations of large macromolecular complexes

Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone.

Function and dynamics of macromolecular complexes explored by integrative structural and computational biology

Three vignettes exemplify the potential of combining EM and X-ray crystallographic data with molecular dynamics (MD) simulation to explore the architecture, dynamics and functional properties of multicomponent, macromolecular complexes. The first two describe how EM and X-ray crystallography were used to solve structures of the ribosome and the Arp2/3-actin complex, which enabled MD simulations that elucidated functional dynamics. The third describes how EM, X-ray crystallography, and microsecond MD simulations of a GPCR:G protein complex were used to explore transmembrane signaling by the β-adrenergic receptor. Recent technical advancements in EM, X-ray crystallography and computational simulation create unprecedented synergies for integrative structural biology to reveal new insights into heretofore intractable biological systems.

Actin binding proteins, spermatid transport and spermiation

The transport of germ cells across the seminiferous epithelium is composed of a series of cellular events during the epithelial cycle essential to the completion of spermatogenesis. Without the timely transport of spermatids during spermiogenesis, spermatozoa that are transformed from step 19 spermatids in the rat testis fail to reach the luminal edge of the apical compartment and enter the tubule lumen at spermiation, thereby arriving the epididymis for further maturation. Step 19 spermatids and/or sperms that remain in the epithelium beyond stage VIII of the epithelial cycle will be removed by the Sertoli cell via phagocytosis to form phagosomes and be degraded by lysosomes, leading to subfertility and/or infertility. However, the biology of spermatid transport, in particular the final events that lead to spermiation remain elusive. Based on recent data in the field, we critically evaluate the biology of spermiation herein by focusing on the actin binding proteins (ABPs) that regulate the organization of actin microfilaments at the Sertoli-spermatid interface, which is crucial for spermatid transport during this event. The hypothesis we put forth herein also highlights some specific areas of research that can be pursued by investigators in the years to come.

The joys and perils of flexible fitting

While performing their functions, biological macromolecules often form large, dynamically changing macromolecular assemblies. Only a relatively small number of such assemblies have been accessible to the atomic-resolution techniques X-ray crystallography and NMR. Electron microscopy in conjunction with image reconstruction has become the preferred alternative for revealing the structures of such macromolecular complexes. However, for most assemblies the achievable resolution is too low to allow accurate atomic modeling directly from the data. Yet, useful models often can be obtained by fitting atomic models of individual components into a low-resolution reconstruction of the entire assembly. Several algorithms for achieving optimal fits in this context were developed recently, many allowing considerable degrees of flexibility to account for binding-induced conformational changes of the assembly components. This chapter describes the advantages and potential pitfalls of these methods and puts them into perspective with alternative approaches such as iterative modular fitting of rigid-body domains.

Actin in action: imaging approaches to study cytoskeleton structure and function

The cytoskeleton plays several fundamental roles in the cell, including organizing the spatial arrangement of subcellular organelles, regulating cell dynamics and motility, providing a platform for interaction with neighboring cells, and ultimately defining overall cell shape. Fluorescence imaging has proved to be vital in furthering our understanding of the cytoskeleton, and is now a mainstay technique used widely by cell biologists. In this review we provide an introduction to various imaging modalities used to study focal adhesions and the actin cytoskeleton, and using specific examples we highlight a number of recent studies in animal cells that have advanced our knowledge of cytoskeletal behavior.

MICAL, the flavoenzyme participating in cytoskeleton dynamics

MICAL (from the Molecule Interacting with CasL) indicates a family of recently discovered cytosolic, multidomain proteins, which uniquely couple an N-terminal FAD-containing monooxygenase-like domain to typical calponine homology, LIM and coiled-coil protein-interaction modules. Genetic and cell biology approaches have demonstrated an essential role of the catalytic activity of the monooxygenase-like domain in transducing the signal initiated by semaphorins interaction with their plexin receptors, which results in local actin cytoskeleton disassembly as part of fundamental processes that include differentiation, migration and cell-cell contacts in neuronal and non-neuronal cell types. This review focuses on the structure-function relations of the MICAL monooxygenase-like domain as they are emerging from the available in vitro studies on mouse, human and Drosophila MICAL forms that demonstrated a NADPH-dependent actin depolymerizing activity of MICAL. With Drosophila MICAL forms, actin depolymerization was demonstrated to be associated to conversion of Met44 to methionine sulfone through a postulated hydroxylating reaction. Arguments supporting the concept that MICAL effect on F-actin may be reversible will be discussed.

Insertions within the actin core of actin-related protein 3 (Arp3) modulate branching nucleation by Arp2/3 complex

The Arp2/3 (actin-related protein 2/3) complex nucleates branched actin filaments involved in multiple cellular functions, including endocytosis and cellular motility. Two subunits (Arp2 and Arp3) in this seven-subunit assembly are closely related to actin and upon activation of the complex form a "cryptic dimer" that stably mimics an actin dimer to nucleate a new filament. Both Arps contain a shared actin core structure, and each Arp contains multiple insertions of unknown function at conserved positions within the core. Here we characterize three key insertions within the actin core of Arp3 and show that each one plays a distinct role in modulating Arp2/3 function. The β4/β5 insert mediates interactions of Arp2/3 complex with actin filaments and "dampers" the nucleation activity of the complex. The Arp3 hydrophobic plug plays an important role in maintaining the integrity of the complex but is not absolutely required for formation of the daughter filament nucleus. Deletion of the αK/β15 insert did not constitutively activate the complex, as previously hypothesized. Instead, it abolished in vitro nucleation activity and caused defects in endocytic actin patch assembly in fission yeast, indicating a role for the αK/β15 insert in the activated state of the complex. Biochemical characterization of each mutant revealed steps in the nucleation pathway influenced by each Arp3-specific insert to provide new insights into the structural basis of activation of the complex.

Two deafness-causing (DFNA20/26) actin mutations affect Arp2/3-dependent actin regulation

Hearing requires proper function of the auditory hair cell, which is critically dependent upon its actin-based cytoskeletal structure. Currently, ten point mutations in nonmuscle γ-actin have been identified as causing progressive autosomal dominant nonsyndromic hearing loss (DFNA20/26), highlighting these ten residues as functionally important to actin structure and/or regulation. Two of the mutations, K118M and K118N, are located near the putative binding site for the ubiquitously expressed Arp2/3 complex. We therefore hypothesized that these mutations may affect Arp2/3-dependent regulation of the actin cytoskeleton. Using in vitro bulk polymerization assays, we show that the Lys-118 mutations notably reduce actin + Arp2/3 polymerization rates compared with WT. Further in vitro analysis of the K118M mutant using TIRF microscopy indicates the actual number of branches formed per filament is reduced compared with WT and, surprisingly, branch location is altered such that the majority of K118M branches form near the pointed end of the filament. These results highlight a previously unknown role for the Lys-118 residue in the actin-Arp2/3 interaction and also further suggest that Lys-118 may play a more significant role in intra- and intermonomer interactions than was initially hypothesized.

Optimal number of coarse-grained sites in different components of large biomolecular complexes

The computational study of large biomolecular complexes (molecular machines, cytoskeletal filaments, etc.) is a formidable challenge facing computational biophysics and biology. To achieve biologically relevant length and time scales, coarse-grained (CG) models of such complexes usually must be built and employed. One of the important early stages in this approach is to determine an optimal number of CG sites in different constituents of a complex. This work presents a systematic approach to this problem. First, a universal scaling law is derived and numerically corroborated for the intensity of the intrasite (intradomain) thermal fluctuations as a function of the number of CG sites. Second, this result is used for derivation of the criterion for the optimal number of CG sites in different parts of a large multibiomolecule complex. In the zeroth-order approximation, this approach validates the empirical rule of taking one CG site per fixed number of atoms or residues in each biomolecule, previously widely used for smaller systems (e.g., individual biomolecules). The first-order corrections to this rule are derived and numerically checked by the case studies of the Escherichia coli ribosome and Arp2/3 actin filament junction. In different ribosomal proteins, the optimal number of amino acids per CG site is shown to differ by a factor of 3.5, and an even wider spread may exist in other large biomolecular complexes. Therefore, the method proposed in this paper is valuable for the optimal construction of CG models of such complexes.