Conformational and helicoidal analysis of the molecular dynamics of proteins: "curves," dials and windows for a 50 psec dynamic trajectory of BPTI

A new procedure for the graphic analysis of molecular dynamics (MD) simulations on proteins is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the conformational and helicoidal parameters for each structure entering the analysis via the method "Curves," developed for proteins by Sklenar, Etchebest, and Lavery (Proteins: Structure, Function Genet. 6:46-60, 1989) followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels ("dials") for each torsional parameter and "windows" on the range values assumed by the linear and angular helicoidal parameters, and is present in a form isomorphous with the primary structure per se. The complete time evolution of dynamic structure can then be depicted in a set of four composite figures. Dynamic aspects of secondary and tertiary structure are also provided. The procedure is illustrated with an analysis of a 50 psec in vacuo simulation on the 58 residue protein, bovine pancreatic trypsin inhibitor (BPTI), in the vicinity of the local minimum on the energy surface corresponding to a high resolution crystal structure. The time evolution of 272 conformational and 788 helicoidal parameters for BPTI is analyzed. A number of interesting features can be discerned in the analysis, including the dynamic range of conformational and helicoidal motions, the dynamic extent of 2 degrees structure motifs, and the calculated fluctuations in the helix axis. This approach is expected to be useful for a critical analysis of the effects of various assumptions about force field parameters, truncation of potentials, solvation, and electrostatic effects, and can thus contribute to the development of more reliable simulation protocols for proteins. Extensions of the analysis to present differential changes in conformational and helicoidal parameters is expected to be valuable in MD studies of protein complexes with substrates, inhibitors, and effectors and in determining the nature of structural changes in protein-protein interactions.     

DNA sequence analyses of the ribosomal spacer regions in theTriticeae

Two regions of the ribosomal DNA (rDNA) were sequenced from a range of species from the tribeTriticeae. One region, the central spacer, was found to be more divergent in sequence than the other, the 18 S-spacer junction. Both regions contained sequences 20–30 bp long which were more highly conserved than the remainder of the region and their possible significance in rDNA expression is discussed. Phenetic relationships based on the sequence data were generally consistent with the relationships based on other criteria. Species possessing the S, E, J₁J₂, D, and B genomes clustered together, with the H genome species being the most distinct of those examined. The R, P, and V genome species occupy an intermediate position in the overall pattern of relationships. Some relationships differed in detail from those established by other parameters, for example the position of the N genome species, and explanations for discrepancies of this type are discussed.     

Structural similarities among the high molecular weight protein fractions of oilseeds

Data on the physico-chemical properties of proteins from soybean, groundnut, sesame seed, sunflower seed, safflower seed, mustard seed, rapeseed and cotton seed are fairly extensive. An examination of the available data on high molecular weight proteins suggests that there are similarities in many of their properties. In this report the similarity in amino acid composition, size and shape, molecular weight, secondary structure, subunit composition, association-dissociation at high and low pH, stability towards denaturants, hydrolysis by enzymes and quaternary structure of the high molecular weight proteins is discussed. Based on these similarities a model has been proposed for the associationdissociation, denaturation and reassociation behaviour of the high molecular weight proteins of oilseeds.     

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.     

转铁蛋白结构与功能的研究——骆驼血清转铁蛋自的分离纯化及其含单一铁结合部位的结构域的制备和鉴定

分离纯化获得的骆驼血清转铁蛋白由分子量为73,000和63,000两个组分组成。两者至少N-端五肽顺序相同(Met-Pro-Asp-Lys-Thr)。骆驼血清转铁蛋白在生理pH下不能与人胎盘转铁蛋白受体结合。用胰蛋白酶酶解骆驼转铁蛋白可以同时得到两个合单一铁结合部位的结构域,分别来自转铁蛋白分子的N-端称N-端结构域(分子量34,700和40,700)和C-端称C-端结构域(分子量35,100)。在上述结果的基础上指出并讨论了反刍动物转铁蛋白在结构和功能上存在更多的共同性,而与其它哺乳动物的转铁蛋白有着明显的区别。     

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K+-ATPase

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center ofRhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na⁺, K⁺-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na⁺, K⁺-ATPase is discussed in the final part of the work.     

Preferential inhibition of acetylcholinesterase molecular forms in rat brain

