Flexibility of the Kir6.2 inward rectifier K+ channel pore

Interactions of sulfhydryl reagents with introduced cysteines in the pore-forming (Kir6.2) subunits of the KATP channel were examined. 2-Aminoethyl methanethiosulfonate (MTSEA+) failed to modify Cd2+-insensitive control-Kir6.2 channels, but rapidly and irreversibly modified Kir6.2[L164C] (L164C) channels. Although a single Cd2+ ion is coordinated by L164C, four MTSEA+ “hits” can occur, each sequentially reducing the single-channel current. A dimeric fusion of control-Kir6.2 and L164C subunits generates Cd2+-insensitive channels, confirming that at least three cysteines are required for coordination, but MTSEA+ modification of the dimer occurs in two hits. L164C channels were not modified by bromotrimethyl ammoniumbimane (qBBr+), even though qBBr+ caused voltage-dependent block (as opposed to modification) that was comparable to that of MTSEA+ or 3-(triethylammonium)propyl methanethiosulfonate (MTSPTrEA+), implying that qBBr+ can also enter the inner cavity but does not modify L164C residues. The Kir channel pore structure was modeled by homology with the KcsA crystal structure. A stable conformation optimally places the four L164C side chains for coordination of a single Cd2+ ion. Modification of these cysteines by up to four MTSEA+ (or three MTSPTrEA+, or two qBBr+) does not require widening of the cavity to accommodate the derivatives within it. However, like the KcsA crystal structure, the energy-minimized model shows a narrowing at the inner entrance, and in the Kir6.2 model this narrowing excludes all ions. To allow entry of ions as large as MTSPTrEA+ or qBBr+, the entrance must widen to >8 Å, but this widening is readily accomplished by minimal M2 helix motion and side-chain rearrangement.

Crystallographic structures of the ligand-binding domains of the androgen receptor and its T877A mutant complexed with the natural agonist dihydrotestosterone

The structures of the ligand-binding domains (LBD) of the wild-type androgen receptor (AR) and the T877A mutant corresponding to that in LNCaP cells, both bound to dihydrotestosterone, have been refined at 2.0 Å resolution. In contrast to the homodimer seen in the retinoid-X receptor and estrogen receptor LBD structures, the AR LBD is monomeric, possibly because of the extended C terminus of AR, which lies in a groove at the dimerization interface. Binding of the natural ligand dihydrotestosterone by the mutant LBD involves interactions with the same residues as in the wild-type receptor, with the exception of the side chain of threonine 877, which is an alanine residue in the mutant. This structural difference in the binding pocket can explain the ability of the mutant AR found in LNCaP cells (T877A) to accommodate progesterone and other ligands that the wild-type receptor cannot.

Large vortex-like structure of dipole field in computer models of liquid water and dipole-bridge between biomolecules

We propose a framework to describe the cooperative orientational motions of water molecules in liquid water and around solute molecules in water solutions. From molecular dynamics (MD) simulation a new quantity “site-dipole field” is defined as the averaged orientation of water molecules that pass through each spatial position. In the site-dipole field of bulk water we found large vortex-like structures of more than 10 Å in size. Such coherent patterns persist more than 300 ps although the orientational memory of individual molecules is quickly lost. A 1-ns MD simulation of systems consisting of two amino acids shows that the fluctuations of site-dipole field of solvent are pinned around the amino acids, resulting in a stable dipole-bridge between side-chains of amino acids. The dipole-bridge is significantly formed even for the side-chain separation of 14 Å, which corresponds to five layers of water. The way that dipole-bridge forms sensitively depends on the side-chain orientations and thereby explains the specificity in the solvent-mediated interactions between biomolecules.

Ras-catalyzed hydrolysis of GTP: a new perspective from model studies.

