Unraveling the Abnormal Molecular Mechanism of Suicide Inhibition of Cytochrome P450 3A4.

Suicide inhibition of the CYP3A4 enzyme by a drug inactivates the enzyme in the drug biotransformation process and often shows safety concerns about the drug. Despite extensive experimental studies, the abnormal molecular mechanism of a suicide inhibitor that forms a covalent bond with the residue far away from the catalytically active center of CYP3A4 inactivating the enzyme remains elusive. Here, the authors used molecular simulation approaches to study in detail how diquinone methide (DQR), the metabolite product of raloxifene, unbinds from CYP3A4 and inactivates the enzyme at the atomistic level. The results clearly indicate that in one of the intermediate states formed in its unbinding process, DQR covalently binds to Cys239, a residue far away from the catalytically active center of CYP3A4, and hinders the substrate from entering or leaving the enzyme. This work therefore provides an unprecedented way of clarifying the abnormal mechanism of suicide inhibition of the CYP3A4 enzyme.

Mechanistic Studies on the Stereoselectivity of FFAR1 Modulators.

Free fatty acid receptor 1 (FFAR1) is a potential therapeutic target for the treatment of type 2 diabetes (T2D). It has been validated that agonists targeting FFAR1 can achieve the initial therapeutic endpoints of T2D, and the epimer agonists (R,S) AM-8596 can activate FFAR1 differently, with one acting as a partial agonist and the other as a full agonist. Up to now, the origin of the stereoselectivity of FFAR1 agonists remains elusive. In this work, we used molecular simulation methods to elucidate the mechanism of the stereoselectivity of the FFAR1 agonists (R)-AM-8596 and (S)-AM-8596. We found that the full agonist (R)-AM-8596 disrupts the residue interaction network around the receptor binding pocket and promotes the opening of the binding site for the G-protein, thereby resulting in the full activation of FFAR1. In contrast, the partial agonist (S)-AM-8596 forms stable electrostatic interactions with FFAR1, which stabilizes the residue network and hinders the conformational transition of the receptor. Our work thus clarifies the selectivity and underlying molecular activation mechanism of FFAR1 agonists.

Homotropic Cooperativity of Midazolam Metabolism by Cytochrome P450 3A4: Insight from Computational Studies.

Human cytochrome P450 3A4 (CYP3A4) is responsible for the metabolism of ∼50% clinically used drugs. Midazolam (MDZ) is a commonly used sedative drug and serves as a marker substrate for the CYP3A4 activity assessment. MDZ is metabolized by CYP3A4 to two hydroxylation products, 1'-OH-MDZ and 4-OH-MDZ. It has been reported that the ratio of 1'-OH-MDZ and 4-OH-MDZ is dependent on the MDZ concentration, which reflects the homotropic cooperative behavior in MDZ metabolism by CYP3A4. Here, we used quantum chemistry (QC), molecular docking, conventional molecular dynamics (cMD), and Gaussian accelerated molecular dynamics (GaMD) approaches to investigate the mechanism of the interactions between CYP3A4 and MDZ. QC calculations suggest that C1' is less reactive for hydroxylation than C4, which is a pro-chirality carbon. However, the 4-OH-MDZ product is likely to be racemic due to the chirality inversion in the rebound step. The MD simulation results indicate that MDZ at the peripheral allosteric site is not stable and the binding modes of the MDZ molecules at the productive site are in line with the experimental observations.

Defect-Assisted Loading and Docking Conformations of Pharmaceuticals in Metal-Organic Frameworks.

Understanding of drug-carrier interactions is essential for the design and application of metal-organic framework (MOF)-based drug-delivery systems, and such drug-carrier interactions can be fundamentally different for MOFs with or without defects. Herein, we reveal that the defects in MOFs play a key role in the loading of many pharmaceuticals with phosphate or phosphonate groups. The host-guest interaction is dominated by the Coulombic attraction between phosphate/phosphonate groups and defect sites, and it strongly enhances the loading capacity. For similar molecules without a phosphate/phosphonate group or for MOFs without defects, the loading capacity is greatly reduced. We employed solid-state NMR spectroscopy and molecular simulations to elucidate the drug-carrier interaction mechanisms. Through a synergistic combination of experimental and theoretical analyses, the docking conformations of pharmaceuticals at the defects were revealed.

