Molecular docking of balanol to dynamics snapshots of protein kinase A

Even if the structure of a receptor has been determined experimentally, it may not be a conformation to which a ligand would bind when induced fit effects are significant. Molecular docking using such a receptor structure may thus fail to recognize a ligand to which the receptor can bind with reasonable affinity. Here, we examine one way to alleviate this problem by using an ensemble of receptor conformations generated from a molecular dynamics simulation for molecular docking. Two molecular dynamics simulations were conducted to generate snapshots for protein kinase A: one with the ligand bound, the other without. The ligand, balanol, was then docked to conformations of the receptors presented by these trajectories. The Lamarckian genetic algorithm in Autodock [Goodsell et al. J Mol Recognit 1996;9(1):1-5; Morris et al. J Comput Chem 1998;19(14):1639-1662] was used in the docking. Three ligand models were used: rigid, flexible, and flexible with torsional potentials. When the snapshots were taken from the molecular dynamics simulation of the protein-ligand complex, the correct docking structure could be recovered easily by the docking algorithm in all cases. This was an easier case for challenging the docking algorithm because, by using the structure of the protein in a protein-ligand complex, one essentially assumed that the protein already had a pocket to which the ligand can fit well. However, when the snapshots were taken from the ligand-free protein simulation, which is more useful for a practical application when the structure of the protein-ligand complex is not known, several clusters of structures were found. Of the 10 docking runs for each snapshot, at least one structure was close to the correctly docked structure when the flexible-ligand models were used. We found that a useful way to identify the correctly docked structure was to locate the structure that appeared most frequently as the lowest energy structure in the docking experiments to different snapshots.     

Ligand binding and protein relaxation in heme proteins: a room temperature analysis of NO geminate recombination

Ultrafast absorption spectroscopy is used to study heme-NO recombination at room temperature in aqueous buffer on time scales where the ligand cannot leave its cage environment. While a single barrier is observed for the cage recombination of NO with heme in the absence of globin, recombination in hemoglobin and myoglobin is nonexponential. Examination of hemoglobin with and without inositol hexaphosphate points to proximal constraints as important determinants of the geminate rebinding kinetics. Molecular dynamics simulations of myoglobin and heme-imidazole subsequent to ligand dissociation were used to investigate the transient behavior of the Fe-proximal histidine coordinate and its possible involvement in geminate recombination. The calculations, in the context of the absorption measurements, are used to formulate a distinction between nonexponential rebinding that results from multiple protein conformations (substates) present at equilibrium or from nonequilibrium relaxation of the protein triggered by a perturbation such as ligand dissociation. The importance of these two processes is expected to depend on the time scale of rebinding relative to equilibrium fluctuations and nonequilibrium relaxation. Since NO rebinding occurs on the picosecond time scale of the calculated myoglobin relaxation, a time-dependent barrier is likely to be an important factor in the observed nonexponential kinetics. The general implications of the present results for ligand binding in heme proteins and its time and temperature dependence are discussed. It appears likely that, at low temperatures, inhomogeneous protein populations play an important role and that as the temperature is raised, relaxation effects become significant as well.     

Analysis of Delta-Notch interaction by molecular modeling and molecular dynamic simulation studies

Notch is a single-pass transmembrane receptor protein. Delta (member of the DSL protein family), a Notch ligand, is also single-pass transmembrane protein that can interact with Notch to form the Delta-Notch complex. It has been demonstrated that of the 36 Epidermal Growth Factor (EGF) repeats of Notch, 11th and 12th are sufficient to mediate interactions with Delta. Crystal structure of mammalian Notch1 extracellular ligand binding domain shows the presence of 11th and 12th EGF-like repeats. Here a portion of the Drosophila Delta protein, known to interact with Notch extracellular domain, has been modeled using homology modeling. The structure of the Delta-Notch complex was subsequently modeled by protein-docking method using GRAMM. Molecular dynamic simulations of the modeled structures were performed. The probable structures for Delta-Notch complex have been proposed based on interaction energy parameter and planarity studies.     

Exploration of multi-state conformational dynamics and underlying global functional landscape of maltose binding protein

An increasing number of biological machines have been revealed to have more than two macroscopic states. Quantifying the underlying multiple-basin functional landscape is essential for understanding their functions. However, the present models seem to be insufficient to describe such multiple-state systems. To meet this challenge, we have developed a coarse grained triple-basin structure-based model with implicit ligand. Based on our model, the constructed functional landscape is sufficiently sampled by the brute-force molecular dynamics simulation. We explored maltose-binding protein (MBP) which undergoes large-scale domain motion between open, apo-closed (partially closed) and holo-closed (fully closed) states responding to ligand binding. We revealed an underlying mechanism whereby major induced fit and minor population shift pathways co-exist by quantitative flux analysis. We found that the hinge regions play an important role in the functional dynamics as well as that increases in its flexibility promote population shifts. This finding provides a theoretical explanation of the mechanistic discrepancies in PBP protein family. We also found a functional "backtracking" behavior that favors conformational change. We further explored the underlying folding landscape in response to ligand binding. Consistent with earlier experimental findings, the presence of ligand increases the cooperativity and stability of MBP. This work provides the first study to explore the folding dynamics and functional dynamics under the same theoretical framework using our triple-basin functional model.     

