A double catgrip mixed L and D mini protein only 20 residues long

Stereochemistry limits but also defines proteins, as conformational constructs stereospecific for poly-L structure. Employed as a variable in sequence, stereochemistry could make proteins customizable in the letters of L and D amino acid alphabet. In proof of concept, we previously demonstrated stereochemical reengineering of canonical beta-hairpins as bracelet and boat shaped molecules. Illustrating the prospect for functional design, a 20-residue four-stranded mini-beta protein is now customized stereochemically as a canoe shaped molecule. A conformational construct of four side by side hydrogen-bonded strands in alternately (l)beta, (d)beta conformation, joined via Type-II/II' beta-turns, is planned to be preponderantly apolar in beta-sheet favoring residues, interspersing two ion pairs, and suitably L and D in sequence. Synthesis followed by MD, NMR, CD, and MALDI-MS studies established the molecule as a canoe shaped fold in water, demonstrable in affinity of alkali and alkaline-earth metal ions as expected given its catgrip like elements. Another success in accomplishing a synthetic miniprotein complex in stereochemistry and stereospecific in conformation, exceptionally small yet functional in metal ion affinity, affirms the value in combined L and D alphabet in programming molecular shapes and functions stereochemically.     

Synthesis and characterization of polyamides containing unnatural amino acids

We describe the synthesis of several polyamides that retain the secondary structure of proteins and contain derivatizable side chains. The derivatizable side chain allows for further reaction of the polymer chain (e.g., chain cross-linking or addition of pendant groups). Polymers of α-amino acids containing a terminal unsaturated bond on the side chain have been synthesized. Poly-L -pentenyl glycine, poly-L -propargyl glycine, and poly-L -allyl glycine were synthesized chemically via Leuchs' anhydrides and enzymatically using subtilisin Carlsberg. Poly-L -propargyl glycine and poly-D ,L -allyl glycine folded into the β-sheet configuration whereas poly-L -pentenyl glycine assumed a helical conformation. The secondary structure of poly-L -allyl glycine and poly-D ,L -pentenylglycine could not be determined conclusively. Comparison of properties between the polymers obtained chemically and enzymatically is provided. © 1995 John Wiley & Sons, Inc.     

The geometry and efficacy of cation-pi interactions in a diagonal position of a designed beta-hairpin

Cation-pi interactions are common in proteins, but their contribution to the stability and specificity of protein structure has not been well established. In this study, we examined the impact of cation-pi interactions in a diagonal position of a beta-hairpin peptide through comparison of the interaction of Phe or Trp with Lys or Arg. The diagonal interactions ranged from -0.20 to -0.48 kcal/mole. Our experimental values for the diagonal cation-pi interactions are similar to those found in alpha-helical studies. Upfield shifting of the Lys and Arg side chains indicates that the geometries of cation-pi interactions adopted in the 12-residue beta-hairpin are comparable to those found in protein structures. The Lys was found to interact through the polarized Cepsilon, and the Arg is stacked against the aromatic ring of Phe or Trp. Folding of these peptides was found to be enthalpically favorable (DeltaH degrees equals approximately -3 kcal/mole) and entropically unfavorable (DeltaS degrees equals approximately -8 cal mole(-1) K(-1)).     

Stereochemical quality of protein structure coordinates.

Methods have been developed to assess the stereochemical quality of any protein structure both globally and locally using various criteria. Several parameters can be derived from the coordinates of a given structure. Global parameters include the distribution of phi, psi and chi 1 torsion angles, and hydrogen bond energies. There are clear correlations between these parameters and resolution; as the resolution improves, the distribution of the parameters becomes more clustered. These features show a broad distribution about ideal values derived from high-resolution structures. Some structures have tightly clustered distributions even at relatively low resolutions, while others show abnormal scatter though the data go to high resolution. Additional indicators of local irregularity include proline phi angles, peptide bond planarities, disulfide bond lengths, and their chi 3 torsion angles. These stereochemical parameters have been used to generate measures of stereochemical quality which provide a simple guide as to the reliability of a structure, in addition to the most important measures, resolution and R-factor. The parameters used in this evaluation are not novel, and are easily calculated from structure coordinates. A program suite is currently being developed which will quickly check a given structure, highlighting unusual stereochemistry and possible errors.     

The interplay of sequence and stereochemistry in defining conformation in proteins and polypeptides