Mechanistic Insights into the Regio- and Stereoselectivities of Testosterone and Dihydrotestosterone Hydroxylation Catalyzed by CYP3A4 and CYP19A1.

The hydroxylation of nonreactive C-H bonds can be easily catalyzed by a variety of metalloenzymes, especially cytochrome P450s (P450s). The mechanism of P450 mediated hydroxylation has been intensively studied, both experimentally and theoretically. However, understanding the regio- and stereoselectivities of substrates hydroxylated by P450s remains a great challenge. Herein, we use a multi-scale modeling approach to investigate the selectivity of testosterone (TES) and dihydrotestosterone (DHT) hydroxylation catalyzed by two important P450s, CYP3A4 and CYP19A1. For CYP3A4, two distinct binding modes for TES/DHT were predicted by dockings and molecular dynamics simulations, in which the experimentally identified sites of metabolism of TES/DHT can access to the catalytic center. The regio- and stereoselectivities of TES/DHT hydroxylation were further evaluated by quantum mechanical and ONIOM calculations. For CYP19A1, we found that sites 1ß, 2ß and 19 can access the catalytic center, with the intrinsic reactivity 2ß>1ß>19. However, our ONIOM calculations indicate that the hydroxylation is favored at site 19 for both TES and DHT, which is consistent with the experiments and reflects the importance of the catalytic environment in determining the selectivity. Our study unravels the mechanism underlying the selectivity of TES/DHT hydroxylation mediated by CYP3A4 and CYP19A1 and is helpful for understanding the selectivity of other substrates that are hydroxylated by P450s.

Dissecting the Structural Plasticity and Dynamics of Cytochrome P450 2B4 by Molecular Dynamics Simulations.

The plasticity of cytochromes P450 (P450s) is known to contribute significantly to their catalytic capacity of metabolizing various substrates. Although numerous studies have been performed, factors governing the plasticity and dynamics of P450s are still not fully understood. In this study, taking CYP2B4 as an example, we dissect the protein plasticity and dynamics in different environments. CYP2B4 is featured by a high degree of plasticity, which exhibits open, closed, and intermediate states. By analyzing the CYP2B4 crystal structures, we identified the structural features for the closed, open, and intermediate states. Interestingly, formation of the dimer structure was found in the open and intermediate states. The subsequent molecular dynamics (MD) simulations of the open structure in water confirmed the importance of the dimer form in stabilizing the open conformations. MD simulations of the closed and open structures in the membrane environment and the free energies for opening the F-G cassette obtained from the umbrella sampling calculations indicate that the membrane environment is important for stabilizing the F-G cassette. The dynamical network analysis indicates that Asp105 on the B-C loop plays an important role in transiting the structure from the open to the intermediate state. Our results thus unveil the mechanisms of dimer formation and open-to-intermediate transition for CYP2B4 in the water and membrane environments.

Free Energy Profile and Kinetics of Coupled Folding and Binding of the Intrinsically Disordered Protein p53 with MDM2.

Intrinsically disordered proteins (IDPs) exert their functions by binding to partner proteins via a complex process that includes coupled folding and binding. Because inhibiting the binding of the IDP p53 to its partner MDM2 has become a promising strategy for the design of anticancer drugs, we carried out metadynamics simulations to study the coupled folding and binding process linking the IDP p53 to MDM2 in atomic detail. Using bias-exchange metadynamics (BE-MetaD) and infrequent metadynamics (InMetaD), we estimated the binding free energy, the unbinding rate, and the binding rate. By analyzing the stable intermediates, we uncovered the role non-native interactions played in the p53-MDM2 binding/unbinding process. We used a three-state model to describe the whole binding/unbinding process and to obtain the corresponding rate constants. Our work shows that the binding of p53 favors an induced-fit mechanism which proceeds in a stepwise fashion. Our results can be helpful for gaining an in-depth understanding of the coupled folding and binding process needed for the design of MDM2 inhibitors.

