Key Protein Residues Research Articles

Weighted protein residue networks based on joint recurrences between residues.

BackgroundWeighted and un-weighted protein residue networks can predict key functional residues in proteins based on the closeness centrality C and betweenness centrality B values for each residue. A static snapshot of the protein structure, and a cutoff distance, are used to define edges between the network nodes. In this work we apply the weighted network approach to study the β-Lactamase Inhibitory Protein (BLIP). Joint recurrences extracted from molecular dynamics MD trajectory positions of the protein residue carbon alpha atoms are used to define edge weights between nodes, and no cutoff distance is used. The results for B and C from our approach are compared with those extracted from an un-weighted network, and a weighted network that uses interatomic contacts to define edge weights between nodes, respectively.ResultsThe joint recurrence weighted network approach performs well in pointing out key protein residues. Furthermore, it seems to emphasize residues with medium to high relative solvent accessibility that lie in loop regions between secondary structure elements of the protein.ConclusionsProtein residue networks that use joint recurrences extracted from molecular dynamics simulations of a solvated protein perform well in pointing to hotspot residues and hotspot clusters. This approach uses no distance cutoff threshold, and does not exclude any interactions between the residues, including water-mediated interactions.

Theory uncovers an unusual mechanism of DNA repair of a lesioned adenine by AlkB enzymes.

DNA-base lesions cause cancer and propagate into the genome. We use in-protein QM/MM calculations to study the repair of etheno-bridged adenine (εA) by the iron(IV)-oxo species of AlkB enzymes. Recent experimental investigations, using mass-spectrometry and in crystallo isolation, suggested that εA was repaired by formation of an epoxide (εA-ep) that further transforms to a glycol (εA-gl), ending finally in adenine and glyoxal. Theory reproduces the experimentally observed barrier for the rate-determining step and its pH dependence. However, as we show, the mass-spectrometrically identified species are side-byproducts unassociated with the repair mechanism. The repair is mediated by a zwitterionic species, of the same molecular mass as the epoxide, which transforms to an intermediate that matches the in crystallo trapped species in structure and mass, but is NOT the assumed εA-gl iron-glycol complex. Verifiable/falsifiable predictions, regarding the key protein residues, follow. The paper underscores the indispensable role of theory by providing atomistic descriptions of this vital mechanism, and guiding further experimental investigations.

Molecular simulations of β-lactoglobulin complexed with fatty acids reveal the structural basis of ligand affinity to internal and possible external binding sites.

The interaction of saturated fatty acids of different length (C8:0 to C18:0) with β-lactoglobulin (βLG) was investigated by molecular dynamics simulation and docking approaches. The results show that the presence of such ligands in the hydrophobic central cavity of βLG, known as the protein calyx, determines an enhancement of atomic fluctuations compared with the unliganded form, especially for loops at the entrance of the binding site. Concerted motions are evidenced for protein regions that could favor the binding of ligands. The mechanism of anchoring of fatty acids of different length is similar for the carboxylate head-group, through electrostatic interactions with the side chains of Lys60/Lys69. The key protein residues to secure the hydrocarbon chain are Phe105/Met107, which adapt their conformation upon ligand binding. In particular, Phe105 provides an additional hydrophobic clamp only for the tail of the two fatty acids with the longest chains, palmitic, and stearic acid, which are known to bind βLG with a high affinity. The search of additional external binding sites for fatty acids, distinct from the calyx, was also carried out for palmitic acid. Two external sites with a lower affinity were identified as secondary sites, one consisting in a hydrophobic cavity allowing two distinct binding modes for the fatty acid, and the other corresponding to a surface crevice close to the protein α-helix. The overall results provide a comprehensive picture of the dynamical behavior of βLG in complex with fatty acids, and elucidate the structural basis of the binding of these physiological ligands.

Interaction of the Inward Rectifier Potassium Channel Kir 2.2 with Phosphatidylserine

