Molecular dynamics simulations on SDF-1alpha: binding with CXCR4 receptor

Insights into the interacting mode of CXCR4 with SDF-1alpha are crucial in understanding the structural and functional characteristics of CXCR4 receptor. In this paper a computational pipeline, integrating protein structure prediction, molecular dynamics simulations, automated molecular docking, and Brownian dynamics simulations were employed to investigate the dynamic and energetic aspects of CXCR4 associating with SDF-1alpha. The entire simulation revealed the surface distribution feature of electrostatic potentials and conformational "open-close" process of the receptor. The possible binding conformation of CXCR4 was identified, and the CXCR4-SDF-1alpha binding complex was generated. Arg188-Glu277 salt bridge plays an important role for both the extracellular domain conformational change and SDF-1alpha binding. Two binding sites were mapped at the extracellular domain (Site 1) and inside the transmembrane domain (Site 2), which are composed of conserved residues. Sites 1 and 2 contribute approximately 60% and 40% to the binding affinity with SDF-1alpha, respectively. The binding model is in agreement with most of the experimental data. Transmembrane VI has more significant motion in the harmonious conformational transition of CXCR4 during SDF-1alpha binding, which may be possibly associated with signal transduction. Based on the modeling and simulation, a binding mechanism hypothesis between CXCR4 and SDF-1alpha and its relationship to the signal transduction has been proposed.     

Molecular modeling of the GABAC receptor ligand-binding domain

We have constructed a molecular model of the ligand-binding domain of the GABA(C) receptor, which is a member of the Cys-loop ligand-gated ion channel family. The extracellular domains of these receptors share similar sequence homology (20%) with Limnaea acetylcholine-binding protein for which an X-ray crystal structure is available. We used this structure as a template for homology modeling of the GABA(C) receptor extracellular domain using FUGUE and MODELLER software. FlexX was then used to dock GABA into the receptor ligand-binding site, resulting in three alternative energetically favorable orientations. Residues located no more than 5 A from the docked GABA were identified for each model; of these, three were found to be common to all models with 14 others present only in certain models. Using data from experimental studies, we propose that the most likely orientation of GABA is with its amine close to Y198, and its carboxylate close to R104. These studies have therefore provided a model of the ligand-binding domain, which will be useful for both GABA(C) and GABA(A) receptor studies, and have also yielded an experimentally testable hypothesis of the location of GABA in the binding pocket. [Figure: see text].     

Molecular dynamics simulations suggest ligand's binding to nicotinamidase/pyrazinamidase

The research on the binding process of ligand to pyrazinamidase (PncA) is crucial for elucidating the inherent relationship between resistance of Mycobacterium tuberculosis and PncA's activity. In the present study, molecular dynamics (MD) simulation methods were performed to investigate the unbinding process of nicotinamide (NAM) from two PncA enzymes, which is the reverse of the corresponding binding process. The calculated potential of mean force (PMF) based on the steered molecular dynamics (SMD) simulations sheds light on an optimal binding/unbinding pathway of the ligand. The comparative analyses between two PncAs clearly exhibit the consistency of the binding/unbinding pathway in the two enzymes, implying the universality of the pathway in all kinds of PncAs. Several important residues dominating the pathway were also determined by the calculation of interaction energies. The structural change of the proteins induced by NAM's unbinding or binding shows the great extent interior motion in some homologous region adjacent to the active sites of the two PncAs. The structure comparison substantiates that this region should be very important for the ligand's binding in all PncAs. Additionally, MD simulations also show that the coordination position of the ligand is displaced by one water molecule in the unliganded enzymes. These results could provide the more penetrating understanding of drug resistance of M. tuberculosis and be helpful for the development of new antituberculosis drugs.     

Molecular dynamics simulations of ethanol binding to the transmembrane domain of the glycine receptor: implications for the channel potentiation mechanism

The glycine receptor (GlyR) is potentiated by ethanol and other anesthetics. The potentiation mechanism at the molecular level is unknown and remains elusive, but mutagenic studies have shown that ethanol and other volatile anesthetics bind to a pocket between TM1, TM2, and TM3. The present study extends previous studies (Cheng et al., Proteins 2007;68:581-593) wherein we conducted homology modeling and molecular dynamics (MD) simulations to construct models of the homopentameric alpha1 subunits of the GlyR transmembrane domain in open and closed states. To understand the potentiation of GlyR by ethanol we compare the binding of ethanol molecules to the channel in these different states. We observe that ethanol stably resides inside solvent-accessible cavities found in the open state of GlyR that are formed by I229 (of TM1) in one subunit and S267 and A288 (of TM2 and TM3, respectively) in the adjacent subunit. The volume of these putative binding pockets is state-dependent. Selective binding to the open states of receptors has been proposed to explain the potentiating actions of this class of anesthetics. In accordance with this model, our MD simulations suggest that the potentiation of ethanol on GlyR may be effected through preferential binding of ethanol molecules to an inter-subunit binding pocket in the open state.     