Duet of Acetate and Water at the Defects of Metal-Organic Frameworks.

Metal-organic frameworks (MOFs) are porous crystalline materials with promising applications in molecular adsorption, separation, and catalysis. It has been discovered recently that structural defects introduced unintentionally or by design could have a significant impact on their properties. However, the exact chemical composition and structural evolution under different conditions at the defects are still under debate. In this study, we performed multidimensional solid-state nuclear magnetic resonance (SSNMR) coupled with computer simulations to elucidate an important scenario of MOF defects, uncovering the dynamic interplay between residual acetate and water. Acetate, as a defect modulator, and water, as a byproduct, are prevalent defect-associated species, which are among the key factors determining the reactivity and stability of defects. We discovered that acetate molecules coordinate to a single metal site monodentately and pair with water at the neighboring position. The acetates are highly flexible, which undergo fast libration as well as a slow kinetic exchange with water through dynamic hydrogen bonds. The dynamic processes under variable temperatures and different hydration levels have been quantitatively analyzed across a broad time scale from microseconds to seconds. The integration of SSNMR and computer simulations allows a precision probe into defective MOF structures with intrinsic dynamics and disorder.

Enhanced Sampling Simulations of Ligand Unbinding Kinetics Controlled by Protein Conformational Changes.

Understanding unbinding kinetics of protein-ligand systems is of great importance for the design of ligands with desired specificity and safety. In recent years, enhanced sampling techniques have emerged as effective tools for studying unbinding kinetics of protein-ligand systems at the atomistic level. However, in many protein-ligand systems, the ligand unbinding processes are strongly coupled to protein conformational changes and the disclosure of the hidden degrees of freedom closely related to the protein conformational changes so that sampling is enhanced over these degrees of freedom remains a great challenge. Here, we show how potential-scaled molecular dynamics (sMD) and infrequent metadynamics (InMetaD) simulation techniques can be combined to successfully reveal the unbinding mechanism of 3-(1,4-diazabicyclo[3.2.2]nonan-4-yl)-6-[18F]fluo-rodibenzo[b,d]thiophene 5,5-dioxide ([18F]ASEM) from a chimera structure of the α7-nicotinic acetylcholine receptor. By using sMD simulations, we disclosed that the "close" to "open" conformational change of loop C plays a key role in the ASEM unbinding process. By carrying out InMetaD simulations with this conformational change taken into account as an additional collective variable, we further captured the key states in the unbinding process and clarified the unbinding mechanism of ASEM from the protein. Our work indicates that combining sMD and InMetaD simulation techniques can be an effective approach for revealing the unbinding mechanism of a protein-ligand system where protein conformational changes control the unbinding process.

Activity of Antimicrobial Peptide Aggregates Decreases with Increased Cell Membrane Embedding Free Energy Cost.

Antimicrobial peptides (AMPs) are a promising alternative to antibiotics for mitigating bacterial infections, in light of increasing bacterial resistance to antibiotics. However, predicting, understanding, and controlling the antibacterial activity of AMPs remain a significant challenge. While peptide intramolecular interactions are known to modulate AMP antimicrobial activity, peptide intermolecular interactions remain elusive in their impact on peptide bioactivity. Herein, we test the relationship between AMP intermolecular interactions and antibacterial efficacy by controlling AMP intermolecular hydrophobic and hydrogen bonding interactions. Molecular dynamics simulations and Gibbs free energy calculations in concert with experimental assays show that increasing intermolecular interactions via interpeptide aggregation increases the energy cost for the peptide to embed into the bacterial cell membrane, which in turn decreases the AMP antibacterial activity. Our findings provide a route for predicting and controlling the antibacterial activity of AMPs against Gram-negative bacteria via reductions of intermolecular AMP interactions.

