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     

Tackling the conformational sampling of larger flexible compounds and macrocycles in pharmacology and drug discovery

Computational conformational sampling underpins much of molecular modeling and design in pharmaceutical work. The sampling of smaller drug-like compounds has been an active area of research. However, few studies have tested in details the sampling of larger more flexible compounds, which are also relevant to drug discovery, including therapeutic peptides, macrocycles, and inhibitors of protein–protein interactions. Here, we investigate extensively mainstream conformational sampling methods on three carefully curated compound sets, namely the ‘Drug-like’, larger ‘Flexible’, and ‘Macrocycle’ compounds. These test molecules are chemically diverse with reliable X-ray protein-bound bioactive structures. The compared sampling methods include Stochastic Search and the recent LowModeMD from MOE, all the low-mode based approaches from MacroModel, and MD/LLMOD recently developed for macrocycles. In addition to default settings, key parameters of the sampling protocols were explored. The performance of the computational protocols was assessed via (i) the reproduction of the X-ray bioactive structures, (ii) the size, coverage and diversity of the output conformational ensembles, (iii) the compactness/extendedness of the conformers, and (iv) the ability to locate the global energy minimum. The influence of the stochastic nature of the searches on the results was also examined. Much better results were obtained by adopting search parameters enhanced over the default settings, while maintaining computational tractability. In MOE, the recent LowModeMD emerged as the method of choice. Mixed torsional/low-mode from MacroModel performed as well as LowModeMD, and MD/LLMOD performed well for macrocycles. The low-mode based approaches yielded very encouraging results with the flexible and macrocycle sets. Thus, one can productively tackle the computational conformational search of larger flexible compounds for drug discovery, including macrocycles.     

Envisioning the Loop Movements and Rotation of the Two Subdomains of Dihydrofolate Reductase by Elastic Normal Mode Analysis

Abstract

Normal mode analysis, using the elastic network model, has been employed to envision the low frequency normal mode motion trends in the structures of five intermediates and a transition state in the kinetic pathway of E. coli dihydrofolate reductase (DHFR). Five of the reaction pathway analog structures and a crystal structure resembling the transition state, using X-ray analyses determined by Kraut et al., have been adapted as structural models. The motions that poise pathways of the M20 loop transitions from closed to occluded conformations and sub domain rotation to close the substrate cleft, have been predicted and envisioned for the first time by this study. Pathway entries to the movement of the substrate binding cleft helices are also envisioned. These motions play roles in transition structure stabilization and in regulating the release of the product tetrahydrofolate (THF). The motions observed push the ground state conformation of each intermediate towards a higher energy sub state conformation. A set of conserved residues involved in the catalytic reactions and conformational changes, previously studied by kinetic, theoretical and NMR, have been analyzed. The importance of these motions in terms of protein dynamics are revealed and envisioned by the normal mode analysis. Additional residues are proposed as candidates for further study of their potential promotional function.     

Protein-inhibitor flexible docking by a multicanonical sampling: native complex structure with the lowest free energy and a free-energy barrier distinguishing the native complex from the others