Sequential specification of conformation in proteins and polypeptides is a triangular interplay involving the system of linked peptides, the sequences in side chains, and water as solvent. Stereochemistry in side chain linkages is obviously important in the interaction between all of the players, but no specification of its explicit role, if any, in linking sequence with conformation has been made. Flory and coworkers made a puzzling observation in 1967 that, when mutated from poly-L to alternating-L,D stereochemical structure, polypeptides will suffer a reduction in overall dimension or characteristic ratio by an astonishing factor of 10 and to a value even lower than that predicted for free rotation (Miller, W. G.; Brant, D. A.; Flory, P. J. J Mol Biol 1967, 23, 67-80). Enquiring into this longstanding puzzle, Durani and coworkers found that the stereochemical modification will also abolish conformational sensitivity in polypeptide structure to solvent, because electrostatic interactions in the system of linked peptides are transformed from a condition of mutual conflict to one of harmony (Ramakrishnan, V.; Ranbhor, R.; Kumar, A.; Durani, S. J Phys Chem B 2006, 110, 9314-9323). Thus, poly-L stereochemistry could be the fulcrum linking sequences with phi,psis in protein and polypeptide structures, via dielectric arbitrations in a conflicting type of interpeptide electrostatics, in agreement with the electrostatic screening model of Avbelj and Moult (Avbelj, F.; Moult, J. Biochemistry 1995, 34, 755-764).     

A conserved interaction with the chromophore of fluorescent proteins

The chromophore of fluorescent proteins, including the green fluorescent protein (GFP), contains a highly conjugated imidazolidinone ring. In many fluorescent proteins, the carbonyl group of the imidazolidinone ring engages in a hydrogen bond with the side chain of an arginine residue. Prior studies have indicated that such an electrophilic carbonyl group in a protein often accepts electron density from a main-chain oxygen. A survey of high-resolution structures of fluorescent proteins indicates that electron lone pairs of a main-chain oxygen-Thr62 in GFP-donate electron density into an antibonding orbital of the imidazolidinone carbonyl group. This n→π* electron delocalization prevents structural distortion during chromophore excitation that could otherwise lead to fluorescence quenching. In addition, this interaction is present in on-pathway intermediates leading to the chromophore, and thus could direct its biogenesis. Accordingly, this n→π* interaction merits inclusion in computational and photophysical analyses of the chromophore, and in speculations about the molecular evolution of fluorescent proteins.     

Solution structure of the catalytic domain of human stromelysin complexed with a hydrophobic inhibitor.

Stromelysin, a representative matrix metalloproteinase and target of drug development efforts, plays a prominent role in the pathological proteolysis associated with arthritis and secondarily in that of cancer metastasis and invasion. To provide a structural template to aid the development of therapeutic inhibitors, we have determined a medium-resolution structure of a 20-kDa complex of human stromelysin's catalytic domain with a hydrophobic peptidic inhibitor using multinuclear, multidimensional NMR spectroscopy. This domain of this zinc hydrolase contains a mixed beta-sheet comprising one antiparallel strand and four parallel strands, three helices, and a methionine-containing turn near the catalytic center. The ensemble of 20 structures was calculated using, on average, 8 interresidue NOE restraints per residue for the 166-residue protein fragment complexed with a 4-residue substrate analogue. The mean RMS deviation (RMSD) to the average structure for backbone heavy atoms is 0.91 A and for all heavy atoms is 1.42 A. The structure has good stereochemical properties, including its backbone torsion angles. The beta-sheet and alpha-helices of the catalytic domains of human stromelysin (NMR model) and human fibroblast collagenase (X-ray crystallographic model of Lovejoy B et al., 1994b, Biochemistry 33:8207-8217) superimpose well, having a pairwise RMSD for backbone heavy atoms of 2.28 A when three loop segments are disregarded. The hydroxamate-substituted inhibitor binds across the hydrophobic active site of stromelysin in an extended conformation. The first hydrophobic side chain is deeply buried in the principal S'1 subsite, the second hydrophobic side chain is located on the opposite side of the inhibitor backbone in the hydrophobic S'2 surface subsite, and a third hydrophobic side chain (P'3) lies at the surface.     

Evaluation of "click" binaphthyl chiral stationary phases by liquid chromatography

Two "click" binaphthyl chiral stationary phases were synthesized and evaluated by liquid chromatography. Their structures incorporate S-(-)-1,1'-binaphthyl moiety as the chiral selector and 1,2,3-triazole ring as the spacer. These chiral stationary phases (CSPs) allowed the efficient resolution for a wide range of racemic BINOL derivatives, particularly for nonpolar diether derivatives and 3-phenyl indolin-2-one analogs. The chromatographic data showed that the π-π interaction was crucial for enantiorecognition of these CSPs. Loss of enantioselectivity observed on CSP3, which are lacking the triazole ring linkage, indicated that the triazole ring linkage took part in the enantioseparation process, although it was remote from the chiral selector of the CSP. The substitution of the phenyl group at 6 and 6' positions can significantly improve the separation ability of the CSP. The chiral recognition mechanism was also investigated by tracking the elution orders and studying the thermodynamic parameters.     

