Linking structural features of protein complexes and biological function

Protein–protein interaction (PPI) establishes the central basis for complex cellular networks in a biological cell. Association of proteins with other proteins occurs at varying affinities, yet with a high degree of specificity. PPIs lead to diverse functionality such as catalysis, regulation, signaling, immunity, and inhibition, playing a crucial role in functional genomics. The molecular principle of such interactions is often elusive in nature. Therefore, a comprehensive analysis of known protein complexes from the Protein Data Bank (PDB) is essential for the characterization of structural interface features to determine structure–function relationship. Thus, we analyzed a nonredundant dataset of 278 heterodimer protein complexes, categorized into major functional classes, for distinguishing features. Interestingly, our analysis has identified five key features (interface area, interface polar residue abundance, hydrogen bonds, solvation free energy gain from interface formation, and binding energy) that are discriminatory among the functional classes using Kruskal-Wallis rank sum test. Significant correlations between these PPI interface features amongst functional categories are also documented. Salt bridges correlate with interface area in regulator-inhibitors (r = 0.75). These representative features have implications for the prediction of potential function of novel protein complexes. The results provide molecular insights for better understanding of PPIs and their relation to biological functions.     

Effects of single nucleotide polymorphisms on human N-acetyltransferase 2 structure and dynamics by molecular dynamics simulation

Background

Arylamine N-acetyltransferase 2 (NAT2) is an important catalytic enzyme that metabolizes the carcinogenic arylamines, hydrazine drugs and chemicals. This enzyme is highly polymorphic in different human populations. Several polymorphisms of NAT2, including the single amino acid substitutions R64Q, I114T, D122N, L137F, Q145P, R197Q, and G286E, are classified as slow acetylators, whereas the wild-type NAT2 is classified as a fast acetylator. The slow acetylators are often associated with drug toxicity and efficacy as well as cancer susceptibility. The biological functions of these 7 mutations have previously been characterized, but the structural basis behind the reduced catalytic activity and reduced protein level is not clear.

Methodology/Principal Findings

We performed multiple molecular dynamics simulations of these mutants as well as NAT2 to investigate the structural and dynamical effects throughout the protein structure, specifically the catalytic triad, cofactor binding site, and the substrate binding pocket. None of these mutations induced unfolding; instead, their effects were confined to the inter-domain, domain 3 and 17-residue insert region, where the flexibility was significantly reduced relative to the wild-type. Structural effects of these mutations propagate through space and cause a change in catalytic triad conformation, cofactor binding site, substrate binding pocket size/shape and electrostatic potential.

Conclusions/Significance

Our results showed that the dynamical properties of all the mutant structures, especially in inter-domain, domain 3 and 17-residue insert region were affected in the same manner. Similarly, the electrostatic potential of all the mutants were altered and also the functionally important regions such as catalytic triad, cofactor binding site, and substrate binding pocket adopted different orientation and/or conformation relative to the wild-type that may affect the functions of the mutants. Overall, our study may provide the structural basis for reduced catalytic activity and protein level, as was experimentally observed for these polymorphisms.     

Electronic structure simulations of 2,6-dimethyl-2,5-heptadien-4-one by FTIR,FT-Raman,NMR, UV–vis,NBO and density functional theory

Density functional theory (DFT) (B3LYP and B3PW91) calculations have been carried out for 2,6-dimethyl-2,5-heptadien-4-one (DMHD4O) using 6–311++ G** basis set. Complete vibrational assignment and analysis of the fundamental modes of the compound were carried out from the FTIR and FT-Raman spectral data. The theoretical electronic absorption has been calculated by using time-dependent DFT (TD-DFT) methods and compared with the experimental spectra. The theoretically computed Frontier energy gaps and TD-DFT calculations are in good agreement with the experimental UV–vis spectral absorption. The chemical hardness measured from the Frontier molecular orbital energies of DMHD4O is 0.0693 eV. Electronic stability of the compound arising from hyperconjugative interactions and charge delocalisation were also investigated based on the natural bond orbital (NBO) analysis. Effective stabilisation energy E ⁽²⁾ associated with the interactions of the π and the lone pair of electrons was determined by the NBO analysis. ¹³C and ¹H NMR chemical shifts of the compound have been calculated by means of Gauge-Invariant Atomic Orbital using B3LYP/6–311++ G** method. The partial ionic character of the carbonyl group due to resonance render a partially positive charge to the carbonyl carbon, and thus C4 chemical shift lie in the very downfield 191.6 ppm. Comparison between the experimental and the theoretical results indicates that B3LYP method is able to provide satisfactory results for predicting vibrational, electronic and NMR properties.     

Expression of industrially relevant laccases: prokaryotic style

Laccases are a class of multi-copper oxidases (MCOs) that catalyze the one-electron oxidation of four equivalents of a reducing substrate, with the concomitant four-electron reduction of dioxygen to water. They can catalyze a multitude of reactions, including the degradation of polymers and oxidative coupling of phenolic compounds, positioning them as significant industrial enzymes. Although fungal laccases are well known and well characterized, only recently has in silico biology led to rapid advances in the discovery, characterization and engineered expression of prokaryotic laccases. We describe the recent burgeoning of prokaryotic laccases, their catalytic properties, structural features and molecular evolution, vis-à-vis fungal laccases where possible. Special focus is given to the application of laccases to the emerging cellulosic biofuel industry.     

Hot spots for allosteric regulation on protein surfaces

Recent work indicates a general architecture for proteins in which sparse networks of physically contiguous and coevolving amino acids underlie basic aspects of structure and function. These networks, termed sectors, are spatially organized such that active sites are linked to many surface sites distributed throughout the structure. Using the metabolic enzyme dihydrofolate reductase as a model system, we show that: (1) the sector is strongly correlated to a network of residues undergoing millisecond conformational fluctuations associated with enzyme catalysis, and (2) sector-connected surface sites are statistically preferred locations for the emergence of allosteric control in vivo. Thus, sectors represent an evolutionarily conserved "wiring" mechanism that can enable perturbations at specific surface positions to rapidly initiate conformational control over protein function. These findings suggest that sectors enable the evolution of intermolecular communication and regulation.