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

Crystal structure and mechanistic studies of the PPM1D serine/threonine phosphatase catalytic domain

Protein phosphatase 1D (PPM1D, Wip1) is induced by the tumor suppressor p53 during DNA damage response signaling and acts as an oncoprotein in several human cancers. Although PPM1D is a potential therapeutic target, insights into its atomic structure were challenging due to flexible regions unique to this family member. Here, we report the first crystal structure of the PPM1D catalytic domain to 1.8 Å resolution. The structure reveals the active site with two Mg2+ ions bound, similar to other structures. The flap subdomain and B-loop, which are crucial for substrate recognition and catalysis, were also resolved, with the flap forming two short helices and three short β-strands that are followed by an irregular loop. Unexpectedly, a nitrogen-oxygen-sulfur bridge was identified in the catalytic domain. Molecular dynamics simulations and kinetic studies provided further mechanistic insights into the regulation of PPM1D catalytic activity. In particular, the kinetic experiments demonstrated a magnesium concentration-dependent lag in PPM1D attaining steady-state velocity, a feature of hysteretic enzymes that show slow transitions compared with catalytic turnover. All combined, these results advance the understanding of PPM1D function and will support the development of PPM1D-targeted therapeutics.

Engineering Human Pancreatic RNase 1 as an Immunotherapeutic Agent for Cancer Therapy Through Computational and Experimental Studies

Most plant and bacterial toxins are highly immunogenic with non-specific toxic effects. Human ribonucleases are thought to provide a promising basis for reducing the toxic agent's immunogenic properties, which are candidates for cancer therapy. In the cell, the ribonuclease inhibitor (RI) protein binds to the ribonuclease enzyme and forms a tight complex. This study aimed to engineer and provide a gene construct encoding an improved version of Human Pancreatic RNase 1 (HP-RNase 1) to reduce connection to RI and modulate the immunogenic effects of immunotoxins. To further characterize the interaction complex of HP-RNase 1 and RI, we established various in silico and in vitro approaches. These methods allowed us to specifically monitor interactions within native and engineered HP-RNase 1/RI complexes. In silico research involved molecular dynamics (MD) simulations of native and mutant HP-RNase 1 in their free form and when bound to RI. For HP-RNase 1 engineering, we designed five mutations (K8A/N72A/N89A/R92D/E112/A) based on literature studies, as this combination proved effective for the intended investigation. Then, the cDNA encoding HP-RNase 1 was generated by RT-PCR from blood and cloned into the pSYN2 expression vector. Consequently, wild-type and the engineered HP-RNase 1 were over-expressed in E. coli TG1 and purified using an IMAC column directed against a poly-his tag. The protein products were detected by SDS-PAGE and Western blot analysis. HP-RNase 1 catalytic activity, in the presence of various concentrations of RI, demonstrated that the mutated version of the protein is able to escape the ribonuclease inhibitor and target the RNA substrate 2.5 folds more than that of the wild type. From these data, we tend to suggest the engineered recombinant HP-RNase 1 potentially as a new immunotherapeutic agent for application in human cancer therapy.

Evaluation of mutations on O6-methylguanine methyl transferase structure and its interactions: molecular dynamics simulation study

O6-methylguanine DNA methyl transferase (MGMT) is a significant vehicle for the cellular clearance of alkyl lesions, particularly the methyl group of the O-6 and O-4 positions of guanine and thymine, respectively. Many publications have studied the correlation between polymorphisms in MGMT and susceptibility to various cancers. In the present study, we investigated the consequence of L84F, common single-nucleotide polymorphism, K125E, site-specific mutagenesis, and L84F/K125E on conformation, stability, and behavior of MGMT in the free form and interaction with proliferating cell nuclear antigen (PCNA) and DNA as partners in the biochemical network by using molecular dynamics simulation method. Our results showed that all free variants of MGMT differed from the native form. However, among all free variants of MGMT, the L84F/K125E variant exhibited similar properties compared with the wild-type. In contrast, in complex modes, only amino acid residues of the L84F variant are involved in the interactions with PCNA and DNA somewhat differently relative to the wild-type. Furthermore, L84F SNP showed the highest binding free energy compared to other variants and native forms. These alterations in the amino acids and binding free energy of L84F relative to the native are the reasons for changing its region connection compared to the native form. Therefore, we propose conducting further investigations into the impact of inhibitors or chemotherapeutic agents to assess their effectiveness on MGMT variants compared to the wild-type, aiming to reduce the cost of cancer treatment that will depend on inhibiting native MGMT protein.Communicated by Ramaswamy H. Sarma.