Stably folded de novo proteins from a designed combinatorial library

Binary patterning of polar and nonpolar amino acids has been used as the key design feature for constructing large combinatorial libraries of de novo proteins. Each position in a binary patterned sequence is designed explicitly to be either polar or nonpolar; however, the precise identities of these amino acids are varied extensively. The combinatorial underpinnings of the "binary code" strategy preclude explicit design of particular side chains at specified positions. Therefore, packing interactions cannot be specified a priori. To assess whether the binary code strategy can nonetheless produce well-folded de novo proteins, we constructed a second-generation library based upon a new structural scaffold designed to fold into 102-residue four-helix bundles. Characterization of five proteins chosen arbitrarily from this new library revealed that (1) all are alpha-helical and quite stable; (2) four of the five contain an abundance of tertiary interactions indicative of well-ordered structures; and (3) one protein forms a well-folded structure with native-like features. The proteins from this new 102-residue library are substantially more stable and dramatically more native-like than those from an earlier binary patterned library of 74-residue sequences. These findings demonstrate that chain length is a crucial determinant of structural order in libraries of de novo four-helix bundles. Moreover, these results show that the binary code strategy--if applied to an appropriately designed structural scaffold--can generate large collections of stably folded and/or native-like proteins.     

Polymorphism and conformational flexibility of DNA: right and left handed duplexes

Left handed duplexes are shown to be in agreement with the X-ray intensity data of A-, B- and D-forms of DNA. The structures are stereochemically satisfactory because they were obtained following a stereochemical guideline derived from theory and single crystal structure data of nucleic acid components. The same stereochemical guideline also led to right handed duplexes for B- and D-forms of DNA which have stereochemically preferred conformation and hence are superior to those given by Arnott and coworkers^10,11.     

Zinc-finger hydrolase: Computational selection of a linker and a sequence towards metal activation with a synthetic αββ protein

The zinc-finger protein is targeted for computational redesign as a hydrolase enzyme. Successful in having zinc activated for hydrolase function, the study validates the stepwise approach to having the protein tuned in main-chain structure stereochemically and over side chains chemically. A leucine homopolypeptide, harboring histidines to tri coordinate zinc and d-amino-acid-nucleated α-helix and β-hairpin building blocks of an αββ protein, is taken up for modeling, first with cyana, in a mixed-chirality linker between the building blocks, and then with IDeAS, in a sequence over side chains. The designed mixed-chirality polypeptide structure is proven to order as an intended αββ fold and capture zinc to activate its role as a hydrolase catalyst. The design approach to have protein folds defined stereochemically and receptor and catalysis functions defined chemically is presented, and illustrates L- and D-α-amino-acid structures as the alphabet integrating chemical- and stereochemical-structure variables as its letters.     

Engineering of betabellin 15D: Copper(II)-induced folding of a fibrillar β-sandwich protein

The inverse protein-folding problem has been explored by designing de novo the betabellin target structure (a 64-residue -sandwich protein), synthesizing a 32-residue peptide chain (HSLTAKIpkLTFSIAphTYTCAVpkYTAKVSH, where p = DPro, k = DLys, and h = DHis) that might fold into this structure, and studying how its disulfide-bridged form (betabellin 15D) folds in 10 mM ammonium acetate with and without Cu2+. Circular dichroic spectropolarimetry indicated that at pH 5.8, 6.4, or 6.7 betabellin 15D exhibited -sheet structure in the presence of Cu2+ but not in its absence. Electrospray mass spectrometry demonstrated that at pH 6.3 each molecule of betabellin 15D bound one or two Cu(II) ions. Electron microscopy showed that at pH 6.7 betabellin 15D formed short broad fibrils in the presence of Cu2+ but not in its absence. The observed width of the fibrils (7 ± 2 nm) was consistent with the width (6.8 nm) of a structural model of a fibril that contained two adjacent rows of betabellin 15D -sandwiches joined lengthwise by multiple intersheet hydrogen bonds and widthwise by multiple Cu(II)-imidazole bonds. Electron paramagnetic resonance spectrometry revealed that some pairs of Cu(II) ions in a Cu(II)/betabellin 15D complex were magnetically coupled, which is consistent with the structural model of the Cu(II)/betabellin 15D fibril.     

