Three-dimensional view of the surface motif associated with the P-loop structure: cis and trans cases of convergent evolution

Here we identify the determinants of the nucleotide-binding ability associated with the P-loop-containing proteins, inferring their functional importance from their structural convergence to a unique three- dimensional (3D) motif. (1) A new surface 3D pattern is identified for the P-loop nucleotide-binding region, which is more selective than the corresponding sequence pattern; (2) the signature displays one residue that we propose is the determinant for the guanine-binding ability (the residues aligned to ras D119; this residue is known to be important only in the G-proteins, we extend the prediction to all the other P-loop- containing proteins); and (3) two cases of convergent evolution at the molecular level are highlighted in the analysis of the active site: the positive charge aligned to ras K117 and the arginine residues aligned to the GAP arginine finger.The analysis of the residues conserved on protein surfaces allows one to identify new functional or evolutionary relationships among protein structures that would not be detectable by conventional sequence or structure comparison methods.     

The exposed N-terminal tail of the D1 subunit is required for rapid D1 degradation during photosystem II repair in Synechocystis sp PCC 6803

The selective replacement of photodamaged D1 protein within the multisubunit photosystem II (PSII) complex is an important photoprotective mechanism in chloroplasts and cyanobacteria. FtsH proteases are involved at an early stage of D1 degradation, but it remains unclear how the damaged D1 subunit is recognized, degraded, and replaced. To test the role of the N-terminal region of D1 in PSII biogenesis and repair, we have constructed mutants of the cyanobacterium Synechocystis sp PCC 6803 that are truncated at the exposed N terminus. Removal of 5 or 10 residues blocked D1 synthesis, as assessed in radiolabeling experiments, whereas removal of 20 residues restored the ability to assemble oxygen-evolving dimeric PSII complexes but inhibited PSII repair at the level of D1 degradation. Overall, our results identify an important physiological role for the exposed N-terminal tail of D1 at an early step in selective D1 degradation. This finding has important implications for the recognition of damaged D1 and its synchronized replacement by a newly synthesized subunit.     

Exploring the binding mechanism of HDAC8 selective inhibitors: Lessons from the modification of Cap group

The abnormal expression of histone deacetylase 8 (HDAC8) has been reported to associate with various cancer entities (colon, breast cancer, pancreas, etc.) as well as parasitic diseases, making HDAC8 gradually develop into an attractive and potential therapeutic target. Among the various design strategies of selective HDAC8 inhibitors (modification of Cap, Linker, or zinc binding group regions), the optimization of Cap region has aroused great interest among the researchers. However, the detailed information underlying how the modification of Cap region influences the inhibitory activities is still unclear, and in this study, compounds 2c, 3g, and 3n were selected to explore the differences in binding mechanisms brought by Cap modifications via various computational approaches at the atomic level. Five residues (Y293, H167, D254, D165, and M261) have a large difference in energy contributions to the constructed systems, and the subpocket formed by Y293 and M261 could interact with Cap groups, triggering the differences in the energy contributions of the residues (H167, D254, and D165) located in metal-catalytic center. In summary, the compounds 2c, 3g, and 3n were selected as molecular probes to explore the binding mechanism, and the residues (Y293 and M261) forming the subpocket should be paid special attention in the design and synthesis of novel selective HDAC8 inhibitors.     

Structural Basis of Drug Recognition by the Multidrug Transporter ABCG2