Ultrafast dynamics of ligands within heme proteins

Physiological bond formation and bond breaking events between proteins and ligands and their immediate consequences are difficult to synchronize and study in general. However, diatomic ligands can be photodissociated from heme, and thus in heme proteins ligand release and rebinding dynamics and trajectories have been studied on timescales of the internal vibrations of the protein that drive many biochemical reactions, and longer. The rapidly expanding number of characterized heme proteins involved in a large variety of functions allows comparative dynamics-structure-function studies. In this review, an overview is given of recent progress in this field, and in particular on initial sensing processes in signaling proteins, and on ligand and electron transfer dynamics in oxidases and cytochromes.     

Enhanced sampling of protein conformational states for dynamic cross-docking within the protein-protein docking server SwarmDock

The formation of specific protein-protein interactions is often a key to a protein's function. During complex formation, each protein component will undergo a change in the conformational state, for some these changes are relatively small and reside primarily at the sidechain level; however, others may display notable backbone adjustments. One of the classic problems in the protein-docking field is to be able to a priori predict the extent of such conformational changes. In this work, we investigated three protocols to find the most suitable input structure conformations for cross-docking, including a robust sampling approach in normal mode space. Counterintuitively, knowledge of the theoretically best combination of normal modes for unbound-bound transitions does not always lead to the best results. We used a novel spatial partitioning library, Aether Engine (see Supplementary Materials ), to efficiently search the conformational states of 56 receptor/ligand pairs, including a recent CAPRI target, in a systematic manner and selected diverse conformations as input to our automated docking server, SwarmDock, a server that allows moderate conformational adjustments during the docking process. In essence, here we present a dynamic cross-docking protocol, which when benchmarked against the simpler approach of just docking the unbound components shows a 10% uplift in the quality of the top docking pose.     

Dynamical differences of hemoglobin and the ionotropic glutamate receptor in different states revealed by a new dynamics alignment method

A new algorithm for comparison of protein dynamics is presented. Compared protein structures are superposed and their modes of motions are calculated using the anisotropic network model. The obtained modes are aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. Dynamical comparison of hemoglobin in the T and R2 states reveals that the dynamics of the allosteric effector 2,3‐bisphosphoglycerate binding site is different in the two states. These differences can contribute to the selectivity of the effector to the T state. Similar comparison of the ionotropic glutamate receptor in the kainate+(R,R)‐2b and ZK bound states reveals that the kainate+(R,R)‐2b bound states slow modes describe upward motions of ligand binding domain and the transmembrane domain regions. Such motions may lead to the opening of the receptor. The upper lobes of the LBDs of the ZK bound state have a smaller interface with the amino terminal domains above them and have a better ability to move together. The present study exemplifies the use of dynamics comparison as a tool to study protein function. Proteins 2017; 85:1507–1517. © 2014 Wiley Periodicals, Inc.     

Local segmental dynamics of cis-1,4-polybutadiene,polypropylene and polyethylene terephthalate via molecular dynamics simulations

The local segmental dynamics of cis-1,4-polybutadiene, polypropylene and polyethylene terephthalate have been investigated via isothermal-isobaric molecular dynamics simulations. The simulation pressure was 1 atm for all systems, with all simulation temperatures being at least 150 K above the polymer's glass transition temperature. The trajectories have been analysed via autocorrelation functions (ACFs) of chord vectors spanning different number of chain backbone bonds. Inverse Laplace transformations of these ACFs using the CONTIN algorithm afforded the corresponding distribution of relaxation times (DRTs) for the simulated dynamics. All DRTs illustrated a peak on fast timescales corresponding to short length scale segmental motion and a peak at longer timescales corresponding to longer length scale relaxations. A third peak, intermediate between the fast and slow processes, appears as the relaxation of chord vectors spanning increasing number of backbone bonds is considered. The temperature dependence of the relaxation dynamics is also investigated.     

Computer simulation of protein molecular dynamics