Residues from Homologous Transmembrane Helices 4 and 10 Are Critical for P-Glycoprotein (ABCB1)-Mediated Drug Transport

P-glycoprotein (P-gp, ABCB1) transports structurally dissimilar hydrophobic and amphipathic compounds, including anticancer drugs, thus contributing to multidrug-resistant cancer. Cryo-EM structures of human P-gp revealed that TMHs 4 and 10 contribute to the formation of the drug-binding cavity and undergo conformational changes during drug transport. To assess the role of the conformational changes in TMH4 and TMH10 during drug transport, we generated two mutants (TMH4-7A and TMH10-7A), each containing seven alanine substitutions. Analysis of the drug efflux function of these mutants using 15 fluorescent substrates revealed that most of the substrates were transported, indicating that even seven mutations in an individual helix have no significant effect on transport function. We then designed the TMH4,10-14A mutant combining seven mutations in both TMHs 4 and 10. Interestingly, when the TMH4,10-14A mutant was tested with 15 substrates, there was no efflux observed for fourteen. The basal ATPase activity of the TMH4,10-14A mutant, similar to that of the WT protein, was inhibited by zosuquidar but was not stimulated by verapamil or rhodamine 6G. Molecular dynamics simulations indicated that the mutations cause TMHs 4 and 10 to pack tighter to their proximal helices, reducing their independent mobility. In aggregate, our findings demonstrate the critical role of the residues of homologous TMHs 4 and 10 for substrate transport, consistent with conformational changes observed in the structure of P-gp.

Peonidin-3-O-Glucoside from Purple Corncob Ameliorates Nonalcoholic Fatty Liver Disease by Regulating Mitochondrial and Lysosome Functions to Reduce Oxidative Stress and Inflammation

A frequent chronic liver condition across the world is nonalcoholic fatty liver disease (NAFLD). Oxidative stress caused by lipid accumulation is generally considered to be the main cause of NAFLD. Anthocyanins can effectively inhibit the production of reactive oxygen species and improve oxidative stress. In this work, six major anthocyanins were separated from purple corncob by semi-preparative liquid chromatography. The effects of the 6 kinds of anthocyanins against NAFLD were investigated using a free fatty acid (FFA)-induced cell model. The results showed that peonidin 3-O-glucoside (P3G) can significantly reduce lipid accumulation in the NAFLD cell model. The treatment with P3G also inhibited oxidative stress via inhibiting the excessive production of reactive oxygen species and superoxide anion, increasing glutathione levels, and enhancing the activities of SOD, GPX, and CAT. Further studies unveiled that treatment with P3G not only alleviated inflammation but also improved the depletion of mitochondrial content and damage of the mitochondrial electron transfer chain developed concomitantly in the cell model. P3G upregulated transcription factor EB (TFEB)-mediated lysosomal function and activated the peroxisome proliferator-activated receptor alpha (PPARα)-mediated peroxisomal lipid oxidation by interacting with PPARα possibly. Overall, this study added to our understanding of the protective effects of purple corn anthocyanins against NAFLD and offered suggestions for developing functional foods containing these anthocyanins.

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 (3CL^pro), 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 3CL^pro 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 (3CL^pro) and P3_Val (peptide), followed by binding to the oxyanion hole and completed by the sequence-specific recognition at P1_Gln.

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

