Ligand-insensitive State of Cardiac ATP-sensitive K+ Channels : Basis for Channel Opening

The mechanism by which ATP-sensitive K⁺ (K_ATP) channels open in the presence of inhibitory concentrations of ATP remains unknown. Herein, using a four-state kinetic model, we found that the nucleotide diphosphate UDP directed cardiac K_ATP channels to operate within intraburst transitions. These transitions are not targeted by ATP, nor the structurally unrelated sulfonylurea glyburide, which inhibit channel opening by acting on interburst transitions. Therefore, the channel remained insensitive to ATP and glyburide in the presence of UDP. “Rundown” of channel activity decreased the efficacy with which UDP could direct and maintain the channel to operate within intraburst transitions. Under this condition, the channel was sensitive to inhibition by ATP and glyburide despite the presence of UDP. This behavior of the K_ATP channel could be accounted for by an allosteric model of ligand-channel interaction. Thus, the response of cardiac K_ATP channels towards inhibitory ligands is determined by the relative lifetime the channel spends in a ligand-sensitive versus -insensitive state. Interconversion between these two conformational states represents a novel basis for K_ATP channel opening in the presence of inhibitory concentrations of ATP in a cardiac cell.     

Crystal structures of the psychrophilic alpha-amylase from Alteromonas haloplanctis in its native form and complexed with an inhibitor.

Alteromonas haloplanctis is a bacterium that flourishes in Antarctic sea-water and it is considered as an extreme psychrophile. We have determined the crystal structures of the alpha-amylase (AHA) secreted by this bacterium, in its native state to 2.0 angstroms resolution as well as in complex with Tris to 1.85 angstroms resolution. The structure of AHA, which is the first experimentally determined three-dimensional structure of a psychrophilic enzyme, resembles those of other known alpha-amylases of various origins with a surprisingly greatest similarity to mammalian alpha-amylases. AHA contains a chloride ion which activates the hydrolytic cleavage of substrate alpha-1,4-glycosidic bonds. The chloride binding site is situated approximately 5 angstroms from the active site which is characterized by a triad of acid residues (Asp 174, Glu 200, Asp 264). These are all involved in firm binding of the Tris moiety. A reaction mechanism for substrate hydrolysis is proposed on the basis of the Tris inhibitor binding and the chloride activation. A trio of residues (Ser 303, His 337, Glu 19) having a striking spatial resemblance with serine-protease like catalytic triads was found approximately 22 angstroms from the active site. We found that this triad is equally present in other chloride dependent alpha-amylases, and suggest that it could be responsible for autoproteolytic events observed in solution for this cold adapted alpha-amylase.     

Active and inactive state structures of unliganded Lactobacillus casei allosteric L‐lactate dehydrogenase

Lactobacillus casei L ‐lactate dehydrogenase (LCLDH) is activated through the homotropic and heterotropic activation effects of pyruvate and fructose 1,6‐bisphosphate (FBP), respectively, and exhibits unusually high pH‐dependence in the allosteric effects of these ligands. The active (R) and inactive (T) state structures of unliganded LCLDH were determined at 2.5 and 2.6 Å resolution, respectively. In the catalytic site, the structural rearrangements are concerned mostly in switching of the orientation of Arg171 through the flexible intersubunit contact at the Q‐axis subunit interface. The distorted orientation of Arg171 in the T state is stabilized by a unique intra‐helix salt bridge between Arg171 and Glu178, which is in striking contrast to the multiple intersubunit salt bridges in Lactobacillus pentosus nonallosteric L ‐lactate dehydrogenase. In the backbone structure, major structural rearrangements of LCLDH are focused in two mobile regions of the catalytic domain. The two regions form an intersubunit linkage through contact at the P‐axis subunit interface involving Arg185, replacement of which with Gln severely decreases the homotropic and hetertropic activation effects on the enzyme. These two regions form another intersubunit linkage in the Q‐axis related dimer through the rigid NAD‐binding domain, and thus constitute a pivotal frame of the intersubunit linkage for the allosteric motion, which is coupled with the concerted structural change of the four subunits in a tetramer, and of the binding sites for pyruvate and FBP. The unique intersubunit salt bridges, which are observed only in the R state structure, are likely involved in the pH‐dependent allosteric equilibrium. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

The stereodynamics of 1,2‐dipropyldiaziridines

