PubMed HealthSearch

Biomedical subjects

G M Maggiora

Publications and source records attributed to G M Maggiora.

18 recordsLinked to original sources

Solitary wave dynamics as a mechanism for explaining the internal motion during microtubule growth.

Microtubules, which play many diverse and important roles in biological systems, are usually made up of 13 nearly axial protofilaments formed from individual tubulin molecules. In this paper, a nonlinear dynamic model has been developed to elucidate the mechanism of the internal motion occurring during the assembly of microtubules. The results derived from the model indicate that such internal motion is associated with a solitary wave, or kink, excited by the energy released from the hydrolysis of GTP-->GDP in microtubular solutions. As the kink moves forward, the individual tubulin molecules involved in the kink undergo motions that can be likened to the dislocation of atoms within the crystal lattice. Thus, the dynamic instability of microtubules may be characterized by a series of dislocation motions of the tubulin molecules. An energy estimate shows that a kink in the system possesses about 0.36-0.44 eV, which is quite close to but smaller than the 0.49 eV of energy released from the hydrolysis of GTP. Therefore, the relevant energy derived from our model is fully consistent with experimental observations; this finding also suggests that the hydrolysis energy may be responsible for exciting the solitary wave, or kink, leading to tubulin dislocation in microtubules. Our model, and its intrinsic properties, i.e., dynamic nonlinearity, thermodynamic irreversibility, as well as an energy input from a sustained source, implies that the growth of microtubules is a typical dissipative process and that their structure in vivo is typical of dissipative structures.

Energy Metabolism

A consensus procedure for predicting the location of alpha-helical transmembrane segments in proteins.

To aid in the development of three-dimensional models of membrane-bound proteins, a consensus procedure for predicting alpha-helical transmembrane segments from amino acid sequence is presented. The algorithm combines the results of six individual prediction methods and some basic properties of membrane-spanning helices to obtain a final consensus prediction. Comparison with experiment and several other recently developed methods shows that the consensus procedure performs quite well in comparison to other recent methods. A FORTRAN program has been developed which takes an input file containing an amino acid sequence in one-letter code and outputs a list of the alpha-helical transmembrane segments predicted by the consensus algorithm.

Algorithms

Role of loop-helix interactions in stabilizing four-helix bundle proteins.

One of the critical issues regarding proteins with a four-helix bundle motif is which interactions play the major role in stabilizing this type of folded structure: the interaction among the four alpha-helices or the interaction between the loop and helix segments. To answer this question, an energetic analysis has been carried out for three proteins with a four-helix bundle--namely, methemerythrin, cytochrome b-562, and cytochrome c'. The structures on which the analysis has been made were derived from their respective crystallographic coordinates. All three proteins have long helices (16-26 residues) and most of their loops are short (3-5 residues). However, it was found in all three proteins that loop-helix interactions were stronger than helix-helix interactions. Moreover, not only the nonbonded component but also the electrostatic component of the interaction energy were dominated by loop-helix interactions rather than by interhelix interactions, although the latter involve favorable helix-dipole interactions due to the antiparallel arrangement of neighboring helices. The results of the energetic analysis indicate that the loop segments, whether they are in a theoretical model or in real proteins, play a significant role in stabilizing proteins with four-helix bundles.

Cytochrome b Group

An energy-based approach to packing the 7-helix bundle of bacteriorhodopsin.

Based on the heavy-atom coordinates determined by the electron microscopy for the seven main helical regions of bacteriorhodopsin with the all-trans retinal isomer, energy optimizations were carried out for helix bundles containing the all-trans retinal and 13-cis retinal chromophores, respectively. A combination of simulated annealing and energy minimization was utilized during the process of energy optimization. It was found that the 7-helix bundle containing the all-trans isomer is about 10 kcal/mol lower in conformational energy than that containing the 13-cis isomer. An energetic analysis indicates that such a difference in energy is consistent with the observation that absorption of a 570-nm proton is required for the conversion of a bacteriorhodopsin from its all-trans to 13-cis form. It was also found that the above conversion process is accompanied by a significant conformational perturbation around the chromophore, as reflected by the fact that the beta-ionone ring of retinal moves about 5.6 A along the direction perpendicular to the membrane plane. This is consistent with the observation by Fodor et al. (Fodor, S.P.A., Ames, J.B., Gebhard, R., van der Berg, E.M.M., Stoeckenius, W., Lugtenburg, J., & Mathies, R.A., 1988, Biochemistry 27, 7097-7101). Furthermore, it is interesting to observe that although the retinal chromophore undergoes a significant change in its spatial position, the orientation of its transition dipole changes only slightly, in accord with experimental observations. In other words, even though orientation of the retinal transition dipole is very restricted, there is sufficient room, and degrees of freedom, for the retinal chromophore to readjust its position considerably. This finding provides new insight into the subtle change of the retinal microenvironment, which may be important for revealing the proton-pumping mechanism of bacteriorhodopsin. The importance of electrostatic and nonbonded interactions in stabilizing the 7-helix bundle structure has also been analyzed. Electrostatic interactions favor an antiparallel arrangement among adjacent helices. Nonbonded interactions, however, drive most of the closely packed helices into an arrangement in which the packing angles lie around -160 degrees, a value very near the -154 degrees value computed earlier as the most favorable packing arrangement of two poly(Ala) alpha-helices (Chou, K.-C., Némethy, G., & Scheraga, H.A., 1983, J. Phys. Chem. 87, 2869-2881). The structural features of the 7-helix bundle and their relationship to those found in typical 4-helix bundle proteins are also discussed.(ABSTRACT TRUNCATED AT 400 WORDS)

