B56δ long-disordered arms form a dynamic PP2A regulation interface coupled with global allostery and Jordan's syndrome mutations.

Intrinsically disordered regions (IDR) and short linear motifs (SLiMs) play pivotal roles in the intricate signaling networks governed by phosphatases and kinases. B56δ (encoded by PPP2R5D) is a regulatory subunit of protein phosphatase 2A (PP2A) with long IDRs that harbor a substrate-mimicking SLiM and multiple phosphorylation sites. De novo missense mutations in PPP2R5D cause intellectual disabilities (ID), macrocephaly, Parkinsonism, and a broad range of neurological symptoms. Our single-particle cryo-EM structures of the PP2A-B56δ holoenzyme reveal that the long, disordered arms at the B56δ termini fold against each other and the holoenzyme core. This architecture suppresses both the phosphatase active site and the substrate-binding protein groove, thereby stabilizing the enzyme in a closed latent form with dual autoinhibition. The resulting interface spans over 190 Šand harbors unfavorable contacts, activation phosphorylation sites, and nearly all residues with ID-associated mutations. Our studies suggest that this dynamic interface is coupled to an allosteric network responsive to phosphorylation and altered globally by mutations. Furthermore, we found that ID mutations increase the holoenzyme activity and perturb the phosphorylation rates, and the severe variants significantly increase the mitotic duration and error rates compared to the normal variant.

Molecular mechanism of GTP binding- and dimerization-induced enhancement of Sar1-mediated membrane remodeling.

The Sar1 GTPase initiates coat protein II (COPII)-mediated protein transport by generating membrane curvature at subdomains on the endoplasmic reticulum, where it is activated by the guanine nucleotide exchange factor (GEF) Sec12. Crystal structures of GDP- and GTP-bound forms of Sar1 suggest that it undergoes a conformational switch in which GTP binding enhances the exposure of an amino-terminal amphipathic helix necessary for efficient membrane penetration. However, key residues in the amino terminus were not resolved in crystal structures, and experimental studies have suggested that the amino terminus of Sar1 is solvent-exposed in the absence of a membrane, even in the GDP-bound state. Therefore, the molecular mechanism by which GTP binding activates the membrane-remodeling activity of Sar1 remains unclear. Using atomistic molecular dynamics simulations, we compare the membrane-binding and curvature generation activities of Sar1 in its GDP- and GTP-bound states. We show that in the GTP-bound state, Sar1 inserts into the membrane with its complete (residues 1 to 23) amphipathic amino-terminal helix, while Sar1-GDP binds to the membrane only through its first 12 residues. Such differential membrane-binding modes translate into significant differences in the protein volume inserted into the membrane. As a result, Sar1-GTP generates positive membrane curvature 10 to 20 times higher than Sar1-GDP. Dimerization of the GTP-bound form of Sar1 further amplifies curvature generation. Taken together, our results present a detailed molecular mechanism for how the nucleotide-bound state of Sar1 regulates its membrane-binding and remodeling activities in a concentration-dependent manner, paving the way toward a better understanding COPII-mediated membrane transport.

Sensitive and selective polymer condensation at membrane surface driven by positive co-operativity.

Biomolecular phase separation has emerged as an essential mechanism for cellular organization. How cells respond to environmental stimuli in a robust and sensitive manner to build functional condensates at the proper time and location is only starting to be understood. Recently, lipid membranes have been recognized as an important regulatory center for biomolecular condensation. However, how the interplay between the phase behaviors of cellular membranes and surface biopolymers may contribute to the regulation of surface condensation remains to be elucidated. Using simulations and a mean-field theoretical model, we show that two key factors are the membrane's tendency to phase-separate and the surface polymer's ability to reorganize local membrane composition. Surface condensate forms with high sensitivity and selectivity in response to features of biopolymer when positive co-operativity is established between coupled growth of the condensate and local lipid domains. This effect relating the degree of membrane-surface polymer co-operativity and condensate property regulation is shown to be robust by different ways of tuning the co-operativity, such as varying membrane protein obstacle concentration, lipid composition, and the affinity between lipid and polymer. The general physical principle emerged from the current analysis may have implications in other biological processes and beyond.