N‐alkylated trans‐diaziridines are an intriguing class of compounds with two stereogenic nitrogen atoms which easily interconvert. In the course of our investigations of the nature of the interconversion process via nitrogen inversion or electrocyclic ring opening ring closure, we synthesized and characterized the three constitutionally isomeric diaziridines 1,2‐di‐n‐propyldiaziridine 1 , 1‐isopropyl‐2‐n‐propyldiaziridine 2 , and 1,2‐diisopropyldiaziridine 3 to study the influence of the substituents on the interconversion barriers. Enantiomer separation was achieved by enantioselective gas chromatography on the chiral stationary phase Chirasil‐β‐Dex with high separation factors α (1‐isopropyl‐2‐n‐propyldiaziridine: 1.18; 1, 2‐diisopropyldiaziridine: 1.24; 100°C 50 kPa He) for the isopropyl substituted diaziridines. These compounds showed pronounced plateau formation between 100 and 150°C, and peak coalescence at elevated temperatures. The enantiomerization barriers ΔG? and activation parameters ΔH? and ΔS? were determined by enantioselective dynamic gas chromatography (DGC) and direct evaluation of the elution profiles using the unified equation implemented in the software DCXplorer. Interestingly, 1‐isopropyl‐2‐n‐propyldiaziridine and 1,2‐diisopropyldiaziridine exhibit similar high interconversion barriers ΔG? (100°C) of 128.3 ± 0.4 kJ mol^?1 and 129.8 ± 0.4 kJ mol^?1, respectively, which indicates that two sterically demanding substituents do not substantially increase the barrier as expected for a distinct nitrogen inversion process. Chirality, 2010. © 2009 Wiley‐Liss, Inc.     

Protein folding pathways and state transitions described by classical equations of motion of an elastic network model

Protein topology defined by the matrix of residue contacts has proved to be a fruitful basis for the study of protein dynamics. The widely implemented coarse-grained elastic network model of backbone fluctuations has been used to describe crystallographic temperature factors, allosteric couplings, and some aspects of the folding pathway. In the present study, we develop a model of protein dynamics based on the classical equations of motion of a damped network model (DNM) that describes the folding path from a completely unfolded state to the native conformation through a single-well potential derived purely from the native conformation. The kinetic energy gained through the collapse of the protein chain is dissipated through a friction term in the equations of motion that models the water bath. This approach is completely general and sufficiently fast that it can be applied to large proteins. Folding pathways for various proteins of different classes are described and shown to correlate with experimental observations and molecular dynamics and Monte Carlo simulations. Allosteric transitions between alternative protein structures are also modeled within the DNM through an asymmetric double-well potential.     

A recombinant allosteric lectin antagonist of HIV-1 envelope gp120 interactions

The first, critical stage of HIV-1 infection is fusion of viral and host cellular membranes initiated by a viral envelope glycoprotein gp120. We evaluated the potential to form a chimeric protein entry inhibitor that combines the action of two gp120-targeting molecules, an allosteric peptide inhibitor 12p1 and a higher affinity carbohydrate-binding protein cyanovirin (CVN). In initial mixing experiments, we demonstrated that the inhibitors do not interfere with each other and instead show functional synergy in inhibiting viral cell infection. Based on this, we created a chimera, termed L5, with 12p1 fused to the C-terminal domain of CVN through a linker of five penta-peptide repeats. L5 revealed the same broad specificity as CVN for gp120 from a variety of clades and tropisms. By comparison to CVN, the L5 chimera exhibited substantially increased inhibition of gp120 binding to receptor CD4, coreceptor surrogate mAb 17b and gp120 antibody F105. These binding inhibition effects by the chimera reflected both the high affinity of the CVN domain and the allosteric action of the 12p1 domain. The results open up the possibility to form high potency chimeras, as well as noncovalent mixtures, as leads for HIV-1 envelope antagonism that can overcome potency limits and potential virus mutational resistance for either 12p1 or CVN alone.     

Local motions in a benchmark of allosteric proteins

