Binding and Unbinding Pathways of Peptide Substrates on the SARS-CoV-2 3CL Protease.

Based on many crystal structures of ligand complexes, much study has been devoted to understanding the molecular recognition of SARS-CoV-2 3C-like protease (3CLpro), a potent drug target for COVID-19. In this research, to extend this present static view, we examined the kinetic process of binding/unbinding of an eight-residue substrate peptide to/from 3CLpro by evaluating the path ensemble with the weighted ensemble simulation. The path ensemble showed the mechanism of how a highly flexible peptide folded into the bound form. At the early stage, the dominant motion was the diffusion on the protein surface showing a broad distribution, whose center was led into the cleft of the chymotrypsin fold. We observed a definite sequential formation of the hydrogen bonds at the later stage occurring in the cleft, initiated between Glu166 (3CLpro) and P3_Val (peptide), followed by binding to the oxyanion hole and completed by the sequence-specific recognition at P1_Gln.

Functional dynamics of SARS-CoV-2 3C-like protease as a member of clan PA.

SARS-CoV-2 3C-like protease (3CLpro), a potential therapeutic target for COVID-19, consists of a chymotrypsin fold and a C-terminal α-helical domain (domain III), the latter of which mediates dimerization required for catalytic activation. To gain further understanding of the functional dynamics of SARS-CoV-2 3CLpro, this review extends the scope to the comparative study of many crystal structures of proteases having the chymotrypsin fold (clan PA of the MEROPS database). First, the close correspondence between the zymogen-enzyme transformation in chymotrypsin and the allosteric dimerization activation in SARS-CoV-2 3CLpro is illustrated. Then, it is shown that the 3C-like proteases of family Coronaviridae (the protease family C30), which are closely related to SARS-CoV-2 3CLpro, have the same homodimeric structure and common activation mechanism via domain III mediated dimerization. The survey extended to order Nidovirales reveals that all 3C-like proteases belonging to Nidovirales have domain III, but with various chain lengths, and 3CLpro of family Mesoniviridae (family C107) has the same homodimeric structure as that of C30, even though they have no sequence similarity. As a reference, monomeric 3C proteases belonging to the more distant family Picornaviridae (family C3) lacking domain III are compared with C30, and it is shown that the 3C proteases are rigid enough to maintain their structures in the active state. Supplementary Information: The online version contains supplementary material available at 10.1007/s12551-022-01020-x.

Molecular basis of ubiquitin-specific protease 8 autoinhibition by the WW-like domain.

Ubiquitin-specific protease 8 (USP8) is a deubiquitinating enzyme involved in multiple membrane trafficking pathways. The enzyme activity is inhibited by binding to 14-3-3 proteins. Mutations in the 14-3-3-binding motif in USP8 are related to Cushing's disease. However, the molecular basis of USP8 activity regulation remains unclear. This study identified amino acids 645-684 of USP8 as an autoinhibitory region, which might interact with the catalytic USP domain, as per the results of pull-down and single-molecule FRET assays performed in this study. In silico modelling indicated that the region forms a WW-like domain structure, plugs the catalytic cleft, and narrows the entrance to the ubiquitin-binding pocket. Furthermore, 14-3-3 inhibited USP8 activity partly by enhancing the interaction between the WW-like and USP domains. These findings provide the molecular basis of USP8 autoinhibition via the WW-like domain. Moreover, they suggest that the release of autoinhibition may underlie Cushing's disease due to USP8 mutations.

Allosteric Regulation of 3CL Protease of SARS-CoV-2 and SARS-CoV Observed in the Crystal Structure Ensemble.