The third leg: Molecular dynamics simulations of lipid binding proteins

Molecular dynamics computer simulations can provide a third leg which balances the contributions of both structural biology and binding studies performed on the lipid binding protein family. In this context, these calculations help to establish a dialogue between all three communities, by relating experimental observables with details of structure. Working towards this connection is important, since experience has shown the difficulty of inferring thermodynamic properties from a single static conformation. The challenge is exemplified by ongoing attempts to interpret the impact of mutagenesis on structure and function (i.e. binding). A detailed atomic-level understanding of this system could be achieved with the support of all three legs, paving the way towards rational design of proteins with novel specificities. This paper provides an outline of the connections possible between experiment and theory concerning lipid binding proteins.     

Molecular dynamics simulations of human LRH-1: the impact of ligand binding in a constitutively active nuclear receptor

The liver receptor homologue 1 (LRH-1 (NR5A2)) belongs to the orphan nuclear receptor family, indicating that initially no ligand was known. Although recent studies have shown that ligand binding can be obtained, the biological relevance remains elusive. Here, we modify the observed X-ray ligand into a biologically more significant phospholipid (phosphatidylserine, PS) present in human, to study, by molecular dynamics (MD) simulations, the impact of the ligand on the receptor and the interaction with different cofactor peptides. Furthermore, we characterize the interactions between receptor and the cofactor peptides of DAX-1 (NR0B1), Prox1 and SHP LXXLL box 1 and 2 (NR0B2) in terms of specificity. Our MD simulation results show different interaction patterns for the SHP box2 compared to DAX-1, PROX1 and SHP box1. SHP box2 shows specific interactions at its more C-terminal end while the other investigated peptides show specific interactions at several positions but particularly at the +2 site. The peptide +2 side chain interacts with a charged amino acid of the receptor, in hLRH-1 Asp372. Together with the charge clamp residues Arg361 and Glu534, Asp372 forms a triangle shaped charge clamp responsible for peptide orientation and increased affinity. The binding of the PS ligand causes no overall structural changes of the receptor but affects the interactions with cofactor peptides. The cofactor peptides from SHP decrease its interaction with the receptor upon ligand binding while DAX-1 and PROX1 are unchanged or increase. The diverse ligand binding response of the cofactor provides an opportunity for drug design with the possibility to create agonist ligands to modify cofactor interaction.     

Molecular dynamics simulations of the ligand-binding domain of an N-methyl-D-aspartate receptor

The mechanism of partial agonism at N-methyl-D-aspartate receptors is an unresolved issue, especially with respect to the role of protein dynamics. We have performed multiple molecular dynamics simulations (7 x 20 ns) to examine the behavior of the ligand-binding core of the NR1 subunit with a series of ligands. Our results show that water plays an important role in stabilizing different conformations of the core and how a closed cleft conformation of the protein might be stabilized in the absence of ligands. In the case of ligand-bound simulations with both full and partial agonists, we observed that ligands within the binding cleft may undergo distinct conformational changes, without grossly influencing the degree of cleft closure within the ligand-binding domain. In agreement with recently published crystallographic data, we also observe similar changes in backbone torsions corresponding to the hinge region between the two lobes for the partial agonist, D-cycloserine. This observation rationalizes the classification of D-cycloserine as a partial agonist and should provide a basis with which to predict partial agonism in this class of receptor by analyzing the behavior of these torsions with other potential ligands.     

Catalyzed decomposition of urea. Molecular dynamics simulations of the binding of urea to urease

We present the results of molecular dynamics simulations on the urea/urease system. The starting structure was prepared from the 2.0 A crystal structure of Benini et al. [(1999) Struct. Folding Des. 7, 205-216] of DAP-inhibited urease (PDB code ), and the trimeric structure (2479 residues) resulted in 180K atoms after solvation by water. The force field parameters were derived using the bonded model approach described by Hoops et al. [(1991) J. Am. Chem. Soc. 113, 8262-8270]. Three different systems were analyzed, each one modeling a different protonation pattern for the His320 and His219 residues. In each case, the three monomers of urease have been analyzed separately. The time-averaged structures observed in the three monomers suggest that urease could follow two different competitive mechanisms. A "protein-assisted proton transfer" mechanism points to Asp221 as crucial for catalysis. An "Asp-mediated proton transfer" involves the transfer of a proton from the bridging OH to an NH2 moiety of urea, assisted by Asp360 in the active site. The impact of the simulation results on our understanding of urease catalysis is discussed in detail.     