ABCG2 is an ATP-binding cassette (ABC) transporter whose function affects the pharmacokinetics of drugs and contributes to multidrug resistance of cancer cells. While its interaction with the endogenous substrate estrone-3-sulfate (E₁S) has been elucidated at a structural level, the recognition and recruitment of exogenous compounds is not understood at sufficiently high resolution. Here we present three cryo-EM structures of nanodisc-reconstituted, human ABCG2 bound to anticancer drugs tariquidar, topotecan and mitoxantrone. To enable structural insight at high resolution, we used Fab fragments of the ABCG2-specific monoclonal antibody 5D3, which binds to the external side of the transporter but does not interfere with drug-induced stimulation of ATPase activity. We observed that the binding pocket of ABCG2 can accommodate a single tariquidar molecule in a C-shaped conformation, similar to one of the two tariquidar molecules bound to ABCB1, where tariquidar acts as an inhibitor. We also found single copies of topotecan and mitoxantrone bound between key phenylalanine residues. Mutagenesis experiments confirmed the functional importance of two residues in the binding pocket, F439 and N436. Using 3D variability analyses, we found a correlation between substrate binding and reduced dynamics of the nucleotide binding domains (NBDs), suggesting a structural explanation for drug-induced ATPase stimulation. Our findings provide additional insight into how ABCG2 differentiates between inhibitors and substrates and may guide a rational design of new modulators and substrates.     

Functional impacts of non-synonymous single nucleotide polymorphisms: selective constraint and structural environments

In this work, we studied the correlations between selective constraint, structural environments and functional impacts of non-synonymous single nucleotide polymorphisms (nsSNPs). We found that the relation between solvent accessibility and functional impacts of nsSNPs is not as simple as generally thought. Finer structural classifications need to be taken into account to reveal the complex relations between the characteristics of a structure environment and its influence on the functional impacts of nsSNPs. We introduced two parameters for each structural environment, consensus residue percentage and residue distribution distance, to characterize the selective constraint imposed by the environment. Both parameters significantly correlate with the functional bias of nsSNPs across the structural environments. This result shows that selective constraint underlies the bias of a structural environment towards a certain type of nsSNPs (disease-associated or benign).     

Ligand recognition by the vitamin D receptor.

Three-dimensional structure of the ligand binding domain (LBD) of the vitamin D receptor (VDR) docked with the natural ligand 1 alpha,25-dihydroxyvitamin D(3) [1,25-(OH)(2)D(3)] has been mostly solved by the X-ray crystallographic analysis of the deletion mutant (VDR-LBD Delta 165-215). The important focus, from now on, is how the VDR recognizes and interacts with potent synthetic ligands. We now report the docking models of the VDR with three functionally and structurally interesting ligands, 22-oxa-1,25-(OH)(2)D(3) (OCT), 20-epi-1,25-(OH)(2)D(3) and 20-epi-22-oxa-24,26,27-trihomo-1,25-(OH)(2)D(3). In parallel with the computational docking studies, we prepared twelve one-point mutants of amino acid residues lining the ligand binding pocket of the VDR and examined their transactivation potency induced by 1,25-(OH)(2)D(3) and these synthetic ligands. The results indicate that L233, R274, W286, H397 and Y401 are essential for holding the all ligands tested, S278 and Q400 are not important at all, and the importance of S237, V234, S275, C288 and H305 is variable depending on the side-chain structure of the ligands. Based on these studies, we suggested key structural factors to bestow the selective action on OCT and the augmented activities on 20-epi-ligands. Furthermore, the docking models coincided well with our proposed active space-region theory of vitamin D based on the conformational analyses of ligands.     

Molecular modelling of the structures of endothelin antagonists. Identification of a possible structural determinant for ET-A receptor binding.

Computer-aided molecular modelling of the endothelin (ET-A) receptor antagonists, BQ-123 and BE-18257B, shows that they have very similar 3D structures. Parts of their 3D structures are also shown to match closely with that reported for residues 6-8 in endothelin-1. On the basis of these similarities (and with supporting evidence from literature data on endothelin structure-activity relationships) a structural determinant is proposed for ET-A receptor binding, and novel designs of peptide are suggested for providing more potent and selective ET-A receptor antagonists.     

Mapping the mutual information network of enzymatic families in the protein structure to unveil functional features