Flexibility of Binding Site is Essential to the Ca2+ Selectivity in EF-Hand Calcium-Binding Proteins.

High binding affinity and selectivity of metal ions are essential to the function of metalloproteins. Thus, understanding the factors that determine these binding characteristics is of major interest for both fundamental mechanistic investigations and guiding of the design of novel metalloproteins. In this work, we perform QM cluster model calculations and quantum mechanics/molecular mechanics (QM/MM) free energy simulations to understand the binding selectivity of Ca2+ and Mg2+ in the wild-type carp parvalbumin and its mutant. While a nonpolarizable MM model (CHARMM36) does not lead to the correct experimental trend, treatment of the metal binding site with the DFTB3 model in a QM/MM framework leads to relative binding free energies (ΔΔGbind) comparable with experimental data. For the wild-type (WT) protein, the calculated ΔΔGbind is ∼6.6 kcal/mol in comparison with the experimental value of 5.6 kcal/mol. The good agreement highlights the value of a QM description of the metal binding site and supports the role of electronic polarization and charge transfer to metal binding selectivity. For the D51A/E101D/F102W mutant, different binding site models lead to considerable variations in computed binding affinities. With a coordination number of seven for Ca2+, which is shown by QM/MM metadynamics simulations to be the dominant coordination number for the mutant, the calculated relative binding affinity is ∼4.8 kcal/mol, in fair agreement with the experimental value of 1.6 kcal/mol. The WT protein is observed to feature a flexible binding site that accommodates a range of coordination numbers for Ca2+, which is essential to the high binding selectivity for Ca2+ over Mg2+. In the mutant, the E101D mutation reduces the flexibility of the binding site and limits the dominant coordination number of Ca2+ to be seven, thereby leading to reduced binding selectivity against Mg2+. Our results highlight that the binding selectivity of metal ions depends on both the structural and dynamical properties of the protein binding site.

Modulation of Allostery with Multiple Mechanisms by Hotspot Mutations in TetR.

Modulating allosteric coupling offers unique opportunities for biomedical applications. Such efforts can benefit from efficient prediction and evaluation of allostery hotspot residues that dictate the degree of cooperativity between distant sites. We demonstrate that effects of allostery hotspot mutations can be evaluated qualitatively and semiquantitatively by molecular dynamics simulations in a bacterial tetracycline repressor (TetR). The simulations recapitulate the effects of these mutations on abolishing the induction function of TetR and provide a rationale for the different rescuabilities observed to restore allosteric coupling of the hotspot mutations. We demonstrate that the same noninducible phenotype could be the result of perturbations in distinct structural and energetic properties of TetR. Our work underscores the value of explicitly computing the functional free energy landscapes to effectively evaluate and rank hotspot mutations despite the prevalence of compensatory interactions and therefore provides quantitative guidance to allostery modulation for therapeutic and engineering applications.

Dynamics and number of trans-SNARE complexes determine nascent fusion pore properties.

The fusion pore is the first crucial intermediate formed during exocytosis, yet little is known about the mechanisms that determine the size and kinetic properties of these transient structures. Here, we reduced the number of available SNAREs (proteins that mediate vesicle fusion) in neurons and observed changes in transmitter release that are suggestive of alterations in fusion pores. To investigate these changes, we employed reconstituted fusion assays using nanodiscs to trap pores in their initial open state. Optical measurements revealed that increasing the number of SNARE complexes enhanced the rate of release from single pores and enabled the escape of larger cargoes. To determine whether this effect was due to changes in nascent pore size or to changes in stability, we developed an approach that uses nanodiscs and planar lipid bilayer electrophysiology to afford microsecond resolution at the single event level. Both pore size and stability were affected by SNARE copy number. Increasing the number of vesicle (v)-SNAREs per nanodisc from three to five caused a twofold increase in pore size and decreased the rate of pore closure by more than three orders of magnitude. Moreover, pairing of v-SNAREs and target (t)-SNAREs to form trans-SNARE complexes was highly dynamic: flickering nascent pores closed upon addition of a v-SNARE fragment, revealing that the fully assembled, stable SNARE complex does not form at this stage of exocytosis. Finally, a deletion at the base of the SNARE complex, which mimics the action of botulinum neurotoxin A, markedly reduced fusion pore stability. In summary, trans-SNARE complexes are dynamic, and the number of SNAREs recruited to drive fusion determines fundamental properties of individual pores.