Despite the biological and medical importance of signal transduction via Ras proteins and despite considerable kinetic and structural studies of wild-type and mutant Ras proteins, the mechanism of Ras-catalyzed GTP hydrolysis remains controversial. We take a different approach to this problem: the uncatalyzed hydrolysis of GTP is analyzed, and the understanding derived is applied to the Ras-catalyzed reaction. Evaluation of previous mechanistic proposals from this chemical perspective suggests that proton abstraction from the attacking water by a general base and stabilization of charge development on the gamma-phosphoryl oxygen atoms would not be catalytic. Rather, this analysis focuses attention on the GDP leaving group, including the beta-gamma bridge oxygen of GTP, the atom that undergoes the largest change in charge in going from the ground state to the transition state. This leads to a new catalytic proposal in which a hydrogen bond from the backbone amide of Gly-13 to this bridge oxygen is strengthened in the transition state relative to the ground state, within an active site that provides a template complementary to the transition state. Strengthened transition state interactions of the active site lysine, Lys-16, with the beta-nonbridging phosphoryl oxygens and a network of interactions that positions the nucleophilic water molecule and gamma-phosphoryl group with respect to one another may also contribute to catalysis. It is speculated that a significant fraction of the GAP-activated GTPase activity of Ras arises from an additional interaction of the beta-gamma bridge oxygen with an Arg side chain that is provided in trans by GAP. The conclusions for Ras and related G proteins are expected to apply more widely to other enzymes that catalyze phosphoryl (-PO(3)2-) transfer, including kinases and phosphatases.

A new isoleucine substitution of Val-20 in transthyretin tetramers selectively impairs dimer-dimer contacts and causes systemic amyloidosis.

The most frequent form of inherited amyloidoses is associated with mutations in the transthyretin (TTR) gene coding for 127-amino acid residues of four identical, noncovalently linked subunits that form a pair of dimers in the plasma protein complex. Amyloid fibrils containing the variant and to a lesser extent the wild-type form of the TTR molecule are deposited in various organs, including peripheral nerves and the myocardium, with polyneuropathy and cardiomyopathy as major clinical manifestations. So far, more than 40 distinct amino acid substitutions distributed throughout the TTR sequence over 30 positions have been found to be correlated with an increased amyloidogenicity of TTR. Most of these amyloidogenic amino acid substitutions are suspected to alter the conformation and stability of the monomer. Here we identify and characterize by protein and DNA analysis a novel amyloidogenic Val-20 to Ile mutation in a German three-generation family. The index patient suffered from severe amyloid cardiomyopathy at the age of 60. Conformational stability and unfolding behavior of the Ile-20 monomer in urea gradients was found to be almost indistinguishable from that of wild-type TTR. In contrast, tetramer stability was significantly reduced in agreement with the expected change in the interactions between the two opposing dimers via the side chain of Ile-20. Our observations provide strong evidence for the view that amyloidogenic amino acid substitutions in TTR facilitate the conversion of tetrameric TTR complexes into those conformational intermediates of the TTR folding pathway that have an intrinsic amyloidogenic potential.

Local densities orthogonal to beta-sheet amide planes: patterns of packing in globular proteins.

We have investigated the efficiency of packing by calculating intramolecular packing density above and below peptide planes of internal beta-pleated sheet residues in five globular proteins. The orientation of interest was chosen to allow study of regions that are approximately perpendicular to the faces of beta-pleated sheets. In these locations, nonbonded van der Waals packing interactions predominate over hydrogen bonding and solvent interactions. We observed considerable variability in packing densities within these regions, confirming that the interior packing of a protein does not result in uniform occupation of the available space. Patterns of fluctuation in packing density suggest that the regular backbone-to-backbone network of hydrogen bonds is not likely to be interrupted to maximize van der Waals interactions. However, high-density packing tends to occur toward the ends of beta-structure strands where hydrogen bonds are more likely to involve nonpolar side-chain groups or solvent molecules. These features result in internal protein folding with a central low-density core surrounded by a higher-density subsurface shell, consistent with our previous calculations regarding overall protein packing density.

Crystallographic analysis of endogenous peptides associated with HLA-DR1 suggests a common, polyproline II-like conformation for bound peptides.

The structure of the human major histocompatibility complex (MHC) class II molecule HLA-DR1 derived from the human lymphoblastoid cell line LG-2 has been determined in a complex with the Staphylococcus aureus enterotoxin B superantigen. The HLA-DR1 molecule contains a mixture of endogenous peptides derived from cellular or serum proteins bound in the antigen-binding site, which copurify with the class II molecule. Continuous electron density for 13 amino acid residues is observed in the MHC peptide-binding site, suggesting that this is the core length of peptide that forms common interactions with the MHC molecule. Electron density is also observed for side chains of the endogenous peptides. The electron density corresponding to peptide side chains that interact with the DR1-binding site is more clearly defined than the electron density that extends out of the binding site. The regions of the endogenous peptides that interact with DRI are therefore either more restricted in conformation or sequence than the peptide side chains or amino acids that project out of the peptide-binding site. The hydrogen-bond interactions and conformation of a peptide model built into the electron density are similar to other HLA-DR-peptide structures. The bound peptides assume a regular conformation that is similar to a polyproline type II helix. The side-chain pockets and conserved asparagine residues of the DR1 molecule are well-positioned to interact with peptides in the polyproline type II conformation and may restrict the range of acceptable peptide conformations.