The 3C-like protease (3CL^pro) 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 3CL^pro of SARS-CoV (101 PDB entries), we constructed the crystal structure ensemble of 3CL^pro to analyze the dynamic regulation of its catalytic function. The structural dynamics of the 3CL^pro 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 3CL^pro to SARS-CoV-2 3CL^pro 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 3CL^pro seems to have increased activity relative to that of SARS-CoV 3CL^pro.

Multiscale Enhanced Sampling Using Machine Learning

Multiscale enhanced sampling (MSES) allows for an enhanced sampling of all-atom protein structures by coupling with the accelerated dynamics of the associated coarse-grained (CG) model. In this paper, we propose an MSES extension to replace the CG model with the dynamics on the reduced subspace generated by a machine learning approach, the variational autoencoder (VAE). The molecular dynamic (MD) trajectories of the ribose-binding protein (RBP) in both the closed and open forms were used as the input by extracting the inter-residue distances as the structural features in order to train the VAE model, allowing the encoded latent layer to characterize the difference in the structural dynamics of the closed and open forms. The interpolated data characterizing the RBP structural change in between the closed and open forms were thus efficiently generated in the low-dimensional latent space of the VAE, which was then decoded into the time-series data of the inter-residue distances and was useful for driving the structural sampling at an atomistic resolution via the MSES scheme. The free energy surfaces on the latent space demonstrated the refinement of the generated data that had a single basin into the simulated data containing two closed and open basins, thus illustrating the usefulness of the MD simulation together with the molecular mechanics force field in recovering the correct structural ensemble.

Coarse-Grained Modeling of Multiple Pathways in Conformational Transitions of Multi-Domain Proteins

Large-scale conformational transitions in multi-domain proteins are often essential for their functions. To investigate the transitions, it is necessary to explore multiple potential pathways, which involve different intermediate structures. Here, we present a multi-basin (MB) coarse-grained (CG) structure-based Go̅ model for describing transitions in proteins with more than two moving domains. This model is an extension of our dual-basin Go̅ model in which system-dependent parameters are determined systematically using the multistate Bennett acceptance ratio method. In the MB Go̅ model for multi-domain proteins, we assume that intermediate structures may have partial inter-domain native contacts. This approach allows us to search multiple transition pathways that involve distinct intermediate structures using the CG molecular dynamics (MD) simulations. We apply this scheme to an enzyme, adenylate kinase (AdK), which has three major domains and can move along two different pathways. Using the optimized mixing parameters for each pathway, AdK shows frequent transitions between the Open, Closed, and the intermediate basins and samples a wide variety of conformations within each basin. The explored multiple transition pathways could be compared with experimental data and examined in more detail by atomistic MD simulations.

Exploring Large Domain Motions in Proteins Using Atomistic Molecular Dynamics with Enhanced Conformational Sampling

Conformational transitions in multidomain proteins are essential for biological functions. The Apo conformations are typically open and flexible, while the Holo states form more compact conformations stabilized by protein-ligand interactions. Unfortunately, the atomically detailed mechanisms for such open-closed conformational changes are difficult to be accessed experimentally as well as computationally. To simulate the transitions using atomistic molecular dynamics (MD) simulations, efficient conformational sampling algorithms are required. In this work, we propose a new approach based on generalized replica-exchange with solute tempering (gREST) for exploring the open-closed conformational changes in multidomain proteins. Wherein, selected surface charged residues in a target protein are defined as the solute region in gREST simulation and the solute temperatures are different in replicas and exchanged between them to enhance the domain motions. This approach is called gREST selected surface charged residues (gREST_SSCR) and is applied to the Apo and Holo states of ribose binding protein (RBP) in solution. The conformational spaces sampled with gREST_SSCR are much wider than those with the conventional MD, sampling open-closed conformational changes while maintaining RBP domains’ stability. The free-energy landscapes of RBP in the Apo and Holo states are drawn along with twist and hinge angles of the two moving domains. The inter-domain salt-bridges that are not observed in the experimental structures are also important in the intermediate states during the conformational changes.

