Residues remote from the binding pocket control the antagonist selectivity towards the corticotropin-releasing factor receptor-1

Optimization of degradation behavior and conditions for the protease K of polylactic acid films by simulation

With global enforcement of plastic bans and restrictions, the biodegradable plastic, e.g., polylactic acid (PLA), has been extensively employed as a primary substitute for traditional petroleum-based plastics. However, the growing problem associated with PLA waste accumulation is posing grand environmental challenges. In addition, although PLA has the degrading property under natural conditions, the degradation process takes too long and the degradation products cannot be recycled. In this context, enzymatic degradation of PLA arouses great attention in scientific communities. This study aims at selecting the most cost-effective protease from various enzymes and optimizing the enzymolysis conditions towards the degradation of PLA. We will demonstrate that under an optimal temperature of 45 °C, a pH vale of 11, and an enzyme concentration of 0.6 mg mL-1, the protease K would achieve a remarkable degradation efficiency (> 99 %) for PLA films within just 50 min. Finally, molecular dynamics (MD) simulation and molecular docking studies reveal the mechanism behind the protease-induced PLA degradation, providing a promising direction for waste treatment and resource utilization for future biodegradable plastics.

Deciphering the mechanism of action of 5FDQD and the design of new neutral analogues for the FMN riboswitch: a well-tempered metadynamics simulation study

The FMN riboswitch is a novel drug target for the design of new antibiotics, and efforts have been made to design new charged and uncharged ligands. Uncharged ligands have shown advantages of not requiring any transporter for intracellular transport or proteins for their phosphorylation. 5FDQD (5-(3-(4-fluorophenyl)butyl)-7,8-dimethylpyrido(3,4-b)quinoxaline-1,3(2H,5H)-dione) is a recently reported neutral ligand for the FMN riboswitch active against Clostridium difficile infection in mice. However, the crystal structure of the 5FDQD bound FMN riboswitch is not available, and the mechanism of ligand binding and triggering the function of the riboswitch is not well understood. We have examined 5FDQD for its binding affinity with the FMN riboswitch using the well-tempered metadynamics (WT-MtD) simulation technique. The crystal structure of the FMN riboswitch shows that the FMN interacts with the J4/5 region through the phosphate group with G62; however, the uncharged ligands take advantage of π-π stacking interactions with the same residue of the riboswitch observed from the WT-MtD simulation results. The simulation results show that the presence of fluorine on the phenyl ring in 5FDQD is important to enhance the binding affinity of the neutral ligands with the FMN riboswitch. The WT-MtD results showed that the 1,2-difluoro substitution on the phenyl ring in 5FDQD (FMN-difluoro2) and the 1,3 positions in the phenyl ring (FMN-difluoro1) showed weaker binding energy with the FMN riboswitch compared to 5FDQD. The substitution of another fluorine atom at the 5-position of the phenyl ring (FMN-trifluoro) showed a comparable binding affinity (∼-31.4 kcal mol^-1) to 5FDQD. Electron-donating substitution on the phenyl ring such as the amino group also lowered the binding affinity (-28.8 kcal mol^-1) with the riboswitch compared to 5FDQD. The computed results suggest that the position and nature of substitution in the phenyl ring of the uncharged ligands affect the overall binding and such a delicate balance is important to achieve superior binding affinity with the FMN riboswitch.

Activation Mechanism of Corticotrophin Releasing Factor Receptor Type 1 Elucidated Using Molecular Dynamics Simulations

The corticotropin-releasing factor receptor type 1 (CRF1R), a member of class B G-protein-coupled receptors (GPCRs), is a good drug target for treating depression, anxiety, and other stress-related neurodisorders. However, there is no approved drug targeting the CRF1R to date, partly due to inadequate structural information and its elusive activation mechanism. Here, by use of the crystal structures of its transmembrane domain (TMD) and the N-terminal extracellular domain (ECD) as a template, a full-length homology model of CRF1R was built and its complexes with peptide agonist urocortin 1 or small molecule antagonist CP-376395 were subjected to all-atom molecular dynamics simulations. We observed well preserved helical contents in the TMD through simulations, while the transmembrane (TM) helices showed clear rearrangements. The TM rearrangement is especially pronounced for the TM6 in the agonist-bound CRF1R system. The observed conformational changes are likely due to breakage of interhelical/inter-regional hydrogen bonds in the TMD. Dynamical network analysis identifies communities with high connections to TM6. Simulations reveal three key residues, Y356^6.53, Q384^7.49, and L395^7.60, which corroborate experimental mutagenesis data, implying the important roles in the receptor activation. The observed large-scale conformational changes are related to CRF1R activation by agonist binding, providing guidance for ligand design.

