Prediction of the three-dimensional structure of the rap-1A protein from its homology to theras-gene-encoded p21 protein

rap-1A, an anti-oncogene-encoded protein, is aras-p21-like protein whose sequence is over 80% homologous to p21 and which interacts with the same intracellular target proteins and is activated by the same mechanisms as p21, e.g., by binding GTP in place of GDP. Both interact with effector proteins in the same region, involving residues 32–47. However, activated rap-1A blocks the mitogenic signal transducing effects of p21. Optimal sequence alignment of p21 and rap-1A shows two insertions of rap-1A atras positions 120 and 138. We have constructed the three-dimensional structure of rap-1A bound to GTP by using the energy-minimized three-dimensional structure ofras-p21 as the basis for the modeling using a stepwise procedure in which identical and homologous amino acid residues in rap-1A are assumed to adopt the same conformation as the corresponding residues in p21. Side-chain conformations for homologous and nonhomologous residues are generated in conformations that are as close as possible to those of the corresponding side chains in p21. The entire structure has been subjected to a nested series of energy minimizations. The final predicted structure has an overall backbone deviation of 0.7 å from that ofras-p21. The effector binding domains from residues 32–47 are identical in both proteins (except for different side chains of different residues at position 45). A major difference occurs in the insertion region at residue 120. This region is in the middle of another effector loop of the p21 protein involving residues 115–126. Differences in sequence and structure in this region may contribute to the differences in cellular functions of these two proteins.     

Studies on evolutionary and selective properties of hypercycles using a Monte Carlo method

Summary The most relevant properties of hypercycles were previously studied mainly from a theoretical point of view. We have developed a Monte Carlo method simulating hypercyclic organization to obtain information about the dynamics of this prebiotic organization. Nucleation, growth, and selective properties have been tested and the results obtained are in good agreement with those of the theoretical predictions. The influence of hypercyclic organization of the error threshold has also been studied. As a consequence of the emergence of a hypercycle, the value of this threshold decreases. The amount of this decrease depends on the population size. Moreover, for some interval of quality factor values, either the hypercycle organization or an error catastrophe can be produced, depending on the initial conditions. The influence of these phenomena on both the dynamic behavior and evolutionary advantages of the hypercycle, as well as their decisive roles on genome size, are discussed.Presented at the FEBS Symposium on Genome Organization and Evolution, held in Crete, Greece, September 1–5, 1986     

Bias of heritability estimates from twin data in presence of errors of zygosity determination

Simple heritability estimators of continuous as well as discrete traits from twin data are known to overestimate the degree of genetic determination of the measured traits for several reasons. Errors of zygosity determination will, however, underestimate the true heritability. The bias due to wrong assignment of dizygous twin pairs into monozygous type is evaluated here, and the results indicate that this negative bias has a compensatory effect on the estimate of the degree of genetic determination when other factors of similarity between twin pairs are taken into account. It is shown that when an estimate of zygosity error is available, the bias due to this factor can be evaluated quantitatively, and hence the adjustment for zygosity error can be incorporated in the estimation of the degree of genetic determination of a trait. Although this theory is explicitly developed here for twin studies, the general principle also applies for other types of errors of determining the degree of biological relationships for estimation of heritability, in which case this type of error may be more important than the simple zygosity error.     

Generation of potential structures for the G-domain of chloroplast EF-Tu using comparative molecular modeling

Comparative molecular modeling has been used to generate several possible structures for the G-domain of chloroplast elongation factor Tu (EF-Tu(chl)) based on the crystallographic data of the homologous E. coli protein. EF-Tu(chl) contains a 10 amino acid insertion not present in the E. coli protein and this region has been modeled based on its predicted secondary structure. The insertion appears to lie on the surface of the protein. Its orientation could not be determined unequivocally but several likely structures for the nucleotide binding domain of EF-Tu(chl) have been developed. The effects of the presence of water in the Mg2+ coordination sphere and of the protonation state of the GDP ligand on the conformation of the guanine nucleotide binding site have been examined. Relative binding constants of several guanine nucleotide analogs for EF-Tu(chl) have been obtained. The interactions between EF-Tu(chl) and GDP predicted to be important by the models that have been developed are discussed in relation to the nucleotide binding properties of this factor and to the interactions proposed to be important in the binding of guanine nucleotides to related proteins.     

Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison

A method for the prediction of hydrogen positions in proteins is presented. The method is based on the knowledge of the heavy atom positions obtained, for instance, from X-ray crystallography. It employs an energy minimization limited to the environment of the hydrogen atoms bound to a common heavy atom or to a single water molecule. The method is not restricted to proteins and can be applied without modification to nonpolar hydrogens and to nucleic acids. The method has been applied to the neutron diffraction structures of trypsin ribonuclease A, and bovine pancreatic trypsin inhibitor. A comparison of the constructed and the observed hydrogen positions shows few deviations except in situations in which several energetically similar conformations are possible. Analysis of the potential energy of rotation of Lys amino and Ser, Thr, Tyr hydroxyl groups reveals that the conformations of lowest intrinsic torsion energies are statistically favored in both the crystal and the constructed structures.     

Michaelis complexes of porcine pancreatic elastase with 7-[(alkylcarbamoyl)amino]-4-chloro-3-ethoxyisocoumarins: translational sampling of inhibitor position and kinetic measurements.