Mechanistic insights into peptide and ligand binding of the ATAD2-bromodomain via atomistic simulations disclosing a role of induced fit and conformational selection.

ATAD2 has emerged as a promising bromodomain (BRD)-containing therapeutic drug target in multiple human cancers. However, recent druggability assessment studies predicted ATAD2's BRD as a target 'difficult to drug' because its binding pocket possesses structural features that are unfeasible for ligand binding. Here, by using all-atom molecular dynamics simulations and an advanced metadynamics method, we demonstrate a dynamic view of the binding pocket features which can hardly be obtained from the "static" crystal data. The most important features disclosed from our simulation data, include: (1) a distinct 'open-to-closed' conformational switch of the ZA loop region in the context of peptide or ligand binding, akin to the induced fit mechanism of molecular recognition, (2) a dynamic equilibrium of the BC loop "in" and "out" conformations, highlighting a role in the conformational selection mechanism for ligand binding, and (3) a new binding region identified distal to the histone-binding pocket that might have implications in bromodomain biology and in inhibitor development. Moreover, based on our simulation results, we propose a model for an "auto-regulatory" mechanism of ATAD2's BRD for histone binding. Overall, the results of this study will not only have implications in bromodomain biology but also provide a theoretical basis for the discovery of new ATAD2's BRD inhibitors.

Organelle-Specific Triggered Release of Immunostimulatory Oligonucleotides from Intrinsically Coordinated DNA-Metal-Organic Frameworks with Soluble Exoskeleton.

DNA has proven of high utility to modulate the surface functionality of metal-organic frameworks (MOFs) for various biomedical applications. Nevertheless, current methods for preparing DNA-MOF nanoparticles rely on either inefficient covalent conjugation or specific modification of oligonucleotides. In this work, we report that unmodified oligonucleotides can be loaded on MOFs with high density (∼2500 strands/particle) via intrinsic, multivalent coordination between DNA backbone phosphate and unsaturated zirconium sites on MOFs. More significantly, surface-bound DNA can be efficiently released in either bulk solution or specific organelles in live cells when free phosphate ions are present. As a proof-of-concept for using this novel type of DNA-MOFs in immunotherapy, we prepared a construct of immunostimulatory DNA-MOFs (isMOFs) by intrinsically coordinating cytosine-phosphate-guanosine (CpG) oligonucleotides on biocompatible zirconium MOF nanoparticles, which was further armed by a protection shell of calcium phosphate (CaP) exoskeleton. We demonstrated that isMOFs exhibited high cellular uptake, organelle specificity, and spatiotemporal control of Toll-like receptors (TLR)-triggered immune responses. When isMOF reached endolysosomes via microtubule-mediated trafficking, the CaP exoskeleton dissolved in the acidic environment and in situ generated free phosphate ions. As a result, CpG was released from isMOFs and stimulated potent immunostimulation in living macrophage cells. Compared with naked CpG-MOF, isMOFs exhibited 83-fold up-regulation in stimulated secretion of cytokines. We thus expect this isMOF design with soluble CaP exoskeleton and an embedded sequential "protect-release" program provides a highly generic approach for intracellular delivery of therapeutic nucleic acids.

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 Trp2746.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.

Insight into the adsorption profiles of the Saprolegnia monoica chitin synthase MIT domain on POPA and POPC membranes by molecular dynamics simulation studies.