The 3C-like protease (3CLpro) of SARS-CoV-2 is a potential therapeutic target for COVID-19. Importantly, it has an abundance of structural information solved as a complex with various drug candidate compounds. Collecting these crystal structures (83 Protein Data Bank (PDB) entries) together with those of the highly homologous 3CLpro of SARS-CoV (101 PDB entries), we constructed the crystal structure ensemble of 3CLpro to analyze the dynamic regulation of its catalytic function. The structural dynamics of the 3CLpro dimer observed in the ensemble were characterized by the motions of four separate loops (the C-loop, E-loop, H-loop, and Linker) and the C-terminal domain III on the rigid core of the chymotrypsin fold. Among the four moving loops, the C-loop (also known as the oxyanion binding loop) causes the order (active)-disorder (collapsed) transition, which is regulated cooperatively by five hydrogen bonds made with the surrounding residues. The C-loop, E-loop, and Linker constitute the major ligand binding sites, which consist of a limited variety of binding residues including the substrate binding subsites. Ligand binding causes a ligand size dependent conformational change to the E-loop and Linker, which further stabilize the C-loop via the hydrogen bond between the C-loop and E-loop. The T285A mutation from SARS-CoV 3CLpro to SARS-CoV-2 3CLpro significantly closes the interface of the domain III dimer and allosterically stabilizes the active conformation of the C-loop via hydrogen bonds with Ser1 and Gly2; thus, SARS-CoV-2 3CLpro seems to have increased activity relative to that of SARS-CoV 3CLpro.

Inter-lobe Motions Allosterically Regulate the Structure and Function of EGFR Kinase.

Protein kinases play important roles in cellular signaling and have been one of the best-studied drug targets. The kinase domain of epidermal growth factor receptor (EGFR) is a receptor tyrosine kinase that has been extensively studied for cancer drug discovery and for understanding the unique activation mechanism by dimerization. Here, we analyzed all available 206 crystal structures of the EGFR kinase and the dynamics observed in molecular simulations to identify how these structures are determined. It was found that the arrangement between the N- and C-terminal lobes plays a key role in regulating the kinase structure by sensitively responding to the intermolecular interactions, or the crystal environment. A whole variety of crystal forms in the database is thus reflected in the broad distribution of the inter-lobe arrangement. The configuration of the catalytically important motifs as well as the bound ATP is closely coupled with the inter-lobe motion. When the intermolecular interactions are those of the activating asymmetric dimer, EGFR kinase takes the open-lobe arrangement that constructs the catalytically active configuration.

Allosteric response to ligand binding: Molecular dynamics study of the N-terminal domains in IP 3 receptor.

Inositol 1,4,5-trisphosphate (IP3) receptor (IP3R) is a huge tetrameric intracellular Ca2+ channel that mediates cytoplasmic Ca2+ signaling. The structural basis of the gating in IP3R has been studied by X-ray crystallography and cryo-electron microscopy, focusing on the domain rearrangements triggered by IP3 binding. Here, we conducted molecular dynamics (MD) simulations of the three N-terminal domains of IP3R responsible for IP3 binding (IBC/SD; two domains of the IP3 binding core, IBCß and IBCα, and suppressor domain, SD) as a model system to study the initial gating stage. The response upon removal of IP3 from the IP3-bound form of IBC/SD was traced in MD trajectories. The two IBC domains showed an immediate response of opening after removal of IP3, and SD showed a simultaneous opening motion indicating a tight dynamic coupling with IBC. However, when IBC remained in a more closed form, the dynamic coupling broke and SD exhibited a more amplified closing motion independently of IBC. This amplified SD motion was caused by the break of connection between SD and IBCß at the hinge region, but was suppressed in the native tetrameric state. The analyses using Motion Tree and the linear response theory clarified that in the open form, SD and IBCα moved collectively relative to IBCß with a response upon IP3 binding within the linear regime, whereas in the closed form, such collectiveness disappeared. These results suggest that the regulation of dynamics via the domain arrangement and multimerization is requisite for large-scale allosteric communication in IP3R gating machinery.

Dynamic recognition and linkage specificity in K63 di-ubiquitin and TAB2 NZF domain complex.

