Sigmoid Accelerated Molecular Dynamics: An Efficient Enhanced Sampling Method for Biosystems.

Gaussian accelerated molecular dynamics (GaMD) is recognized as a popular enhanced sampling method for tackling long-standing challenges in biomolecular simulations. Inspired by GaMD, Sigmoid accelerated molecular dynamics (SaMD) is proposed in this work by adding a Sigmoid boost potential to improve the balance between the highest acceleration and accurate reweighting. Compared with GaMD, SaMD extends the accessible time scale and improves the computational efficiency as tested in three tasks. In the alanine dipeptide task, SaMD can produce the free energy landscape with better accuracy and efficiency. In the chignolin folding task, the estimated Gibbs free energy difference can converge to the experimental value ∼30% faster. In the protein-ligand binding task, the bound conformations are closer to the crystal structure with a minimal ligand root-mean-square deviation of 1.7 Å. The binding of the ligand XK263 to the HIV protease is reproduced by SaMD in ∼60% less simulation time.

A combined computational and experimental approach predicts thrombin adsorption to zeolites.

Robust protein-nanomaterial surface analysis is important, but also a challenge. Thrombin plays an important role in the coagulant activity of protein corona mediated by Ca2+ ion exchanged zeolites. However, the mechanism for this modulation remains unresolved. In this study, we proposed a combined computational and experimental approach to determine the adsorbed sites and orientations of thrombin binding to Ca2+ -exchanged LTA-type (CaA) zeolite. Specifically, fourteen ensembles of simulated annealing molecular dynamics (SAMD) simulations and experimental surface residues microenvironment analysis were used to reduce the starting orientations needed for further molecular dynamics (MD) simulations. The combined MD simulations and procoagulant activity characterization also reveal the consequent corresponding deactivation of thrombin on CaA zeolite. It is mainly caused by two aspects: (1) the secondary structure of thrombin can change after its adsorption on the CaA zeolite. (2) The positively charged area of thrombin mediates the preferential interaction between thrombin and CaA zeolite. Some thrombin substrate sites are thus blocked by zeolite after its adsorption. This study not only provides a promising method for characterizing the protein-nanoparticle interaction, but also gives an insight into the design and application of zeolite with high procoagulant activity.

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.

Binding modes of prothrombin cleavage site sequences to the factor Xa catalytic triad: Insights from atomistic simulations.

Prothrombin is a key zymogen of the coagulation process and can be converted to thrombin by the prothrombinase complex, which consists of factor Xa (FXa), cofactor Va (FVa), and phospholipids. Prothrombin can be activated at two cleavage sites, R271 and R320, which generates two intermediates: prethrombin-2 via the initial cleavage at R271, and meizothrombin via the first cleavage at R320. Several mechanisms have been proposed to explain this activation preference, but the role of cleavage site sequences in prothrombin activation has not been thoroughly investigated. Here, we used an advanced sampling technique, parallel tempering metadynamics with a well-tempered ensemble (PTMetaD-WTE), to study the binding modes of prothrombin cleavage site sequences R266AIEGRTATSEY277 (denoted as Pep271) and S315YIDGRIVEGSD326 (denoted as Pep320) to the FXa catalytic triad. Our study indicates that there exist three binding modes for Pep271 to the FXa catalytic triad but only one binding mode for Pep320 to the FXa catalytic triad. Further molecular dynamics simulations revealed that due to the strong electrostatic interactions, especially the H-bond interactions and salt bridges formed between Pep320 and FXa, the binding mode in the Pep320-FXa system is more stable than the binding modes in the Pep271-FXa system. In view of experimental observations and our results that there exists only one binding mode for Pep320 to the FXa catalytic triad and especially R320 in Pep320 can stably bind to the FXa catalytic triad, we believe that the first cleavage at R320 is favored.

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.

Development of 11C-Labeled ASEM Analogues for the Detection of Neuronal Nicotinic Acetylcholine Receptors (α7-nAChR).

The homo-pentameric alpha 7 receptor is one of the major types of neuronal nicotinic acetylcholine receptors (α7-nAChRs) related to cognition, memory formation, and attention processing. The mapping of α7-nAChRs by PET pulls a lot of attention to realize the mechanism and development of CNS diseases such as AD, PD, and schizophrenia. Several PET radioligands have been explored for the detection of the α7-nAChR. 18F-ASEM is the most functional for in vivo quantification of α7-nAChRs in the human brain. The first aim of this study was to initially use results from in silico and machine learning techniques to prescreen and predict the binding energy and other properties of ASEM analogues and to interpret these properties in terms of atomic structures using 18F-ASEM as a lead structure, and second, to label some selected candidates with carbon-11/hydrogen-3 (11C/3H) and to evaluate the binding properties in vitro and in vivo using the labeled candidates. In silico predictions are obtained from perturbation free-energy calculations preceded by molecular docking, molecular dynamics, and metadynamics simulations. Machine learning techniques have been applied for the BBB and P-gp-binding properties. Six analogues of ASEM were labeled with 11C, and three of them were additionally labeled with 3H. Binding properties were further evaluated using autoradiography (ARG) and PET measurements in non-human primates (NHPs). Radiometabolites were measured in NHP plasma. All six compounds were successfully synthesized. Evaluation with ARG showed that 11C-Kln83 was preferably binding to the α7-nAChR. Competition studies showed that 80% of the total binding was displaced. Further ARG studies using 3H-KIn-83 replicated the preliminary results. In the NHP PET study, the distribution pattern of 11C-KIn-83 was similar to other α7 nAChR PET tracers. The brain uptake was relatively low and increased by the administration of tariquidar, indicating a substrate of P-gp. The ASEM blocking study showed that 11C-KIn-83 specifically binds to α7 nAChRs. Preliminary in vitro evaluation of KIn-83 by ARG with both 11C and 3H and in vivo evaluation in NHP showed favorable properties for selectively imaging α7-nAChRs, despite a relatively low brain uptake.

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.