Multiple deprotonation paths of the nucleophile 3'-OH in the DNA synthesis reaction.

DNA synthesis by polymerases is essential for life. Deprotonation of the nucleophile 3'-OH is thought to be the obligatory first step in the DNA synthesis reaction. We have examined each entity surrounding the nucleophile 3'-OH in the reaction catalyzed by human DNA polymerase (Pol) η and delineated the deprotonation process by combining mutagenesis with steady-state kinetics, high-resolution structures of in crystallo reactions, and molecular dynamics simulations. The conserved S113 residue, which forms a hydrogen bond with the primer 3'-OH in the ground state, stabilizes the primer end in the active site. Mutation of S113 to alanine destabilizes primer binding and reduces the catalytic efficiency. Displacement of a water molecule that is hydrogen bonded to the 3'-OH using the 2'-OH of a ribonucleotide or 2'-F has little effect on catalysis. Moreover, combining the S113A mutation with 2'-F replacement, which removes two potential hydrogen acceptors of the 3'-OH, does not reduce the catalytic efficiency. We conclude that the proton can leave the O3' via alternative paths, supporting the hypothesis that binding of the third Mg2+ initiates the reaction by breaking the α-ß phosphodiester bond of an incoming deoxyribonucleoside triphosphate (dNTP).

Mapping temperature-dependent conformational change in the voltage-sensing domain of an engineered heat-activated K+ channel.

Temperature-dependent regulation of ion channel activity is critical for a variety of physiological processes ranging from immune response to perception of noxious stimuli. Our understanding of the structural mechanisms that underlie temperature sensing remains limited, in part due to the difficulty of combining high-resolution structural analysis with temperature stimulus. Here, we use NMR to compare the temperature-dependent behavior of Shaker potassium channel voltage sensor domain (WT-VSD) to its engineered temperature sensitive (TS-VSD) variant. Further insight into the molecular basis for temperature-dependent behavior is obtained by analyzing the experimental results together with molecular dynamics simulations. Our studies reveal that the overall secondary structure of the engineered TS-VSD is identical to the wild-type channels except for local changes in backbone torsion angles near the site of substitution (V369S and F370S). Remarkably however, these structural differences result in increased hydration of the voltage-sensing arginines and the S4-S5 linker helix in the TS-VSD at higher temperatures, in contrast to the WT-VSD. These findings highlight how subtle differences in the primary structure can result in large-scale changes in solvation and thereby confer increased temperature-dependent activity beyond that predicted by linear summation of solvation energies of individual substituents.

O to bR transition in bacteriorhodopsin occurs through a proton hole mechanism.

Extensive classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations are used to establish the structural features of the O state in bacteriorhodopsin (bR) and its conversion back to the bR ground state. The computed free energy surface is consistent with available experimental data for the kinetics and thermodynamics of the O to bR transition. The simulation results highlight the importance of the proton release group (PRG, consisting of Glu194/204) and the conserved arginine 82 in modulating the hydration level of the protein cavity. In particular, in the O state, deprotonation of the PRG and downward rotation of Arg82 lead to elevated hydration level and a continuous water network that connects the PRG to the protonated Asp85. Proton exchange through this water network is shown by ∼0.1-µs semiempirical QM/MM free energy simulations to occur through the generation and propagation of a proton hole, which is relayed by Asp212 and stabilized by Arg82. This mechanism provides an explanation for the observation that the D85S mutant of bacteriorhodopsin pumps chloride ions. The electrostatics-hydration coupling mechanism and the involvement of all titration states of water are likely applicable to many biomolecules involved in bioenergetic transduction.

The impact of Chinese airport infrastructure on airline pollutant emissions: A hybrid stochastic-neural network approach based on utility functions.