The critical role of chitin synthases in oomycete hyphal tip growth has been established. A microtubule interacting and trafficking (MIT) domain was discovered in the chitin synthases of the oomycete model organism, Saprolegnia monoica. MIT domains have been identified in diverse proteins and may play a role in intracellular trafficking. The structure of the Saprolegnia monoica chitin synthase 1 (SmChs1) MIT domain has been recently determined by our group. However, although our in vitro assay identified increased strength in interactions between the MIT domain and phosphatidic acid (PA) relative to other phospholipids including phosphatidylcholine (PC), the mechanism used by the MIT domain remains unknown. In this work, the adsorption behavior of the SmChs1 MIT domain on POPA and POPC membranes was systematically investigated by molecular dynamics simulations. Our results indicate that the MIT domain can adsorb onto the tested membranes in varying orientations. Interestingly, due to the specific interactions between MIT residues and lipid molecules, the binding affinity to the POPA membrane is much higher than that to the POPC membrane. A binding hotspot, which is critical for the adsorption of the MIT domain onto the POPA membrane, was also identified. The lower binding affinity to the POPC membrane can be attributed to the self-saturated membrane surface, which is unfavorable for hydrogen-bond and electrostatic interactions. The present study provides insight into the adsorption profile of SmChs1 and additionally has the potential to improve our understanding of other proteins containing MIT domains.

Theoretical study of the binding profile of an allosteric modulator NS-1738 with a chimera structure of the α7 nicotinic acetylcholine receptor.

Potentiation of the function of the α7 nicotinic acetylcholine receptor (α7-nAChR) is believed to provide a possible way for the treatment of cholinergic system dysfunctions such as Alzheimer's disease and schizophrenia. Positive allosteric modulators (PAMs) are able to augment the peak current response of the endogenous agonist of α7-nAChR by binding to some allosteric sites. In this study, the binding profile of a potent type I PAM, NS-1738, with a chimera structure (termed α7-AChBP) constructed from the extracellular domain of α7-nAChR and an acetylcholine binding protein was investigated with molecular docking, molecular dynamics simulation, and free energy calculation methods. We found that NS-1738 could bind to three allosteric sites of α7-AChBP, namely, the top pocket, the vestibule pocket and the agonist sub-pocket. NS-1738 has moderate binding affinities (-6.76 to -9.15 kcal mol-1) at each allosteric site. The urea group is critical for binding and can form hydrogen-bond interactions with the protein. The bulky trifluoromethyl group also has a great impact on the binding modes and binding affinities. We believe that our study provides valuable insight into the binding profiles of type I PAMs with α7-nAChR and is helpful for the development of novel PAMs.

Molecular adhesion at clay nanocomposite interfaces depends on counterion hydration-molecular dynamics simulation of montmorillonite/xyloglucan.

Nacre-mimetic clay/polymer nanocomposites with clay platelet orientation parallel to the film surface show interesting gas barrier and mechanical properties. In moist conditions, interfacial adhesion is lowered and mechanical properties are reduced. Molecular dynamic simulations (MD) have been performed to investigate the effects of counterions on molecular adhesion at montmorillonite clay (Mnt)-xyloglucan (XG) interfaces. We focus on the role of monovalent cations K(+), Na(+), and Li(+) and the divalent cation Ca(2+) for mediating and stabilizing the Mnt/XG complex formation. The conformation of adsorbed XG is strongly influenced by the choice of counterion and so is the simulated work of adhesion. Free energy profiles that are used to estimate molecular adhesion show stronger interaction between XG and clay in the monovalent cation system than in divalent cation system, following a decreasing order of K-Mnt, Na-Mnt, Li-Mnt, and Ca-Mnt. The Mnt clay hydrates differently in the presence of different counterions, leading to a chemical potential of water that is highest in the case of K-Mnt, followed by Na-Mnt and Li-Mnt, and lowest in the case of Ca-Mnt. This means that water is most easily displaced from the interface in the case of K-Mnt, which contributes to the relatively high work of adhesion. In all systems, the penalty of replacing polymer with water at the interface gives a positive contribution to the work of adhesion of between 19 and 35%. Our work confirms the important role of counterions in mediating the adsorption of biopolymer XG to Mnt clays and predicts potassium or sodium as the best choice of counterions for a Mnt-based biocomposite design.

Origin of ion selectivity at the air/water interface.