Poly-ubiquitin (poly-Ub) is involved in various cellular processes through the linkage-specific recognition of Ub-binding domains (UBD). In this study, using molecular dynamics (MD) simulation together with an enhanced sampling method, we demonstrated that K63-linked di-Ub recognizes the NZF domain of TAB2, a zinc finger UBD, in an ensemble of highly dynamic structures that form from the weak interactions between UBD and the flexible linker connecting the two Ubs. However, the K63 di-Ub/TAB2 NZF complex showed a much more compact and stable ensemble than the non-native complexes, linear di-Ub/TAB2 NZF and K33 di-Ub/TAB2 NZF, that were modeled from linear di-Ub/HOIL-1L NZF and K33 di-Ub/TRABID NZF1, respectively. We further demonstrated the importance of the length and position of the Ub-Ub linker in the results of MD simulations of K63 di-Ub/TAB2 NZF by changing the Ub linkage from the native K63 to four different non-native linkages, linear, K6, K11, and K48, while maintaining inter-molecular contacts in the native complex. No systems with non-native linkage maintained the native binding configuration. These simulation results provide an atomistic picture of the linkage specific recognition of poly-Ubs leading to the biological functions such as cellular colocalization of various component proteins in the signal transduction pathways.

Crystal structure of the Ube2K/E2-25K and K48-linked di-ubiquitin complex provides structural insight into the mechanism of K48-specific ubiquitin chain synthesis.

Ubiquitin-conjugating enzymes (E2) form thioester bonds with ubiquitin (Ub), which are subsequently transferred to target proteins for cellular progress. Ube2K/E2-25K (a class II E2 enzyme) contains a C-terminal ubiquitin-associated (UBA) domain that has been suggested to control ubiquitin recognition, dimerization, or poly-ubiquitin chain formation. Ube2K is a special E2 because it synthesizes K48-linked poly-ubiquitin chains without E3 ubiquitin ligase. We found that a novel interaction between the acceptor di-Ub (Ub2) and the auxiliary Ube2K promotes the discharging reaction and production of tri-Ub (Ub3), probably by guiding and positioning the K48 (in the distal Ub) of the acceptor Ub2 in the active site. We also determined the crystal structure of Ube2K-Ub2 at 2.47 Šresolution. Based on our structural and biochemical data, we proposed a structural model of Ub3 synthesis by Ube2K without E3.

Multiscale enhanced sampling of glucokinase: Regulation of the enzymatic reaction via a large scale domain motion.

Enhanced sampling yields a comprehensive structural ensemble or a free energy landscape, which is beyond the capability of a conventional molecular dynamics simulation. Our recently developed multiscale enhanced sampling (MSES) method employs a coarse-grained model coupled with the target physical system for the efficient acceleration of the dynamics. MSES has demonstrated applicability to large protein systems in solution, such as intrinsically disordered proteins and protein-protein and protein-ligand interactions. Here, we applied the MSES simulation to an important drug discovery target, glucokinase (GCK), to elucidate the structural basis of the positive cooperativity of the enzymatic reaction at an atomistic resolution. MSES enabled us to compare two sets of the free energy landscapes of GCK, for the glucose-bound and glucose-unbound forms, and thus demonstrated the drastic change of the free energy surface depending on the glucose concentration. In the glucose-bound form, we found two distinct basins separated by a high energy barrier originating from the domain motion and the folding/unfolding of the α13 helix. By contrast, in the glucose-unbound form, a single flat basin extended to the open and super-open states. These features illustrated the two distinct phases achieving the cooperativity, the fast reaction cycle staying in the closed state at a high glucose concentration and the slow cycle primarily in the open/super-open state at a low concentration. The weighted ensemble simulations revealed the kinetics of the structural changes in GCK with the synergetic use of the MSES results; the rate constant of the transition between the closed state and the open/super-open states, kC/O = 1.1 ms-1, is on the same order as the experimental catalytic rate, kcat = 0.22 ms-1. Finally, we discuss the pharmacological activities of GCK activators (small molecular drugs modulating the GCK activity) in terms of the slight changes in the domain motion, depending on their chemical structures as regulators. The present study demonstrated the capability of the enhanced sampling and the associated kinetic calculations for understanding the atomistic structural dynamics of protein systems in physiological environments.

Ionic scattering factors of atoms that compose biological molecules.

