The role of carbon-donor hydrogen bonds in stabilizing tryptophan conformations

Although most side-chain torsion angles correspond to low-energy rotameric positions, deviations occur with significant frequency. One striking example arises in Trp residues, which have an important role in stabilizing protein structures because of their size and mixed hydrophobic/hydrophilic character. Ten percent of Trp side-chains have unexplained conformations with chi(2) near 0 degrees instead of the expected 90 degrees. The current work is a structural and energetic analysis of these conformations. It is shown that many Trp residues with these orientations are stabilized by three-center carbon-donor hydrogen bonds of the form C-H...X...H-C, where X is a polar hydrogen-bond acceptor in the environment of the side-chain. The bridging hydrogen bonds occur both within the Trp side-chain and between the side-chain and the local protein backbone. Free energy maps of an isolated Trp residue in an explicit water environment show a minimum corresponding to the off-rotamer peak observed in the crystallographic data. Bridging carbon-donor hydrogen bonds are also shown to stabilize on-rotamer Trp conformations, and similar bridging hydrogen bonds also stabilize some off-rotamer Asp conformations. The present results suggest a previously unrecognized role for three-center carbon-donor hydrogen bonds in protein structures and support the view that the off-rotamer Trp side-chain orientations are real rather than artifacts of crystallographic refinements. Certain of the off-rotamer Trp conformations appear to have a functional role.     

Structural Aspects of Plant Ferredoxin : NADP+ Oxidoreductases

Ferredoxin reductase (FNR) is ubiquitous among photosynthetic organisms as the enzyme directly responsible for the generation of NADPH. Structural studies over the last 15 years have generated over 30 crystal structures of wild-type and mutant FNRs that have yielded a great deal of insight into its structure-function relations. These insights are summarized and combined to propose a structurally informed cycle for FNR catalysis in vivo.     

Product-assisted catalysis in base-excision DNA repair

Most spontaneous damage to bases in DNA is corrected through the action of the base-excision DNA repair pathway. Base excision repair is initiated by DNA glycosylases, lesion-specific enzymes that intercept aberrant bases in DNA and catalyze their excision. How such proteins accomplish the feat of catalyzing no fewer than five sequential reaction steps using a single active site has been unknown. To help answer this, we report the structure of a trapped catalytic intermediate in DNA repair by human 8-oxoguanine DNA glycosylase. This structure and supporting biochemical results reveal that the enzyme sequesters the excised lesion base and exploits it as a cofactor to participate in catalysis. To our knowledge, the present example represents the first documented case of product-assisted catalysis in an enzyme-catalyzed reaction.     

NO rebinding to myoglobin: a reactive molecular dynamics study

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.     

The role of glutamine 114 in old yellow enzyme

Glutamine 114 of OYE1 is a well conserved residue in the active site of the Old Yellow Enzyme family. It forms hydrogen bonds to the O2 and N3 of the flavoprotein prosthetic group, FMN. Glutamine 114 was mutated to asparagine, introducing an R-group that is one methylene group shorter. The resultant enzyme was characterized to determine the effect of the mutation on the mechanistic behavior of the enzyme, and the crystal structure was solved to determine the effect of the mutation on the structure of the protein. The Q114N mutation results in little change in the protein structure, moving the amide group of residue 114 out of H-bonding distance, allowing repositioning of the FMN prosthetic group to form new interactions that replace the lost H-bonds. The mutation decreases the ability to bind ligands, as all dissociation constants for substituted phenols are larger than for the wild type enzyme. The rate constant for the reductive half-reaction with beta-NADPH is slightly greater, whereas that for the oxidative half-reaction with 2-cyclohexenone is smaller than for the wild type enzyme. Oxidation with molecular oxygen is biphasic and involves formation and reaction with O(2), a phenomenon that is more pronounced with this mutation than with wild type enzyme. When superoxide dismutase is added to the reaction, we observe a single-phase reaction typical of the wild type enzyme. Turnover reactions using beta-NADPH with 2-cyclohexenone and molecular oxygen were studied to further characterize the mutant enzyme.     

Mechanical unfolding of a titin Ig domain: structure of unfolding intermediate revealed by combining AFM,molecular dynamics simulations,NMR and protein engineering

The mechanical unfolding of an immunoglobulin domain from the human muscle protein titin (TI I27) has been shown to proceed via a metastable intermediate in which the A-strand is detached. The structure and properties of this intermediate are characterised in this study. A conservative destabilising mutation in the A-strand has no effect on the unfolding force, nor the dependence of the unfolding force on the pulling speed, indicating that the unfolding forces measured in an AFM experiment are those required for the unfolding of the intermediate and not the native state. A mutant of TI I27 with the A-strand deleted (TI I27-A) is studied by NMR and standard biophysical techniques, combined with protein engineering. Molecular dynamics simulations show TI I27-A to be a good model for the intermediate. It has a structure very similar to the native state, and is surprisingly stable. Comparison with a Phi-value analysis of the unfolding pathway clearly shows that the protein unfolds by a different pathway under an applied force than on addition of denaturant.     