How are close residues of protein structures distributed in primary sequence?

Structurally neighboring residues are categorized according to their separation in the primary sequence as proximal (1-4 positions apart) and otherwise distal, which in turn is divided into near (5-20 positions), far (21-50 positions), very far ( > 50 positions), and interchain (from different chains of the same structure). These categories describe the linear distance histogram (LDH) for three-dimensional neighboring residue types. Among the main results are the following: (i) nearest-neighbor hydrophobic residues tend to be increasingly distally separated in the linear sequence, thus most often connecting distinct secondary structure units. (ii) The LDHs of oppositely charged nearest-neighbors emphasize proximal positions with a subsidiary maximum for very far positions. (iii) Cysteine-cysteine structural interactions rarely involve proximal positions. (iv) The greatest numbers of interchain specific nearest-neighbors in protein structures are composed of oppositely charged residues. (v) The largest fraction of side-chain neighboring residues from beta-strands involves near positions, emphasizing associations between consecutive strands. (vi) Exposed residue pairs are predominantly located in proximal linear positions, while buried residue pairs principally correspond to far or very far distal positions. The results are principally invariant to protein sizes, amino acid usages, linear distance normalizations, and over- and underrepresentations among nearest-neighbor types. Interpretations and hypotheses concerning the LDHs, particularly those of hydrophobic and charged pairings, are discussed with respect to protein stability and functionality. The pronounced occurrence of oppositely charged interchain contacts is consistent with many observations on protein complexes where multichain stabilization is facilitated by electrostatic interactions.

Crystal and molecular structure of paclitaxel (taxol).

Paclitaxel (formerly called taxol), an important anticancer drug, inhibits cell replication by binding to and stabilizing microtubule polymers. As drug-receptor interactions are governed by the three-dimensional stereochemistries of both participants, we have determined the crystal structure of paclitaxel to identify its conformational preferences that may be related to biological activity. The monoclinic crystals contain two independent paclitaxel molecules in the asymmetric unit plus several water and dioxane solvent molecules. Taxane ring conformation is very similar in both paclitaxel molecules and is similar to the taxane ring conformation found in the crystal structure of the paclitaxel analogue docetaxel (formerly called taxotere). The two paclitaxel molecules have carbon-13 side-chain conformations that differ from each other and from that of the corresponding side chain in the docetaxel crystal structure. The carbon-13 side-chain conformation of one paclitaxel molecule is similar to what was proposed from NMR studies done in polar solvents, while that of the other paclitaxel molecule is different and hitherto unobserved. The paclitaxel molecules interact with each other and with solvent atoms through an extensive network of hydrogen bonds. Analysis of the hydrogen-bonding network together with structure-activity studies may suggest which atoms of paclitaxel are important for binding to microtubule receptors.

Efeito do macrodipolo sobre a estabilidade térmica de derivados de 1,3,5-tricarboxamida-ciclo-hexano