Inward-rectifier K+ (Kir) channels are ion channels that transport potassium into the cell. They are essential to maintain the resting membrane potential and to regulate the action potential duration in excitable cells [1]. As a consequence, Kir mutations result in several diseases, such as periodic paralysis or cardiac arrhythmia [1, 2].The activity of Kir channels is regulated by phosphatidylinositol-(4,5)-bisphosphate (PIP2), a negatively charged phospholipid that has been recognized as one of the major regulators of membrane excitability [2]. Moreover, it has been suggested that Kir activation is not only regulated by PIP2, but it has a secondary, non-specific requirement for other anionic phospholipids [3].Here we have investigated the interaction of phosphatidylserine (PS) with Kir 2.2 by means of all-atom molecular dynamics (MD). Simulations have been performed on both the apo and PIP2-bound states of the channel [4], embedded in a POPS membrane in the presence of KCl. These trajectories reveal the key protein residues interacting with POPS, complementing previous docking [3] and coarse-grain [5] studies. Furthermore, they suggest that the Kir 2.2-POPS interactions drive the cytoplasmic domain closer to the membrane and help to pre-assemble the PIP2 binding site. In other words, our simulations provide a molecular picture of the sensitized state proposed by Nichols and coworkers to explain the synergistic effect of anionic phospholipids [3].[1] Physiol. Rev. 90:291-366 (2010).[2] Pflugers Arch. 460:321-341 (2010).[3] (a) Biophys. J. 100:620-628 (2011). (b) J. Biol. Chem. 288:16726-16737 (2013). (c) Biophys. J., 104:433a (2013).[4] (a) Science, 326:1668-1674 (2009). (b) Nature (Letter), 477:495-499 (2011)[5] Biochemistry, 52:279-281 (2013).

A Conserved Mechanism of GABA Binding and Antagonism Is Revealed by Structure-Function Analysis of the Periplasmic Binding Protein Atu2422 in Agrobacterium tumefaciens

Bacterial periplasmic binding proteins (PBPs) and eukaryotic PBP-like domains (also called as Venus flytrap modules) of G-protein-coupled receptors are involved in extracellular GABA perception. We investigated the structural and functional basis of ligand specificity of the PBP Atu2422, which is implicated in virulence and transport of GABA in the plant pathogen Agrobacterium tumefaciens. Five high-resolution x-ray structures of Atu2422 liganded to GABA, Pro, Ala, and Val and of point mutant Atu2422-F77A liganded to Leu were determined. Structural analysis of the ligand-binding site revealed two essential residues, Phe(77) and Tyr(275), the implication of which in GABA signaling and virulence was confirmed using A. tumefaciens cells expressing corresponding Atu2422 mutants. Phe(77) restricts ligand specificity to α-amino acids with a short lateral chain, which act as antagonists of GABA signaling in A. tumefaciens. Tyr(275) specifically interacts with the GABA γ-amino group. Conservation of these two key residues in proteins phylogenetically related to Atu2422 brought to light a subfamily of PBPs in which all members could bind GABA and short α-amino acids. This work led to the identification of a fingerprint sequence and structural features for defining PBPs that bind GABA and its competitors and revealed their occurrence among host-interacting proteobacteria.

Advanced in silico Approaches in Antiviral Research

Computer-aided drug design has seen constantly increasing application over the past two decades in every area of drug discovery. It can offer significant advantages over conventional approaches, being far less expensive and faster than conventional methods, or offering the possibility to predict molecular behaviours that cannot be elucidated in any other way. Recent developments in software and hardware make it possible to simulate increasingly complex molecular environments, widening the applicability of in silico studies from the interactions of small molecules with key protein residues, to the simulation of the dynamic evolution of complex biological systems with atomic resolution. Antiviral research offers several open challenges, from a biological, biochemical and pharmaceutical point of view. Computational approaches are already providing some answers and will undoubtedly give more in the near future. Here, we present a brief overview of the cutting-edge computational methods that play a major role in present and future antiviral research.

Identifying stabilizing key residues in proteins using interresidue interaction energy matrix

We are proposing an interresidue interaction energy map (IEM)--a new tool for protein structure analysis and protein bioinformatics. This approach employs the sum of pair-wise interaction energies of a particular residue as a measure of its structural importance. We will show that the IEM can serve as a means for identifying key residues responsible for the stability of a protein. Our method can be compared with the interresidue contact map but has the advantage of weighting the contacts by the stabilization energy content which they bring to the protein structure. For the theoretical adjustment of the proposed method, we chose the Trp-cage mini protein as a model system to compare a spectrum of computational methods ranging from the ab initio MP2 level through the DFT method to empirical force-field methods. The IEM method correctly identifies Tryptophane 6 as the key residue in the Trp-cage. The other residues with the highest stabilizing contributions correspond to the structurally important positions in the protein. We have further tested our method on the Trp2Cage miniprotein--a P12W mutant of the Trp-cage and on two proteins from the rubredoxin family that differ in their thermostability. Our method correctly identified the thermodynamically more stable variants in both cases and therefore can also be used as a tool for the relative measurement of protein stability. Finally, we will point out the important role played by dispersion energy, which contributes significantly to the total stabilization energy and whose role in aromatic pairs is clearly dominant. Surprisingly, the dispersion energy plays an even more important role in the interaction of prolines with aromatic systems.