Amino Acid Sequence

A preliminary 3-D model of the tertiary fold of the polymerase domain of HIV-1 reverse transcriptase.

Development of a 3-D model of the reverse transcriptase from type 1 human immunodeficiency virus (HIV-1 RT), a key enzyme in the pathogenesis of the virus, presents a significant challenge. Three-dimensional structural information is not available for any close homolog, the only 3-D structural data being that of the Klenow fragment (KF) of Escherichia coli DNA polymerase I, for which coordinates of only the alpha-carbons are available. A recently published study of the sequences of a large number of polymerases led to the identification of three common sequence patterns, nominally motif A, motif B and motif C, and to the hypothesis that the various DNA and RNA polymerases including E. coli DNA polymerase I and HIV-1 RT share a common structural motif around their respective polymerase active sites. The preliminary results of recent structural studies on two other polymerases also support this hypothesis. Based on the assumption of structural homology in the active site regions of their polymerase domains, the HIV-1 RT and KF sequences were aligned using pattern-based secondary structure predictions as a guide and motifs A, B and C as 'anchor points'. However, as suggested by the results of chemical modification experiments, it was assumed that the order of the motifs in KF, viz. A, B and C, differed from that of the related motifs A, C and B' in HIV-1 RT, a rearrangement that could have been brought about by an exon shuffling type of mechanism.(ABSTRACT TRUNCATED AT 250 WORDS)

Amino Acid Sequence

A heuristic approach to predicting the tertiary structure of bovine somatotropin.

A combination of a heuristic approach and energy minimization was used to predict the three-dimensional structure of bovine somatotropin (bSt), also known as bovine growth hormone, a protein of 191 amino acids. The starting points for energy minimizations were generated from the following two types of inputs: (a) the amino acid sequence and (b) the heuristic inputs, which were derived according to physical, chemical, and biological principles by piecing together all useful information available. The predicted 3-D structure of the bSt molecule has all the features observed in four-helix bundle proteins. The four alpha-helices in bSt are intimately packed to form an assembly with an approximately square cross section. All the adjacent alpha-helices are antiparallel, with a somewhat tilted angle between each of the adjacent pairs so that the assembly of the four helices looks like a left-handed twisted bundle. There are two disulfide bonds in the bSt structure: one "hooking" the middle of a long loop with helix 4 so as to pull the long loop onto the surface of the helix bundle and the other "hooking" the C-terminal segment with the same helix so as to force the C-terminal segment to bend toward the helix bundle. As a consequence, a considerable part of the surface of the four-helix bundle is closely packed or intimately embraced by the loop segments. The predicted bSt structure has a hydrophobic core and a hydrophilic exterior surface. The energetic analysis of the predicted bSt structure indicates that the interaction between helices and loops plays a dominant role in stabilizing the four-helix bundle structure from the viewpoint of both electrostatic and nonbonded interactions. A technique called FOLD was meanwhile developed, by which one can fold a polypeptide chain into any shape as desired. This tool proved to be very useful during the heuristic model-building process.

Algorithms

Mass-weighted molecular dynamics simulation of cyclic polypeptides.

A modified molecular dynamics (MD) method in which atomic masses are weighted was developed previously for studying the conformational flexibility of neuroregulating tetrapeptide Phe-Met-Arg-Phe-amide (FMRF-amide). The method has now been applied to longer and constrained molecules, namely a disulfide-linked cyclic hexapeptide, c[CYFQNC], and its linear and "pseudo-cyclic" analogues. The sampling of dehedral conformational space of teh linear hexapeptide in mass-weighted MD simulations was found to be improved significantly over conventional MD simulations, as in the case of the shorter FMRF-amide molecule studied previously. In the cyclic hexapeptide, the internal constraint of the molecule due to the intramolecular disulfide bond (hence the absence of free terminals in the molecule) does not adversely affect the significant improvement of conformational sampling in mass-weighted MD simulations over normal MD simulations. The pseudo-cyclic polypeptide is identical to the linear CYFQNC molecule in amino acid sequence (i.e., side chains of the cysteine residues are reduced), but the positions of its two terminal heavy atoms were held fixed in space such that the molecule has a nearly cyclic conformation. For this molecule, the mass-weighted MD simulation generated a wide range of polypeptide backbone conformations covering the internal dihedral degrees of freedom; moreover, the physical space of the pseudo-cyclic structure was also sampled in a complete revolution of the entire molecular fragment about the two fixed termini during the simulation. These characteristics suggest that mass-weighted MD can also be an extremely useful method for conformational analyses of constrained molecules and, in particular, for modeling loops on protein surfaces.