The effect of eight different acetylcholinesterase inhibitors (AChEIs) on the activity of acetylcholinesterase (AChE) molecular forms was investigated. Aqueous-soluble and detergent-soluble AChE molecular forms were separated from rat brain homogenate by sucrose density sedimentation. The bulk of soluble AChE corresponds to globular tetrameric (G₄), and monomeric (G₁) forms. Heptylphysostigmine (HEP) and diisopropylfluorophosphate were more selective for the G₁ than for the G₄ form in aqueous-soluble extract. Neostigmine showed slightly more selectivity for the G₁ form both in aqueous- and detergent-soluble extracts. Other drugs such as physostigmine, echothiophate, BW284C51, tetrahydroaminoacridine, and metrifonate inhibited both aqueous- and detergent-soluble AChE molecular forms with similar potency. Inhibition of aqueous-soluble AChE by HEP was highly competitive with Triton X-100 in a gradient, indicating that HEP may bind to a detergent-sensitive non-catalytic site of AChE. These results suggest a differential sensitivity among AChE molecular forms to inhibition by drugs through an allosteric mechanism. The application of these properties in developing AChEIs for treatment of Alzheimer disease is considered.Special issue dedicated to Dr. Morris H. Aprison.     

Molecular mimicry and the autoimmune response to the peripheral nerve myelin PO glycoprotein

In the Lewis rat immunisation with the myelin PO glycoprotein can induce an inflammatory demyelinating disease of the peripheral nervous system, experimental allergic neuritis (EAN), which has many clinical and histopathological parallels with the human disease the Guillain-Barre syndrome. In view of the reported association of GBS with a number of infectious agents we have investigated whether molecular mimicry may occur between microbial antigens and the PO protein that could possibly trigger a similar pathogenic autoimmune response in man. A computer search of the available protein sequence data bases identified several absolute sequence homologies between PO and viral proteins that involve five or more consecutive amino acid residues. Four of these sequence homologies involved viral pathogens previously associated with the Guillain-Barre syndrome, namely Epstein-Barr virus (EBV), cytomegalovirus (CMV), Varicella zoster virus (VZV) and human immunodeficiency virus I (HIV I). Although, sequence homologies were also found between viral peptides and the neuritogenic determinants of PO, residues 56–71 and 180–199, these homologies proved incapable of eliciting EAN in the Lewis rat. These observations are discussed with reference to the role that molecular mimicry between T cell epitopes on pathogen derived antigens and the PO protein may play in the pathogenesis of the Guillain-Barre syndrome.Abbreviations EAN Experimental allergic neuritis - EAE experimental allergic encephalomyelitis - PNS peripheral nervous system - CNS central nervous system - MBP myelin basic protein - GBS Guillain Barre syndrome - CFA complete Freund's Adjuvant - LPC lysophosphatidyl choline - VZV Varicella zoster virus - CMV cytomegalovirus - EBV Epstein Barr virus - HIV I human immunodeficiency virus I Special issue dedicated to Dr. Alan N. Davison     

Molecular modeling of the amphipathic helices of the plasma apolipoproteins.

In this paper we propose a classification of the amphipathic helical repeats occurring in the plasma apolipoprotein sequences. It is based upon the calculation of the molecular hydrophobicity potential around the helical segments. The repeats were identified using a new autocorrelation matrix, based upon similarities of hydrophobic and hydrophilic properties of the amino acid residues within the apolipoprotein sequences. The helices were constructed by molecular modeling, the molecular hydrophobicity potential was calculated, and isopotential contour lines drawn around the helices yielded a three-dimensional visualization of the hydrophobicity potential. Two classes of apolipoproteins could be differentiated by comparing the hydrophobic angles obtained by projection of the isopotential contour lines on a plane perpendicular to the long axis of the helix. The isopotential contour lines around apo AI, AIV, and E are more hydrophilic than hydrophobic, whereas they are of similar intensity for apo AII, CI, and CIII. In both cases discoidal lipid-protein complexes are generated, with the amphipathic helices around the edge of the lipid core. The long axis of the helices is oriented parallel to the phospholipid acyl chains and the hydrophilic side of the helix toward the aqueous phase. As a result of the differences in hydrophobicity potential, the contact between the hydrophobic side of the helices and the phospholipid acyl chains is larger for apo AII, CI, and CIII than for the other apolipoproteins. This might account for the greater stability of the discoidal complexes generated between phospholipids and these apoproteins.     

Crystal structure of human immunoglobulin fragment Fab New refined at 2.0 A resolution.

The three-dimensional structure of the human immunoglobulin fragment Fab New (IgG1, lambda) has been refined to a crystallographic R-factor of 16.9% to 2 A resolution. Rms deviations of the final model from ideal geometry are 0.014 A for bond distances and 3.03 degrees for bond angles. Refinement was based on a new X-ray data set including 28,301 reflections with F > 2.5 sigma(F) from 6.0 to 2.0 A resolution. The starting model for the refinement procedure reported here is from the Brookhaven Protein Data Bank entry 3FAB (rev. 1981). Differences between the initial and final models include modified polypeptide-chain folding in the third complementarity-determining region (CDR3) and the third framework region (FR3) of VH and in some exposed loops of CL and CH1. Amino acid sequence changes were determined at a number of positions by inspection of difference electron density maps. The incorporation of amino acid sequence changes results in an improved VH framework model for the "humanization" of monoclonal antibodies.