A Single Second Shell Amino Acid Determines Affinity and Kinetics of Linagliptin Binding to Type 4 Dipeptidyl Peptidase and Fibroblast Activation Protein

Drugs targeting type 4 dipeptidyl peptidase (DPP‐4) are beneficial for glycemic control, whereas fibroblast activation protein alpha (FAP‐α) is a potential target for cancer therapies. Unlike other gliptins, linagliptin displays FAP inhibition. We compared biophysical and structural characteristics of linagliptin binding to DPP‐4 and FAP to better understand what differentiates linagliptin from other gliptins. Linagliptin exhibited high binding affinity (K_D) and a slow off‐rate (k_off) when dissociating from DPP‐4 (K_D 6.6 pM; k_off 5.1×10⁻⁵ s⁻¹), and weaker inhibitory potency to FAP (K_D 301 nM; k_off>1 s⁻¹). Co‐structures of linagliptin with DPP‐4 or FAP were similar except for one second shell amino acid difference: Asp663 (DPP‐4) and Ala657 (FAP). pH dependence of enzymatic activities and binding of linagliptin for DPP‐4 and FAP are dependent on this single amino acid difference. While linagliptin may not display any anticancer activity at therapeutic doses, our findings may guide future studies for the development of optimized inhibitors.

Computational studies on horseshoe shape pocket of human orexin receptor type 2 and boat conformation of suvorexant by molecular dynamics simulations

The FDA approved drug suvorexant binds to the horseshoe shape pocket of OX₂ R with the boat conformation. The horseshoe shape pocket plays an important role on the biological activity of OX₂ R in the cell membrane. To study the binding mechanism between the horseshoe shape pocket of OX₂ R and boat conformation of suvorexant, the crystal structures of wild type and N324A mutant of OX₂ R in complex with antagonist suvorexant are chosen to perform molecular dynamics (MD) simulations, QM/MM, and MMGBSA calculations. By comparison with the wild type of OX₂ R, the results show the 1,2,3-triazole and p-toluamide groups of suvorexant are changed in the N324A mutant of OX₂ R during 200 ns MD simulations. The QM/MM and weak interaction analysis are employed to calculate the non-covalent bonds interaction between suvorexant and key residues in the wild type and N324A mutant of OX₂ R. The MMGBSA calculations indicate the entropy energy is an important influence factor for suvorexant affinity in the distorted horseshoe shape pocket of OX₂ R. Our results not only show the horseshoe shape pocket of OX₂ R is the necessary conformation for the binding of antagonist suvorexant, but also give the important sites and structural features for antagonist design of OX₂ R.

Current understanding of the structure and function of family B GPCRs to design novel drugs

Family B of G-protein-coupled receptors (GPCRs) and their ligands play a central role in a number of homeostatic mechanisms in the endocrine, gastrointestinal, skeletal, immune, cardiovascular and central nervous systems. Alterations in family B GPCR-regulated homeostatic mechanisms may cause a variety of potentially life-threatening conditions, signifying the necessity to develop novel ligands targeting these receptors. Obtaining structural and functional information on family B GPCRs will accelerate the development of novel drugs to target these receptors. Family B GPCRs are proteins that span the plasma membrane seven times, thus forming seven transmembrane domains (TM1-TM7) which are connected to each other by three extracellular (EL) and three intracellular (IL) loops. In addition, these receptors have a long extracellular N-domain and an intracellular C-tail. The upper parts of the TMs and ELs form the J-domain of receptors. The C-terminal region of peptides first binds to the N-domain of receptors. This 'first-step' interaction orients the N-terminal region of peptides towards the J-domain of receptors, thus resulting in a 'second-step' of ligand-receptor interaction that activates the receptor. Activation-associated structural changes of receptors are transmitted through TMs to their intracellular regions and are responsible for their interaction with the G proteins and activation of the latter, thus resulting in a biological effect. This review summarizes the current information regarding the structure and function of family B GPCRs and their physiological and pathophysiological roles.