De novo proteins from designed combinatorial libraries

Combinatorial libraries of de novo amino acid sequences can provide a rich source of diversity for the discovery of novel proteins with interesting and important activities. Randomly generated sequences, however, rarely fold into well-ordered proteinlike structures. To enhance the quality of a library, features of rational design must be used to focus sequence diversity into those regions of sequence space that are most likely to yield folded structures. This review describes how focused libraries can be constructed by designing the binary pattern of polar and nonpolar amino acids to favor proteins that contain abundant secondary structure, while simultaneously burying hydrophobic side chains and exposing hydrophilic side chains to solvent. The "binary code" for protein design was used to construct several libraries of de novo proteins, including both alpha-helical and beta-sheet structures. The recently determined solution structure of a binary patterned four-helix bundle is well ordered, thereby demonstrating that sequences that have neither been selected by evolution (in vivo or in vitro) nor designed by computer can form nativelike proteins. Examples are presented demonstrating how binary patterned libraries have successfully produced well-ordered structures, cofactor binding, catalytic activity, self-assembled monolayers, amyloid-like nanofibrils, and protein-based biomaterials.     

Structure of the 21-30 fragment of amyloid beta-protein

Folding and self-assembly of the 42-residue amyloid beta-protein (Abeta) are linked to Alzheimer's disease (AD). The 21-30 region of Abeta, Abeta(21-30), is resistant to proteolysis and is believed to nucleate the folding of full-length Abeta. The conformational space accessible to the Abeta(21-30) peptide is investigated by using replica exchange molecular dynamics simulations in explicit solvent. Conformations belonging to the global free energy minimum (the "native" state) from simulation are in good agreement with reported NMR structures. These conformations possess a bend motif spanning the central residues V24-K28. This bend is stabilized by a network of hydrogen bonds involving the side chain of residue D23 and the amide hydrogens of adjacent residues G25, S26, N27, and K28, as well as by a salt bridge formed between side chains of K28 and E22. The non-native states of this peptide are compact and retain a native-like bend topology. The persistence of structure in the denatured state may account for the resistance of this peptide to protease degradation and aggregation, even at elevated temperatures.     

Cation-pi interactions in protein-protein interfaces

Arginine is an abundant residue in protein-protein interfaces. The importance of this residue relates to the versatility of its side chain in intermolecular interactions. Different classes of protein-protein interfaces were surveyed for cation-pi interactions. Approximately half of the protein complexes and one-third of the homodimers analyzed were found to contain at least one intermolecular cation-pi pair. Interactions between arginine and tyrosine were found to be the most abundant. The electrostatic interaction energy was calculated to be approximately 3 kcal/mol, on average. A distance-based search of guanidinium:aromatic interactions was also performed using the Macromolecular Structure Database (MSD). This search revealed that half of the guanidinium:aromatic pairs pack in a coplanar manner. Furthermore, it was found that the cationic group of the cation-pi pair is frequently involved in intermolecular hydrogen bonds. In this manner the arginine side chain can participate in multiple interactions, providing a mechanism for inter-protein specificity. Thus, the cation-pi interaction is established as an important contributor to protein-protein interfaces.     

Biophysical characterization of betabellin 16D: a beta-sandwich protein that forms narrow fibrils which associate into broad ribbons.

The betabellin structure is a de novo designed beta-sandwich protein consisting of two 32-residue beta sheets packed against one another by hydrophobic interactions. Betabellin 16S (B16S), a 32-residue peptide chain (HSLTAKIakLTFSIAahTYTCAVakYTAKVSH, where a is DAla, h is DHis, and k is DLys), did not have beta structure in water at pH 6.5. Air oxidation of B16S furnished betabellin 16D (B16D), a 64-residue disulfide-bridged two-chain protein, which also did not fold in water at pH 6.5. However, the extent of beta structure observed for B16D increased with pH and ionic strength of the solution and the B16D concentration as observed by circular dichroism spectropolarimetry. Transmission electron microscopy showed that B16D formed narrow fibrils that associated into broad ribbons in 5.0 mM Mops and 0.25 M NaCl at pH 6.9.     

Exploratory studies of ab initio protein structure prediction: multiple copy simulated annealing, AMBER energy functions, and a generalized born/solvent accessibility solvation model.