Amino Acid Sequence

Conformational and geometrical properties of idealized beta-barrels in proteins.

An equation for calculating the distances between the atoms involved in forming an idealized hydrogen bond in a parallel or antiparallel beta-barrel has been derived by adjusting the corresponding data given by Pauling and Corey for a beta-sheet. Based on these distances, a geometrical optimization method was developed, by which one can generate various idealized beta-barrels: parallel or antiparallel, tilted or non-tilted, right-tilted or left-tilted. For each type of idealized beta-barrel thus obtained, the corresponding conformation and characteristic geometric parameters as well as their relationship are analyzed and discussed. Since the strand in a tilted beta-barrel traces a curve rather than a straight line on a cylinder-like surface, a regular chain in which the dihedral angles of each residue are the same cannot form a tilted beta-barrel but only a non-tilted beta-barrel. As observed, the strands of a right-tilted beta-barrel possess a very strong right-handed twist. The radii of the idealized tilted parallel and antiparallel beta-barrels are greater than those of the corresponding non-tilted ones by approximately 1 A and approximately 1.5 A, respectively. Consequently, there is relatively more room for a tilted beta-barrel to accommodate the internal side-chains, suggesting that a conformational change from a non-tilted beta-barrel to a tilted one would ease the repulsion among the crowded internal side-chains so as to make the structure more stable. The values of root-mean-square fits indicate that the idealized right-tilted beta-barrels coincide quite well with the observed beta-barrels in both parallel and antiparallel cases.

Chemical Phenomena

Topological analysis of hydrogen bonding in protein structure.

A recent study has shown that topological stereoisomers exist for the polypeptide chain in disulfide-containing proteins that are represented by non-planar graphs. This topological stereochemistry is demonstrated in the structure of variant 3 toxin in the venom of the North American scorpion Centruroides sculpturatus Ewing and the structure of toxin II from the North African scorpion Androctonus australis Hector. In this report, we found that a similar topological analysis can be applied to the hydrogen bonding in alpha-helices and beta-sheets within protein molecules, and we described the topological characteristics of chiral properties of protein secondary structure elements. Specifically, a closed right-handed alpha-helix of more than six residues long is shown formally to be non-planar and has the L topology. Antiparallel beta-sheets are planar. Two parallel beta-strands each of at least three residues in length, however, constitute a non-planar structural element and can have either L or D topology. The favored right-handed crossover for parallel beta-sheets has the L form, the same as the right-handed alpha-helix. This topological description of the hydrogen bonding in secondary structures may be extended to higher levels of protein structure and may provide a conceptual framework for studying complex protein architecture in general.

Chemical Phenomena

Quasi-continuum models of twist-like and accordion-like low-frequency motions in DNA.

Formulae for calculating low-frequency twist-like and accordion-like modes of DNA molecules have been derived using a quasi-continuum model. The formulae can be employed in essentially all (viz. A, B, C, D, E, and Z) forms of DNA. Calculated results indicate that the experimentally observed low-frequency modes at 22 cm-1 for the A-form octanucleotide (d[CCCCGGGG]) and at 18 cm-1 for the B-form dodecanucleotide (d[CGCAA ATTTGCG]) may result from accordion-like motions, while those observed at 12 cm-1 and 15 cm-1 may result from combinations of twist-like oscillations excited in the intact segments of B- and A-DNA's, respectively. Frequency shifts in the low-frequency modes observed when DNA molecules undergo conformational changes among different forms are also discussed in terms of the current model.

DNA

Energetics of the structure of the four-alpha-helix bundle in proteins.

The main features of the four-alpha-helix bundle, one of the characteristic structural elements of many proteins, can be explained in terms of noncovalent interactions between the constituent helices. Conformational energy computations have been carried out on four types of four-alpha-helix bundles, each consisting of four CH3CO-(L-Ala)10-NHCH3 polypeptide chains, with various combinations of parallel and antiparallel orientations of the helices. In the bundle with the most favorable energy, all pairs of neighboring helices are oriented antiparallel--i.e., in the orientation that is favored by electrostatic interactions between the helices. In this structure, the orientation angle between neighboring helix axes is -168 degrees, within +/- 7 degrees, in close agreement with the orientation angles observed in proteins and with the value that we computed earlier for the most favorable packing of pairs of interacting alpha-helices. This orientation corresponds to a left-handed twisting of the helical bundle. The preferred handedness of this twisting arises as a result of favorable nonbonded interactions between the alpha-helices.

Macromolecular Substances