Steroid nuclear hormone receptors: The allosteric conversation

AbstractNuclear hormone receptors (NHRs) regulate gene expression by forming complexes with small molecule ligands, other proteins, and DNA. Specific ligands, which include steroids, thyroids, and retinoids, are responsible for the regulation of complex processes in cellular differentiation, homeostasis, reproduction, and development. NHRs responding to steroid ligands (e.g., androgens, estrogens, glucocorticoids, mineralocorticoids, and progesterones) are the focus of this review, in particular, highlighting the subtle relationships between ligand structure and receptor function. In the light of x‐ray crystallographic information obtained in recent years, the structural organization of the ligand‐ and DNA‐binding domains of steroid nuclear hormone receptors is well understood. Less clear is the means by which a steroid ligand affects the three‐dimensional structure of a receptor, and how subtle differences in ligand structure can lead to major differences in functional activity (agonism, partial agonism, antagonism). Potential mechanisms by which ligands elicit such effects, statistical and mutational techniques used to identify key protein residues likely involved in the transmission of structural information, and novel sites of ligand–protein interaction are discussed in this review. Drug Dev Res 68:95–106, 2007. © 2007 Wiley‐Liss, Inc.

Identification of key residues in proteins by using their physical characters

Key residues in proteins are important to their stability, folding, and functions. They usually are highly conserved and can be identified by sequence or structure alignments. However, these methods can only determine the locations of key residues in sequences and structures and give less information about their physical characters. In this paper, we try to identify key residues by analyzing their inter-residue interactions. The model we study is the protein domain from transducin. We show that the usual Gaussian network analysis and distance-based contact analysis have difficulty identifying the key residues in this protein, but the contact energies can do it well. We find that most key residues can be located by the lowest contact energies. This enables us to predict and analyze the key residues in other proteins. Our results suggest that contact energy analysis may provide an alternative approach to investigating the folding and stability of proteins.

First Crystal Structures of Human Carbonic Anhydrase II in Complex with Dual Aromatase−Steroid Sulfatase Inhibitors,

Carbonic anhydrase (CA) catalyzes the reversible hydration of carbon dioxide to hydrogen carbonate. The role of CA in maintaining pH balance has made it an attractive drug target for the treatment of cancer, and it has recently been implicated in the delivery of sulfamate-containing drugs. With the acceptance of steroid sulfatase as a target for hormone-dependent cancer, novel dual aromatase-steroid sulfatase inhibitors (DASIs) containing a sulfamate group are now being developed. In this study, we show that CA II is potently inhibited by several members of this class of inhibitor. The structures of CA II complexed with 4-[(4-O-sulfamoylbenzyl)(4-cyanophenyl)amino]-4H-[1,2,4]triazole (K(D) = 84 +/- 5 nM) and 4-[(3-bromo-4-O-sulfamoylbenzyl)(4-cyanophenyl)amino]-4H-[1,2,4]triazole (K(D) = 454 +/- 29 nM) are reported to 2.02 and 1.76 A, respectively. Both inhibitors ligate to the active site zinc(II) atom via their sulfamate nitrogen, while the rest of the molecule is contained within the hydrophobic binding pocket. Key protein residues include Val-121, Phe-131, Val-135, Val-143, Leu-141, Leu-198, Pro-202, and Leu-204. Despite being structurally similar, the two ligands experience different types of binding particularly in the sulfamate-containing aromatic ring and the opposite geometric arrangement of the triazole and cyanophenyl groups around the configurationally invertible central nitrogen atom. Small changes in inhibitor structure can cause large changes in binding to CA II, and this underlines the importance of structure-based drug design with this enzyme and other isoforms relevant to potential anticancer therapy. Moreover, these results underpin the idea that binding to erythrocyte CA II may be a general method of stabilizing and delivering sulfamate-based drugs in vivo.

