Structure prediction and molecular dynamics simulations of a G-protein coupled receptor: human CCR2 receptor

CC chemokine receptor type-2 (CCR2) is a member of G-protein coupled receptors superfamily, expressed on the cell surface of monocytes and macrophages. It binds to the monocyte chemoattractant protein-1, a CC chemokine, produced at the sites of inflammation and infection. A homology model of human CCR2 receptor based on the recently available C-X-C chemokine recepor-4 crystal structure has been reported. Ligand information was used as an essential element in the homology modeling process. Six known CCR2 antagonists were docked into the model using simple and induced fit docking procedure. Docked complexes were then subjected to visual inspection to check their suitability to explain the experimental data obtained from site directed mutagenesis and structure-activity relationship studies. The homology model was refined, validated, and assessed for its performance in docking-based virtual screening on a set of CCR2 antagonists and decoys. The docked complexes of CCR2 with the known antagonists, TAK779, a dual CCR2/CCR5 antagonist, and Teijin-comp1, a CCR2 specific antagonist were subjected to molecular dynamics (MD) simulations, which further validated the binding modes of these antagonists. B-factor analysis of 20?ns MD simulations demonstrated that Cys190 is helpful in providing structural rigidity to the extracellular loop (EL2). Residues important for CCR2 antagonism were recognized using free energy decomposition studies. The acidic residue Glu291 from TM7, a conserved residue in chemokine receptors, is favorable for the binding of Teijin-comp1 with CCR2 by ΔG of ?11.4?kcal/mol. Its contribution arises more from the side chains than the backbone atoms. In addition, Tyr193 from EL2 contributes ?0.9?kcal/mol towards the binding of the CCR2 specific antagonist with the receptor. Here, the homology modeling and subsequent molecular modeling studies proved successful in probing the structure of human CCR2 chemokine receptor for the structure-based virtual screening and predicting the binding modes of CCR2 antagonists.     

Insight into the key active sites of afChiA1 based on molecular dynamics simulations and free energy calculations

In this study, the binding of the enzyme chitinase A1 (afChiA1) from the plant-type Aspergillus fumigatus with four potent inhibitors, allosamidin (ASM), acetazolamide (AZM), 8-chloro-theophylline (CTP) and kinetin (KIT) is investigated by molecular docking, molecular dynamics simulation and binding free energy calculation. The results reveal that the electrostatic interactions play an important role in the stabilisation of the binding of afChiA1 with inhibitors. Based on the binding energy of afChiA1-ligands, the key residues (Gln37 and Trp312) in the active binding pocket of the complex systems are confirmed by molecular mechanics/Poisson–Boltzmann surface area method, and the active inhibitors, ASM and AZM, both could form strong interaction with Gln37 and Trp312, and the non-active ligands, CTP and KIT, could not interact with these two residues, which is consistent with the result of experimental report. Then, it is identified that Gln37 and Trp312 should be one of the important active site residues of afChiA1.     

All-atom molecular dynamics simulations reveal how kinesin transits from one-head-bound to two-heads-bound state

Kinesin dimer walks processively along a microtubule (MT) protofilament in a hand-over-hand manner, transiting alternately between one-head-bound (1HB) and two-heads-bound (2HB) states. In 1HB state, one head bound by adenosine diphosphate (ADP) is detached from MT and the other head is bound to MT. Here, using all-atom molecular dynamics simulations we determined the position and orientation of the detached ADP-head relative to the MT-bound head in 1HB state. We showed that in 1HB state when the MT-bound head is in ADP or nucleotide-free state, with its neck linker being undocked, the detached ADP-head and the MT-bound head have the high binding energy, and after adenosine triphosphate (ATP) binds to the MT-bound head, with its neck linker being docked, the binding energy between the two heads is reduced greatly. These results reveal how the kinesin dimer retains 1HB state before ATP binding and how the dimer transits from 1HB to 2HB state after ATP binding. Key residues involved in the head-head interaction in 1HB state were identified.     

Functional roles of magnesium binding to extracellular signal-regulated kinase 2 explored by molecular dynamics simulations and principal component analysis

Molecular dynamics (MD) simulations coupled with principal component (PC) analysis were carried out to study functional roles of Mg²⁺ binding to extracellular signal-regulated kinase 2 (ERK2). The results suggest that Mg²⁺ binding heavily decreases eigenvalue of the first principal component and totally inhibits motion strength of ERK2, which favors stabilization of ERK2 structure. Binding free energy predictions indicate that Mg²⁺ binding produces an important effect on binding ability of adenosine triphosphate (ATP) to ERK2 and strengthens the ATP binding. The calculations of residue-based free energy decomposition show that lack of Mg²⁺ weakens interactions between the hydrophobic rings of ATP and five residues I29, V37, A50, L105, and L154. Hydrogen bond analyses also prove that Mg²⁺ binding increases occupancies of hydrogen bonds formed between ATP and residues K52, Q103, D104, and M106. We expect that this study can provide a significant theoretical hint for designs of anticancer drugs targeting ERK2.     

