Prediction of electrostatic effects of engineering of protein charges

Accurate prediction of electrostatic effects on catalytic activity is an essential component of protein design. Site-directed mutagenesis of charged groups in subtilisin of Bacillus amyloliquefaciens has provided experimental measurements of electrostatic interactions which may be used to test such theoretical methods. The pKa of the histidine of the active site has been perturbed by +0.08 to -1.0 units by modifying one or two residues. Electrostatic effects in proteins can be modelled by the algorithm of Warwicker and Watson, which uses classical electrostatics and considers both the charge position and the shape of the molecule. Here we report that the algorithm can model several pKa shifts in subtilisin to fair accuracy.     

Capping and alpha-helix stability

The first and last four residues of alpha-helices differ from the rest by not being able to make the intrehelical hydrogen bonds between the backbone greater than C=O groups of one turn and the greater than NH groups of the next. Physico-chemical arguments and statistical analysis suggest that there is a preference for certain residues at the C and N termini (The C- and N-caps) that can fulfil the hydrogen bonding requirements. We have tested this hypothesis by constructing a series of mutations in the two N-caps of barnase (Bacillus amyloliquefaciens ribonuclease, positions Thr 6 and Thr 26) and determining the change in their stability. The N-cap is found to stabilize the protein by up to approximately 2.5 kcal mol(-1). The presence of a negative charge of the N-cap adds some 1.6 kcal mol(-1) of stabilization energy because of the interaction with the macroscopic electrostatic dipole of the helix.     

Long-range structural changes in proteinase K triggered by calcium ion removal

The X-ray crystal structure of the subtilisin-type enzyme proteinase K at 1.5 A resolution shows that is has two binding sites for Ca2+. Scatchard analysis indicates that one Ca2+ binds tightly, with pK 7.6 x 10(-8) M-1, and the other only weakly. Although Ca2+ is not directly involved in the catalytic mechanism and is 16.6 A away from the alpha-carbon atoms of the catalytic triad Asp 39-His 69-Ser 224, the activity of proteinase K towards the synthetic substrate succinyl-Ala-Ala-Ala-p-nitroanilide drops slowly to approximately 20% of its original value when it is depleted of Ca2+. This is not due to autolysis of the enzyme. The X-ray crystal structure of Ca2+-free proteinase K shows that removal of Ca2+ from the tight binding site triggers a concerted domino-like movement of five peripheral loops and of two alpha-helices. At a distance of 25 A from this calcium-binding site, the geometry of both the secondary substrate binding site and of the catalytic triad is affected by this movement thereby reducing the activity of the enzyme.     

Free energy perturbation calculations on binding and catalysis after mutating Asn 155 in subtilisin

Site-directed mutagenesis is a very powerful approach to altering the biological functions of proteins, the structural stability of proteins and the interactions of proteins with other molecules. Several experimental studies in recent years have been directed at estimating the changes in catalytic properties, (rates of binding and catalysis) in site-directed mutants of enzymes compared to the native enzymes. Simulation approaches to the study of complex molecules have also become more powerful, in no small measure owing to the increase in computer power. These simulations have often allowed results of experiments to be rationalized and understood mechanistically. A new approach called the free-energy pertubation method, which uses statistical mechanics and molecular dynamics can often be used for quantitative calculation of free energy differences. We have applied such a technique to calculate the differential free energy of binding and free energy of activation for catalysis of a tripeptide substrate by native subtilisin and a subtilisin mutant (Asn 155----Ala 155). Our studies lead to a calculated difference in free energy of binding which is relatively small, but a calculated change in free energy of catalysis which is substantial. These energies are very close to those determined experimentally (J. A. Wells and D. A. Estell, personal communication), which were not known to us until the simulations were completed. This demonstrates the predictive power and utility of theoretical simulation methods in studies of the effects of site-specific mutagenesis on both enzyme binding and catalysis.     

Identification of human CD4 residues affecting class II MHC versus HIV-1 gp120 binding