A theoretical and computational approach to ab initio structure prediction for polypeptides in water is described and applied to selected amino acid sequences for testing and preliminary validation. The method builds systematically on the extensive efforts applied to parameterization of molecular dynamics (MD) force fields, employs an empirically well-validated continuum dielectric model for solvation, and an eminently parallelizable approach to conformational search. The effective free energy of polypeptide chains is estimated from AMBER united atom potential functions, with internal degrees of freedom for both backbone and amino acid side chains explicitly treated. The hydration free energy of each structure is determined using the Generalized Born/Solvent Accessibility (GBSA) method, modified and reparameterized to include atom types consistent with the AMBER force field. The conformational search procedure employs a multiple copy, Monte Carlo simulated annealing (MCSA) protocol in full torsion angle space, applied iteratively on sets of structures of progressively lower free energy until a prediction of a structure with lowest effective free energy is obtained. Calibration tests for the effective energy function and search algorithm are performed on the alanine dipeptide, selected protein crystal structures, and united atom decoys on barnase, crambin, and six examples from the Rosetta set. Specific demonstration cases of the method are provided for the 8-mer sequence of Ala residues, a 12-residue peptide with longer side chains QLLKKLLQQLKQ, a de novo designed 16 residue peptide of sequence (AAQAA)3Y, a 15-residue sequence with a beta sheet motif, GEWTWDATKTFTVTE, and a 36 residue small protein, Villin headpiece. The Ala 8-mer readily formed an alpha-helix. An alpha-helix structure was predicted for the 16-mer, consistent with observed results from IR and CD spectroscopy and with the pattern in psi/straight phi angles of known protein structures. The predicted structure for the 12-mer, composed of a mix of helix and less regular elements of secondary structure, lies 2.65 A RMS from the observed crystal structure. Structure prediction for the 8-mer beta-motif resulted in form 4.50 A RMS from the crystal geometry. For Villin, the predicted native form is very close to the crystal structure, RMS values of 3.5 A (including sidechains), and 1.01 A (main chain only). The methodology permits a detailed analysis of the molecular forces which dominate various segments of the predicted folding trajectory. Analysis of the results in terms of internal torsional, electrostatic and van der Waals and the electrostatic and non-electrostatic contributions to hydration, including the hydrophobic effect, is presented.     

Molecular dynamics (MD) investigations of preformed structures of the transmembrane domain of the oncogenic Neu receptor dimer in a DMPC bilayer

The critical Val/Glu mutation in the membrane spanning domain of the rat Neu receptor confers the ability for ligand-independent signaling and leads to increased dimerization and transforming ability. There is evidence that the two transmembrane interacting helices play a role in receptor activation by imposing orientation constraints to the intracellular tyrosine kinase domains. By using MD simulations we have attempted to discriminate between correct and improper helix-helix packing by examining the structural and energetic properties of preformed left-handed and right-handed structures in a fully hydrated DMPC bilayer. The best energetic balance between the residues at the helix-helix interface and the residues exposed to the lipids is obtained for helices in symmetrical left-handed interactions packed together via Glu side chain/Ala backbone interhelical hydrogen bonds. Analyses demonstrate the importance of the ATVEG motif in helix-helix packing and point to additional contacting residues necessary for association. Our findings, all consistent with experimental data, suggest that a symmetrical left-handed structure of the helices could be the transmembrane domain configuration that promotes receptor activation and transformation. The present study may provide further insight into signal transduction mechanisms of the ErbB/Neu receptors.     

Design and characterization of peptides with amphiphilic beta-strand structures

To extend our studies on peptides and proteins with amphiphilic secondary structures, a series of peptides designed to form amphiphilic beta-strand structures was designed, synthesized, and characterized by circular dichroism and infrared spectroscopy. Amphiphilic beta-strand conformations may be likely to appear in a variety of surface-active proteins, including apolipoprotein B and fibronectin. In a beta-strand conformation, the synthetic peptides will possess a hydrophobic face composed of valine side chains and a hydrophilic face composed of alternating acidic (glutamic acid) and basic (ornithine or lysine) residues. The peptides studied had a variety of chain lengths (5, 9, and 13 residues), and had the amino groups either free or protected with the trifluoroacetyl group. While the peptides did not possess a high potential for beta-sheet formation based on the Chou Fasman parameters, they possessed significant beta-sheet content, with up to 90% beta-sheet calculated for the 13-residue protected peptide. The driving force for beta-sheet formation is the potential amphiphilicity of this conformation. The beta-strand conformation of the 13-residue deprotected peptide was stable in 50% trifluoroethanol, 6 M guanidine hydrochloride, and octanol. The peptides are strongly self-associating in water, which would reduce the unfavorable contacts of the hydrophobic residues with water. It is clear that small peptides can be designed to form stable beta-strand conformations.