Amino acids committed to a particular function correlate tightly along evolution and tend to form clusters in the 3D structure of the protein. Consequently, a protein can be seen as a network of co-evolving clusters of residues. The goal of this work is two-fold: first, we have combined mutual information and structural data to describe the amino acid networks within a protein and their interactions. Second, we have investigated how this information can be used to improve methods of prediction of functional residues by reducing the search space. As a main result, we found that clusters of co-evolving residues related to the catalytic site of an enzyme have distinguishable topological properties in the network. We also observed that these clusters usually evolve independently, which could be related to a fail-safe mechanism. Finally, we discovered a significant enrichment of functional residues (e.g. metal binding, susceptibility to detrimental mutations) in the clusters, which could be the foundation of new prediction tools.     

Structural prediction of the beta-domain of metallothionein-3 by molecular dynamics simulation

The beta-domain of metallothionein-3 (MT3) has been reported to be crucial to the neuron growth inhibitory bioactivity. Little detailed three-dimensional structural information is available to present a reliable basis for elucidation on structure-property-function relationships of this unique protein by experimental techniques. So, molecular dynamics simulation is adopted to study the structure of beta-domain of MT3. In this article, a 3D structural model of beta-domain of MT3 was generated. The molecular simulations provide detailed protein structural information of MT3. As compared with MT2, we found a characteristic conformation formed in the fragment (residue 1-13) at the N-terminus of MT3 owing to the constraint induced by 5TCPCP9, in which Pro7 and Pro9 residues are on the same side of the protein, both facing outward and the two 5-member rings of prolines are arranged almost in parallel, while Thr5 is on the opposite side. Thr5 in MT3 is also found to make the first four residues relatively far from the fragment (residue 23-26) as compared with MT2. The simulated structure of beta-domain of MT3 is looser than that of MT2. The higher energy of MT3 than that of MT2 calculated supports these conclusions. Simulation on the four isomer arising from the cis- or trans-configuration of 6CPCP9 show that the trans-/trans-isomer is energetic favorable. The partially unfolding structure of beta-domain of MT3 is also simulated and the results show the influence of 6CPCP9 sequence on the correct folding of this domain. The correlations between the bioactivity of MT3 and the simulated structure as well as the folding of beta-domain of MT3 are discussed based on our simulation and previous results.     

Structural mapping of protein interactions reveals differences in evolutionary pressures correlated to mRNA level and protein abundance

Genome-wide studies in Saccharomyces cerevisiae concluded that the dominant determinant of protein evolutionary rates is expression level: highly expressed proteins generally evolve most slowly. To determine how this constraint affects the evolution of protein interactions, we directly measure evolutionary rates of protein interface, surface, and core residues by structurally mapping domain interactions to yeast genomes. We find that mRNA level and protein abundance, though correlated, report on pressures affecting regions of proteins differently. Pressures proportional to mRNA level slow evolutionary rates of all structural regions and reduce the variability in rate differences between interfaces and other surfaces. In contrast, the evolutionary rate variation within a domain is much less correlated to protein abundance. Distinct pressures may be associated primarily with the cost (mRNA level) and functional (protein abundance) benefit of protein production. Interfaces of proteins with low mRNA levels may have higher evolutionary flexibility and could constitute the raw material for new functions.     

Phosphorylation of the LIR Domain of SCOC Modulates ATG8 Binding Affinity and Specificity

Autophagy is a highly conserved degradative pathway, essential for cellular homeostasis and implicated in diseases including cancer and neurodegeneration. Autophagy-related 8 (ATG8) proteins play a central role in autophagosome formation and selective delivery of cytoplasmic cargo to lysosomes by recruiting autophagy adaptors and receptors. The LC3-interacting region (LIR) docking site (LDS) of ATG8 proteins binds to LIR motifs present in autophagy adaptors and receptors. LIR-ATG8 interactions can be highly selective for specific mammalian ATG8 family members (LC3A-C, GABARAP, and GABARAPL1-2) and how this specificity is generated and regulated is incompletely understood.We have identified a LIR motif in the Golgi protein SCOC (short coiled-coil protein) exhibiting strong binding to GABARAP, GABARAPL1, LC3A and LC3C. The residues within and surrounding the core LIR motif of the SCOC LIR domain were phosphorylated by autophagy-related kinases (ULK1-3, TBK1) increasing specifically LC3 family binding. More distant flanking residues also contributed to ATG8 binding. Loss of these residues was compensated by phosphorylation of serine residues immediately adjacent to the core LIR motif, indicating that the interactions of the flanking LIR regions with the LDS are important and highly dynamic.Our comprehensive structural, biophysical and biochemical analyses support and provide novel mechanistic insights into how phosphorylation of LIR domain residues regulates the affinity and binding specificity of ATG8 proteins towards autophagy adaptors and receptors.     