Interactions of CD4 with the class II major histocompatibility complex (MHC) are crucial during thymic ontogeny and subsequently for helper and cytotoxic functions of CD4+CD8- T lymphocytes. CD4 is the receptor for the T-lymphotropic human immunodeficiency virus and binds its envelope glycoprotein, gp120. The residues involved in gp120 binding have been localized to a region within the immunoglobulin-like domain I of CD4, which corresponds to CDR2 of an immunoglobulin variable region, but the CD4 residues important in MHC class II interaction have not been characterized. Here, using a cell-binding assay dependent specifically on the CD4-MHC class II association, we analyse the effects of mutations in CD4 on class II versus gp120 binding. Mutations in CDR2 that destroy gp120 binding affect CD4-MHC class II binding similarly. In addition, binding of soluble gp120 to CD4-transfected cells abrogates their ability to interact with class II-bearing B lymphocytes. In contrast, other mutations within domains I or II that have no effect on gp120 binding eliminate or substantially decrease class II interaction. Thus, the CD4 binding site for class II MHC is more complex than the gp120 binding site, possibly reflecting a broader area of contact with the former ligand and a requirement for appropriate juxtaposition of the two N-terminal domains. The ability of gp120 to inhibit the binding of class II MHC to CD4 could be important in disrupting normal T-cell physiology, acting both to inhibit immune responses and to prevent differentiation of CD4+CD8+ thymocytes into CD4+CD8- T lymphocytes.     

Hydrolysis of carboxylate ester catalyzed by a new artificial abzyme based on molecularly imprinted polymer

A new kind of artificial abzyme MIP-3, which contains N-phenyl benzonic amide group and “nanochannel” formed by embedded ZnO nano materials, and is imprinted by a transition-state analogue of p-nitrophenyl methyphosphonate in the hydrolysis of p-nitrophenyl acetate, was prepared by radical co-polymerization. Hydrolytic kinetics of p-nitrophenyl carboxylate catalyzed by MIP-3 was investigated. The results showed that the artificial abzyme exhibited notable substructure selectivity and strong catalytic ability in hydrolysis of p-nitrophenyl acetate     

Cleavage of pre-tRNAs by the splicing endonuclease requires a composite active site

Splicing is required for the removal of introns from a subset of transfer RNAs in all eukaryotic organisms. The first step of splicing, intron recognition and cleavage, is performed by the tRNA-splicing endonuclease, a tetrameric enzyme composed of the protein subunits Sen54, Sen2, Sen34 and Sen15. It has previously been demonstrated that the active sites for cleavage at the 5' and 3' splice sites of precursor tRNA are contained within Sen2 and Sen34, respectively. A recent structure of an archaeal endonuclease complexed with a bulge-helix-bulge RNA has led to the unexpected hypothesis that catalysis requires a critical 'cation-pi sandwich' composed of two arginine residues that serve to position the RNA substrate within the active site. This motif is derived from a cross-subunit interaction between the two catalytic subunits. Here we test the role of this interaction within the eukaryotic endonuclease and show that catalysis at the 5' splice site requires the conserved cation-pi sandwich derived from the Sen34 subunit in addition to the catalytic triad of Sen2. The catalysis of pre-tRNA by the eukaryotic tRNA-splicing endonuclease therefore requires a previously unrecognized composite active site.     

邻苯二甲酸二丁酯与胰蛋白酶的相互作用

在模拟人体生理条件下(pH值为7.4),应用荧光光谱法、紫外光谱法和圆二色谱法(CD)并结合原子力显微镜(AFM)和分子模拟技术,研究了邻苯二甲酸二丁酯(DBP)对胰蛋白酶(Trypsin)的光谱性质、结构及催化活性的影响.结果表明,DBP通过氢键和范德华力与胰蛋白酶形成基态复合物而猝灭胰蛋白酶的内源荧光,DBP在胰蛋白酶上只有1个结合位点.同步荧光、紫外和CD光谱研究发现,DBP与胰蛋白酶的结合诱导了酶的α-螺旋、β-折叠和β-转角含量的减少,而增加了无规卷曲的含量.原子力显微镜图像显示,DBP的存在引起胰蛋白酶的表面形态发生变化,蛋白质发生了聚集.分子模拟结果表明,DBP结合于胰蛋白酶S1疏水空腔附近,与氨基酸His 57,Ser 195和Gly 193形成氢键.酶活测定结果显示,DBP的存在导致胰蛋白酶活性被抑制.     