Molecular dynamics simulations of matrix metalloproteinase 2: role of the structural metal ions

Herein we investigate the role played by the so-called "structural metal ions" in the catalytic domain of the matrix metalloproteinase 2 enzyme (MMP-2 or gelatinase A). We performed seven molecular dynamics simulations that differ in the number and position of the noncatalytic zinc and calcium ions bound to the MMP-2 catalytic domain. An additional simulation including the three fibronectin-type modules inserted into the catalytic domain was also carried out. The analysis of the trajectories confirms that the binding/removal of the structural ions does not perturb the secondary structure elements but influences the position of several solvent-exposed loop regions that are placed near the active site cleft. The position of these loops modulates the accessibility of important anchorage points for substrate binding that have been identified in the active site groove. On the basis of semiempirical quantum chemical calculations, we estimated the relative free energies of the MMP-2 models, obtaining thus that the binding of two zinc and two calcium ions to the MMP-2 catalytic domain is energetically favored. In this MMP-2 model, which shows the most compact structure, all of the substrate binding sites are readily accessible. Globally, our results help to rationalize at the atomic level the calcium and zinc dependence of the hydrolytic activity catalyzed by the MMPs.     

Molecular dynamics simulations of N-terminal peptides from a nucleotide binding protein

Molecular dynamics (MD) simulations of N-terminal peptides from lactate dehydrogenase (LDH) with increasing length and individual secondary structure elements were used to study their stability in relation to folding. Ten simulations of 1–2 ns of different peptides in water starting from the coordinates of the crystal structure were performed. The stability of the peptides was compared qualitatively by analyzing the root mean square deviation (RMSD) from the crystal structure, radius of gyration, secondary and tertiary structure, and solvent accessible surface area. In agreement with earlier MD studies, relatively short (< 15 amino acids) peptides containing individual secondary structure elements were generally found to be unstable; the hydrophobic α₁-helix of the nucleotide binding fold displayed a significantly higher stability, however. Our simulations further showed that the first βαβ supersecondary unit of the characteristic dinucleotide binding fold (Rossmann fold) of LDH is somewhat more stable than other units of similar length and that the α₂-helix, which unfolds by itself, is stabilized by binding to this unit. This finding suggests that the first βαβ unit could function as an N-terminal folding nucleus, upon which the remainder of the polypeptide chain can be assembled. Indeed, simulations with longer units (βαβα and βαβαββ) showed that all structural elements of these units are rather stable. The outcome of our studies is in line with suggestions that folding of the N-terminal portion of LDH in vivo can be a cotranslational process that takes place during the ribosomal peptide synthesis.     

Molecular dynamics simulations of the bacterial periplasmic heme binding proteins ShuT and PhuT

ShuT and PhuT are two periplasmic heme binding proteins that shuttle heme between the outer and inner membranes of the Gram-negative bacteria. Periplasmic binding proteins (PBPs) generally exhibit considerable conformational changes during the ligand binding process, whereas ShuT and PhuT belong to a class of PBPs that do not show such behavior based on their apo and holo crystal structures. By employing a series of molecular dynamic simulations on the ShuT and the PhuT, the dynamics and functions of the two PBPs were investigated. Through monitoring the distance changes between the two conserved glutamates of ShuT and PhuT, it was found the two PBPs were more flexible than previously assumed, exhibiting obvious opening-closing motions which were more remarkable in the apo runs of ShuT. Based on the results of the domain motion analysis, large scale conformational transitions were found in all apo runs of ShuT and PhuT, hinting that the domain motions of the two PBPs may be intrinsic. On the basis of the results of the principle component analysis, distinct opening-closing and twisting motion tendencies were observed not only in the apo, but also in the holo simulations of the two PBPs. The Gaussian network model was applied in order to analyze the hinge bending regions. The most important bending regions of ShuT and PhuT are located around the midpoints of their respective connecting helixes. Finally, the flexibilities and the details of the simulations of ShuT and PhuT were discussed. Characterized by the remarkably large flexibilities, the loop constituted by Ala 169, Gly170 and Gly171 of ShuT and the beta-turn constituted by Ala176, Gly177 and Gly178 of PhuT may be important for the functions of the two PBPs. Furthermore, the Asn254 of ShuT and the Arg228 of PhuT may be indispensable for the binding or unbinding of heme, since it is involved in the important hydrogen bonding to the propionate side-chains of heme.     

Molecular dynamics simulations of the amyloid-beta binding alcohol dehydrogenase (ABAD) enzyme