In silico studies of ASEM analogues targeting α7-nAChR and experimental verification.

The α7 nicotinic acetylcholine receptor (α7-nAChR) is implicated in a variety of neurodegenerative and neuropsychiatric disorders, such as Alzheimer's disease (AD) and schizophrenia. The progress of these disorders can be studied using positron emission tomography (PET) with radiotracers for α7-nAChR. [18F]ASEM and [18F] para-ASEM (also referred to as [18F]DBT-10) are novel and potent α7-nAChR PET radiotracers which have successfully been used in human subjects and nonhuman primates, though further improvement of them is still a pressing task in the community of neurodegeneration research. In this work, we demonstrate the use of modern in silico techniques to predict the binding modes, binding strengths, and residence times for molecular PET tracers binding to proteins, using ASEM and DBT-10 as a showcase of the predictive and interpretational power of such techniques, in particular free energy perturbation theory. The corresponding compounds were synthesized and further tested by in vitro binding experiment for validation. Encouragingly, our in silico modeling can correctly predict the binding affinities of the ASEM analogues. The structure-activity relationships for the ortho- and para-substitutions are well explained at the atomistic level and provide structure-based guiding for the future development of PET tracers for α7-nAChR. A discussion is presented on the complementary use of in silico rational methods based on atomic and electronic principles for in vitro characterization of PET tracers.

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.

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.

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.

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.

Free Energy Profile for Penetration of Pittsburgh Compound-B into the Amyloid ß Fibril.

The amyloid ß (Aß) fibril is a hallmark of Alzheimer's disease (AD) and has therefore served as an important target for early diagnosis of AD. The Pittsburgh Compound-B (PiB) is one of the most famous positron emission tomography (PET) tracers commonly used for in vivo detection of Aß fibrils. Many theoretical studies have predicted the existence of various core binding sites with different microenvironments for probes binding to the Aß fibril. However, little attention has been devoted to how the probes actually penetrate into the different core binding sites. In this study, an integrated molecular modeling scheme is used to study the penetration of PiB into the core binding sites of the Aß1-42 fibril structure recently obtained by cryogenic electron microscopy. We find that there are two core binding sites for PiB with dramatic differences in cavity size and microenvironment properties, and furthermore that the penetration of PiB into site-1 is energetically prohibitive, whereas the penetration into site-2 is much more favorable. Therefore, the binding capacity at site-2 may be larger than that at site-1 despite its lower binding affinity. Our results thus suggest that site-2 may be a major binding site for PiB binding to Aß fibril and emphasize the importance to adopt a full dynamical picture when studying tracer-fibril binding problems in general, something that in turn can be used to guide the development of tracers with higher affinity and selectivity for the Aß fibril.

Computational Insight Into Vitamin K1 ω-Hydroxylation by Cytochrome P450 4F2.

Vitamin K1 (VK1) plays an important role in the modulation of bleeding disorders. It has been reported that ω-hydroxylation on the VK1 aliphatic chain is catalyzed by cytochrome P450 4F2 (CYP4F2), an enzyme responsible for the metabolism of eicosanoids. However, the mechanism of VK1 ω-hydroxylation by CYP4F2 has not been disclosed. In this study, we employed a combination of quantum mechanism (QM) calculations, homology modeling, molecular docking, molecular dynamics (MD) simulations, and combined quantum mechanism/molecular mechanism (QM/MM) calculations to investigate the metabolism profile of VK1 ω-hydroxylation. QM calculations based on the truncated VK1 model show that the energy barrier for ω-hydroxylation is about 6-25 kJ/mol higher than those at other potential sites of metabolism. However, results from the MD simulations indicate that hydroxylation at the ω-site is more favorable than at the other potential sites, which is in accordance with the experimental observation. The evaluation of MD simulations was further endorsed by the QM/MM calculation results. Our studies thus suggest that the active site residues of CYP4F2 play a determinant role in the ω-hydroxylation. Our results provide structural insights into the mechanism of VK1 ω-hydroxylation by CYP4F2 at the atomistic level and are helpful not only for characterizing the CYP4F2 functions but also for looking into the ω-hydroxylation mediated by other CYP4 enzymes.

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.

Discovery of Highly Potent Pinanamine-Based Inhibitors against Amantadine- and Oseltamivir-Resistant Influenza A Viruses.

Influenza pandemic is a constant major threat to public health caused by influenza A viruses (IAVs). IAVs are subcategorized by the surface proteins hemagglutinin (HA) and neuraminidase (NA), in which they are both essential targets for drug discovery. While it is of great concern that NA inhibitor oseltamivir resistant strains are frequently identified from human or avian influenza virus, structural and functional characterization of influenza HA has raised hopes for new antiviral therapies. In this study, we explored a structure-activity relationship (SAR) of pinanamine-based antivirals and discovered a potent inhibitor M090 against amantadine-resistant viruses, including the 2009 H1N1 pandemic strains, and oseltamivir-resistant viruses. Mechanism of action studies, particularly hemolysis inhibition, indicated that M090 targets influenza HA and it occupied a highly conserved pocket of the HA2 domain and inhibited virus-mediated membrane fusion by "locking" the bending state of HA2 during the conformational rearrangement process. This work provides new binding sites within the HA protein and indicates that this pocket may be a promising target for broad-spectrum anti-influenza A drug design and development.