Molecular mechanisms of serpin deficiencies

The study of serpin deficiency is currently one of the most active areas in basic medical research. Recently, three hypotheses concerning serpin deficiency have been proposed, which are referred to as the conformational disturbance hypothesis (CDH) , loop-sheet polymerisation hypothesis (LSPH) and multiple binding site hypothesis (MB-SH) . CDH was put forward to explicit serpin deficiency due to conformational change of reactive loop of serpins as a result of mutations occurring away from the reactive site residues and LSPH was to explain deficient serpins due to the formation of polymers. MBSH was proposed to explain the mechanism of the formation of stable enzyme-serpin complex via more than one binding site and blockage or mutation in any of the sites resulting in serpin deficiency. A combination of these mechanisms may be critical in understanding the roles of the many documented mutations and autoimmunities which result in qualitative and quantitative serpin deficiency.     

Synthesis and structural analysis of 2-quinuclidonium tetrafluoroborate

The amide functional group is one of the most fundamental motifs found in chemistry and biology, and it has been studied extensively for the past century. Typical acyclic amides are planar. But the amide groups of bicyclic bridgehead lactams are highly twisted, and this distortion from planarity can dramatically affect the stability and reactivity of these amides; it also increases the basicity of the nitrogen so that it often behaves more like an amine than a typical planar amide. As a result, the structures and reactivity profiles of these 'anti-Bredt' amides differ significantly from those of planar amides. It is possible that this twisting phenomenon is not exclusive to cyclic systems-non-planarity may also be a critical biological design element that leads to amide ground-state destabilization and alters the reactivity, selectivity and mechanism of various protein and enzymatic processes (such as amide hydrolysis). The intriguing qualities of these twisted amides were first recognized in 1938 (ref. 11), wherein one of the simplest families was introduced--molecules containing the 1-azabicyclo[2.2.2]octan-2-one system. But the parent member of this group, 2-quinuclidone (molecule 1 in this paper), has not yet been unambiguously synthesized. Here, we report the chemical synthesis, isolation and full characterization of the HBF4 salt of 1. Critical to the success of the synthesis and isolation was the decision to generate 1 by a route other than classical amide bond formation. We anticipate that these results will provide a greater understanding of the properties of amide bonds.     

The structure and catalytic mechanism of a poly(ADP-ribose) glycohydrolase

Post-translational modification of proteins by poly(ADP-ribosyl)ation regulates many cellular pathways that are critical for genome stability, including DNA repair, chromatin structure, mitosis and apoptosis. Poly(ADP-ribose) (PAR) is composed of repeating ADP-ribose units linked via a unique glycosidic ribose-ribose bond, and is synthesized from NAD by PAR polymerases. PAR glycohydrolase (PARG) is the only protein capable of specific hydrolysis of the ribose-ribose bonds present in PAR chains; its deficiency leads to cell death. Here we show that filamentous fungi and a number of bacteria possess a divergent form of PARG that has all the main characteristics of the human PARG enzyme. We present the first PARG crystal structure (derived from the bacterium Thermomonospora curvata), which reveals that the PARG catalytic domain is a distant member of the ubiquitous ADP-ribose-binding macrodomain family. High-resolution structures of T. curvata PARG in complexes with ADP-ribose and the PARG inhibitor ADP-HPD, complemented by biochemical studies, allow us to propose a model for PAR binding and catalysis by PARG. The insights into the PARG structure and catalytic mechanism should greatly improve our understanding of how PARG activity controls reversible protein poly(ADP-ribosyl)ation and potentially of how the defects in this regulation are linked to human disease.     

捕食线虫真菌丝氨酸蛋白酶PII结构的同源模建