Flexible docking between a protein (lysozyme) and an inhibitor (tri-N-acetyl-D-glucosamine, tri-NAG) was carried out by an enhanced conformational sampling method, multicanonical molecular dynamics simulation. We used a flexible all-atom model to express lysozyme, tri-NAG, and water molecules surrounding the two bio-molecules. The advantages of this sampling method are as follows: the conformation of system is widely sampled without trapping at energy minima, a thermally equilibrated conformational ensemble at an arbitrary temperature can be reconstructed from the simulation trajectory, and the thermodynamic weight can be assigned to each sampled conformation. During the simulation, exchanges between the binding and free (i.e., unbinding) states of the protein and the inhibitor were repeatedly observed. The conformational ensemble reconstructed at 300 K involved various conformational clusters. The main outcome of the current study is that the most populated conformational cluster (i.e., the cluster of the lowest free energy) was assigned to the native complex structure (i.e., the X-ray complex structure). The simulation also produced non-native complex structures, where the protein and the inhibitor bound with different modes from that of the native complex structure, as well as the unbinding structures. A free-energy barrier (i.e., activation free energy) was clearly detected between the native complex structures and the other structures. The thermal fluctuations of tri-NAG in the lowest free-energy complex correlated well with the X-ray B-factors of tri-NAG in the X-ray complex structure. The existence of the free-energy barrier ensures that the lowest free-energy structure can be discriminated naturally from the other structures. In other words, the multicanonical molecular dynamics simulation can predict the native complex structure without any empirical objective function. The current study also manifested that the flexible all-atom model and the physico-chemically defined atomic-level force field can reproduce the native complex structure. A drawback of the current method is that it requires a time consuming computation due to the exhaustive conformational sampling. We discussed a possibility for combining the current method with conventional docking methods.     

Conformational equilibria and free energy profiles for the allosteric transition of the ribose-binding protein

The ribose-binding protein (RBP) is a sugar-binding bacterial periplasmic protein whose function is associated with a large allosteric conformational change from an open to a closed conformation upon binding to ribose. The crystal structures of RBP in open and closed conformations have been solved. It has been hypothesized that the open and closed conformations exist in a dynamic equilibrium in solution, and that sugar binding shifts the population from open conformations to closed conformations. Here, we study by computer simulations the thermodynamic changes that accompany this conformational change, and model the structural changes that accompany the allosteric transition, using umbrella sampling molecular dynamics and the weighted histogram analysis method. The open state is comprised of a diverse ensemble of conformations; the open ribose-free X-ray crystal conformations being representative of this ensemble. The unligated open form of RBP is stabilized by conformational entropy. The simulations predict detectable populations of closed ribose-free conformations in solution. Additional interdomain hydrogen bonds stabilize this state. The predicted shift in equilibrium from the open to the closed state on binding to ribose is in agreement with experiments. This is driven by the energetic stabilization of the closed conformation due to ribose-protein interactions. We also observe a significant population of a hitherto unobserved ribose-bound partially open state. We believe that this state is the one that has been suggested to play a role in the transfer of ribose to the membrane-bound permease complex.     

Structural variability of C3larvin toxin. Intrinsic dynamics of the α/β fold of the C3-like group of mono-ADP-ribosyltransferase toxins

C3larvin toxin is a new member of the C3 class of the mono-ADP-ribosyltransferase toxin family. The C3 toxins are known to covalently modify small G-proteins, e.g. RhoA, impairing their function, and serving as virulence factors for an offending pathogen. A full-length X-ray structure of C3larvin (2.3 Å) revealed that the characteristic mixed α/β fold consists of a central β-core flanked by two helical regions. Topologically, the protein can be separated into N and C lobes, each formed by a β-sheet and an α-motif, and connected by exposed loops involved in the recognition, binding, and catalysis of the toxin/enzyme, i.e. the ADP-ribosylation turn–turn and phosphate–nicotinamide PN loops. Herein, we provide two new C3larvin X-ray structures and present a systematic study of the toxin dynamics by first analyzing the experimental variability of the X-ray data-set followed by contrasting those results with theoretical predictions based on Elastic Network Models (GNM and ANM). We identify residues that participate in the stability of the N-lobe, putative hinges at loop residues, and energy-favored deformation vectors compatible with conformational changes of the key loops and 3D-subdomains (N/C-lobes), among the X-ray structures. We analyze a larger ensemble of known C3bot1 conformations and conclude that the characteristic ‘crab-claw’ movement may be driven by the main intrinsic modes of motion. Finally, via computational simulations, we identify harmonic and anharmonic fluctuations that might define the C3larvin ‘native state.’ Implications for docking protocols are derived.     

Prediction of the native conformation of a polypeptide by a statistical-mechanical procedure. I. Backbone structure of enkephalin