Among many characteristics of ions, their capability to accumulate at air/water interfaces is a particular issue that has been the subject of much research attention. For example, the accumulation of halide anions (Cl(-), Br(-), I(-)) at the water surface is of great importance to heterogeneous reactions that are of environmental concern. However, the actual mechanism that drives anions towards the air/water interface remains unclear. In this work, we have performed atomistic simulations using polarizable models to mimic ionic behavior under atmospheric conditions. We find that larger anions are abundant at the water surface and that the cations are pulled closer to the surface by the counterions. We propose that polarization effects stabilize the anions with large radii when approaching the surface. This energetically more favorable situation is caused by the fact that the more polarized anions at the surface attract water molecules more strongly. Of relevance is also the ordering of the surface water molecules with their hydrogen atoms pointing outwards which induce an external electronic field that leads to a different surface behavior of anions and cations. The water-water interaction is weakened by the distinct water-ion attraction, a point contradicting the proposition that F(-) is a kosmotrope. The simulation results thus allow us to obtain a more holistic understanding of the interfacial properties of ionic solutions and atmospheric aerosols.

A theoretical study on the molecular determinants of the affibody protein Z(Aß3) bound to an amyloid ß peptide.

Amyloid beta (Aß) peptides are small cleavage products of the amyloid precursor protein. Aggregates of Aß peptides are thought to be linked with Alzheimer's and other neurodegenerative diseases. Strategies aimed at inhibiting amyloid formation and promoting Aß clearance have been proposed and investigated in in vitro experiments and in vivo therapies. A recent study indicated that a novel affibody protein ZAß3, which binds to an Aß40 monomer with a binding affinity of 17 nM, is able to prevent the aggregation of Aß40. However, little is known about the energetic contribution of each residue in ZAß3 to the formation of the (ZAß3)2:Aß complex. To address this issue, we carried out unbiased molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area calculations. Through the per-residue decomposition scheme, we identified that the van der Waals interactions between the hydrophobic residues of (ZAß3)2 and those at the exterior and interior faces of Aß are the main contributors to the binding of (ZAß3)2 to Aß. Computational alanine scanning identified 5 hot spots, all residing in the binding interface and contributing to the binding of (ZAß3)2 to Aß through the hydrophobic effect. In addition, the amide hydrogen bonds in the 4-strand ß-sheet and the π-π stacking were also analyzed. Overall, our study provides a theoretical basis for future experimental improvement of the ZAß3 peptide binding to Aß.

Molecular dynamics simulations reveal the assembly mechanism of polysaccharides in marine aerosols.

The high Arctic marine environment has recently detected polymer gels in atmospheric aerosol particles and cloud water originating from the surface microlayer of the open leads within the pack ice area. These polysaccharide molecules are water insoluble but water solvated, highly surface-active and highly hydrated (99% water). In order to add to the understanding and to complement missing laboratory characterization of marine polymer gels we have in this work performed an atomistic study of the assembly process and interfacial properties of polysaccharides. Our study reveals a number of salient features of the microscopic process behind polysaccharide assembly into nanogels. With three- and four-repeating units the polysaccharides assemble into a cluster in 50 ns. The aggregates grow quicker by absorbing one or two polymers each time, depending on the unit length and the type of inter-bridging cation. Although both the hydrophobic and hydrophilic domains are contracted, the latter dominates distinctly upon the contraction of solvent accessible surface areas. The establishment of inter-chain hydrogen-bonds is the key to the assembly while ionic bridges can further promote aggregation. During the assembly of the more bent four-unit polymers, intra-chain hydrogen bonds are significantly diminished by Ca(2+). Meanwhile, the percentage of Ca(2+) acting as an ionic bridge is more eminent, highlighting the significance of Ca(2+) ions for longer-chain polysaccharides. The aggregates are able to enhance surface tension more in the presence of Ca(2+) than in the presence of Na(+) owing to their more compact structure. These conclusions all demonstrate that studies of the present kind provide insight into the self-assembly process and interfacial properties of marine gels. We hope this understanding will keep up the interest in the complex and the fascinating relationship between marine microbiology, atmospheric aerosols, clouds and climate.