Sequence-dependent structural variation in DNA undergoing intrahelical inspection by the DNA glycosylase MutM

MutM, a bacterial DNA-glycosylase, plays a critical role in maintaining genome integrity by catalyzing glycosidic bond cleavage of 8-oxoguanine (oxoG) lesions to initiate base excision DNA repair. The task faced by MutM of locating rare oxoG residues embedded in an overwhelming excess of undamaged bases is especially challenging given the close structural similarity between oxoG and its normal progenitor, guanine (G). MutM actively interrogates the DNA to detect the presence of an intrahelical, fully base-paired oxoG, whereupon the enzyme promotes extrusion of the target nucleobase from the DNA duplex and insertion into the extrahelical active site. Recent structural studies have begun to provide the first glimpse into the protein-DNA interactions that enable MutM to distinguish an intrahelical oxoG from G; however, these initial studies left open the important question of how MutM can recognize oxoG residues embedded in 16 different neighboring sequence contexts (considering only the 5'- and 3'-neighboring base pairs). In this study we set out to understand the manner and extent to which intrahelical lesion recognition varies as a function of the 5'-neighbor. Here we report a comprehensive, systematic structural analysis of the effect of the 5'-neighboring base pair on recognition of an intrahelical oxoG lesion. These structures reveal that MutM imposes the same extrusion-prone ("extrudogenic") backbone conformation on the oxoG lesion irrespective of its 5'-neighbor while leaving the rest of the DNA relatively free to adjust to the particular demands of individual sequences.     

Effect of mutations of N- and C-terminal charged residues on the activity of LCAT

On the basis of structural homology calculations, we previously showed that lecithin:cholesterol acyltransferase (LCAT), like lipases, belongs to the alpha/beta hydrolase fold family. As there is higher sequence conservation in the N-terminal region of LCAT, we investigated the contribution of the N- and C-terminal conserved basic residues to the catalytic activity of this enzyme. Most basic, and some acidic residues, conserved among LCAT proteins from different species, were mutated in the N-terminal (residues 1;-210) and C-terminal (residues 211;-416) regions of LCAT. Measurements of LCAT-specific activity on a monomeric substrate, on low density lipoprotein (LDL), and on reconstituted high density lipoprotein (rHDL) showed that mutations of N-terminal conserved basic residues affect LCAT activity more than those in the C-terminal region. This agrees with the highest conservation of the alpha/beta hydrolase fold and structural homology with pancreatic lipase observed for the N-terminal region, and with the location of most of the natural mutants reported for human LCAT. The structural homology between LCAT and pancreatic lipase further suggests that residues R80, R147, and D145 of LCAT might correspond to residues R37, K107, and D105 of pancreatic lipase, which form the salt bridges D105-K107 and D105-R37. Natural and engineered mutations at residues R80, D145, and R147 of LCAT are accompanied by a substantial decrease or loss of activity, suggesting that salt bridges between these residues might contribute to the structural stability of the enzyme.     

The 1.7 A crystal structure of BPI: a study of how two dissimilar amino acid sequences can adopt the same fold

