Rate limiting step of the allosteric activation of the bacterial adhesin FimH investigated by molecular dynamics simulations

The bacterial adhesin FimH is a model for the study of protein allostery because its structure has been resolved in multiple configurations, including the active and the inactive state. FimH consists of a pilin domain (PD) that anchors it to the rest of the fimbria and an allosterically regulated lectin domain (LD) that binds mannose on the surface of infected cells. Under normal conditions, the two domains are docked to each other and LD binds mannose weakly. However, in the presence of tensile force generated by shear the domains separate and conformational changes propagate across LD resulting in a stronger bond to mannose. Recently, the crystallographic structure of a variant of FimH has been resolved, called ${\text{FimH}}^{\text{FocH}}$, where PD contains 10 mutations near the inter-domain interface. Although the X-ray structures of FimH and ${\text{FimH}}^{\text{FocH}}$ are almost identical, experimental evidence shows that ${\text{FimH}}^{\text{FocH}}$ is activated even in the absence of shear. Here, molecular dynamics simulations combined with the Jarzynski equality were used to investigate the discrepancy between the crystallographic structures and the functional assays. The results indicate that the free energy barrier of the unbinding process between LD and PD is drastically reduced in ${\text{FimH}}^{\text{FocH}}$. Rupture of inter-domain hydrogen bonds involving R166 constitutes a rate limiting step of the domain separation process and occurs more readily in ${\text{FimH}}^{\text{FocH}}$ than FimH. In conclusion, the mutations in ${\text{FimH}}^{\text{FocH}}$ shift the equilibrium toward an equal occupancy of bound and unbound states for LD and PD by reducing a rate limiting step.     

Crystal structure of inorganic pyrophosphatase from Thermus thermophilus.

The 3-dimensional structure of inorganic pyrophosphatase from Thermus thermophilus (T-PPase) has been determined by X-ray diffraction at 2.0 A resolution and refined to R = 15.3%. The structure consists of an antiparallel closed beta-sheet and 2 alpha-helices and resembles that of the yeast enzyme in spite of the large difference in size (174 and 286 residues, respectively), little sequence similarity beyond the active center (about 20%), and different oligomeric organization (hexameric and dimeric, respectively). The similarity of the polypeptide folding in the 2 PPases provides a very strong argument in favor of an evolutionary relationship between the yeast and bacterial enzymes. The same Greek-key topology of the 5-stranded beta-barrel was found in the OB-fold proteins, the bacteriophage gene-5 DNA-binding protein, toxic-shock syndrome toxin-1, and the major cold-shock protein of Bacillus subtilis. Moreover, all known nucleotide-binding sites in these proteins are located on the same side of the beta-barrel as the active center in T-PPase. Analysis of the active center of T-PPase revealed 17 residues of potential functional importance, 16 of which are strictly conserved in all sequences of soluble PPases. Their possible role in the catalytic mechanism is discussed on the basis of the present crystal structure and with respect to site-directed mutagenesis studies on the Escherichia coli enzyme. The observed oligomeric organization of T-PPase allows us to suggest a possible mechanism for the allosteric regulation of hexameric PPases.     

Control of K+ (Rb+) influx and transport with changed internal K+ concentration and age in two varieties of barley