Propagation of the Allosteric Modulation Induced by Sodium in the δ-Opioid Receptor

Allosteric sodium in the helix bundle of a G protein-coupled receptor (GPCR) can modulate the receptor activation on the intracellular side. This phenomenon has confounded the GPCR community for decades. In this work, we present a theoretical model that reveals the mechanism of the allosteric modulation induced by sodium in the δ-opioid receptor. We found that the allosteric sodium ion exploits a distinct conformation of the key residue Trp274^6.48 to propagate the modulation to helices 5 and 6, which further transmits along the helices and regulates their positions on the intracellular side. This mechanism is supported by subsequent functional assays. Remarkably, our results highlight the contrast between the allosteric effects towards two GPCR partners, the G protein and β-arrestin, as indicated by the fact that the allosteric modulation initiated by the sodium ion significantly affects the β-arrestin recruitment, while it alters the G protein signaling only moderately. We believe that the mechanism revealed in this work can be used to explain allosteric effects initiated by sodium in other GPCRs since the allosteric sodium is highly conserved across GPCRs.

Structure and Dynamics of PD-L1 and an Ultra-High-Affinity PD-1 Receptor Mutant

The immune checkpoint receptor PD-1 and its ligand, PD-L1, have emerged as key regulators of anti-tumor immunity in humans. Recently, we reported an ultra-high-affinity PD-1 mutant, termed high-affinity consensus (HAC) PD-1, which shows superior therapeutic efficacy in mice compared with antibodies. However, the molecular details underlying the action of this agent remain incompletely understood, and a molecular view of PD-1/PD-L1 interactions in general is only beginning to emerge. Here, we report the structure of HAC PD-1 in complex with PD-L1, showing that it binds PD-L1 using a unique set of polar interactions. Biophysical studies and long-timescale molecular dynamics experiments reveal the mechanisms by which ten point mutations confer a 35,000-fold enhancement in binding affinity, and offer atomic-scale views of the role of conformational dynamics in PD-1/PD-L1 interactions. Finally, we show that the HAC PD-1 exhibits pH-dependent affinity, with pseudo-irreversible binding in a low pH setting akin to the tumor microenvironment.

Investigation of allosteric modulation mechanism of metabotropic glutamate receptor 1 by molecular dynamics simulations, free energy and weak interaction analysis

Metabotropic glutamate receptor 1 (mGlu₁), which belongs to class C G protein-coupled receptors (GPCRs), can be coupled with G protein to transfer extracellular signal by dimerization and allosteric regulation. Unraveling the dimer packing and allosteric mechanism can be of great help for understanding specific regulatory mechanism and designing more potential negative allosteric modulator (NAM). Here, we report molecular dynamics simulation studies of the modulation mechanism of FITM on the wild type, T815M and Y805A mutants of mGlu₁ through weak interaction analysis and free energy calculation. The weak interaction analysis demonstrates that van der Waals (vdW) and hydrogen bonding play an important role on the dimer packing between six cholesterol molecules and mGlu₁ as well as the interaction between allosteric sites T815, Y805 and FITM in wild type, T815M and Y805A mutants of mGlu₁. Besides, the results of free energy calculations indicate that secondary binding pocket is mainly formed by the residues Thr748, Cys746, Lys811 and Ser735 except for FITM-bound pocket in crystal structure. Our results can not only reveal the dimer packing and allosteric regulation mechanism, but also can supply useful information for the design of potential NAM of mGlu₁.

Computational studies of the binding profile of phosphoinositide PtdIns (3,4,5) P₃ with the pleckstrin homology domain of an oomycete cellulose synthase