A new methodology for theoretically predicting the native, three-dimensional structure of a polypeptide is presented. Based on equilibrium statistical mechanics, an algorithm has been designed to determine the probable conformation of a polypeptide by calculating conditional free-energy maps for each residue of the macromolecule. The conditional free-energy map of each residue is computed from a set of probability integrals, obtained by summing over the interaction energies of all pairs of nonbonded atoms of the whole molecule. By locating the region(s) of lowest free energy for each map, the probable conformation for each residue can be identified. The native structure of the polypeptide is assumed to be the combination of the probable conformations of the individual residues. All multidimensional probability integrals are evaluated by an adaptive Monte Carlo algorithm (SMAPPS —Statistical-Mechanical Algorithm for Predicting Protein Structure). The Monte Carlo algorithm searches the entire conformational space, adjusting itself automatically to concentrate its sampling in regions where the magnitude of the integrand is largest (“importance sampling”). No assumptions are made about the native conformation. The only prior knowledge necessary for the prediction of the native conformation is the amino acid sequence of the polypeptide. To test the effectiveness of the algorithm, SMAPPS was applied to the prediction of the native conformation of the backbone of Met-enkephalin, a pentapeptide. In the calculations, only the backbone dihedral angles (? and ψ) were allowed to vary; all side-chain (χ) and peptide-bond (ω) dihedral angles were kept fixed at the values corresponding to the alleged global minimum energy previously determined by direct energy minimization. For each conformation generated randomly by the Monte Carlo algorithm, the total conformational energy of the polypeptide was obtained from established empirical potential energy functions. Solvent effects were not included in the computations. With this initial application of SMAPPS , three distinct low-free-energy β-bend structures of Met-enkephalin were found. In particular, one of the structures has a conformation remarkably similar to the one associated with the previously alleged global minimum energy. The two additional structures of the pentapeptide have conformational energies lower than the previously computed low-energy structure. However, the Monte Carlo results are in agreement with an improved energy-minimization procedure. These initial results on the backbone structure of Met-enkephalin indicate that an equilibrium statistical-mechanical procedure, coupled with an adaptive Monte Carlo algorithm, can overcome many of the problems associated with the standard methods of direct energy minimization.     