Position-resolved free energy of solvation for amino acids in lipid membranes from molecular dynamics simulations

Studies of insertion and interactions of amino acids in lipid membranes are pivotal to our understanding of membrane protein structure and function. Calculating the insertion cost as a function of transmembrane helix sequence is thus an important step towards improved membrane protein prediction and eventually drug design. Here, we present position-dependent free energies of solvation for all amino acid analogs along the membrane normal. The profiles cover the entire region from bulk water to hydrophobic core, and were produced from all-atom molecular dynamics simulations. Experimental differences corresponding to mutations and costs for entire segments match experimental data well, and in addition the profiles provide the spatial resolution currently not available from experiments. Polar side-chains largely maintain their hydration and assume quite ordered conformations, which indicates the solvation cost is mainly entropic. The cost of solvating charged side-chains is not only significantly lower than for implicit solvation models, but also close to experiments, meaning these could well maintain their protonation states inside the membrane. The single notable exception to the experimental agreement is proline, which is quite expensive to introduce in vivo despite its hydrophobicity--a difference possibly explained by kinks making it harder to insert helices in the translocon.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

Molecular mechanism study of several inhibitors binding to BRD9 bromodomain based on molecular simulations

Bromodomain-containing protein 9 (BRD9) has been employed as a potential target for anticancer drugs in recent years. In this work, molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and per residue energy decomposition approaches were performed to elucidate the different binding modes between four pyridinone-like scaffold inhibitors and BRD9 bromodomain. Analysis results indicate that non-polar contribution mainly deriving from van der Waals energy is a critical impact on binding affinity of inhibitors against BRD9. Some key residues Phe44, Phe47, Val49, and Ile53 (at ZA loop) enhance the binding energy of inhibitors in BRD9 by means of providing hydrophobic interactions. Moreover, it is observed that BRD9 is anchored by the formation of a stable hydrogen bond between the carbonyl of the inhibitors and the residue Asn100 (at BC loop), and a strong π–π stacking interaction formed between the residue Tyr106 (at BC loop) and the inhibitors. The existence of dimethoxyphenyl structure and the aromatic ring merged to pyridinone scaffold are useful to enhance the BRD9 binding affinity. These findings should guide the rational design of more prospective inhibitors targeting BRD9.

Communicated by Ramaswamy H. Sarma     

Reproducible polypeptide folding and structure prediction using molecular dynamics simulations

The folding of a polypeptide from an extended state to a well-defined conformation is studied using microsecond classical molecular dynamics (MD) simulations and replica exchange molecular dynamics (REMD) simulations in explicit solvent and in vacuo. It is shown that the solvated peptide folds many times in the REMD simulations but only a few times in the conventional simulations. From the folding events in the classical simulations we estimate an approximate folding time of 1-2 micros. The REMD simulations allow enough sampling to deduce a detailed Gibbs free energy landscape in three dimensions. The global minimum of the energy landscape corresponds to the native state of the peptide as determined previously by nuclear magnetic resonance (NMR) experiments. Starting from an extended state it takes about 50 ns before the native structure appears in the REMD simulations, about an order of magnitude faster than conventional MD. The calculated melting curve is in good qualitative agreement with experiment. In vacuo, the peptide collapses rapidly to a conformation that is substantially different from the native state in solvent.     

Partition of protein solvation into group contributions from molecular dynamics simulations

Linear response theory coupled to molecular dynamics simulations with an explicit solvent representation is used to derive fractional contributions of amino acid residues to the solvation of proteins. The new fractional methods developed here are compared with standard approaches based on empirical 1D and 3D statistical potentials, as well as with estimates obtained from the analysis of classical molecular interaction potentials. The new fractional methods, which have a clear physical basis and explicitly account for the effects due to protein structure and flexibility, provide an accurate picture of the contribution to solvation of different regions of the protein.     

Exploration of the chlorpyrifos escape pathway from acylpeptide hydrolases using steered molecular dynamics simulations

Acylpeptide hydrolases (APH) catalyze the removal of an N-acylated amino acid from blocked peptides. APH is significantly more sensitive than acetylcholinesterase, a target of Alzheimer’s disease, to inhibition by organophosphorus (OP) compounds. Thus, OP compounds can be used as a tool to probe the physiological functions of APH. Here, we report the results of a computational study of molecular dynamics simulations of APH bound to the OP compounds and an exploration of the chlorpyrifos escape pathway using steered molecular dynamics (SMD) simulations. In addition, we apply SMD simulations to identify potential escape routes of chlorpyrifos from hydrolase hydrophobic cavities in the APH-inhibitor complex. Two previously proposed APH pathways were reliably identified by CAVER 3.0, with the estimated relative importance of P1 > P2 for its size. We identify the major pathway, P2, using SMD simulations, and Arg526, Glu88, Gly86, and Asn65 are identified as important residues for the ligand leaving via P2. These results may help in the design of APH-targeting drugs with improved efficacy, as well as in understanding APH selectivity of the inhibitor binding in the prolyl oligopeptidase family.     