1,3,5-tricarboxamida-ciclo-hexano são compostos capazes de se autoagregarem formando colunas supramoleculares as quais se mantêm unidas não só devido às interações das cadeias laterais mas também devido às ligações de hidrogênio de cada um dos três grupos amida por monômero. Cada monômero possui momento de dipolo elétrico associado aos grupos amida. Quando as amidas dos vários monômeros dentro da mesma coluna estão apontadas para a mesma direção, os momentos de dipolo individuais de todas as amidas se somam formando elevado dipolo ao longo do eixo da coluna, chamado de macrodipolo, o qual influencia as interações intercolunares. Neste trabalho foram investigadas quatro conformações as quais diferem entre si em relação à orientação dos grupos carbonila: a conformação Up-Up contém grupos carbonilas paralelos dentro das colunas e colunas paralela, a conformação Up-Down possui grupos carbonilas paralelos dentro das colunas e colunas antiparalelas, a conformação Intra-Up-Up contém grupos carbonilas antiparalelos dentro das colunas e colunas paralelas e a conformação Intra-Up-Down possui grupos carbonilas antiparalelos dentro das colunas e colunas antiparalelas. Foi usado Dinâmica Molecular Clássica para investigar o efeito das interações macrodipolo-macrodipolo das quatro diferentes conformações sobre a estabilidade térmica de três diferentes compostos derivados de 1,3,5-tricarboxamida-ciclo-hexano. Foi verificado que as conformações com colunas antiparalelas tendem a ser ligeiramente mais estáveis do que as conformações com orientação paralela. O efeito da orientação dos grupos carbonila dentro das colunas sob a estabilidade do material está relacionado a vários fatores, tais como cargas atômicas parciais, arranjo colunar ou natureza das cadeias laterais, e os resultados não são tão diretos como quando se compara as orientações entre colunas. Outro tópico investigado foi o comportamento do material durante a transição da fase colunar para a fase desordenada. As colunas podem se desmontar em três diferentes formas: elas podem completamente se desintegrar rapidamente, podem primeiro se desintegrar lentamente e então perder a ordem colunar ou primeiro perdem a ordem colunar e então se desmontam em um processo demorado. Tais comportamentos estão associados com as interações dentro e entre colunas.

Countering cooperative effects in protease inhibitors using constrained b-strand-mimicking templates in focused combinatorial libraries

A major problem in de novo design of enzyme inhibitors is the unpredictability of the induced fit, with the shape of both ligand and enzyme changing cooperatively and unpredictably in response to subtle structural changes within a ligand. We have investigated the possibility of dampening the induced fit by using a constrained template as a replacement for adjoining segments of a ligand. The template preorganizes the ligand structure, thereby organizing the local enzyme environment. To test this approach, we used templates consisting of constrained cyclic tripeptides, formed through side chain to main chain linkages, as structural mimics of the protease-bound extended beta-strand conformation of three adjoining amino acid residues at the N- or C-terminal sides of the scissile bond of substrates. The macrocyclic templates were derivatized to a range of 30 structurally diverse molecules via focused combinatorial variation of nonpeptidic appendages incorporating a hydroxyethylamine transition-state isostere. Most compounds in the library were potent inhibitors of the test protease (HIV-1 protease). Comparison of crystal structures for five protease-inhibitor complexes containing an N-terminal macrocycle and three protease-inhibitor complexes containing a C-terminal macrocycle establishes that the macrocycles fix their surrounding enzyme environment, thereby permitting independent variation of acyclic inhibitor components with only local disturbances to the protease. In this way, the location in the protease of various acyclic fragments on either side of the macrocyclic template can be accurately predicted. This type of templating strategy minimizes the problem of induced fit, reducing unpredictable cooperative effects in one inhibitor region caused by changes to adjacent enzyme-inhibitor interactions. This idea might be exploited in template-based approaches to inhibitors of other proteases, where a beta-strand mimetic is also required for recognition, and also other protein-binding ligands where different templates may be more appropriate.

Comparative agonist/antagonist responses in mutant human C5a receptors define the ligand binding site

The C terminus is responsible for all of the agonist activity of C5a at human C5a receptors (C5aRs). In this report we have mapped the ligand binding site on the C5aR using a series of agonist and antagonist peptide mimics of the C terminus of C5a as well as receptors mutated at putative interaction sites ( Ile(116), Arg(175), Arg(206), Glu(199), Asp(282), and Val(286)). Agonist peptide 1 (Phe-Lys-Pro-D-cyclohexylalanine-cyclohexylalanine-D-Arg) can be converted to an antagonist by substituting the bulkier Trp for cyclohexylalanine at position 5 ( peptide 2). Conversely, mutation of C5aR transmembrane residue Ile(116) to the smaller Ala (I116A) makes the receptor respond to peptide 2 as an agonist (Gerber, B. O., Meng, E. C., Dotsch, V., Baranski, T. J., and Bourne, H. R. (2001) J. Biol. Chem. 276, 3394 - 3400). However, a potent cyclic hexapeptide antagonist, Phe-cyclo-[Orn-Pro-D-cyclohexylalanine-Trp-Arg] ( peptide 3), derived from peptide 2 and which binds to the same receptor site, remains a full antagonist at I116AC5aR. This suggests that although the residue at position 5 might bind near to Ile(116), the latter is not essential for either activation or antagonism. Arg(206) and Arg(175) both appear to interact with the C-terminal carboxylate of C5a agonist peptides, suggesting a dynamic binding mechanism that may be a part of a receptor activation switch. Asp(282) has been previously shown to interact with the side chain of the C-terminal Arg residue, and Glu(199) may also interact with this side chain in both C5a and peptide mimics. Using these interactions to orient NMR-derived ligand structures in the binding site of C5aR, a new model of the interaction between peptide antagonists and the C5aR is presented.