Solid and solution state conformations of (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-allo-inositol and (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-6-O-methyl-allo-inositol

The synthesis and conformational studies of (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-allo-inositol and (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-6-O-methyl-allo-inositol are described. Solid state conformations of the title compounds have been studied by solving their X-ray crystal structures. The inositol ring in both the compounds deviate considerably from the ideal chair conformation to flattened chair conformation in the solid state. Their conformations in solution were studied by the use of 1H NMR spectroscopy. These conformational analyses revealed that the title compounds adopt similar conformations in solid and solution states irrespective of the solvent polarity.     

Reaction path and free energy calculations of the transition between alternate conformations of HIV-1 protease

Two different structures of ligand-free HIV protease have been determined by X-ray crystallography. These structures differ in the position of two 12 residue, β-hairpin regions (or “flaps”) which cap the active site. The movements of the flaps must be involved in the binding of substrates since, in either conformation, the flaps block the binding site. One of these structures is similar to structures of the ligand-bound enzyme; however, the importance of both structures to enzyme function is unclear. This transformation takes place on a time scale too long for conventional molecular dynamics simulations, so the process was studied by first identifying a reaction path between the two structures and then calculating the free energy along this path using umbrella sampling. For the ligand-free enzyme, it is found that the two structures are nearly equally stable, with the ligand-bound-type structure being less stable, consistent with X-ray crystallography data. The more stable open structure does not have a lower potential energy, but is stabilized by entropy. The transition occurs through a collapse and reformation of the β-sheet structure of the conformationally flexible, glycine-rich flap ends. Additionally, some problems in studying conformational changes in proteins through the use of a single reaction path are addressed. Proteins 32:7–16, 1998. © 1998 Wiley-Liss, Inc.     

Flexible Gates Generate Occluded Intermediates in the Transport Cycle of LacY

The major facilitator superfamily (MFS) transporter lactose permease (LacY) alternates between cytoplasmic and periplasmic open conformations to co-transport a sugar molecule together with a proton across the plasma membrane. Indirect experimental evidence suggested the existence of an occluded transition intermediate of LacY, which would prevent leaking of the proton gradient. As no experimental structure is known, the conformational transition is not fully understood in atomic detail. We simulated transition events from a cytoplasmic open conformation to a periplasmic open conformation with the dynamic importance sampling molecular dynamics method and observed occluded intermediates. Analysis of water permeation pathways and the electrostatic free-energy landscape of a solvated proton indicated that the occluded state contains a solvated central cavity inaccessible from either side of the membrane. We propose a pair of geometric order parameters that capture the state of the pathway through the MFS transporters as shown by a survey of available crystal structures and models. We present a model for the occluded state of apo-LacY, which is similar to the occluded crystal structures of the MFS transporters EmrD, PepT_So, NarU, PiPT and XylE. Our simulations are consistent with experimental double electron spin–spin distance measurements that have been interpreted to show occluded conformations. During the simulations, a salt bridge that has been postulated to be involved in driving the conformational transition formed. Our results argue against a simple rigid-body domain motion as implied by a strict “rocker-switch mechanism” and instead hint at an intricate coupling between two flexible gates.     

The c‐ring ion binding site of the ATP synthase from Bacillus pseudofirmus OF4 is adapted to alkaliphilic lifestyle

In the c‐ring rotor of ATP synthases ions are shuttled across the membrane during ATP synthesis by a unique rotary mechanism. We investigated characteristics of the c‐ring from the alkaliphile Bacillus pseudofirmus OF4 with respect to evolutionary adaptations to operate with protons at high environmental pH. The X‐ray structures of the wild‐type c₁₃ ring at pH 9.0 and a ‘neutralophile‐like’ mutant (P51A) at pH 4.4, at 2.4 and 2.8 Å resolution, respectively, reveal a dependency of the conformation and protonation state of the proton‐binding glutamate (E⁵⁴) on environmental hydrophobicity. Faster labelling kinetics with the inhibitor dicyclohexylcarbodiimide (DCCD) demonstrate a greater flexibility of E⁵⁴ in the mutant due to reduced water occupancy within the H⁺ binding site. A second ‘neutralophile‐like’ mutant (V21N) shows reduced growth at high pH, which is explained by restricted conformational freedom of the mutant's E⁵⁴ carboxylate. The study directly connects subtle structural adaptations of the c‐ring ion binding site to in vivo effects of alkaliphile cell physiology.     

Intrinsic nature of the three-dimensional structure of proteins as determined by distance geometry with good sampling properties

Summary A protocol for distance geometry calculation is shown to have excellent sampling properties in the determination of three-dimensional structures of proteins from nuclear magnetic resonance (NMR) data. This protocol uses a simulated annealing optimization employing mass-weighted molecular dynamics in four-dimensional space (Havel, T.F. (1991) Prog. Biophys. Mol. Biol., 56, 43–78). It attains an extremely large radius of convergence, allowing a random coil conformation to be used as the initial estimate for the succeeding optimization process. Computations are performed with four systems of simulated distance data as tests of the protocol, using an unconstrained l-alanine 30mer and three different types of proteins, bovine pancreatic trypsin inhibitor, the -amylase inhibitor Tendamistat, and the N-terminal domain of the 434-repressor. The test of the unconstrained polypeptide confirms that the sampled conformational space is that of the statistical random coil. In the larger and more complicated systems of the three proteins, the protocol gives complete convergence of the optimization without any trace of initial structure dependence. As a result of an exhaustive conformational sampling by the protocol, the intrinsic nature of the structures generated with distance restraints derived from NMR data has been revealed. When the sampled structures are compared with the corresponding X-ray structures, we find that the averages of the sampled structures always show a certain pattern of discrepancy from the X-ray structure. This discrepancy is due to the short distance nature of the distance restraints, and correlates with the characteristic shape of the protein molecule.Abbreviations r.m.s.d. root-mean-square deviation - MD molecular dynamics - NMR nuclear magnetic resonance - NOE nuclear Overhauser enhancement - BPTI bovine pancreatic trypsin inhibitor     

Identification and function of conformational dynamics in the multidomain GTPase dynamin

Vesicle release upon endocytosis requires membrane fission, catalyzed by the large GTPase dynamin. Dynamin contains five domains that together orchestrate its mechanochemical activity. Hydrogen–deuterium exchange coupled with mass spectrometry revealed global nucleotide‐ and membrane‐binding‐dependent conformational changes, as well as the existence of an allosteric relay element in the α2^S helix of the dynamin stalk domain. As predicted from structural studies, FRET analyses detect large movements of the pleckstrin homology domain (PHD) from a ‘closed’ conformation docked near the stalk to an ‘open’ conformation able to interact with membranes. We engineered dynamin constructs locked in either the closed or open state by chemical cross‐linking or deletion mutagenesis and showed that PHD movements function as a conformational switch to regulate dynamin self‐assembly, membrane binding, and fission. This PHD conformational switch is impaired by a centronuclear myopathy‐causing disease mutation, S619L, highlighting the physiological significance of its role in regulating dynamin function. Together, these data provide new insight into coordinated conformational changes that regulate dynamin function and couple membrane binding, oligomerization, and GTPase activity during dynamin‐catalyzed membrane fission.     

Stereochemistry of nucleic acids and polynucleotides. VI. Minimum energy conformations of dimethyl phosphate

As part of a study on the conformation of polynucleotides and nucleic acids the preferred conformations of the model conpound dimethyl phosphate are worked out using potential energy functions. In calculating the total potential energy associated with the conformation, nonbonded, torsional, and electrostatic terms have been considered. The variation of the total conformational energy is represented as a function of two torsion angles ? and ψ which are the rotations about the two phosphoester bonds. The most stable conformations are found to be the gauche–gauche conformations about these bonds. The conformations observed for phosphodiesters in the solid state and in the proposed structures of polynucleotides and nucleic acids cluster around the minimum. Also, regions of minimum energy correspond well with the typical allowed regions of a representative dinucleotide.     

Functional conformational changes of endo-1,4-xylanase II from Trichoderma reesei: A molecular dynamics study

Recent crystallographic studies have revealed a range of structural changes in the three-dimensional structure of endo-1,4-xylanase (XYNII) from Trichoderma reesei. The observed conformational changes can be described as snapshots of an open-close movement of the active site of XYNII. These structures were further analyzed in this study. In addition, a total of four 1 ns molecular dynamics (MD) simulations were performed representing different states of the enzyme. A comparison of the global and local changes found in the X-ray structures and the MD runs suggested that the simulations reproduced a similar kind of active site opening and closing as predicted by the crystal structures. The open-close movement was characterized by the use of distance difference matrixes and the Hingefind program (Wriggers and Schulten, Proteins 29:1–14, 1997) to be a ‘hinge-bending’ motion involving two large rigidly-moving regions and an extended hinge. This conformational feature is probably inherent to this molecular architecture and probably plays a role in the function of XYNII. Proteins 31:434–444, 1998. © 1998 Wiley-Liss, Inc.     

Mass-weighted molecular dynamics simulation and conformational analysis of polypeptide.

Atomic motions in protein molecules have been studied by molecular dynamics (MD) simulations; dynamics simulation methods have also been employed in conformational studies of polypeptide molecules. It was found that when atomic masses are weighted, the molecular dynamics method can significantly increase the sampling of dihedral conformation space in such studies, compared to a conventional MD simulation of the same total simulation time length. Herein the theoretical study of molecular conformation sampling by the molecular dynamics-based simulation method in which atomic masses are weighted is reported in detail; moreover, a numerical scheme for analyzing the extensive conformational sampling in the simulation of a tetrapeptide amide molecule is presented. From numerical analyses of the mass-weighted molecular dynamics trajectories of backbone dihedral angles, low-resolution structures covering the entire backbone dihedral conformation space of the molecule were determined, and the distribution of rotationally stable conformations in this space were analyzed quantitatively. The theoretical analyses based on the computer simulation and numerical analytical methods suggest that distinctive regimes in the conformational space of the peptide molecule can be identified.     

Prediction of the folding of short polypeptide segments by uniform conformational sampling

A procedure, CONGEN, for uniformly sampling the conformational spaceof short polypeptide segments in proteins has been implemented. Because thetime required for this sampling grows exponentially with the number of residues, parameters are introduced to limit the conformational space that has to be explored. This is done by the use of the empirical energy function ofCHARMM [B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan and M. Karplus (1983) J. Comput. Chem. 4 , 187-217] and truncating the search when conformations of grossly unfavorable energy are sampled. Tests are made to determine control parameters that optimize the search without excluding important configurations. When applied to known protein structures, the resulting procedure is generally capable of generating conformations where the lowest energy conformation matches the known structure within a rms deviation of 1 Å.     

From signal perception to signal transduction: ligand‐induced dimeric switch of DctB sensory domain in solution

Sinorhizobium meliloti DctB is a typical transmembrane sensory histidine kinase, which senses C₄‐dicarboxylic acids (DCA) and regulates the expression of DctA, the DCA transporter. We previously reported the crystal structures of its periplasmic sensory domain (DctBp) in apo and succinate‐bound states, and these structures showed dramatic conformational changes at dimeric level. Here we show a ligand‐induced dimeric switch in solution and a strong correlation between DctBp's dimerization states and the in vivo activities of DctB. Using site‐directed mutagenesis, we identify important determinants for signal perception and transduction. Specifically, we show that the ligand‐binding pocket is essential for DCA‐induced ‘on’ activity of DctB. Mutations at different sections of DctBp's dimerization interface can lock full‐length DctB at either ‘on’ or ‘off’ state, independent of ligand binding. Taken together, these results suggest that DctBp's signal perception and transduction occur through a ‘ligand‐induced dimeric switch’, in which the changes in the dimeric conformations upon ligand binding are responsible for the signal transduction in DctB.     

Effect of Sec62 on the conformation of the Sec61 channel in yeast

Most eukaryotic secretory and membrane proteins are funneled by the Sec61 complex into the secretory pathway. Furthermore, some substrate peptides rely on two essential accessory proteins, Sec62 and Sec63, being present to assist with their translocation via the Sec61 channel in post-translational translocation. Cryo-electron microscopy (cryo-EM) recently succeeded in determining atomistic structures of unbound and signal sequence-engaged Sec complexes from Saccharomyces cerevisiae, involving the Sec61 channel and the proteins Sec62, Sec63, Sec71 and Sec72. In this study, we investigated the conformational effects of Sec62 on Sec61. Indeed, we observed in molecular dynamics simulations that the conformational dynamics of lateral gate, plug and pore region of Sec61 are altered by the presence/absence of Sec62. In molecular dynamics simulations that were started from the cryo-EM structures of Sec61 coordinated to Sec62 or of apo Sec61, we observed that the luminal side of the lateral gate gradually adopts a closed conformation similar to the apo state during unbound state simulations. In contrast, it adopts a wider conformation in the bound state. Furthermore, we demonstrate that the conformation of the active (substrate-bound) state of the Sec61 channel shifts toward an alternative conformation in the absence of the substrate. We suggest that the signal peptide holds/stabilizes the active state conformation of Sec61 during post-translational translocation. Thus, our study explains the effect of Sec62 on the conformation of the Sec61 channel and describes the conformational transitions of Sec61 channel.