Ionic scattering factors of atoms that compose biological molecules have been computed by the multi-configuration Dirac-Fock method. These ions are chemically unstable and their scattering factors had not been reported except for O-. Yet these factors are required for the estimation of partial charges in protein molecules and nucleic acids. The electron scattering factors of these ions are particularly important as the electron scattering curves vary considerably between neutral and charged atoms in the spatial-resolution range explored in structural biology. The calculated X-ray and electron scattering factors have then been parameterized for the major scattering curve models used in X-ray and electron protein crystallography and single-particle cryo-EM. The X-ray and electron scattering factors and the fitting parameters are presented for future reference.

Energetics and conformational pathways of functional rotation in the multidrug transporter AcrB.

The multidrug transporter AcrB transports a broad range of drugs out of the cell by means of the proton-motive force. The asymmetric crystal structure of trimeric AcrB suggests a functionally rotating mechanism for drug transport. Despite various supportive forms of evidence from biochemical and simulation studies for this mechanism, the link between the functional rotation and proton translocation across the membrane remains elusive. Here, calculating the minimum free energy pathway of the functional rotation for the complete AcrB trimer, we describe the structural and energetic basis behind the coupling between the functional rotation and the proton translocation at atomic resolution. Free energy calculations show that protonation of Asp408 in the transmembrane portion of the drug-bound protomer drives the functional rotation. The conformational pathway identifies vertical shear motions among several transmembrane helices, which regulate alternate access of water in the transmembrane as well as peristaltic motions that pump drugs in the periplasm.

Structure of the Dnmt1 Reader Module Complexed with a Unique Two-Mono-Ubiquitin Mark on Histone H3 Reveals the Basis for DNA Methylation Maintenance.

The proper location and timing of Dnmt1 activation are essential for DNA methylation maintenance. We demonstrate here that Dnmt1 utilizes two-mono-ubiquitylated histone H3 as a unique ubiquitin mark for its recruitment to and activation at DNA methylation sites. The crystal structure of the replication foci targeting sequence (RFTS) of Dnmt1 in complex with H3-K18Ub/23Ub reveals striking differences to the known ubiquitin-recognition structures. The two ubiquitins are simultaneously bound to the RFTS with a combination of canonical hydrophobic and atypical hydrophilic interactions. The C-lobe of RFTS, together with the K23Ub surface, also recognizes the N-terminal tail of H3. The binding of H3-K18Ub/23Ub results in spatial rearrangement of two lobes in the RFTS, suggesting the opening of its active site. Actually, incubation of Dnmt1 with H3-K18Ub/23Ub increases its catalytic activity in vitro. Our results therefore shed light on the essential role of a unique ubiquitin-binding module in DNA methylation maintenance.

Free-Energy Landscape of Protein-Ligand Interactions Coupled with Protein Structural Changes.

Protein-ligand interactions are frequently coupled with protein structural changes. Focusing on the coupling, we present the free-energy surface (FES) of the ligand-binding process for glutamine-binding protein (GlnBP) and its ligand, glutamine, in which glutamine binding accompanies large-scale domain closure. All-atom simulations were performed in explicit solvents by multiscale enhanced sampling (MSES), which adopts a multicopy and multiscale scheme to achieve enhanced sampling of systems with a large number of degrees of freedom. The structural ensemble derived from the MSES simulation yielded the FES of the coupling, described in terms of both the ligand's and protein's degrees of freedom at atomic resolution, and revealed the tight coupling between the two degrees of freedom. The derived FES led to the determination of definite structural states, which suggested the dominant pathways of glutamine binding to GlnBP: first, glutamine migrates via diffusion to form a dominant encounter complex with Arg75 on the large domain of GlnBP, through strong polar interactions. Subsequently, the closing motion of GlnBP occurs to form ligand interactions with the small domain, finally completing the native-specific complex structure. The formation of hydrogen bonds between glutamine and the small domain is considered to be a rate-limiting step, inducing desolvation of the protein-ligand interface to form the specific native complex. The key interactions to attain high specificity for glutamine, the "door keeper" existing between the two domains (Asp10-Lys115) and the "hydrophobic sandwich" formed between the ligand glutamine and Phe13/Phe50, have been successfully mapped on the pathway derived from the FES.

Extended string-like binding of the phosphorylated HP1α N-terminal tail to the lysine 9-methylated histone H3 tail.