Saprolegnia monoica is a model organism to investigate Saprolegnia parasitica, an important oomycete which causes considerable loss in aquaculture every year. S. monoica contains cellulose synthases vital for oomycete growth. However, the molecular mechanism of the cellulose biosynthesis process in the oomycete growth is still poorly understood. Some cellulose synthases of S. monoica, such as SmCesA2, are found to contain a plecsktrin homology (PH) domain, which is a protein module widely found in nature and known to bind to phosphoinositides, a class of signaling compounds involved in many biological processes. Understanding the molecular interactions between the PH domain and phosphoinositides would help to unravel the cellulose biosynthesis process of oomycetes. In this work, the binding profile of PtdIns (3,4,5) P₃, a typical phosphoinositide, with SmCesA2-PH was studied by molecular docking, molecular dynamics and metadynamics simulations. PtdIns (3,4,5) P₃ is found to bind at a specific site located at β1, β2 and β1-β2 loop of SmCesA2-PH. The high affinity of PtdIns (3,4,5) P₃ to SmCesA2-PH is contributed by the free phosphate groups, which have electrostatic and hydrogen-bond interactions with Lys88, Lys100 and Arg102 in the binding site.

Molecular switches of the κ opioid receptor triggered by 6'-GNTI and 5'-GNTI

The κ opioid receptor (κOR) is a member of G-protein-coupled receptors, and is considered as a promising drug target for treating neurological diseases. κOR selective 6'-GNTI was proved to be a G-protein biased agonist, whereas 5'-GNTI acts as an antagonist. To investigate the molecular mechanism of how these two ligands induce different behaviors of the receptor, we built two systems containing the 5'-GNTI-κOR complex and the 6'-GNTI-κOR complex, respectively, and performed molecular dynamics simulations of the two systems. We observe that transmembrane (TM) helix 6 of the κOR rotates about 4.6(o) on average in the κOR-6'-GNTI complex. Detailed analyses of the simulation results indicate that E297(6.58) and I294(6.55) play crucial roles in the rotation of TM6. In the simulation of the κOR-5'-GNTI system, it is revealed that 5'-GNTI can stabilize TM6 in the inactive state form. In addition, the kink of TM7 is stabilized by a hydrogen bond between S324(7.47) and the residue V69(1.42) on TM1.

Josephin Domain Structural Conformations Explored by Metadynamics in Essential Coordinates

The Josephin Domain (JD), i.e. the N-terminal domain of Ataxin 3 (At3) protein, is an interesting example of competition between physiological function and aggregation risk. In fact, the fibrillogenesis of Ataxin 3, responsible for the spinocerebbellar ataxia 3, is strictly related to the JD thermodynamic stability. Whereas recent NMR studies have demonstrated that different JD conformations exist, the likelihood of JD achievable conformational states in solution is still an open issue. Marked differences in the available NMR models are located in the hairpin region, supporting the idea that JD has a flexible hairpin in dynamic equilibrium between open and closed states. In this work we have carried out an investigation on the JD conformational arrangement by means of both classical molecular dynamics (MD) and Metadynamics employing essential coordinates as collective variables. We provide a representation of the free energy landscape characterizing the transition pathway from a JD open-like structure to a closed-like conformation. Findings of our in silico study strongly point to the closed-like conformation as the most likely for a Josephin Domain in water.

Computational insight into conformational states of glucagon-like peptide-1 receptor (GLP-1R) and its binding mode with GLP-1

It was observed that the apo-GLP-1R stabilized in the ‘closed’ state, while GLP-1R coupled with GLP-1 stabilized in the ‘open’ state.

Computational investigation of the interaction mechanism between the estrogen related receptor α and its agonists

The agonists may promote the binding of coactivator PGC-1α to ERRα by stabilizing the conformation and the site of H12.

New insight into the difference of PC lipase-catalyzed degradation on poly(butylene succinate)-based copolymers from molecular levels

In the present work, the difference of enzymatic degradation of PBS-based polyesters was investigated from the molecular level with molecular modeling.

Expanding the horizons of G protein-coupled receptor structure-based ligand discovery and optimization using homology models

We show the key role of structural homology models in GPCR structure-based lead discovery and optimization, highlighting methodological aspects, recent progress and future directions.