Comparison of the binding of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines with their isosteric sulfonamides to the active site of phenylethanolamine N-methyltransferase

3-Fluoromethyl-7-(N-substituted aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines (14, 16, and 18-22) are highly potent and selective inhibitors of phenylethanolamine N-methyltransferase (PNMT). Molecular modeling studies with 3-fluoromethyl-7-(N-alkyl aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines, such as 16, suggested that the sulfonamide -NH-could form a hydrogen bond with the side chain of Lys57. However, SAR studies and analysis of the crystal structure of human PNMT (hPNMT) in complex with 7 indicated that the sulfonamide oxygens, and not the sulfonamide -NH-, formed favorable interactions with the enzyme. Thus, we hypothesized that replacement of the sulfonamide -NH-with a methylene group could result in compounds that would retain potency at PNMT and that would have increased lipophilicity, thus increasing the likelihood they will cross the blood brain barrier. A series of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines (23-30) were synthesized and evaluated for their PNMT inhibitory potency and affinity for the R2-adrenoceptor. A comparison of these compounds with their isosteric sulfonamides (14, 16, and 18-22) showed that the sulfones were more lipophilic but less potent than their corresponding sulfonamides. Sulfone 24 (hPNMT K-i = 1.3 mu M) is the most potent compound in this series and is quite selective for PNMT versus the R2-adrenoceptor, but 24 is less potent than the corresponding sulfonamide, 16 (hPNMT K-i = 0.13 mu M). We also report the crystal structure of hPNMT in complex with sulfonamide 15, from which a potential hydrogen bond acceptor within the hPNMT active site has been identified, the main chain carbonyl oxygen of Asn39. The interaction of this residue with the sulfonamide -NH-is likely responsible for much of the enhanced inhibitory potency of the sulfonamides versus the sulfones.

Modular alpha-helical mimetics with antiviral activity against respiratory syncitial virus

A 13-residue peptide sequence from a respiratory syncitial virus fusion protein was constrained in an alpha-helical conformation by fusing two back-to-back cyclic alpha-turn mimetics. The resulting peptide, Ac-(3 -> 7; 8 -> 12)-bicyclo-FP[KDEFD][KSIRD]V-NH2, was highly alpha-helical in water by CD and NMR spectroscopy, correctly positioning crucial binding residues (F488, I491, V493) on one face of the helix and side chain-side chain linkers on a noninteracting face of the helix. This compound displayed potent activity in both a recombinant fusion assay and an RSV antiviral assay (IC50 = 36 nM) and demonstrates for the first time that back-to-back modular alpha-helix mimetics can produce functional antagonists of important protein-protein interactions.

Structural and functional characterization of the conserved salt bridge in mammalian Paneth cell alpha-defensins - Solution structures of mouse cryptdin-4 AND (E15D)-cryptdin-4

Defensins are mediators of mammalian innate immunity, and knowledge of their structure-function relationships is essential for understanding their mechanisms of action. We report here the NMR solution structures of the mouse Paneth cell α-defensin cryptdin-4 (Crp4) and a mutant (E15D)-Crp4 peptide, in which a conserved Glu15 residue was replaced by Asp. Structural analysis of the two peptides confirms the involvement of this Glu in a conserved salt bridge that is removed in the mutant because of the shortened side chain. Despite disruption of this structural feature, the peptide variant retains a well defined native fold because of a rearrangement of side chains, which result in compensating favorable interactions. Furthermore, salt bridge-deficient Crp4 mutants were tested for bactericidal effects and resistance to proteolytic degradation, and all of the variants had similar bactericidal activities and stability to proteolysis. These findings support the conclusion that the function of the conserved salt bridge in Crp4 is not linked to bactericidal activity or proteolytic stability of the mature peptide.