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.

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.

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.

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.

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.

PSCDB: a database for protein structural change upon ligand binding.

Proteins are flexible molecules that undergo structural changes to function. The Protein Data Bank contains multiple entries for identical proteins determined under different conditions, e.g. with and without a ligand molecule, which provides important information for understanding the structural changes related to protein functions. We gathered 839 protein structural pairs of ligand-free and ligand-bound states from monomeric or homo-dimeric proteins, and constructed the Protein Structural Change DataBase (PSCDB). In the database, we focused on whether the motions were coupled with ligand binding. As a result, the protein structural changes were classified into seven classes, i.e. coupled domain motion (59 structural changes), independent domain motion (70), coupled local motion (125), independent local motion (135), burying ligand motion (104), no significant motion (311) and other type motion (35). PSCDB provides lists of each class. On each entry page, users can view detailed information about the motion, accompanied by a morphing animation of the structural changes. PSCDB is available at http://idp1.force.cs.is.nagoya-u.ac.jp/pscdb/.

The conserved Arg241-Glu439 salt bridge determines flexibility of the inositol 1,4,5-trisphosphate receptor binding core in the ligand-free state.

Inositol 1,4,5-trisphosphate receptor (InsP3 R) is an intracellular Ca(2+) -release channel activated by binding of inositol 1,4,5-trisphosphate (InsP3 ) to the InsP3 binding core (IBC). Structural change in the IBC upon InsP3 binding is the key process in channel pore opening. In this study, we performed molecular dynamics (MD) simulations of the InsP3 -free form of the IBC, starting with removal of InsP3 from the InsP3 -bound crystal structure, and obtained the structural ensemble of the InsP3 -free form of the IBC. The simulation revealed that the two domains of the IBC largely fluctuate around the average structure with the hinge angle opened 17° more than in the InsP3 -bound form, and the twist angle rotated by 45°, forming interdomain contacts that are different from those in the bound form. The InsP3 binding loop was disordered. The InsP3 -free form thus obtained was reproduced four times in simulations started from a fully extended configuration of the two domains. Simulations beginning with the fully extended form indicated that formation of a salt bridge between Arg241 and Glu439 is crucial for stabilizing the closed form of the two domains. Mutation of Arg241 to Gln prevented formation of the compact structure by the two domains, but the fully flexible domain arrangement was maintained. Thus, the Arg241-Glu439 salt bridge determines the flexibility of the InsP3 -free form of the IBC.

Minimum free energy path of ligand-induced transition in adenylate kinase.

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.

Multiscale enhanced path sampling based on the Onsager-Machlup action: application to a model polymer.

We propose a novel path sampling method based on the Onsager-Machlup (OM) action by generalizing the multiscale enhanced sampling technique suggested by Moritsugu and co-workers [J. Chem. Phys. 133, 224105 (2010)]. The basic idea of this method is that the system we want to study (for example, some molecular system described by molecular mechanics) is coupled to a coarse-grained (CG) system, which can move more quickly and can be computed more efficiently than the original system. We simulate this combined system (original + CG system) using Langevin dynamics where different heat baths are coupled to the two systems. When the coupling is strong enough, the original system is guided by the CG system, and is able to sample the configuration and path space with more efficiency. We need to correct the bias caused by the coupling, however, by employing the Hamiltonian replica exchange, where we prepare many path replicas with different coupling strengths. As a result, an unbiased path ensemble for the original system can be found in the weakest coupling path ensemble. This strategy is easily implemented because a weight for a path calculated by the OM action is formally the same as the Boltzmann weight if we properly define the path "Hamiltonian." We apply this method to a model polymer with Asakura-Oosawa interaction, and compare the results with the conventional transition path sampling method.

SAHG, a comprehensive database of predicted structures of all human proteins.