Allosteric proteins have been studied extensively in the last 40 years, but so far, no systematic analysis of conformational changes between allosteric structures has been carried out. Here, we compile a set of 51 pairs of known inactive and active allosteric protein structures from the Protein Data Bank. We calculate local conformational differences between the two structures of each protein using simple metrics, such as backbone and side-chain Cartesian displacement, and torsion angle change and rearrangement in residue-residue contacts. Thresholds for each metric arise from distributions of motions in two control sets of pairs of protein structures in the same biochemical state. Statistical analysis of motions in allosteric proteins quantifies the magnitude of allosteric effects and reveals simple structural principles about allostery. For example, allosteric proteins exhibit substantial conformational changes comprising about 20% of the residues. In addition, motions in allosteric proteins show strong bias toward weakly constrained regions such as loops and the protein surface. Correlation functions show that motions communicate through protein structures over distances averaging 10-20 residues in sequence space and 10-20 A in Cartesian space. Comparison of motions in the allosteric set and a set of 21 nonallosteric ligand-binding proteins shows that nonallosteric proteins also exhibit bias of motion toward weakly constrained regions and local correlation of motion. However, allosteric proteins exhibit twice as much percent motion on average as nonallosteric proteins with ligand-induced motion. These observations may guide efforts to design flexibility and allostery into proteins.     

Binding hotspots on K‐ras: Consensus ligand binding sites and other reactive regions from probe‐based molecular dynamics analysis

We have used probe‐based molecular dynamics (pMD) simulations to search for interaction hotspots on the surface of the therapeutically highly relevant oncogenic K‐Ras G12D. Combining the probe‐based query with an ensemble‐based pocket identification scheme and an analysis of existing Ras‐ligand complexes, we show that (i) pMD is a robust and cost‐effective strategy for binding site identification, (ii) all four of the previously reported ligand binding sites are suitable for structure‐based ligand design, and (iii) in some cases probe binding and expanded sampling of configurational space enable pocket expansion and increase the likelihood of site identification. Furthermore, by comparing the distribution of hotspots in nonpocket‐like regions with known protein‐ and membrane‐interacting interfaces, we propose that pMD has the potential to predict surface patches responsible for protein‐biomolecule interactions. These observations have important implications for future drug design efforts and will facilitate the search for potential interfaces responsible for the proposed transient oligomerization or interaction of Ras with other biomolecules in the cellular milieu. Proteins 2015; 83:898–909. © 2015 Wiley Periodicals, Inc.     

Nucleotide‐dependent structural fluctuations and regulation of microtubule‐binding affinity of KIF1A

Molecular motors such as kinesin regulate affinity to a rail protein during the ATP hydrolysis cycle. The regulation mechanism, however, is yet to be determined. To understand this mechanism, we investigated the structural fluctuations of the motor head of the single‐headed kinesin called KIF1A in different nucleotide states using molecular dynamics simulations of a Gō‐like model. We found that the helix at the microtubule (MT) binding site intermittently exhibits a large structural fluctuation when MT is absent. Frequency of this fluctuation changes systematically according to the nucleotide states and correlates strongly with the experimentally observed binding affinity to MT. We also showed that thermal fluctuation enhances the correlation and the interaction with the nucleotide suppresses the fluctuation of the helix . These results suggest that KIF1A regulates affinity to MT by changing the flexibility of the helix during the ATP hydrolysis process: the binding site becomes more flexible in the strong binding state than in the weak binding state. Proteins 2015; 83:809–819. © 2015 Wiley Periodicals, Inc.     

Allosteric sites can be identified based on the residue–residue interaction energy difference

Allosteric drugs act at a distance to regulate protein functions. They have several advantages over conventional orthosteric drugs, including diverse regulation types and fewer side effects. However, the rational design of allosteric ligands remains a challenge, especially when it comes to the identification allosteric binding sites. As the binding of allosteric ligands may induce changes in the pattern of residue–residue interactions, we calculated the residue–residue interaction energies within the allosteric site based on the molecular mechanics generalized Born surface area energy decomposition scheme. Using a dataset of 17 allosteric proteins with structural data for both the apo and the ligand‐bound state available, we used conformational ensembles generated by molecular dynamics simulations to compute the differences in the residue–residue interaction energies in known allosteric sites from both states. For all the known sites, distinct interaction energy differences (>25%) were observed. We then used CAVITY, a binding site detection program to identify novel putative allosteric sites in the same proteins. This yielded a total of 31 “druggable binding sites,” of which 21 exhibited >25% difference in residue interaction energies, and were hence predicted as novel allosteric sites. Three of the predicted allosteric sites were supported by recent experimental studies. All the predicted sites may serve as novel allosteric sites for allosteric ligand design. Our study provides a computational method for identifying novel allosteric sites for allosteric drug design. Proteins 2015; 83:1375–1384. © 2014 Wiley Periodicals, Inc.