Building a macro-mixing dual-basin Gō model using the Multistate Bennett Acceptance Ratio

The dual-basin Gō-model is a structural-based coarsegrained model for simulating a conformational transition between two known structures of a protein. Two parameters are required to produce a dual-basin potential mixed using two single-basin potentials, although the determination of mixing parameters is usually not straightforward. Here, we have developed an efficient scheme to determine the mixing parameters using the Multistate Bennett Acceptance Ratio (MBAR) method after short simulations with a set of parameters. In the scheme, MBAR allows us to predict observables at various unsimulated conditions, which are useful to improve the mixing parameters in the next round of iterative simulations. The number of iterations that are necessary for obtaining the converged mixing parameters are significantly reduced in the scheme. We applied the scheme to two proteins, the glutamine binding protein and the ribose binding protein, for showing the effectiveness in the parameter determination. After obtaining the converged parameters, both proteins show frequent conformational transitions between open and closed states, providing the theoretical basis to investigate structure-dynamics-function relationships of the proteins.

Analysis of molecular dynamics simulations of 10-residue peptide, chignolin, using statistical mechanics: Relaxation mode analysis and three-dimensional reference interaction site model theory

Molecular dynamics simulation is a fruitful tool for investigating the structural stability, dynamics, and functions of biopolymers at an atomic level. In recent years, simulations can be performed on time scales of the order of milliseconds using special purpose systems. Since the most stable structure, as well as meta-stable structures and intermediate structures, is included in trajectories in long simulations, it is necessary to develop analysis methods for extracting them from trajectories of simulations. For these structures, methods for evaluating the stabilities, including the solvent effect, are also needed. We have developed relaxation mode analysis to investigate dynamics and kinetics of simulations based on statistical mechanics. We have also applied the three-dimensional reference interaction site model theory to investigate stabilities with solvent effects. In this paper, we review the results for designing amino-acid substitution of the 10-residue peptide, chignolin, to stabilize the misfolded structure using these developed analysis methods.

Allosteric response to ligand binding: Molecular dynamics study of the N-terminal domains in IP 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.

All Atom Motion Tree detects side chain-related motions and their coupling with domain motion in proteins

Structural changes of proteins are closely related with their molecular function. We previously developed a computational tool, Motion Tree (MT), to compare protein structures and describe structural changes using solely the Cα atoms. Here, we have extended MT to incorporate all heavy atoms to analyze side chain-related (SCR) motions. All Atom Motion Tree (AAMT) was applied to 76 proteins that exhibited a simple domain motion identified by MT. AAMT also detected 921 SCR motions. We examined the coupling of domain and SCR motions and classified the structural changes in terms of coupling. The statistical results indicated that it is common for coupled SCR motions to also couple with the domain motion. The classification correlates properties of domain motions and SCR motions. The AAMT results suggest that a large domain motion with a sizable domain boundary is accompanied by SCR motions composed of more than a single residue, which induces further couplings of SCR motions.

Relaxation mode analysis for molecular dynamics simulations of proteins

Molecular dynamics simulation is a powerful method for investigating the structural stability, dynamics, and function of biopolymers at the atomic level. In recent years, it has become possible to perform simulations on time scales of the order of milliseconds using special hardware. However, it is necessary to derive the important factors contributing to structural change or function from the complicated movements of biopolymers obtained from long simulations. Although some analysis methods for protein systems have been developed using increasing simulation times, many of these methods are static in nature (i.e., no information on time). In recent years, dynamic analysis methods have been developed, such as the Markov state model and relaxation mode analysis (RMA), which was introduced based on spin and homopolymer systems. The RMA method approximately extracts slow relaxation modes and rates from trajectories and decomposes the structural fluctuations into slow relaxation modes, which characterize the slow relaxation dynamics of the system. Recently, this method has been applied to biomolecular systems. In this article, we review RMA and its improved versions for protein systems.