Most proteins from higher organisms are known to be multi-domain proteins and contain substantial numbers of intrinsically disordered (ID) regions. To analyse such protein sequences, those from human for instance, we developed a special protein-structure-prediction pipeline and accumulated the products in the Structure Atlas of Human Genome (SAHG) database at http://bird.cbrc.jp/sahg. With the pipeline, human proteins were examined by local alignment methods (BLAST, PSI-BLAST and Smith-Waterman profile-profile alignment), global-local alignment methods (FORTE) and prediction tools for ID regions (POODLE-S) and homology modeling (MODELLER). Conformational changes of protein models upon ligand-binding were predicted by simultaneous modeling using templates of apo and holo forms. When there were no suitable templates for holo forms and the apo models were accurate, we prepared holo models using prediction methods for ligand-binding (eF-seek) and conformational change (the elastic network model and the linear response theory). Models are displayed as animated images. As of July 2010, SAHG contains 42,581 protein-domain models in approximately 24,900 unique human protein sequences from the RefSeq database. Annotation of models with functional information and links to other databases such as EzCatDB, InterPro or HPRD are also provided to facilitate understanding the protein structure-function relationships.

Comparative simulations of the ground state and the M-intermediate state of the sensory rhodopsin II-transducer complex with a HAMP domain model.

The complex of sensory rhodopsin II (SRII) and its cognate transducer HtrII (2:2 SRII-HtrII complex) consists of a photoreceptor and its signal transducer, respectively, associated with negative phototaxis in extreme halophiles. In this study to investigate how photoexcitation in SRII affects the structures of the complex, we conducted two series of molecular dynamics simulations of the complex of SRII and truncated HtrII (residues 1-136) of Natronomonas pharaonis linked with a modeled HAMP domain in the lipid bilayer using the two crystal structures of the ground state and the M-intermediate state as the starting structures. The simulation results showed significant enhancements of the structural differences observed between the two crystal structures. Helix F of SRII showed an outward motion, and the C-terminal end of transmembrane domain 2 (TM2) in HtrII rotated by ∼10°. The most significant structural changes were observed in the overall orientations of the two SRII molecules, closed in the ground state and open in the M-state. This change was attributed to substantial differences in the structure of the four-helix bundle of the HtrII dimer causing the apparent rotation of TM2. These simulation results established the structural basis for the various experimental observations explaining the structural differences between the ground state and the M-intermediate state.

Disorder-to-order transition of an intrinsically disordered region of sortase revealed by multiscale enhanced sampling.

Molecular functions of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs), such as molecular recognition and cellular signaling, are ascribed to dynamic changes in the conformational space in response to binding of target molecules. Sortase, a transpeptitase in Gram-positive bacteria, has an IDR in a loop which undergoes a disordered-to-ordered transition (called "disordered loop"), accompanying a tilt of another loop ("dynamic loop"), upon binding of a signal peptide and a calcium ion. In this study, all-atom conformational ensembles of sortase were calculated for the four different binding states (with/without the peptide and with/without a calcium ion) by the multiscale enhanced sampling (MSES) simulation to examine how the binding of the peptide and/or calcium influences the conformational ensemble. The MSES is a multiscale and multicopy simulation method that allows an enhanced sampling of the all-atom model of large proteins including explicit solvent. A 100 ns MSES simulation of the ligand-free sortase using 20 replicas (in total 2 µs) demonstrated large flexibility in both the disordered and dynamic loops; however, their distributions were not random but had a clear preference which populates the N-terminal part of the disordered loop near the bound form. The MSES simulations of the three binding states clarified the allosteric mechanism of sortase: the N- and C-terminal parts of the disordered loop undergo a disorder-to-order transition independently of each other upon binding of the peptide and a calcium ion, respectively; however, upon binding of both ligands, the two parts work cooperatively to stabilize the bound peptide.

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.

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.

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.

Side-chain conformational changes of transcription factor PhoB upon DNA binding: a population-shift mechanism.