A review of the works on the computer simulation of the globular protein dynamics is given. Methodological aspects of the simulation procedure are outlined briefly. Main peculiarities of protein dynamics revealed in the course of simulation of pancreatic trypsin inhibitor and cytochrome c are presented. The causes of "anomalous" processes, inherent in the simulated behaviour of model proteins are discussed. These "anomalous" processes are: continuous drift of the structure and its deviation from the experimental one, determined by X-ray analysis. Both processes are supposed to be the consequence of the reduced conformational rigidity of the model protein in comparison to the real one. Among the possible reasons for this reduced rigidity absence of the water molecules, hydrating peptide groups in the real protein, may be mentioned. Analogy between "anomalous" processes in the simulated protein dynamics and some phenomena observed in the real proteins during their functioning is drawn.     

Order changes within receptor systems upon ligand binding: receptor tightening/oligomerisation and the interpretation of binding parameters

Recent hydrogen-deuterium exchange experiments have highlighted tightening and loosening of protein structures upon ligand binding, with changes in bonding (DeltaH) and order (DeltaS) which contribute to the overall thermodynamics of ligand binding. Tightening and loosening show that ligand binding respectively stabilises or destabilises the internal structure of the protein, i.e. it shows positive or negative cooperativity between ligand binding and the receptor structure. In the case of membrane-bound receptors, such as G protein-coupled receptors (GPCRs) and ligand gated ion channel receptors (LGICRs), most binding studies have focussed on association/dissociation constants. Where these have been broken down into enthalpic and entropic contributions, the phenomenon of "thermodynamic discrimination" between antagonists and agonists has often been noted; e.g. for a receptor where agonist binding is predominantly enthalpy driven, antagonist binding is predominantly entropy driven and vice versa. These data have not previously been considered in terms of the tightening, or loosening, of receptor structures that respectively occurs upon positively, or negatively, cooperative binding of ligand. Nor have they been considered in light of the homo- and hetero-oligomerisation of GPCRs and the possibility of ligand-induced changes in oligomerisation. Here, we argue that analysis of the DeltaH and DeltaS of ligand binding may give useful information on ligand-induced changes in membrane-bound receptor oligomers, relevant to the differing effects of agonists and antagonists.     

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.     

Thermodynamic interpretation of protein dynamics from NMR relaxation measurements

Protein dynamics and thermodynamics can be characterized through measurements of relaxation rates of side chain (2)H and (13)C, and backbone (15)N nuclei using NMR spectroscopy. The rates reflect protein motions on timescales from picoseconds to milliseconds. Backbone and methyl side chain NMR relaxation measurements for several proteins are beginning to reveal the role of protein dynamics in protein stability and ligand binding.     

Millisecond Timescale Dynamics of Human Liver Fatty Acid Binding Protein: Testing of Its Relevance to the Ligand Entry Process

For over a decade, scientists have been attempting to know more about the conformational dynamics of fatty acid binding proteins (FABPs), to answer the puzzling question of how ligands could access the internalized binding site(s). Conformational exchange of FABPs on the microsecond to millisecond timescales has been found in many FABPs and offers an important hypothesis for the ligand entry mechanism. Despite the potential significance, the validity of this hypothesis has not been verified yet. In this study, the slow dynamics of human liver fatty acid binding protein (hLFABP) that was shown previously to be highly flexible on millisecond timescales was quantitatively characterized in detail. In addition, the interaction between hLFABP and 1,8-ANS was studied using NMR spectroscopy, and the kinetic rate of ANS association to hLFABP was measured. We believe the current result excludes the possibility that the intrinsic millisecond dynamics of hLFABP represents a critical conformational reorganization process required for ligand entry, but implies that it may represent the exchange between the apo-state and a state resembling the singly-bound conformation. Furthermore, we suggest these results show that the ligand-entry related functional dynamics could occur on the microsecond/submicrosecond timescales, highly encouraging future computational studies on this topic.     

Backbone dynamics of the human CC chemokine eotaxin: fast motions, slow motions, and implications for receptor binding.

Eotaxin is a member of the chemokine family of about 40 proteins that induce cell migration. Eotaxin binds the CC chemokine receptor CCR3 that is highly expressed by eosinophils, and it is considered important in the pathology of chronic respiratory disorders such as asthma. The high resolution structure of eotaxin is known. The 74 amino acid protein has two disulfide bridges and shows a typical chemokine fold comprised of a core of three antiparallel beta-strands and an overlying alpha-helix. In this paper, we report the backbone dynamics of eotaxin determined through 15N-T1, T2, and [1H]-15N nuclear Overhauser effect heteronuclear multidimensional NMR experiments. This is the first extensive study of the dynamics of a chemokine derived from 600, 500, and 300 MHz NMR field strengths. From the T1, T2, and NOE relaxation data, parameters that describe the internal motions of eotaxin were derived using the Lipari-Szabo model free analysis. The most ordered regions of the protein correspond to the known secondary structure elements. However, surrounding the core, the regions known to be functionally important in chemokines show a range of motions on varying timescales. These include extensive subnanosecond to picosecond motions in the N-terminus, C-terminus, and the N-loop succeeding the disulfides. Analysis of rotational diffusion anisotropy of eotaxin and chemical exchange terms at multiple fields also allowed the confident identification of slow conformational exchange through the "30s" loop, disulfides, and adjacent residues. In addition, we show that these motions may be attenuated in the dimeric form of a synthetic eotaxin. The structure and dynamical basis for eotaxin receptor binding is discussed in light of the dynamics data.     