With China being the world's largest emitter of greenhouse gases and its aviation sector burgeoning, the environmental performance of Chinese airlines has global significance. Amidst rising demands for eco-friendly practices from both customers and regulators, the interplay between airport infrastructure and environmental performance becomes pivotal. This research offers an innovative methodology to gauge the environmental performance of Chinese airlines, emphasizing the distance traveled between airports using weighted additive utility functions. Leveraging neural networks, the study investigates the impact of various airport infrastructural characteristics on environmental performance. Noteworthy findings indicate that ground control measures, automatic information services at origin airports, surface concrete on runways at both ends, and a centerline lighting system in destination airports positively influence environmental performance. In contrast, longer and wider runways at origin airports, increased distances to control towers, and asphalt runways at destination airports adversely affect it. These insights not only underscore the importance of strategic infrastructure enhancements for reducing carbon footprints but also hold profound policy implications. As global climate change remains at the forefront, fostering sustainable airport infrastructure in China can significantly contribute to worldwide mitigation efforts.

Thermodynamic analysis of the GASright transmembrane motif supports energetic model of dimerization.

The GASright motif, best known as the fold of the glycophorin A transmembrane dimer, is one of the most common dimerization motifs in membrane proteins, characterized by its hallmark GxxxG-like sequence motifs (GxxxG, AxxxG, GxxxS, and similar). Structurally, GASright displays a right-handed crossing angle and short interhelical distance. Contact between the helical backbones favors the formation of networks of weak hydrogen bonds between Cα-H carbon donors and carbonyl acceptors on opposing helices (Cα-H···O=C). To understand the factors that modulate the stability of GASright, we previously presented a computational and experimental structure-based analysis of 26 predicted dimers. We found that the contributions of van der Waals packing and Cα-H hydrogen bonding to stability, as inferred from the structural models, correlated well with relative dimerization propensities estimated experimentally with the in vivo assay TOXCAT. Here we test this model with a quantitative thermodynamic analysis. We used Förster resonance energy transfer (FRET) to determine the free energy of dimerization of a representative subset of seven of the 26 original TOXCAT dimers using FRET. To overcome the technical issue arising from limited sampling of the dimerization isotherm, we introduced a globally fitting strategy across a set of constructs comprising a wide range of stabilities. This strategy yielded precise thermodynamic data that show strikingly good agreement between the original propensities and ΔG° of association in detergent, suggesting that TOXCAT is a thermodynamically driven process. From the correlation between TOXCAT and thermodynamic stability, the predicted free energy for all the 26 GASright dimers was calculated. These energies correlate with the in silico ΔE scores of dimerization that were computed on the basis of their predicted structure. These findings corroborate our original model with quantitative thermodynamic evidence, strengthening the hypothesis that van der Waals and Cα-H hydrogen bond interactions are the key modulators of GASright stability.

The coiled-coil domain of Escherichia coli FtsLB is a structurally detuned element critical for modulating its activation in bacterial cell division.

The FtsLB complex is a key regulator of bacterial cell division, existing in either an off state or an on state, which supports the activation of septal peptidoglycan synthesis. In Escherichia coli, residues known to be critical for this activation are located in a region near the C-terminal end of the periplasmic coiled-coil domain of FtsLB, raising questions about the precise role of this conserved domain in the activation mechanism. Here, we investigate an unusual cluster of polar amino acids found within the core of the FtsLB coiled coil. We hypothesized that these amino acids likely reduce the structural stability of the domain and thus may be important for governing conformational changes. We found that mutating these positions to hydrophobic residues increased the thermal stability of FtsLB but caused cell division defects, suggesting that the coiled-coil domain is a "detuned" structural element. In addition, we identified suppressor mutations within the polar cluster, indicating that the precise identity of the polar amino acids is important for fine-tuning the structural balance between the off and on states. We propose a revised structural model of the tetrameric FtsLB (named the "Y-model") in which the periplasmic domain splits into a pair of coiled-coil branches. In this configuration, the hydrophilic terminal moieties of the polar amino acids remain more favorably exposed to water than in the original four-helix bundle model ("I-model"). We propose that a shift in this architecture, dependent on its marginal stability, is involved in activating the FtsLB complex and triggering septal cell wall reconstruction.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA.

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

Untangling the Interactions between Anionic Polystyrene Nanoparticles and Lipid Membranes Using Laurdan Fluorescence Spectroscopy and Molecular Simulations.