Using molecular dynamics (MD) simulations and analyses of NMR relaxation order parameters, we investigated conformational changes of side chains in hydrophobic cores upon DNA binding for the DNA binding/transactivation domain of the transcription factor PhoB, in which backbone conformational changes upon DNA binding are small. The simulation results correlated well with experimental order parameters for the backbone and side-chain methyl groups, showing that the order parameters generally represent positional fluctuations of the backbone and side-chain methyl groups. However, topological effects of the side chains on the order parameters were also found and could be eliminated using normalized order parameters for each amino acid type. Consistent with the NMR experiments, the normalized order parameters from the MD simulations showed that the side chains in one of the two hydrophobic cores (the soft core) were highly flexible in comparison with those in the other hydrophobic core (the hard core) before DNA binding and that the flexibility of the hydrophobic cores, particularly of the soft core, was reduced upon DNA binding. Principal component analysis of methyl group configurations revealed strikingly different side-chain dynamics for the soft and hard cores. In the hard core, side-chain configurations were simply distributed around one or two average configurations. In contrast, the side chains in the soft core dynamically varied their configurations in an equilibrium ensemble that included binding configurations as minor components before DNA binding. DNA binding led to a restriction of the side-chain dynamics and a shift in the equilibrium toward binding configurations, in clear correspondence with a population-shift model.

Onsager-Machlup action-based path sampling and its combination with replica exchange for diffusive and multiple pathways.

For sampling multiple pathways in a rugged energy landscape, we propose a novel action-based path sampling method using the Onsager-Machlup action functional. Inspired by the Fourier-path integral simulation of a quantum mechanical system, a path in Cartesian space is transformed into that in Fourier space, and an overdamped Langevin equation is derived for the Fourier components to achieve a canonical ensemble of the path at a finite temperature. To avoid "path trapping" around an initially guessed path, the path sampling method is further combined with a powerful sampling technique, the replica exchange method. The principle and algorithm of our method is numerically demonstrated for a model two-dimensional system with a bifurcated potential landscape. The results are compared with those of conventional transition path sampling and the equilibrium theory, and the error due to path discretization is also discussed.

Latent dynamics of a protein molecule observed in dihedral angle space.

Dihedral angles are alternative set of variables to Cartesian coordinates for representing protein dynamics. The two sets of variables exhibit extremely different behavior. Motions in dihedral angle space are characterized by latent dynamics, in which motion induced in each dihedral angle is always compensated for by motions of many other dihedral angles, in order to maintain a rigid globular shape. Using molecular dynamics simulations, we propose a molecular mechanism for the latent dynamics in dihedral angle space. It was found that, due to the unique structure of dihedral principal components originating in the globular shape of the protein, the dihedral principal components with large (small) amplitudes are highly correlated with the eigenvectors of the metric matrix with small (large) eigenvalues. Such an anticorrelation in the eigenmode structures minimizes the mean square displacement of Cartesian coordinates upon rotation of dihedral angles. In contrast, a short peptide, deca-alanine in this study, does not show such behavior of the latent dynamics in the dihedral principal components, but shows similar behaviors to those of the Cartesian principal components, due to the absence of constraints to maintain a rigid globular shape.

Normal mode analysis of protein dynamics in a non-Eckart frame.

Normal mode analysis, with the all-atom or coarse-grained elastic network model, represents the equilibrium fluctuation of protein molecule in the Eckart frame, where contributions from external motions (translation and rotation) of the entire protein molecule are eliminated. On the other hand, domain motion is frequently exhibited by the relative motion of one domain to the other. Such a representation of fluctuations in the non-Eckart frame cannot be achieved by conventional normal mode analysis. Here, we propose normal mode analysis in a non-Eckart frame, where the external degrees of freedom are fixed for any portion of the system. In this analysis, the covariance matrix in the Eckart frame is transformed into one in the non-Eckart frame. Using a molecular dynamics simulation, we have confirmed the validity of the transformation formula and discussed the physical implication of the formula.