Myoglobin cavities provide interior ligand pathway

The myoglobin protein binds oxygen and catalyzes NO oxidation. As a key model protein, its dynamics have been well studied by spectroscopy and by crystallography as well as by simulation. Nonetheless, visualization of the mechanism of movement of ligands within myoglobin has been difficult. Coordinates of the A1 and A3 taxonomic spectral states of myoglobin from the 1 A crystal structure (1a6g) are generated as consistent sets of correlated clusters of residues with A or B crystal alternates. Analysis of cavities in these A1 and A3 conformations clarifies the pathway of ligand motion from distal entry through interior movement to the proximal side of the heme. Cavities opened up by buried alternate conformations link the distal to the proximal side of the heme. Structural conservation highlights the relevance of this pathway to human neuroglobin. Cavity migration via myoglobin crystal alternates provides a specific link of protein structure to protein dynamics and protein function and demonstrates the relevance of substates (discrete disorder) to function for all proteins.     

The leghemoglobin proximal heme pocket directs oxygen dissociation and stabilizes bound heme

Sperm whale myoglobin (Mb) and soybean leghemoglobin (Lba) are two small, monomeric hemoglobins that share a common globin fold but differ widely in many other aspects. Lba has a much higher affinity for most ligands, and the two proteins use different distal and proximal heme pocket regulatory mechanisms to control ligand binding. Removal of the constraint provided by covalent attachment of the proximal histidine to the F-helices of these proteins decreases oxygen affinity in Lba and increases oxygen affinity in Mb, mainly because of changes in oxygen dissociation rate constants. Hence, Mb and Lba use covalent constraints in opposite ways to regulate ligand binding. Swapping the F-helices of the two proteins brings about similar effects, highlighting the importance of this helix in proximal heme pocket regulation of ligand binding. The F7 residue in Mb is capable of weaving a hydrogen-bonding network that holds the proximal histidine in a fixed orientation. On the contrary, the F7 residue in Lba lacks this property and allows the proximal histidine to assume a conformation favorable for higher ligand binding affinity. Geminate recombination studies indicate that heme iron reactivity on picosecond timescales is not the dominant cause for the effects observed in each mutation. Results also indicate that in Lba the proximal and distal pocket mutations probably influence ligand binding independently. These results are discussed in the context of current hypotheses for proximal heme pocket structure and function.     

Dynamic instabilities in the microtubule cytoskeleton: A state diagram

The dynamics of microtubule growth and disassembly is considered in the framework of the theory of nonequilibrium reaction-diffusion systems. The phase diagram contains regions corresponding to stable stationary and nonstationary solutions. Dynamic instabilities can arise from nonequilibrium kinetic transitions. Agents affecting the microtubule dynamics are classed into four types, and the interplay of their effects is analyzed.     

Characterizing and controlling the inherent dynamics of cyclophilin-A

With the recent advances in NMR relaxation techniques, protein motions on functionally important timescales can be studied at atomic resolution. Here, we have used NMR-based relaxation experiments at several temperatures and both 600 and 900 MHz to characterize the inherent dynamics of the enzyme cyclophilin-A (CypA). We have discovered multiple chemical exchange processes within the enzyme that form a “dynamic continuum” that spans 20–30 Å comprising active site residues and residues proximal to the active site. By combining mutagenesis with these NMR relaxation techniques, a simple method of counting the dynamically sampled conformations has been developed. Surprisingly, a combination of point mutations has allowed for the specific regulation of many of the exchange processes that occur within CypA, suggesting that the dynamics of an enzyme may be engineered.     

Quantitative analysis of ligand migration from transition networks

In this work we use transition network analysis for the first time to investigate ligand migration in truncated hemoglobin (trHbN) and obtain kinetic information about the docking-site dynamics in the protein. A comparison with explicit water molecular dynamics simulations (100 ns in total) shows that the rate constants derived from the network analysis are realistic. The transition network analysis provides 1) The time-resolved connectivity network in the protein; 2) The half-lives of the docking sites; 3) The transition timescales between two given docking sites; and 4) The extent of population transfer among different docking sites of the protein as a function of lag time. We investigate the role of the Tyr³³ and Gln⁵⁸ residues in ligand migration by studying ligand migration in four mutants of trHbN. The mutation study suggests that residues Tyr³³ and Gln⁵⁸ stabilize the NO ligand in the Xe2 docking site of trHbN, thus facilitating the efficiency of the NO detoxification reaction.