The full activation mechanism of the adenosine A1 receptor revealed by GaMD and Su-GaMD simulations.

The full activation process of G protein-coupled receptor (GPCR) plays an important role in cellular signal transduction. However, it remains challenging to simulate the whole process in which the GPCR is recognized and activated by a ligand and then couples to the G protein on a reasonable simulation timescale. Here, we developed a molecular dynamics (MD) approach named supervised (Su) Gaussian accelerated MD (GaMD) by incorporating a tabu-like supervision algorithm into a standard GaMD simulation. By using this Su-GaMD method, from the active and inactive structure of adenosine A1 receptor (A1R), we successfully revealed the full activation mechanism of A1R, including adenosine (Ado)-A1R recognition, preactivation of A1R, and A1R-G protein recognition, in hundreds of nanoseconds of simulations. The binding of Ado to the extracellular side of A1R initiates conformational changes and the preactivation of A1R. In turn, the binding of Gi2 to the intracellular side of A1R causes a decrease in the volume of the extracellular orthosteric site and stabilizes the binding of Ado to A1R. Su-GaMD could be a useful tool to reconstruct or even predict ligand-protein and protein-protein recognition pathways on a short timescale. The intermediate states revealed in this study could provide more detailed complementary structural characterizations to facilitate the drug design of A1R in the future.

Mechanism of the Conformational Change of the Protein Methyltransferase SMYD3: A Molecular Dynamics Simulation Study.

SMYD3 is a SET-domain-containing methyltransferase that catalyzes the transfer of methyl groups onto lysine residues of substrate proteins. Methylation of MAP3K2 by SMYD3 has been implicated in Ras-driven tumorigenesis, which makes SMYD3 a potential target for cancer therapy. Of all SMYD family proteins, SMYD3 adopt a closed conformation in a crystal structure. Several studies have suggested that the conformational changes between the open and closed forms may regulate the catalytic activity of SMYD3. In this work, we carried out extensive molecular dynamics simulations on a series of complexes with a total of 21 µs sampling to investigate the conformational changes of SMYD3 and unveil the molecular mechanisms. Based on the C-terminal domain movements, the simulated models could be depicted in three different conformational states: the closed, intermediate and open states. Only in the case that both the methyl donor binding pocket and the target lysine-binding channel had bound species did the simulations show SMYD3 maintaining its conformation in the closed state, indicative of a synergetic effect of the cofactors and target lysine on regulating the conformational change of SMYD3. In addition, we performed analyses in terms of structure and energy to shed light on how the two regions might regulate the C-terminal domain movement. This mechanistic study provided insights into the relationship between the conformational change and the methyltransferase activity of SMYD3. The more complete understanding of the conformational dynamics developed here together with further work may lay a foundation for the rational drug design of SMYD3 inhibitors.

Comprehensive Insights into the Catalytic Mechanism of Middle East Respiratory Syndrome 3C-Like Protease and Severe Acute Respiratory Syndrome 3C-Like Protease.

Coronavirus 3C-like protease (3CLPro) is a highly conserved cysteine protease employing a catalytic dyad for its functions. 3CLPro is essential to the viral life cycle and, therefore, is an attractive target for developing antiviral agents. However, the detailed catalytic mechanism of coronavirus 3CLPro remains largely unknown. We took an integrated approach of employing X-ray crystallography, mutational studies, enzyme kinetics study, and inhibitors to gain insights into the mechanism. Such experimental work is supplemented by computational studies, including the prereaction state analysis, the ab initio calculation of the critical catalytic step, and the molecular dynamic simulation of the wild-type and mutant enzymes. Taken together, such studies allowed us to identify a residue pair (Glu-His) and a conserved His as critical for binding; a conserved GSCGS motif as important for the start of catalysis, a partial negative charge cluster (PNCC) formed by Arg-Tyr-Asp as essential for catalysis, and a conserved water molecule mediating the remote interaction between PNCC and catalytic dyad. The data collected and our insights into the detailed mechanism have allowed us to achieve a good understanding of the difference in catalytic efficiency between 3CLPro from SARS and MERS, conduct mutational studies to improve the catalytic activity by 8-fold, optimize existing inhibitors to improve the potency by 4-fold, and identify a potential allosteric site for inhibitor design. All such results reinforce each other to support the overall catalytic mechanism proposed herein.

Crystal Structure of African Swine Fever Virus pS273R Protease and Implications for Inhibitor Design.