A step leading to the formation of the covalent complexes between porcine pancreatic elastase (PPE) and 7-[(alkylcarbamoyl)amino]-4-chloro-3-ethoxyisocoumarins (alkylHNCO-EICs) is the formation of the noncovalent Michaelis complex. No average structures are available for the Michaelis complexes of PPE with alkylHNCO-EICs. We present the results of an initial step in obtaining these structures and have determined kinetic constants as well. The kinetic results indicate that formation of the Michaelis complex is what differentiates the effectiveness of these inhibitors in inactivating PPE. The structural and kinetic results together suggest that the structure of the Michaelis complex is necessary for the design of potent alkylHNCO-EIC inhibitors of PPE. Two novel alkylHNCO-EICs are predicted to be the best inhibitors of this series. An alternate mechanism for serine protease inhibition is also proposed. Evidence for, and studies that may add support to, the hypothesized mechanism are discussed.     

The tertiary structure of Aspergillus saitoi minor ribonuclease (Ms) predicted from the structure of RNase T1

Ribonuclease Ms from Aspergillus saitoi is a small acidic protein (11 714 Da) containing 106 amino acids of known sequence. Unlike other enzymes belonging to the RNase T₁ family this ribonuclease is base-unspecific. Using interactive computer graphics and energy minimisation we predicted the structure of RNase Ms on the basis of sequence homology to RNase T₁ of known structure. In this report the predicted structure of this protein is presented and characterised.     

Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry

Globular proteins have high packing densities as a result of residue side chains in the core achieving a tight, complementary packing. The internal packing is considered the main determinant of native protein structure. From that point of view, we present here a method of energy minimization using an automata network to predict a set of amino acid sequences and their side-chain conformations from a desired backbone geometry for de novo design of proteins. Using discrete side-chain conformations, that is, rotamers, the sequence generation problem from a given backbone geometry becomes one of combinatorial problems. We focused on the residues composing the interior core region and predicted a set of amino acid Sequences and their side-chain conformations only from a given backbone geometry. The kinds of residues were restricted to six hydrophobic amino acids (Ala, Ile, Met, Leu, Phe, and Val) because the core regions are almost always composed of hydrophobic residues. The obtained sequences were well packed as was the native sequence. The method can be used for automated sequence generation in the de novo design of proteins. © 1994 Wiley-Liss, Inc.     

Stereochemical studies on cyclic peptides: Detailed energy minimization studies on hydrogen bonded all-trans cyclic pentapeptide backbones

Conformational studies have been carried out on hydrogenbonded all-trans cyclic pentapeptide backbone. Application of a combination of grid search and energy minimization on this system has resulted in obtaining 23 minimum energy conformations, which are characterized by unique patterns of hydrogen bonding comprising of β- and γ-turns. A study of the minimum energy conformationsvis-a-vis non-planar deviation of the peptide units reveals that non-planarity is an inherent feature in many cases. A study on conformational clustering of minimum energy conformations shows that the minimum energy conformations fall into 6 distinct conformational families. Preliminary comparison with available X-ray structures of cyclic pentapeptide indicates that only some of the minimum energy conformations have formed crystal structures. The set of minimum energy conformations worked out in the present study can form a consolidated database of prototypes for hydrogen bonded backbone and be useful for modelling cyclic pentapeptides both synthetic and bioactive in nature. This is part XV of the series. Part XIV in this series is Ramakrishnanet al 1987.     

High-resolution solution structure of siamycin II: Novel amphipathic character of a 21-residue peptide that inhibits HIV fusion

Summary The 21-amino acid peptides siamycin II (BMY-29303) and siamycin I (BMY-29304), derived from Streptomyces strains AA3891 and AA6532, respectively, have been found to inhibit HIV-1 fusion and viral replication in cell culture. The primary sequence of siamycin II is CLGIGSCNDFAGCGYAIVCFW. Siamycin I differs by only one amino acid; it has a valine residue at position 4. In both peptides, disulfide bonds link Cys¹ with Cys¹³ and Cys⁷ with Cys¹⁹, and the side chain of Asp⁹ forms an amide bond with the N-terminus. Siamycin II, when dissolved in a 50:50 mixture of DMSO and H₂O, yields NOESY spectra with exceptional numbers of cross peaks for a peptide of this size. We have used 335 NOE distance constraints and 13 dihedral angle constraints to generate an ensemble of 30 siamycin II structures; these have average backbone atom and all heavy atom rmsd values to the mean coordinates of 0.24 and 0.52 Å, respectively. The peptide displays an unusual wedge-shaped structure, with one face being predominantly hydrophobic and the other being predominantly hydrophilic. Chemical shift and NOE data show that the siamycin I structure is essentially identical to siamycin II. These peptides may act by preventing oligomerization of the HIV transmembrane glycoprotein gp41, or by interfering with interactions between gp41 and the envelope glycoprotein gp120, the cell membrane or membrane-bound proteins [Frèchet, D. et al. (1994) Biochemistry, 33, 42–50]. The amphipathic nature of siamycin II and siamycin I suggests that a polar (or apolar) site on the target protein may be masked by the apolar (or polar) face of the peptide upon peptide/protein complexation.Abbreviations ABNR adopted basis Newton Raphson - AIDS acquired immunodeficiency syndrome - CW continuous wave - DMSO dimethylsulfoxide - DQF-COSY two-dimensional double-quantum-filtered correlation spectroscopy - HIV human immunodeficiency virus - HSQC heteronuclear single-quantum coherence - NOE nuclear Overhauser enhancement - NOESY two-dimensional nuclear Overhauser enhancement spectroscopy - ppm parts per million - P.E.-COSY two-dimensional primitive exclusive correlation spectroscopy - REDAC redundant dihedral angle constraint - rf radio frequency - rmsd root-mean-square difference - SIV simian immunodeficiency virus - sw spectral width - _m mixing time - TOCSY two-dimensional total correlation spectroscopy - TSP trimethylsilyl-2,2,3,3-²H₄-propionate - 2D two-dimensional