The chromodomain of HP1α binds directly to lysine 9-methylated histone H3 (H3K9me). This interaction is enhanced by phosphorylation of serine residues in the N-terminal tail of HP1α by unknown mechanism. Here we show that phosphorylation modulates flexibility of HP1α's N-terminal tail, which strengthens the interaction with H3. NMR analysis of HP1α's chromodomain with N-terminal tail reveals that phosphorylation does not change the overall tertiary structure, but apparently reduces the tail dynamics. Small angle X-ray scattering confirms that phosphorylation contributes to extending HP1α's N-terminal tail. Systematic analysis using deletion mutants and replica exchange molecular dynamics simulations indicate that the phosphorylated serines and following acidic segment behave like an extended string and dynamically bind to H3 basic residues; without phosphorylation, the most N-terminal basic segment of HP1α inhibits interaction of the acidic segment with H3. Thus, the dynamic string-like behavior of HP1α's N-terminal tail underlies the enhancement in H3 binding due to phosphorylation.

Itinerary profiling to analyze a large number of protein-folding trajectories.

Understanding how proteins fold through a vast number of unfolded states is a major subject in the study of protein folding. Herein, we present itinerary profiling as a simple method to analyze molecular dynamics trajectories, and apply this method to Trp-cage. In itinerary profiling, structural clusters included in a trajectory are represented by a bit sequence, and a number of trajectories, as well as the structural clusters, can be compared and classified. As a consequence, the structural clusters that characterize the foldability of trajectories were able to be identified. The connections between the clusters were then illustrated as a network and the structural features of the clusters were examined. We found that in the true folding funnel, Trp-cage formed a left-handed main-chain topology and the Trp6 side-chain was located at the front of the main-chain ring, even in the initial unfolded states. In contrast, in the false folding funnel of the pseudo-native states, in which the Trp6 side-chain is upside down in the protein core, Trp-cage had a right-handed main-chain topology and the Trp side-chain was at the back. The initial topological partition, as determined by the main-chain handedness and the location of the Trp residue, predetermines Trp-cage foldability and the destination of the trajectory to the native state or the pseudo-native states.

Motion Tree Delineates Hierarchical Structure of Protein Dynamics Observed in Molecular Dynamics Simulation.

Molecular dynamics (MD) simulations of proteins provide important information to understand their functional mechanisms, which are, however, likely to be hidden behind their complicated motions with a wide range of spatial and temporal scales. A straightforward and intuitive analysis of protein dynamics observed in MD simulation trajectories is therefore of growing significance with the large increase in both the simulation time and system size. In this study, we propose a novel description of protein motions based on the hierarchical clustering of fluctuations in the inter-atomic distances calculated from an MD trajectory, which constructs a single tree diagram, named a "Motion Tree", to determine a set of rigid-domain pairs hierarchically along with associated inter-domain fluctuations. The method was first applied to the MD trajectory of substrate-free adenylate kinase to clarify the usefulness of the Motion Tree, which illustrated a clear-cut dynamics picture of the inter-domain motions involving the ATP/AMP lid and the core domain together with the associated amplitudes and correlations. The comparison of two Motion Trees calculated from MD simulations of ligand-free and -bound glutamine binding proteins clarified changes in inherent dynamics upon ligand binding appeared in both large domains and a small loop that stabilized ligand molecule. Another application to a huge protein, a multidrug ATP binding cassette (ABC) transporter, captured significant increases of fluctuations upon binding a drug molecule observed in both large scale inter-subunit motions and a motion localized at a transmembrane helix, which may be a trigger to the subsequent structural change from inward-open to outward-open states to transport the drug molecule. These applications demonstrated the capabilities of Motion Trees to provide an at-a-glance view of various sizes of functional motions inherent in the complicated MD trajectory.

Sequential data assimilation for single-molecule FRET photon-counting data.