African swine fever (ASF) is a highly contagious hemorrhagic viral disease of domestic and wild pigs that is responsible for serious economic and production losses. It is caused by the African swine fever virus (ASFV), a large and complex icosahedral DNA virus of the Asfarviridae family. Currently, there is no effective treatment or approved vaccine against the ASFV. pS273R, a specific SUMO-1 cysteine protease, catalyzes the maturation of the pp220 and pp62 polyprotein precursors into core-shell proteins. Here, we present the crystal structure of the ASFV pS273R protease at a resolution of 2.3 Å. The overall structure of the pS273R protease is represented by two domains named the "core domain" and the N-terminal "arm domain." The "arm domain" contains the residues from M1 to N83, and the "core domain" contains the residues from N84 to A273. A structure analysis reveals that the "core domain" shares a high degree of structural similarity with chlamydial deubiquitinating enzyme, sentrin-specific protease, and adenovirus protease, while the "arm domain" is unique to ASFV. Further, experiments indicated that the "arm domain" plays an important role in maintaining the enzyme activity of ASFV pS273R. Moreover, based on the structural information of pS273R, we designed and synthesized several peptidomimetic aldehyde compounds at a submolar 50% inhibitory concentration, which paves the way for the design of inhibitors to target this severe pathogen.IMPORTANCE African swine fever virus, a large and complex icosahedral DNA virus, causes a deadly infection in domestic pigs. In addition to Africa and Europe, countries in Asia, including China, Vietnam, and Mongolia, were negatively affected by the hazards posed by ASFV outbreaks in 2018 and 2019, at which time more than 30 million pigs were culled. Until now, there has been no vaccine for protection against ASFV infection or effective treatments to cure ASF. Here, we solved the high-resolution crystal structure of the ASFV pS273R protease. The pS273R protease has a two-domain structure that distinguishes it from other members of the SUMO protease family, while the unique "arm domain" has been proven to be essential for its hydrolytic activity. Moreover, the peptidomimetic aldehyde compounds designed to target the substrate binding pocket exert prominent inhibitory effects and can thus be used in a potential lead for anti-ASFV drug development.

Exploration of the Substrate Preference of Lysine Methyltransferase SMYD3 by Molecular Dynamics Simulations.

SMYD3, a SET and MYND domain containing lysine methyltransferase, catalyzes the transfer of the methyl group from a methyl donor onto the Nε group of a lysine residue in the substrate protein. Methylation of MAP3 kinase kinase (MAP3K2) by SMYD3 has been implicated in Ras-driven tumorigenesis. The crystal structure of SMYD3 in complex with MAP3K2 peptide reveals a shallow hydrophobic pocket (P-2), which accommodates the binding of a phenylalanine residue at the -2 position of the substrate (F258) is a crucial determinant of substrate specificity of SMYD3. To better understand the substrate preference of SMYD3 at the -2 position, molecular dynamics (MD) simulations and the MM/GBSA method were performed on the crystal structure of SMYD3-MAP3K2 complex (PDB: 5EX0) after substitution of F258 residue of MAP3K2 to each of the other 19 natural residues, respectively. Binding free energy calculations reveal that the P-2 pocket prefers an aromatic hydrophobic group and none of the substitutions behave better than the wild-type phenylalanine residue does. Furthermore, we investigated the structure-activity relationships (SAR) of a series of non-natural phenylalanine derivative substitutions at the -2 position and found that quite a few modifications on the sidechain of F258 residue could strengthen its binding to the P-2 pocket of SMYD3. These explorations provide insights into developing novel SMYD3 inhibitors with high potency and high selectivity against MAP3K2 and cancer.

Computer-aided design of short peptide ligands targeting tumor necrosis factor-alpha for adsorbent applications.