Several classes of synthetic nanoparticles (NPs) induce rearrangements of cell membranes that can affect membrane function. This paper describes the investigation of the interactions between polystyrene nanoparticles and liposomes, which serve as model cell membranes, using a combination of laurdan fluorescence spectroscopy and coarse-grained molecular dynamics (MD) simulations. The relative intensities of the gel-like and fluid fluorescent peaks of laurdan, which is embedded in the liposome membranes, are quantified from the areas of deconvoluted lognormal laurdan fluorescence peaks. This provides significant advantages in understanding polymer-membrane interactions. Our study reveals that anionic polystyrene NPs, which are not cross-linked, induce significant membrane rearrangement compared to other cationic or anionic NPs. Coarse-grained MD simulations demonstrate that polymer chains from the anionic polystyrene NP penetrate the liposome membrane. The inner leaflet remains intact throughout this process, though both leaflets show a decrease in lipid packing that is indicative of significant local rearrangement of the liposome membrane. These results are attributed to the formation of a hybrid gel made up of a combination of polystyrene (PS) and lipids that forces water molecules away from laurdan. Our study concludes that a combination of negative surface charge to interact electrostatically with positive charges on the membrane, a hydrophobic core to provide a thermodynamic preference for membrane association, and the ability to extend non-cross linked polymer chains into the liposome membrane are necessary for NPs to cause a significant rearrangement in the liposomes.

Pectin mediates the mechanism of host blood glucose regulation through intestinal flora.

Pectin is a complex polysaccharide found in plant cell walls and interlayers. As a food component, pectin is benefit for regulating intestinal flora. Metabolites of intestinal flora, including short-chain fatty acids (SCFAs), bile acids (BAs) and lipopolysaccharides (LPS), are involved in blood glucose regulation. SCFAs promote insulin synthesis through the intestine-GPCRs-derived pathway and hepatic adenosine 5'-monophosphate (AMP)-activated protein kinase (AMPK) pathway to promote hepatic glycogen synthesis. On the one hand, BAs stimulate intestinal L cells and pancreatic α cells to secrete Glucagon-like peptide-1 (GLP-1) and peptide YY (PYY) through receptors G protein-coupled receptor (TGR5) and farnesoid X receptor (FXR). On the other hand, BAs promote hepatic glycogen synthesis through AMPK pathway. LPS inhibits the release of inflammatory cytokines through Toll-like receptors (TLRs)-myeloid differentiation factor 88 (MYD88) pathway and mitogen-activated protein kinase (MAPK) pathway, thereby alleviating insulin resistance (IR). In brief, both SCFAs and BAs promote GLP-1 secretion through different pathways, employing strategies of increasing glucose consumption and decreasing glucose production to maintain normal glucose levels. Notably, pectin can also directly inhibit the release of inflammatory cytokines through the -TLRs-MYD88 pathway. These data provide valuable information for further elucidating the relationship between pectin-intestinal flora-glucose metabolism.

Anionic Phospholipids Stabilize RecA Filament Bundles in Escherichia coli.

We characterize the interaction of RecA with membranes in vivo and in vitro and demonstrate that RecA binds tightly to the anionic phospholipids cardiolipin (CL) and phosphatidylglycerol (PG). Using computational models, we identify two regions of RecA that interact with PG and CL: (1) the N-terminal helix and (2) loop L2. Mutating these regions decreased the affinity of RecA to PG and CL in vitro. Using 3D super-resolution microscopy, we demonstrate that depleting Escherichia coli PG and CL altered the localization of RecA foci and hindered the formation of RecA filament bundles. Consequently, E. coli cells lacking aPLs fail to initiate a robust SOS response after DNA damage, indicating that the membrane acts as a scaffold for nucleating the formation of RecA filament bundles and plays an important role in the SOS response.

Reparameterization of the chemical-potential equalization model with DFTB3: A practical balance between accuracy and transferability.