通过同源模建方法构建捕食线虫真菌丝氨酸蛋白酶PII的三维结构模型,并结合以往研究对其结构进行分析。结果表明,蛋白酶PII的结构具有枯草杆菌丝氨酸蛋白酶特有的a／B脚手架折叠模式,具有保守的催化三聚体和氧负离子孔结构组织,整体表现为坚固的球状折叠构型,不具备二硫键。研究结果将为进一步深入研究PII的结构和功能之间的关系奠定结构基础。     

RanGAP mediates GTP hydrolysis without an arginine finger

GTPase-activating proteins (GAPs) increase the rate of GTP hydrolysis on guanine nucleotide-binding proteins by many orders of magnitude. Studies with Ras and Rho have elucidated the mechanism of GAP action by showing that their catalytic machinery is both stabilized by GAP binding and complemented by the insertion of a so-called 'arginine finger' into the phosphate-binding pocket. This has been proposed as a universal mechanism for GAP-mediated GTP hydrolysis. Ran is a nuclear Ras-related protein that regulates both transport between the nucleus and cytoplasm during interphase, and formation of the mitotic spindle and/or nuclear envelope in dividing cells. Ran-GTP is hydrolysed by the combined action of Ran-binding proteins (RanBPs) and RanGAP. Here we present the three-dimensional structure of a Ran-RanBP1-RanGAP ternary complex in the ground state and in a transition-state mimic. The structure and biochemical experiments show that RanGAP does not act through an arginine finger, that the basic machinery for fast GTP hydrolysis is provided exclusively by Ran and that correct positioning of the catalytic glutamine is essential for catalysis.     

Context-dependent contributions of backbone hydrogen bonding to beta-sheet folding energetics

Backbone hydrogen bonds (H-bonds) are prominent features of protein structures; however, their role in protein folding remains controversial because they cannot be selectively perturbed by traditional methods of protein mutagenesis. Here we have assessed the contribution of backbone H-bonds to the folding kinetics and thermodynamics of the PIN WW domain, a small beta-sheet protein, by individually replacing its backbone amides with esters. Amide-to-ester mutations site-specifically perturb backbone H-bonds in two ways: a H-bond donor is eliminated by replacing an amide NH with an ester oxygen, and a H-bond acceptor is weakened by replacing an amide carbonyl with an ester carbonyl. We perturbed the 11 backbone H-bonds of the PIN WW domain by synthesizing 19 amide-to-ester mutants. Thermodynamic studies on these variants show that the protein is most destabilized when H-bonds that are enveloped by a hydrophobic cluster are perturbed. Kinetic studies indicate that native-like secondary structure forms in one of the protein's loops in the folding transition state, but the backbone is less ordered elsewhere in the sequence. Collectively, our results provide an unusually detailed picture of the folding of a beta-sheet protein.     

Active-site remodelling in the bifunctional fructose-1,6-bisphosphate aldolase/phosphatase

Fructose-1,6-bisphosphate (FBP) aldolase/phosphatase is a bifunctional, thermostable enzyme that catalyses two subsequent steps in gluconeogenesis in most archaea and in deeply branching bacterial lineages. It mediates the aldol condensation of heat-labile dihydroxyacetone phosphate (DHAP) and glyceraldehyde-3-phosphate (GAP) to FBP, as well as the subsequent, irreversible hydrolysis of the product to yield the stable fructose-6-phosphate (F6P) and inorganic phosphate; no reaction intermediates are released. Here we present a series of structural snapshots of the reaction that reveal a substantial remodelling of the active site through the movement of loop regions that create different catalytic functionalities at the same location. We have solved the three-dimensional structures of FBP aldolase/phosphatase from thermophilic Thermoproteus neutrophilus in a ligand-free state as well as in complex with the substrates DHAP and FBP and the product F6P to resolutions up to 1.3??. In conjunction with mutagenesis data, this pinpoints the residues required for the two reaction steps and shows that the sequential binding of additional Mg(2+) cations reversibly facilitates the reaction. FBP aldolase/phosphatase is an ancestral gluconeogenic enzyme optimized for high ambient temperatures, and our work resolves how consecutive structural rearrangements reorganize the catalytic centre of the protein to carry out two canonical reactions in a very non-canonical type of bifunctionality.     