The influx of Rb⁺ into the roots of two barley varieties (Hordeum vulgare L. cv. Salve and cv. Ingrid) from a K⁺-free ⁸⁶Rb-labelled nutrient solution with 2.0 mM Rb⁺, was checked at intervals from day 6 to day 18. The control plants were continuously grown in complete nutrient solution containing 5.0 mM K⁺, while two other groups of plants were grown in K⁺-free nutrient solution starting on day 6 and between day 6 and day 9, respectively. The pattern of Rb⁺ influx was similar for both varieties, although their efficiencies in absorbing Rb⁺ were different. The relationship between Rb⁺ influx and K⁺ concentration of the root could be interpreted in terms of negative feedback through allosteric control of uptake across the plasmalemma of the root cells. Hill plots were bimodal, but in the opposite direction. The Hill coefficients, reflecting the minimum number of interacting allosteric binding sites for K⁺ (Rb⁺), were low (≤–3.0). It is discussed whether the threshold value, that is the breaking point in the Hill plot, is indicative of a changed efficiency of transporting units for K⁺ (Rb⁺) transport to the xylem. Moreover, feedback regulation might be involved in transport of K⁺ between root and shoot. The variation in K⁺ concentrations in the roots and shoots of control plants were cyclic but in phase opposition despite an exponential growth. The average K⁺ concentration varied only slightly with age.     

Structure of D-glyceraldehyde-3-phosphate dehydrogenase fromPalinurus versicolor in a tetragonal crystal form

D-glyceraldehyde-3-phosphate dehydrogenase (holo-GAPDH) from Palinurus versicolor was crystallized in a novel crystal form by the method of sitting-drop vapor diffusion. The crystals have space group P4212, cell parameters a=15.49 nm, c=8.03 nm and two subunits per asymmetric unit. The crystal structure at 0.34 nm was determined by the molecular replacement method. The final model has crystallographic Rfree and R factors of 0.274 and 0.262, and r.m.s. deviations of 0.002 nm for bond lengths and 2.33?for bond angles. The two subunits in asymmetric unit are similar to each other not only in the three-dimensional structure, but also in average temperature factors. This result demonstrates that the obvious difference in average temperature factors for the different subunits in C2 crystal form reported previously may be attributed to the different crystallographic environments of the subunits. This further supports that holo-GAPDH has a good 222 molecular symmetry.     

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.     

开启防御之门: 植物抗病小体

NLR蛋白是存在于植物和动物中的一个免疫受体大家族, 具有核苷酸结合域并富含亮氨酸重复序列。植物NLR通过识别病原菌特异效应子开启免疫信号转导。第1个植物NLR抗性蛋白于25年前克隆, 但其激活机制仍不清楚, 至今仍未获得一个完整的NLR蛋白结构。最近, 柴继杰、周俭民和王宏伟实验室合作解析了第一个植物完整NLR ZAR1激活前后的结构, 研究成果以两篇论文形式发表在“科学”杂志上, 填补了NLR介导的免疫信号转导研究领域的空白。该文简要总结了相关研究进展, 讨论了NLR免疫信号转导研究领域尚需解决的问题。     

Development of strategies for conditional RNA interference

BACKGROUND: RNA interference (RNAi) represents a powerful tool with which to undertake sequence-dependent suppression of gene expression. Synthesized double-stranded RNA (dsRNA) or dsRNA generated endogenously from plasmid or viral vectors can be used for RNAi. For the latter, polymerase III promoters which drive ubiquitous expression in all tissues have typically been adopted. Given that dsRNA molecules must contain few 5' and 3' over-hanging bases to maintain potency, employing polymerase II promoters to drive tissue-specific expression of RNAi may be problematic due to potential inclusion of nucleotides 5' and 3' of siRNA sequences. METHODS: To circumvent this, polymerase II promoters in combination with cis-acting hammerhead ribozymes and short-hairpin RNA sequences have been explored as a means to generate potent dsRNA molecules in tissues defined by the promoter in use. RESULTS: The novel constructs evaluated in this study produced functional siRNA which suppressed the enhanced green fluorescent protein (eGFP) both in vitro and in vivo (in mice). Additionally, the constructs did not appear to elicit a significant type-1 interferon response compared to traditional H1-transcribed shRNA. CONCLUSIONS: Given the potential 'off-target' effects of dsRNAs, it would be preferable in many cases to limit expression of dsRNA to the tissue of interest and moreover would significantly augment the resolution of RNAi technologies. Notably, the system under evaluation in this study could readily be adapted to achieve this objective.     

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.