Tumor necrosis factor alpha (TNF-α) is a pro-inflammatory cytokine active in the bodily immune response and serious inflammatory diseases. Traditional ligands targeting TNF-α focus on antibodies and receptors, which always associate with low efficacy and specificity. In the present study, two peptide ligands (T1: Ac-RKEM-NH2 and T2: Ac-RHCLS-NH2) were designed by computer simulation technology considering the weak interactions between TNF-α and its receptor TNFR1. Calculations of binding free energy (BFE) were made by the Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) method between T1 or T2 and TNF-α (-22.68 and -14.23 kcal mol-1, respectively). To assess the affinity levels, short peptide ligands were fixed on polyvinyl alcohol (PVA) microspheres; adsorption tests showed a stronger affinity of both PVA-T1 and PVA-T2 to TNF-α in PBS buffer than PVA microspheres (79.20 ± 1.32 and 74.27 ± 1.10 vs. 39.03 ± 1.25 pg mg-1, respectively). Moreover, PVA-T1 (74.8%, 17.60 ± 2.98 pg mg-1) and PVA-T2 (63.2%, 15.30 ± 4.81 pg mg-1) exhibit significantly enhanced TNF-α adsorption from the plasma of rats with sepsis to blank PVA and commercial XAD-7 resin. In conclusion, our results show that T1 designed by computer-aided molecular design (CAMD) exhibits a stronger affinity to TNF-α and it can significantly enhance PVA microsphere adsorption efficiency of TNF-α in plasma.

Study of the mechanism of protonated histidine-induced conformational changes in the Zika virus dimeric envelope protein using accelerated molecular dynamic simulations.

The Zika virus has drawn worldwide attention because of the epidemic diseases it causes. It is a flavivirus that has an icosahedral protein shell constituted by an envelope glycoprotein (E-protein) and membrane protein (M-protein) in the mature virion. The multistep process of membrane fusion to infect the host cell is pH-induced. To understand the mechanism of the conformational changes in the (E-M)2 protein homodimer embedded in the membrane, two 200-ns accelerated dynamic simulations were performed under different pH conditions. The low pH condition weakens the interactions and correlations in both E-protein monomers and in the E-M heterodimer. The highly conserved residues, His249, His288, His323 and His446, are protonated under low pH conditions and play key roles in driving the fusion process. The analysis and discussion in this study may provide some insight into the molecular mechanism of Zika virus infection.

A non-canonical binding interface in the crystal structure of HIV-1 gp120 core in complex with CD4.

Numerous crystal structures of HIV gp120 have been reported, alone or with receptor CD4 and cognate antibodies; however, no sole gp120/CD4 complex without stabilization by an antibody is available. Here, we report a crystal structure of the gp120/CD4 complex without the aid of an antibody from HIV-1 CRF07_BC, a strain circulating in China. Interestingly, in addition to the canonical binding surface, a second interacting interface was identified. A mutagenesis study on critical residues revealed that the stability of this interface is important for the efficiency of Env-mediated membrane fusion. Furthermore, we found that a broad neutralizing antibody, ibalizumab, which targets CD4 in the absence of gp120, occupies the same binding surface as the second interface identified here on gp120. Therefore, we identified the possibility of the involvement of a second gp120-CD4 interaction interface during viral entry, and also provided a reasonable explanation for the broad activity of neutralizing antibody ibalizumab.

Biochemical characterization of recombinant Enterovirus 71 3C protease with fluorogenic model peptide substrates and development of a biochemical assay.

Enterovirus 71 (EV71), a primary pathogen of hand, foot, and mouth disease (HFMD), affects primarily infants and children. Currently, there are no effective drugs against HFMD. EV71 3C protease performs multiple tasks in the viral replication, which makes it an ideal antiviral target. We synthesized a small set of fluorogenic model peptides derived from cleavage sites of EV71 polyprotein and examined their efficiencies of cleavage by EV71 3C protease. The novel peptide P08 [(2-(N-methylamino)benzoyl) (NMA)-IEALFQGPPK(DNP)FR] was determined to be the most efficiently cleaved by EV71 3C protease, with a kinetic constant kcat/Km of 11.8 ± 0.82 mM(-1) min(-1). Compared with literature reports, P08 gave significant improvement in the signal/background ratio, which makes it an attractive substrate for assay development. A Molecular dynamics simulation study elaborated the interactions between substrate P08 and EV71 3C protease. Arg39, which is located at the bottom of the S2 pocket of EV71 3C protease, may participate in the proteolysis process of substrates. With an aim to evaluate EV71 3C protease inhibitors, a reliable and robust biochemical assay with a Z' factor of 0.87 ± 0.05 was developed. A novel compound (compound 3) (50% inhibitory concentration [IC50] = 1.89 ± 0.25 µM) was discovered using this assay, which effectively suppressed the proliferation of EV 71 (strain Fuyang) in rhabdomyosarcoma (RD) cells with a highly selective index (50% effective concentration [EC50] = 4.54 ± 0.51 µM; 50% cytotoxic concentration [CC50] > 100 µM). This fast and efficient assay for lead discovery and optimization provides an ideal platform for anti-EV71 drug development targeting 3C protease.