We have extended the resolution of the crystal structure of human bactericidal/permeability-increasing protein (BPI) to 1.7 A. BPI has two domains with the same fold, but with little sequence similarity. To understand the similarity in structure of the two domains, we compare the corresponding residue positions in the two domains by the method of 3D-1D profiles. A 3D-1D profile is a string formed by assigning each position in the 3D structure to one of 18 environment classes. The environment classes are defined by the local secondary structure, the area of the residue which is buried from solvent, and the fraction of the area buried by polar atoms. A structural alignment between the two BPI domains was used to compare the 3D-1D environments of structurally equivalent positions. Greater than 31% of the aligned positions have conserved 3D-1D environments, but only 13% have conserved residue identities. Analysis of the 3D-1D environmentally conserved positions helps to identify pairs of residues likely to be important in conserving the fold, regardless of the residue similarity. We find examples of 3D-1D environmentally conserved positions with dissimilar residues which nevertheless play similar structural roles. To generalize our findings, we analyzed four other proteins with similar structures yet dissimilar sequences. Together, these examples show that aligned pairs of dissimilar residues often share similar structural roles, stabilizing dissimilar sequences in the same fold.     

Molecular Basis of Substrate Polyspecificity of the Candida albicans Mdr1p Multidrug/H+ Antiporter

The molecular basis of polyspecificity of Mdr1p, a major drug/H⁺ antiporter of Candida albicans, is not elucidated. We have probed the nature of the drug-binding pocket by performing systematic mutagenesis of the 12 transmembrane segments. Replacement of the 252 amino acid residues with alanine or glycine yielded 2/3 neutral mutations while 1/3 led to the complete or selective loss of resistance to drugs or substrates transported by the pump. Using the GlpT-based 3D–model of Mdr1p, we roughly categorized these critical residues depending on their type and localization, 1°/ main structural impact (“S” group), 2°/ exposure to the lipid interface (“L” group), 3°/ buried but not facing the main central pocket, inferred as critical for the overall H⁺/drug antiport mechanism (“M” group) and finally 4°/ buried and facing the main central pocket (“B” group). Among “B” category, 13 residues were essential for the large majority of drugs/substrates, while 5 residues were much substrate-specific, suggesting a role in governing polyspecificity (P group). 3D superposition of the substrate-specific MFS Glut1 and XylE with the MDR substrate-polyspecific MdfA and Mdr1p revealed that the B group forms a common substrate interaction core while the P group is only found in the 2 MDR MFS transporters, distributed into 3 areas around the B core. This specific pattern has let us to propose that the structural basis for polyspecificity of MDR MFS transporters is the extended capacity brought by residues located at the periphery of a binding core to accomodate compounds differing in size and type.     

Solution structure of monomeric peptide YY supports the functional significance of the PP-fold

Peptide YY (PYY) belongs to a family of peptides including neuropeptide Y (NPY) and pancreatic peptide (PP) that regulate numerous functions through both central and peripheral receptors. The solution structure of these peptides is hypothesized to be critically important in receptor selectivity and activation, based on prior demonstration of a stable tertiary conformation of PP called the "PP-fold". Circular dichroism (CD) spectra show a pH-dependent structural transition in the pH range 3-4. Thus we describe the tertiary structure of porcine PYY in water at pH 5.5, 25 degrees C, and 150 mM NaCl, as determined from 2D (1)H NMR data recorded at 500 MHz. A constraint set consisting of 396 interproton distances from NOE data was used as input for distance geometry, simulated annealing, and restrained energy minimization calculations in X-PLOR. The RMSDs of the 20 X-PLOR-generated structures were 0.71 +/- 0.14 and 1.16 +/- 0.17 A, respectively, for backbone and heavy atom overlays of residues 1-34. The resulting structure consists of two C-terminal helical segments from residues 17 to 22 and 25 to 33 separated by a kink at residues 23, 24, and 25, a turn centered around residues 12-14, and the N-terminus folded near residues 30 and 31. The well-defined portions of the PYY structure reported here bear a marked similarity to the structure of PP. Our findings strongly support the importance of the stable folded structure of this family of peptides for binding and activation of Y receptor subtypes.     

Structural analysis of a pectic polysaccharide from the leaves of Diospyros kaki

