Stereochemical analysis of the antigenic tip of the V3 loop peptide of HIV-1 gp120

A novel multiple turn conformation has been observed for a segment GPGRAFY in the crystal structure of a complex of HIV-1 gp120 V3 loop peptide with the Fab fragment of a neutralizing antibody [Ghiara ct al. (1994) Science 264, 82-85]. A structural motif has been defined for the peptide segment, employing idealized backbone conformations characterized by ranges of virtual C-alpha torsion angles and bond angles. A search of 122 high-resolution protein crystal structures has permitted identification of 24 examples of similar structural motifs. Two major conformational families have been identified, which differ primarily in the conformation at residue 3. The observed conformation at residue 3 in family 1 is left-handed helical (alpha(L)) and that in family 2 is right-handed helical (alpha(R)). Of the 10 examples in family 1, 9 examples have Gly residues at position 3. Of the 12 examples in family 2, 7 examples have Asn/Asp at position 3. Computer modeling of the V3 loop tip sequence using the two backbone conformational families as starting points leads to minimum-energy conformations in which antigenically important side-chains occupy similar spatial arrangements. This stereochemical analysis of the V3 loop tip sequence suggests a rational basis for the design of synthetic analog peptides for use as viral antagonists or synthetic antigens. (C) Munksgaard 1995.

Stereochemical Analysis of Higher α,α-Dialkylglycine Containing Peptides. Characterization of Local Helical Conformations at Dipropylglycine Residues and Observation of a Novel Hydrated Multiple β-Turn Structure in Crystals of a Glycine Rich Peptide

The peptide Boc-Gly-Dpg-Gly-Gly-Dpg-Gly-NHMe (1) has been synthesized to examine the conformational preferences of Dpg residues in the context of a poor helix promoting sequence. Single crystals of 1 were obtained in the space group P21/c with a = 13.716(2) Å, b = 12.960(2) Å, c = 22.266(4) Å, and β = 98.05(1)°; R = 6.3% for 3660 data with |Fo| > 4σ. The molecular conformation in crystals revealed that the Gly(1)-Dpg(2) segment adopts φ, ψ values distorted from those expected for an ideal type II‘ β-turn (φGly(1) = +72.0°, ψGly(1) = −166.0°; φDpg(2) = −54.0°, ψDpg(2) = −46.0°) with an inserted water molecule between Boc-CO and Gly(3)NH. The Gly(3)-Gly(4) segment adopts φ, ψ values which lie broadly in the right handed helical region (φGly(3) = −78.0°, ψGly(3) = −9.0°; φGly(4) = −80.0°, ψGly(4) = −18.0°). There is a chiral reversal at Dpg(5) which takes up φ, ψ values in the left handed helical region. The Dpg(5)-Gly(6) segment closely resembles an ideal type I‘ β-turn (φDpg(5) = +56.0°, ψDpg(5) = +32.0°; φGly(6) = +85.0°, ψGly(6) = −3.0°). Molecules of both chiral senses are found in the centrosymmetric crystal. The C-terminus forms a hydrated Schellman motif, with water insertion into the potential 6 → 1 hydrogen bond between Gly(1)CO and Gly(6)NH. NMR studies in CDCl3 suggest substantial retention of the multiple turn conformation observed in crystals. In solution the observed NOEs support local helical conformation at the two Dpg residues.

Thermodynamic analysis of three state denaturation of Peanut Agglutinin

Peanut Agglutinin (PNA) is a homotetrameric protein with a very unusual open quaternary structure. During denaturation, it first dissociates into a molten globule like state, which subsequently undergoes complete denaturation. Urea denaturation of PNA at neutral pH has been studied by intrinsic fluorescence spectroscopy and has been fitted to a three state model, A(4) double left right arrow 4I double left right arrow 4U, to get all the relevant thermodynamic parameters. Urea denaturation leads to continuous red shift of wavelength maxima. The molten globule like state is formed in a short range of urea concentration. Refolding of the denatured PNA has been attempted by intrinsic fluorescence study. Refolding by instantaneous dilution shows the occurrence of the formation of an intermediate at a relatively rapid rate, within few seconds. The transition from PNA tetramer to molten globule like state is found to have a Delta G value of similar to 33 kcal/mole while it is similar to 8 kcal/mole for the transition from molten globule like state to a completely denatured state. This in turn indicates that the tetramerization in PNA contributes significantly to the stability of the oligomer.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Molecular analysis of Friedreich's ataxia locus in the Indian population