In this work, we present 10 ns molecular dynamics simulations of the homotetramer of the ABAD enzyme, as well as of the structural units, dimer and monomer, that assemble to form the tetramer, in the presence and absence of a NAD-inhibitor adduct. The aim was to compare the stability of the different structures and to study the effects of the inhibitor binding on the flexibility of the enzyme structure. The results indicate that the tetramer, dimer and monomer show a comparable stability and that tetramerization stabilizes some regions of the protein that when exposed to the solvent in dimer and monomer become more flexible. Binding of the cofactor and inhibitor stabilizes the protein, the main effect being a stabilization of the substrate binding loop. In the absence of the ligand, this region was found to have a much higher flexibility and to adopt an open conformation. An interesting result emerging from this work is the conformational flexibility exhibited by the azepane and benzene rings of the inhibitor moiety of the adduct, which appears to be influenced by the mobility of the substrate binding loop. This highlights the importance of integrate the flexibility of the substrate binding loop into de novo design of inhibitors of ABAD.     

Interdomain dynamics and ligand binding: molecular dynamics simulations of glutamine binding protein

Periplasmic binding proteins from Gram-negative bacteria possess a common architecture, comprised of two domains linked by a hinge region, a fold which they share with the neurotransmitter-binding domains of ionotropic glutamate receptors (GluRs). Glutamine-binding protein (GlnBP) is one such protein, whose crystal structure has been solved in both open and closed forms. Multi-nanosecond molecular dynamics simulations have been used to explore motions about the hinge region and how they are altered by ligand binding. Glutamine binding is seen to significantly reduce inter-domain motions about the hinge region. Essential dynamics analysis of inter-domain motion revealed the presence of both hinge-bending and twisting motions, as has been reported for a related sugar-binding protein. Significantly, the influence of the ligand on GlnBP dynamics is similar to that previously observed in simulations of rat glutamate receptor (GluR2) ligand-binding domain. The essential dynamics analysis of GlnBP also revealed a third class of motion which suggests a mechanism for signal transmission in GluRs.     

Molecular dynamics simulations of xDNA

xDNA is a modified DNA, which contains natural as well as expanded bases. Expanded bases are generated by the addition of a benzene spacer to the natural bases. A set of AMBER force‐field parameters were derived for the expanded bases and the structural dynamics of the xDNA decamer ( xT5 ′ G xT A xC xG C xA xG T3′ ) · ( xA5′ C T xG C G xT A xC A3′) was explored using a 22 ns molecular dynamics simulation in explicit solvent. During the simulation, the duplex retained its Watson‐Crick base‐pairing and double helical structure, with deviations from the starting B‐form geometry towards A‐form; the deviations are mainly in the backbone torsion angles and in the helical parameters. The sugar pucker of the residues were distributed among a variety of modes; C2′ endo, C1′ exo, O4′ endo, C4′ exo, C2′ exo, and C3′ endo. The enhanced stacking interactions on account of the modification in the bases could help to retain the duplex nature of the helix with minor deviations from the ideal geometry. In our simulation, the xDNA showed a reduced minor groove width and an enlarged major groove width in comparison with the NMR structure. Both the grooves are larger than that of standard B‐DNA, but major groove width is larger than that of A‐DNA with almost equal minor groove width. The enlarged groove widths and the possibility of additional hydration in the grooves makes xDNA a potential molecule for various applications. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 351–360, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

Molecular dynamics simulations of carrabiose

Molecular mechanics calculations have been performed for the disaccharide carrabiose, one of the repeat units of β-carrageenan, as a general model for the (1→4)-linkage in the carrageenans. An adiabatic conformational energy map for this unsulfated molecule was prepared by constrained energy minimization and compared to a previously reported rigid-residue energy map for the sulfated molecule and to a similar adiabatic map for neocarrabiose, the related (1→3)-linked dimer repeat unit of β-carrageenan. Molecular dynamics simulations of this molecule in vacuo and in an aqueous (TIP3P) solution were calculated, and the observed motions were found to be generally consistent with the vacuum adiabatic energy map. Unlike the case observed in previous simulations of neocarrabiose, little salvation shift in the molecular conformation was observed for carrabiose. From the dynamics, the linkage was observed to be relatively flexible, as has been inferred from experiment on sulfated carrageenan polymers. © 1997 John Wiley & Sons, Inc.     

Molecular dynamics simulations of metalloproteins

Molecular dynamics simulations are now commonly applied to metalloproteins, despite the challenges introduced by the presence of metal ions. Force field parameters are nowadays available also for these 'exotic' atoms and several biological systems have been successfully studied. Some of the most relevant results and methodological advancements are reviewed.