Data assimilation is a statistical method designed to improve the quality of numerical simulations in combination with real observations. Here, we develop a sequential data assimilation method that incorporates one-dimensional time-series data of smFRET (single-molecule Förster resonance energy transfer) photon-counting into conformational ensembles of biomolecules derived from "replicated" molecular dynamics (MD) simulations. A particle filter using a large number of "replicated" MD simulations with a likelihood function for smFRET photon-counting data is employed to screen the conformational ensembles that match the experimental data. We examine the performance of the method using emulated smFRET data and coarse-grained (CG) MD simulations of a dye-labeled polyproline-20. The method estimates the dynamics of the end-to-end distance from smFRET data as well as revealing that of latent conformational variables. The particle filter is also able to correct model parameter dependence in CG MD simulations. We discuss the applicability of the method to real experimental data for conformational dynamics of biomolecules.

Energy landscape of all-atom protein-protein interactions revealed by multiscale enhanced sampling.

Protein-protein interactions are regulated by a subtle balance of complicated atomic interactions and solvation at the interface. To understand such an elusive phenomenon, it is necessary to thoroughly survey the large configurational space from the stable complex structure to the dissociated states using the all-atom model in explicit solvent and to delineate the energy landscape of protein-protein interactions. In this study, we carried out a multiscale enhanced sampling (MSES) simulation of the formation of a barnase-barstar complex, which is a protein complex characterized by an extraordinary tight and fast binding, to determine the energy landscape of atomistic protein-protein interactions. The MSES adopts a multicopy and multiscale scheme to enable for the enhanced sampling of the all-atom model of large proteins including explicit solvent. During the 100-ns MSES simulation of the barnase-barstar system, we observed the association-dissociation processes of the atomistic protein complex in solution several times, which contained not only the native complex structure but also fully non-native configurations. The sampled distributions suggest that a large variety of non-native states went downhill to the stable complex structure, like a fast folding on a funnel-like potential. This funnel landscape is attributed to dominant configurations in the early stage of the association process characterized by near-native orientations, which will accelerate the native inter-molecular interactions. These configurations are guided mostly by the shape complementarity between barnase and barstar, and lead to the fast formation of the final complex structure along the downhill energy landscape.

Solvent friction effects propagate over the entire protein molecule through low-frequency collective modes.

Protein solvation dynamics has been investigated using atom-dependent Langevin friction coefficients derived directly from molecular dynamics (MD) simulations. To determine the effect of solvation on the atomic friction coefficients, solution and vacuum MD simulations were performed for lysozyme and staphylococcal nuclease and analyzed by Langevin mode analysis. The coefficients thus derived are roughly correlated with the atomic solvent-accessible surface area (ASA), as expected from the fact that friction occurs as the result of collisions with solvent molecules. However, a considerable number of atoms with higher friction coefficients are found inside the core region. Hence, the influence of solvent friction propagates into the protein core. The internal coefficients have large contributions from the low-frequency modes, yielding a simple picture of the surface-to-core long-range damping via solvation governed by collective low-frequency modes. To make use of these findings in implicit-solvent modeling, we compare the all-atom friction results with those obtained using Langevin dynamics (LD) with two empirical representations: the constant-friction and the ASA-dependent (Pastor-Karplus) friction models. The constant-friction model overestimates the core and underestimates the surface damping whereas the ASA-dependent friction model, which damps protein atoms only on the solvent-accessible surface, reproduces well the friction coefficients for both the surface and core regions observed in the explicit-solvent MD simulations. Therefore, in LD simulation, the solvent friction coefficients should be imposed only on the protein surface.

Hierarchical description and extensive classification of protein structural changes by Motion Tree.

The structures of the same protein, determined under different conditions, provide clues toward understanding the role of structural changes in the protein's function. Structural changes are usually identified as rigid-body motions, which are defined using a particular threshold of rigidity, such as domain motions. However, each protein actually undergoes motions with various size and magnitude ranges. In this study, to describe protein structural changes more comprehensively, we propose a method based on hierarchical clustering. This method enables the illustration of a wide range of protein motions in a single tree diagram, named the "Motion Tree". We applied the method to 432 proteins exhibiting large structural changes and classified their Motion Trees in terms of the characteristic indices of the trees. This classification of the Motion Trees revealed clear relationships to their protein functions. Especially, complex structural changes are significantly correlated with multi-step protein functions.