Ser-His-Glu triad forms the catalytic site of the lipase from Geotrichum candidum

The Ser-His-Asp triad is a well known structural feature of the serine proteases. It has also been directly observed in the catalytic sites of two lipases, whose high-resolution three-dimensional structures have been determined 1,2. Lipases show a wide variety of sizes, substrate and positional specificities, and catalytic rates 3. They achieve maximal catalytic rates at oil-water interfaces. The fungus Geotrichum candidum produces several different forms of lipases, two of which have been purified to homogeneity 4,5. Two lipase genes have been identified, cloned and sequenced 6,7. Both code for proteins of 544 amino acids with a total relative molecular mass of about 60,000 (Mr 60K). The two forms are 86% identical. Their isoelectric points differ slightly, being between 4.3 and 4.6. About 7% of the total Mr is carbohydrate. Until now, only a low resolution structure of GCL has been reported 8, but no high resolution structure has followed. We now report the three-dimensional structure of a lipase from G. candidum (GCL) at 2.2 A resolution. Unlike the other lipases and serine proteases, the catalytic triad of GCL is Ser-His-Glu, with glutamic acid replacing the usual aspartate. Although the sequence similarity with the other two lipases is limited to the region near the active-site serine, there is some similarity in their three-dimensional structures. The GCL is also an alpha/beta protein with a central mixed beta sheet whose topology is similar to that of the N-terminal domain of human pancreatic lipase. As in the other lipases 1,2, the catalytic site is buried under surface loops. Sequence comparisons with proteins from the cholinesterase family suggest that they also contain the Ser-His-Glu triad.     

TBC-domain GAPs for Rab GTPases accelerate GTP hydrolysis by a dual-finger mechanism

Rab GTPases regulate membrane trafficking by cycling between inactive (GDP-bound) and active (GTP-bound) conformations. The duration of the active state is limited by GTPase-activating proteins (GAPs), which accelerate the slow intrinsic rate of GTP hydrolysis. Proteins containing TBC (Tre-2, Bub2 and Cdc16) domains are broadly conserved in eukaryotic organisms and function as GAPs for Rab GTPases as well as GTPases that control cytokinesis. An exposed arginine residue is a critical determinant of GAP activity in vitro and in vivo. It has been expected that the catalytic mechanism of TBC domains would parallel that of Ras and Rho family GAPs. Here we report crystallographic, mutational and functional analyses of complexes between Rab GTPases and the TBC domain of Gyp1p. In the crystal structure of a TBC-domain-Rab-GTPase-aluminium fluoride complex, which approximates the transition-state intermediate for GTP hydrolysis, the TBC domain supplies two catalytic residues in trans, an arginine finger analogous to Ras/Rho family GAPs and a glutamine finger that substitutes for the glutamine in the DxxGQ motif of the GTPase. The glutamine from the Rab GTPase does not stabilize the transition state as expected but instead interacts with the TBC domain. Strong conservation of both catalytic fingers indicates that most TBC-domain GAPs may accelerate GTP hydrolysis by a similar dual-finger mechanism.     

Stabilization of charges on isolated ionic groups sequestered in proteins by polarized peptide units

Electrostatic interactions are of considerable importance in protein structure and function, and in a variety of cellular and biochemical processes. Here we report three similar findings from highly refined atomic structures of periplasmic binding proteins. Hydrogen bonds, acting primarily through backbone peptide units, are mainly responsible for the involvement of the positively charged arginine 151 residue in the ligand site of the arabinose-binding protein, for the association between teh sulphate-binding protein and the completely buried sulphate dianion, and for the formation of the complex of the leucine/isoleucine/valine-binding protein with the leucine zwitterion. We propose a general mechanism in which the isolated charges on the various buried, desolvated ionic groups are stabilized by the polarized peptide units. This mechanism also has broad application to processes requiring binding of uncompensated ions and charged ligands and stabilization of enzyme reaction charged intermediates, as well as activation of catalytic residues.