To improve the performance of the third-order density-functional tight-binding method (DFTB3) for non-covalent interactions involving organic and biological molecules, a chemical-potential equalization (CPE) approach was introduced [J. Phys. Chem. A, 116, 9131 (2012)] and parameterized for the H, C, N, O, and S chemical elements [J. Chem. Phys., 143, 084123 (2015)]. Based largely on equilibrium structures, the parameterized DFTB3/CPE models were shown to exhibit improvements in molecular polarizabilities and intermolecular interactions. With more extensive analyses, however, we observe here that the available DFTB3/CPE models have two critical limitations: (1) they lead to sharply varying potential energy surfaces, thus causing numerical instability in molecular dynamics (MD) simulations, and (2) they lead to spurious interactions at short distances for some dimer complexes. These shortcomings are attributed to the employed screening functions and the overfitting of CPE parameters. In this work, we introduce a new strategy to simplify the parameterization procedure and significantly reduce free parameters down to four global (i.e., independent of element type) ones. With this strategy, two new models, DFTB3/CPE(r) and DFTB3/CPE(r†) are parameterized. The new models lead to smooth potential energy surfaces, stable MD simulations, and alleviate the spurious interactions at short distances, thus representing consistent improvements for both neutral and ionic hydrogen bonds.

Functional plasticity and evolutionary adaptation of allosteric regulation.

Allostery is a fundamental regulatory mechanism of protein function. Despite notable advances, understanding the molecular determinants of allostery remains an elusive goal. Our current knowledge of allostery is principally shaped by a structure-centric view, which makes it difficult to understand the decentralized character of allostery. We present a function-centric approach using deep mutational scanning to elucidate the molecular basis and underlying functional landscape of allostery. We show that allosteric signaling exhibits a high degree of functional plasticity and redundancy through myriad mutational pathways. Residues critical for allosteric signaling are surprisingly poorly conserved while those required for structural integrity are highly conserved, suggesting evolutionary pressure to preserve fold over function. Our results suggest multiple solutions to the thermodynamic conditions of cooperativity, in contrast to the common view of a finely tuned allosteric residue network maintained under selection.

Anionic nanoparticle-induced perturbation to phospholipid membranes affects ion channel function.

Understanding the mechanisms of nanoparticle interaction with cell membranes is essential for designing materials for applications such as bioimaging and drug delivery, as well as for assessing engineered nanomaterial safety. Much attention has focused on nanoparticles that bind strongly to biological membranes or induce membrane damage, leading to adverse impacts on cells. More subtle effects on membrane function mediated via changes in biophysical properties of the phospholipid bilayer have received little study. Here, we combine electrophysiology measurements, infrared spectroscopy, and molecular dynamics simulations to obtain insight into a mode of nanoparticle-mediated modulation of membrane protein function that was previously only hinted at in prior work. Electrophysiology measurements on gramicidin A (gA) ion channels embedded in planar suspended lipid bilayers demonstrate that anionic gold nanoparticles (AuNPs) reduce channel activity and extend channel lifetimes without disrupting membrane integrity, in a manner consistent with changes in membrane mechanical properties. Vibrational spectroscopy indicates that AuNP interaction with the bilayer does not perturb the conformation of membrane-embedded gA. Molecular dynamics simulations reinforce the experimental findings, showing that anionic AuNPs do not directly interact with embedded gA channels but perturb the local properties of lipid bilayers. Our results are most consistent with a mechanism in which anionic AuNPs disrupt ion channel function in an indirect manner by altering the mechanical properties of the surrounding bilayer. Alteration of membrane mechanical properties represents a potentially important mechanism by which nanoparticles induce biological effects, as the function of many embedded membrane proteins depends on phospholipid bilayer biophysical properties.

Direct readout of heterochromatic H3K9me3 regulates DNMT1-mediated maintenance DNA methylation.

In mammals, repressive histone modifications such as trimethylation of histone H3 Lys9 (H3K9me3), frequently coexist with DNA methylation, producing a more stable and silenced chromatin state. However, it remains elusive how these epigenetic modifications crosstalk. Here, through structural and biochemical characterizations, we identified the replication foci targeting sequence (RFTS) domain of maintenance DNA methyltransferase DNMT1, a module known to bind the ubiquitylated H3 (H3Ub), as a specific reader for H3K9me3/H3Ub, with the recognition mode distinct from the typical trimethyl-lysine reader. Disruption of the interaction between RFTS and the H3K9me3Ub affects the localization of DNMT1 in stem cells and profoundly impairs the global DNA methylation and genomic stability. Together, this study reveals a previously unappreciated pathway through which H3K9me3 directly reinforces DNMT1-mediated maintenance DNA methylation.