Multiple ligand-binding modes in bacterial R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR), a bacterial plasmid-encoded enzyme associated with resistance to the drug trimethoprim, shows neither sequence nor structural homology with the chromosomal DHFR. It presents a highly symmetrical toroidal structure, where four identical monomers contribute to the unique central active-site pore. Two reactants (dihydrofolate, DHF), two cofactors (NADPH) or one of each (R67*DHF*NADPH) can be found simultaneously within the active site, the last one being the reactive ternary complex. As the positioning of the ligands has proven elusive to empirical determination, we addressed the problem from a theoretical perspective. Several potential structures of the ternary complex were generated using the docking programs AutoDock and FlexX. The variability among the final poses, many of which conformed to experimental data, prompted us to perform a comparative scoring analysis and molecular dynamics simulations to assess the stability of the complexes. Analysis of ligand-ligand and ligand-protein interactions along the 4 ns trajectories of eight different structures allowed us to identify important inter-ligand contacts and key protein residues. Our results, combined with published empirical data, clearly suggest that multipe binding modes of the ligands are possible within R67 DHFR. While the pterin ring of DHF and the nicotinamide ring of NADPH assume a stacked endo-conformation at the centre of the pore, probably assisted by V66, Q67 and I68, the tails of the molecules extend towards opposite ends of the cavity, adopting multiple configurations in a solvent rich-environment where hydrogen-bond interactions with K32 and Y69 may play important roles.

Small‐world network approach to identify key residues in protein–protein interaction

We show that protein complexes can be represented as small-world networks, exhibiting a relatively small number of highly central amino-acid residues occurring frequently at protein-protein interfaces. We further base our analysis on a set of different biological examples of protein-protein interactions with experimentally validated hot spots, and show that 83% of these predicted highly central residues, which are conserved in sequence alignments and nonexposed to the solvent in the protein complex, correspond to or are in direct contact with an experimentally annotated hot spot. The remaining 17% show a general tendency to be close to an annotated hot spot. On the other hand, although there is no available experimental information on their contribution to the binding free energy, detailed analysis of their properties shows that they are good candidates for being hot spots. Thus, highly central residues have a clear tendency to be located in regions that include hot spots. We also show that some of the central residues in the protein complex interfaces are central in the monomeric structures before dimerization and that possible information relating to hot spots of binding free energy could be obtained from the unbound structures.

Network Analysis of Protein Structures Identifies Functional Residues

Identifying active site residues strictly from protein three-dimensional structure is a difficult task, especially for proteins that have few or no homologues. We transformed protein structures into residue interaction graphs (RIGs), where amino acid residues are graph nodes and their interactions with each other are the graph edges. We found that active site, ligand-binding and evolutionary conserved residues, typically have high closeness values. Residues with high closeness values interact directly or by a few intermediates with all other residues of the protein. Combining closeness and surface accessibility identified active site residues in 70% of 178 representative structures. Detailed structural analysis of specific enzymes also located other types of functional residues. These include the substrate binding sites of acetylcholinesterases and subtilisin, and the regions whose structural changes activate MAP kinase and glycogen phosphorylase. Our approach uses single protein structures, and does not rely on sequence conservation, comparison to other similar structures or any prior knowledge. Residue closeness is distinct from various sequence and structure measures and can thus complement them in identifying key protein residues. Closeness integrates the effect of the entire protein on single residues. Such natural structural design may be evolutionary maintained to preserve interaction redundancy and contribute to optimal setting of functional sites.

PH dependent conformational and structural changes of xylanase from an alkalophilic thermophilic Bacillus sp (NCIM 59).

The pH induced conformational and structural changes of Xyl II have been investigated from the alkalophilic thermophilic Bacillus sp. using kinetic, circular dichroism and fluorescence spectroscopy studies. The systematic studies on the folding and stability of cellulase-free xylanases are important, since their biotechnological applications require them to function under extremes of pH and temperature. The Trp fluorescence and the kinetic constants were found dependent on the pH. Above pH 8, the enzyme exhibited unfolding transitions as revealed by a red shift in the emission maximum as well as decreases in the fluorescence intensity. Circular dichroism studies revealed a decrease in the CD ellipticity at 222 nm at pH 9 and 10. The reduced catalytic activity of Xyl II at alkaline pH is correlated to the pH induced unfolding and ionization or protonation of key protein residues. The pH profile of Xyl II showed apparent pK values of 5.5 and 7 for the free enzyme and 5.6 and 6.7 for the enzyme-substrate complex. The abnormally high pK of 6.7 indicated the participation of a carboxyl group present in a non-polar environment. The pH dependence of inactivation kinetics of Xyl II with Woodward's reagent K corroborates evidence for the presence of a catalytically important carboxyl residue. The sequence alignment studies of Xyl II, in combination with kinetic and chemical modification data provide strong evidence for the participation of Asp94 in the catalytic function. The Xyl II produced from an alkalophilic source, was stable at pH 10 with a t(1/2) of 24 h. However, the enzyme exhibited pH optimum at near neutral values, which can be explained by the ionization and microenvironment of the active site residues.