Crystal structures of two mutants that have implications for the folding of bovine pancreatic ribonuclease A.

The Tyr92-Pro93 peptide group of bovine pancreatic ribonuclease A (RNase A) exists in the cis conformation in the native state. From unfolding/refolding kinetic studies of the disulfide-intact wild-type protein and of a variant in which Pro93 had been replaced by Ala, it had been suggested that the Tyr92-Ala93 peptide group also exists in the cis conformation in the native state. Here, we report the crystal structure of the P93A variant. Although there is disorder in the region of residues 92 and 93, the best structural model contains a cis peptide at this position, lending support to the results of the kinetics experiments. We also report the crystal structure of the C[40, 95]A variant, which is an analog of the major rate-determining three-disulfide intermediate in the oxidative folding of RNase A, missing the 40-95 disulfide bond. As had been detected by NMR spectroscopy, the crystal structure of this analog shows disorder in the region surrounding the missing disulfide. However, the global chain fold of the remainder of the protein, including the disulfide bond between Cys65 and Cys72, appears to be unaffected by the mutation.     

New methods for accurate prediction of protein secondary structure.

A primary and a secondary neural network are applied to secondary structure and structural class prediction for a database of 681 non-homologous protein chains. A new method of decoding the outputs of the secondary structure prediction network is used to produce an estimate of the probability of finding each type of secondary structure at every position in the sequence. In addition to providing a reliable estimate of the accuracy of the predictions, this method gives a more accurate Q3 (74.6%) than the cutoff method which is commonly used. Use of these predictions in jury methods improves the Q3 to 74.8%, the best available at present. On a database of 126 proteins commonly used for comparison of prediction methods, the jury predictions are 76.6% accurate. An estimate of the overall Q3 for a given sequence is made by averaging the estimated accuracy of the prediction over all residues in the sequence. As an example, the analysis is applied to the target beta-cryptogein, which was a difficult target for ab initio predictions in the CASP2 study; it shows that the prediction made with the present method (62% of residues correct) is close to the expected accuracy (66%) for this protein. The larger database and use of a new network training protocol also improve structural class prediction accuracy to 86%, relative to 80% obtained previously. Secondary structure content is predicted with accuracy comparable to that obtained with spectroscopic methods, such as vibrational or electronic circular dichroism and Fourier transform infrared spectroscopy.     

Use of quantitative structure-property relationships to predict the folding ability of model proteins

We investigate the folding of a 125-bead heteropolymer model for proteins subject to Monte Carlo dynamics on a simple cubic lattice. Detailed study of a few sequences revealed a folding mechanism consisting of a rapid collapse followed by a slow search for a stable core that served as the transition state for folding to a near-native intermediate. Rearrangement from the intermediate to the native state slowed folding further because it required breaking native-like local structure between surface monomers so that those residues could condense onto the core. We demonstrate here the generality of this mechanism by a statistical analysis of a 200 sequence database using a method that employs a genetic algorithm to pick the sequence attributes that are most important for folding and an artificial neural network to derive the corresponding functional dependence of folding ability on the chosen sequence attributes [quantitative structure-property relationships (QSPRs)]. QSPRs that use three sequence attributes yielded substantially more accurate predictions than those that use only one. The results suggest that efficient search for the core is dependent on both the native state's overall stability and its amount of kinetically accessible, cooperative structure, whereas rearrangement from the intermediate is facilitated by destabilization of contacts between surface monomers. Implications for folding and design are discussed. Proteins 33:177–203, 1998. © 1998 Wiley-Liss, Inc.     

Multiscale affinity maturation simulations to elicit broadly neutralizing antibodies against HIV

The design of vaccines against highly mutable pathogens, such as HIV and influenza, requires a detailed understanding of how the adaptive immune system responds to encountering multiple variant antigens (Ags). Here, we describe a multiscale model of B cell receptor (BCR) affinity maturation that employs actual BCR nucleotide sequences and treats BCR/Ag interactions in atomistic detail. We apply the model to simulate the maturation of a broadly neutralizing Ab (bnAb) against HIV. Starting from a germline precursor sequence of the VRC01 anti-HIV Ab, we simulate BCR evolution in response to different vaccination protocols and different Ags, which were previously designed by us. The simulation results provide qualitative guidelines for future vaccine design and reveal unique insights into bnAb evolution against the CD4 binding site of HIV. Our model makes possible direct comparisons of simulated BCR populations with results of deep sequencing data, which will be explored in future applications.