Dissociation mechanism of GDP from Cdc42 via DOCK9 revealed by molecular dynamics simulations

Cell division control protein 42 homolog (Cdc42) influences a variety of cellular responses such as cell migration and polarity. Deregulation of Cdc42 has been associated with several human diseases and developmental disorders. Over-activation of Cdc42 through guanine nucleotide exchange factor (GEF) is a critical event for Cdc42 involved cancer metastasis. Members of DOCK family of GEF are important activators of Cdc42. However, this activation mechanism is still unknown. Molecular dynamics (MD) simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) calculations were employed to investigate the central step of the activation of Cdc42: the dissociation mechanism of GDP from Cdc42 via DOCK9. Simulation results show that Mg²⁺ ion has a remarkable influence on the conformational change of switch I of Cdc42 through residue Pro34 which functions as a “clasp” to control the flexibility of switch I. In the GDP dissociation process, the Mg²⁺ ion leave first to result in a suitable conformation of Cdc42 for following DOCK9 binding to. When DOCK9 binds to Cdc42, it changes the orientations of residues Lys16, Thr17, Cys18 and Phe28 of Cdc42 to weaken the interactions between Cdc42 and GDP to release GDP. This study first elucidates the dissociation mechanism of GDP from Cdc42 via DOCK9 and identifies the essential residues of Cdc42 in this process. These simulation results are consistent with the recent findings of biochemical and amino acid mutational studies, and the observations are beneficial to understand the activation mechanism of Cdc42 and to provide insights for designing compounds targeting on Cdc42 related cancer metastasis.     

Studying pressure denaturation of a protein by molecular dynamics simulations

Many globular proteins unfold when subjected to several kilobars of hydrostatic pressure. This “unfolding‐up‐on‐squeezing” is counter‐intuitive in that one expects mechanical compression of proteins with increasing pressure. Molecular simulations have the potential to provide fundamental understanding of pressure effects on proteins. However, the slow kinetics of unfolding, especially at high pressures, eliminates the possibility of its direct observation by molecular dynamics (MD) simulations. Motivated by experimental results—that pressure denatured states are water‐swollen, and theoretical results—that water transfer into hydrophobic contacts becomes favorable with increasing pressure, we employ a water insertion method to generate unfolded states of the protein Staphylococcal Nuclease (Snase). Structural characteristics of these unfolded states—their water‐swollen nature, retention of secondary structure, and overall compactness—mimic those observed in experiments. Using conformations of folded and unfolded states, we calculate their partial molar volumes in MD simulations and estimate the pressure‐dependent free energy of unfolding. The volume of unfolding of Snase is negative (approximately ?60 mL/mol at 1 bar) and is relatively insensitive to pressure, leading to its unfolding in the pressure range of 1500–2000 bars. Interestingly, once the protein is sufficiently water swollen, the partial molar volume of the protein appears to be insensitive to further conformational expansion or unfolding. Specifically, water‐swollen structures with relatively low radii of gyration have partial molar volume that are similar to that of significantly more unfolded states. We find that the compressibility change on unfolding is negligible, consistent with experiments. We also analyze hydration shell fluctuations to comment on the hydration contributions to protein compressibility. Our study demonstrates the utility of molecular simulations in estimating volumetric properties and pressure stability of proteins, and can be potentially extended for applications to protein complexes and assemblies. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

The mechanism of high affinity pentasaccharide binding to antithrombin,insights from Gaussian accelerated molecular dynamics simulations

Abstract

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from ‘static’ structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a ‘non-activated’ AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. ‘Generalized correlation’ matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.

Communicated by Ramaswamy H. Sarma     

A metabolomic view of how the human gut microbiota impacts the host metabolome using humanized and gnotobiotic mice

Defining the functional status of host-associated microbial ecosystems has proven challenging owing to the vast number of predicted genes within the microbiome and relatively poor understanding of community dynamics and community–host interaction. Metabolomic approaches, in which a large number of small molecule metabolites can be defined in a biological sample, offer a promising avenue to ‘fingerprint'' microbiota functional status. Here, we examined the effects of the human gut microbiota on the fecal and urinary metabolome of a humanized (HUM) mouse using an optimized ultra performance liquid chromatography–mass spectrometry-based method. Differences between HUM and conventional mouse urine and fecal metabolomic profiles support host-specific aspects of the microbiota''s metabolomic contribution, consistent with distinct microbial compositions. Comparison of microbiota composition and metabolome of mice humanized with different human donors revealed that the vast majority of metabolomic features observed in donor samples are produced in the corresponding HUM mice, and individual-specific features suggest ‘personalized'' aspects of functionality can be reconstituted in mice. Feeding the mice a defined, custom diet resulted in modification of the metabolite signatures, illustrating that host diet provides an avenue for altering gut microbiota functionality, which in turn can be monitored via metabolomics. Using a defined model microbiota consisting of one or two species, we show that simplified communities can drive major changes in the host metabolomic profile. Our results demonstrate that metabolomics constitutes a powerful avenue for functional characterization of the intestinal microbiota and its interaction with the host.     

HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.