Friedreich's ataxia (FRDA) is an autosomal recessive neurodegenerative disorder caused by expansion of GAA repeats in the frataxin gene. We have carried out the first molecular analysis at the Friedreich's ataxia locus in the Indian population. Materials and methods - Three families clinically diagnosed for Friedreich's ataxia were analyzed for GAA expansion at the FRDA locus. The distribution of GAA repeats was also estimated in normal individuals of Indian origin. Results - All patients clinically diagnosed for Friedreich's ataxia were found to be homozygous for GAA repeat expansion. The GAA repeat in the normal population show a bimodal distribution with 94% of alleles ranging from 7-16 repeats. Conclusion - Indian patients with expansion at the FRDA locus showed typical clinical features of Friedreich's ataxia. The low frequency of large normal alleles (6%) could indicate that the prevalence of this disease in the Indian population is likely to be low.

Topological analysis off charge density distribution in concomitant polymorphs of 3-acetylcoumarin, a case off packing polymorphism

Detailed investigation of the charge density distribution in concomitant polymorphs of 3-acetylcoumarin in terms of experimental and theoretical densities shows significant differences in the intermolecular features when analyzed based on the topological properties via the quantum theory of atoms in molecules. The two forms, triclinic and monoclinic (Form A and Form B), pack in the crystal lattice via weak C-H---O and C-H---pi interactions. Form A results in a head-to-head molecular stack, while Form B generates a head-to-tail stack. Form A crystallizes in PI (Z' = 2) and Form B crystallizes in P2(1)/n (Z = 1). The electron density maps of the polymorphs demonstrate the differences in the nature of the charge density distribution in general. The charges derived from experimental and theoretical analysis show significant differences with respect to the polymorphic forms. The molecular dipole moments differ significantly for the two forms. The lattice energies evaluated at the HF and DFT (B3LYP) methods with 6-31G** basis set for the two forms clearly suggest that Form A is the thermodynamically stable form as compared to Form B. Mapping of electrostatic potential over the molecular surface shows dominant variations in the electronegative region, which bring out the differences between the two forms.

Conformational analysis of the right-hand twisted antiparallel β-Structure

The conformational analysis of a pair of two-linked peptide units in the anti-parallel arrangement is reported here with a view to study the effect of association of one chain with the other. The pair of two-linked peptide units were fixed in space through the hydrogen bonds between them, in accordance with certain hydrogen bond criteria. Model building was undertaken to ascertain whether the proximity of the side-chains could be used to eliminate any one of the right-hand twisted, left-hand twisted or regular β-structures. Stereochemically, it was found possible with all of them. The preference for a right-hand twisted β-structure, however, was indicated by the classical energy calculations. The relevance of the results thus obtained is discussed in the context of the preferential right-hand twist of the β-pleated sheets present in globular proteins. The agreement between the minimum energy conformations obtained for the pair of two-linked peptide units and the globular protein data is also indicated.

Vibrational analysis of boldine hydrochloride using QM/MM approach

A complete vibrational analysis was performed on the molecular structure of boldine hydrochloride using QM/MM method. The equilibrium geometry, harmonic vibrational frequencies and infrared intensities were calculated by QM/MM method with B3LYP/6-31G(d) and universal force field (UFF) combination using ONIOM code. We found the geometry obtained by the QM/MM method to be very accurate, and we can use this rapid method in place of time consuming ab initio methods for large molecules. A detailed interpretation of the infrared spectra of boldine hydrochloride is reported. The scaled theoretical wave numbers are in perfect agreement with the experimental values. The FT-IR spectra of boldine hydrochloride in the region 4000-500 cm(-1) were recorded in CsI (solid phase) and in chloroform with concentration 5 and 10 mg/ml.

Conformational analysis and design of cross-strand disulfides in antiparallel beta-sheets

Cross-strand disulfides bridge two cysteines in a registered pair of antiparallel beta-strands. A nonredundant data set comprising 5025 polypeptides containing 2311 disulfides was used to study cross-strand disulfides. Seventy-six cross-strand disulfides were found of which 75 and 1 occurred at non-hydrogen-bonded (NHB) and hydrogen-bonded (HB) registered pairs, respectively. Conformational analysis and modeling studies demonstrated that disulfide formation at HB pairs necessarily requires an extremely rare and positive chi(1) value for at least one of the cysteine residues. Disulfides at HB positions also have more unfavorable steric repulsion with the main chain. Thirteen pairs of disulfides were introduced in NHB and HB pairs in four model proteins: leucine binding protein (LBP), leucine, isoleucine, valine binding protein (LIVBP), maltose binding protein (MBP), and Top7. All mutants LIVBP T247C V331C showed disulfide formation either on purification, or on treatment with oxidants. Protein stability in both oxidized and reduced states of all mutants was measured. Relative to wild type, LBP and MBP mutants were destabilized with respect to chemical denaturation, although the sole exposed NHB LBP mutant showed an increase of 3.1 degrees C in T-m. All Top7 mutants were characterized for stability through guanidinium thiocyanate chemical denaturation. Both exposed and two of the three buried NHB mutants were appreciably stabilized. All four HB Top7 mutants were destabilized (Delta Delta G(0) = -3.3 to -6.7 kcal/mol). The data demonstrate that introduction of cross-strand disulfides at exposed NHB pairs is a robust method of improving protein stability. All four exposed Top7 disulfide mutants showed mild redox activity. Proteins 2011; 79: 244-260. (C) 2010 Wiley-Liss, Inc.

Genome analysis: A new approach for visualization of sequence organization in genomes

In this article we describe and demonstrate the versatility of a computer program, GENOME MAPPING, that uses interactive graphics and runs on an IRIS workstation. The program helps to visualize as well as analyse global and local patterns of genomic DNA sequences. It was developed keeping in mind the requirements of the human genome sequencing programme, which requires rapid analysis of the data. Using GENOME MAPPING one can discern signature patterns of different kinds of sequences and analyse such patterns for repetitive as well as rare sequence strings. Further, one can visualize the extent of global homology between different genomic sequences. An application of our method to the published yeast mitochondrial genome data shows similar sequence organizations in the entire sequence and in smaller subsequences.

Theoretical analysis of x-ray-absorption near-edge fine structure at the O and metal K edges of LaFeO3 and LaCoO3

We present experimental x-ray-absorption spectra at the oxygen and 3d transition-metal K edges of LaFeO3 and LaCoO3. We interpret the experimental results in terms of detailed theoretical calculations based on multiple-scattering theory. Along with providing an understanding of the origin of various experimental features, we investigate the effects of structural distortions and the core-hole potential in determining the experimental spectral shape. The results indicate that the core-hole potential as well as many-body effects within the valence electrons do not have any strong effect on the spectra suggesting that the spectral features can be directly interpreted in terms of the electronic structure of such compounds.

An analysis of substitution, deletion and insertion mutations in cancer genes

Cancer-associated mutations in cancer genes constitute a diverse set of mutations associated with the disease. To gain insight into features of the set, substitution, deletion and insertion mutations were analysed at the nucleotide level, from the COSMIC database. The most frequent substitutions were c -> t, g -> a, g -> t, and the most frequent codon changes were to termination codons. Deletions more than insertions, FS (frameshift) indels more than I-F (in-frame) ones, and single-nucleotide indels, were frequent. FS indels cause loss of significant fractions of proteins. The 5'-cut in FS deletions, and 5'-ligation in FS insertions, often occur between pairs of identical bases. Interestingly, the cut-site and 3'-ligation in insertions, and 3'-cut and join-pair in deletions, were each found to be the same significantly often (p < 0.001). It is suggested that these features aid the incorporation of indel mutations. Tumor suppressors undergo larger numbers of mutations, especially disruptive ones, over the entire protein length, to inactivate two alleles. Proto-oncogenes undergo fewer, less-disruptive mutations, in selected protein regions, to activate a single allele. Finally, catalogues, in ranked order, of genes mutated in each cancer, and cancers in which each gene is mutated, were created. The study highlights the nucleotide level preferences and disruptive nature of cancer mutations.

Analysis of designed beta-hairpin peptides: molecular conformation and packing in crystals

The crystal structures of several designed peptide hairpins have been determined in order to establish features of molecular conformations and modes of aggregation in the crystals. Hairpin formation has been induced using a centrally positioned (D)Pro-Xxx segment (Xxx = (L)Pro, Aib, Ac(6)c, Ala; Aib = alpha-aminoisobutyric acid; Ac(6)c = 1-aminocyclohexane-1-carboxylic acid). Structures of the peptides Boc-Leu-Phe-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (1), Boc-Leu-Tyr-Val-(D)Pro-(L)Pro-Leu-Phe-Val-OMe (2, polymorphic forms labeled as 2a and 2b), Boc-Leu-Val-Val-(D)Pro-(L)Pro-Leu-Val-Val-OMe (3), Boc-Leu-Phe-Val-(D)Pro-Aib-Leu-Phe-Val-OMe (4, polymorphic forms labeled as 4a and 4b), Boc-Leu-Phe-Val-(D)Pro-Ac(6)c-Leu-Phe-Val-OMe (5) and Boc-Leu-Phe-Val-(D)Pro-Ala-Leu-Phe-Val-OMe (6) are described. All the octapeptides adopt type II' beta-turn nucleated hairpins, stabilized by three or four cross-strand intramolecular hydrogen bonds. The angle of twist between the two antiparallel strands lies in the range of -9.8 degrees to -26.7 degrees. A detailed analysis of packing motifs in peptide hairpin crystals is presented, revealing three broad modes of association: parallel packing, antiparallel packing and orthogonal packing. An attempt to correlate aggregation modes in solution with observed packing motifs in crystals has been made by indexing of crystal faces in the case of three of the peptide hairpins. The observed modes of hairpin aggregation may be of relevance in modeling multiple modes of association, which may provide insights into the structure of insoluble polypeptide aggregates.

Topological analysis of electron density and the electrostatic properties of isoniazid: an experimental and theoretical study

Isoniazid (isonicotinohydrazide) is an important first-line antitubercular drug that targets the InhA enzyme which synthesizes the critical component of the mycobacterial cell wall. An experimental charge-density analysis of isoniazid has been performed to understand its structural and electronic properties in the solid state. A high-resolution single-crystal X-ray intensity data has been collected at 90 K. An aspherical multipole refinement was carried out to explore the topological and electrostatic properties of the isoniazid molecule. The experimental results were compared with the theoretical charge-density calculations performed using CRYSTAL09 with the B3LYP/6-31G** method. A topological analysis of the electron density reveals that the Laplacian of electron density of the N-N bond is significantly less negative, which indicates that the charges at the b.c.p. (bond-critical point) of the bond are least accumulated, and so the bond is considered to be weak. As expected, a strong negative electrostatic potential region is present in the vicinity of the O1, N1 and N3 atoms, which are the reactive locations of the molecule. The C-H center dot center dot center dot N, C-H center dot center dot center dot O and N-H center dot center dot center dot N types of intermolecular hydrogen-bonding interactions stabilize the crystal structure. The topological analysis of the electron density on hydrogen bonding shows the strength of intermolecular interactions.