A pectic polysaccharide DL-2A with a molar mass of 8.5 x 10(5), was obtained from the boiling water extract of Diospyros kaki leaves. It had [alpha]20D -21.8 degrees (c 0.22, H2O) and consisted of rhamnose, arabinose, galactose, xylose and galacturonic acid units in the molar ratio of 0.4:3.4:2.4:1.0:0.8, along with traces of glucuronic acid. About 16.7% of galacturonic acid existed as the methyl ester. A combination of linkage analyses, periodate oxidation, partial acid hydrolysis, selective lithium-degraded reaction, ESIMS, 1H- and 13C- NMR spectral analyses revealed its structural features. It was found that DL-2A possessed an alpha-(1-->4)-galacturonan backbone with some insertions of alpha-1,2-Rhap residues. The side-chains of arabino-3,6-galactan were attached to the backbone via O-4 of Rhap residues and O-3 of GalAp residues, while 4-linked xylose residues (forming short linear chains) were directly linked to O-4 of rhamnose residues, not as part of the xylogalacturonan. These novel structural features enlarge the knowledge on the fine structure of pectic substances in the plant kingdom.     

Regulation of Protein Function and Signaling by Reversible Cysteine S-Nitrosylation

NO is a versatile free radical that mediates numerous biological functions within every major organ system. A molecular pathway by which NO accomplishes functional diversity is the selective modification of protein cysteine residues to form S-nitrosocysteine. This post-translational modification, S-nitrosylation, impacts protein function, stability, and location. Despite considerable advances with individual proteins, the in vivo biological chemistry, the structural elements that govern the selective S-nitrosylation of cysteine residues, and the potential overlap with other redox modifications are unknown. In this minireview, we explore the functional features of S-nitrosylation at the proteome level and the structural diversity of endogenously modified residues, and we discuss the potential overlap and complementation that may exist with other cysteine modifications.     

NMR structural studies of interactions of a small, nonpeptidyl Tpo mimic with the thrombopoietin receptor extracellular juxtamembrane and transmembrane domains

Thrombopoietin (Tpo) is a glycoprotein growth factor that supports hematopoietic stem cell survival and expansion and is the principal regulator of megakaryocyte growth and differentiation. Several small, nonpeptidyl molecules have been identified as selective human Tpo receptor (hTpoR) agonists. To understand how the small molecule Tpo mimic SB394725 interacts and activates hTpoR, we performed receptor domain swap and mutagenesis studies. The results suggest that SB394725 interacts specifically with the extracellular juxtamembrane region (JMR) and the transmembrane (TM) domain of hTpoR. Solution and solid-state NMR structural studies using a peptide containing the JMR-TM sequences showed that this region of hTpoR, unexpectedly, consists of two alpha-helices separated by a few nonhelical residues. SB394725 interacts specifically with His-499 in the TM domain and a few distinct residues in the JMR-TM region and affects several specific C-terminal TM domain residues. The unique structural information provided by these studies both sheds light on the distinctive mechanism of action of SB394725 and provides valuable insight into the mechanism of ligand-induced cytokine receptor activation.     

Exploring the Roles of Proline in Three-Dimensional Domain Swapping from Structure Analysis and Molecular Dynamics Simulations

Three-dimensional (3D) domain swapping is a mechanism to form protein oligomers. It has been proposed that several factors, including proline residues in the hinge region, may affect the occurrence of 3D domain swapping. Although introducing prolines into the hinge region has been found to promote domain swapping for some proteins, the opposite effect has also been observed in several studies. So far, how proline affects 3D domain swapping remains elusive. In this work, based on a large set of 3D domain-swapped structures, we performed a systematic analysis to explore the correlation between the presence of proline in the hinge region and the occurrence of 3D domain swapping. We further analyzed the conformations of proline and pre-proline residues to investigate the roles of proline in 3D domain swapping. We found that more than 40% of the domain-swapped structures contained proline residues in the hinge region. Unexpectedly, conformational transitions of proline residues were rarely observed upon domain swapping. Our analyses showed that hinge regions containing proline residues preferred more extended conformations, which may be beneficial for the occurrence of domain swapping by facilitating opening of the exchanged segments.