Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety-six of these alignment positions correspond to conserved Gly residue with (phi, ) values allowed for non-glycyl residues. Reasons for this observation were investigated by in-silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (phi, ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (phi, ) values which are disallowed for Ala. In-silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.

NON-GLYCINE RESIDUES IN THE DISALLOWED REGIONS OF THE RAMACHANDRAN MAP AND THEIR NEIGHBORHOOD PREFERENCES

The allowed and the ``disallowed'' regions in the celebrated Ramachandran map (phi-psi] map) was elegantly deduced by Ramachandran, Ramakrishnan and Sasisekharan even before the protein crystal structures became available. This powerful map was derived based on rigid geometry of the peptide group and later several investigations on protein crystal structures reported the occurrence of a small fraction of the phi-psi] torsion angles in the disallowed region. The question is what factors make these residues adopt disallowed conformations? Is it driven by the necessity to maintain the overall topology or is it associated with function or is it just that the disallowed conformations are extreme limits of the allowed conformations? Today, with the availability of a large number of high resolution crystal structures, we have revisited this problem. Apart from validating some of the earlier findings such as residue propensities, preferred location in the secondary structure, we have explored their spatial neighborhood preferences using the protein structure network PSN] approach developed in our lab. Finally, the structural and functional implications of the disallowed conformations are examined.

PROPENSITIES OF AMINO ACID RESIDUES IN PROTEINS FOR DIFFERENT REGIONS OF THE RAMACHANDRAN MAP

Using a dataset of 1164 crystal structures of largely non-homologous proteins defined at a resolution of 1.5 angstrom or better, we have investigated the (phi,psi) preferences of 20 residue types by considering the residues which occur in loops. Propensities of residue types to occur in the loops with (phi,psi) values in the aa region of the Ramachandran map has a poor correlation coefficient of 0.48 to the Chou-Fasman propensities of the residue types to occur in the a-helical segments. However the correlation coefficient between propensities of residues in loops to adopt beta conformations and those in beta-sheet is much higher (0.95). These observations suggest that a-helix formation is well influenced by the local amino acid sequence while intrinsic preference of residue types for beta-sheet plays a major role in the formation of beta-sheet. The main chain polar groups of residues in loops, that can affect the (phi,psi) values, can be involved in intra-molecular hydrogen bonding. Therefore we investigated further by considering subset of residues in loops with low (0 to 2) number of intra-molecular hydrogen bonds per residue involving main chain polar atoms. For this subset, the correlation coefficients between propensities for alpha-helix and alpha(R) region and between beta-sheet and beta-region are 0.26 and 0.64 respectively. This reiterates higher intrinsic tendency of beta-region favouring residues to adopt beta-sheet than alpha(R) region favouring residues to adopt alpha-helical structure.

INDUCTION OF FOLDED STRUCTURES IN DESIGNED PEPTIDES USING CONFORMATIONALLY CONSTRAINED RESIDUES

Folding into compact globular structures, with well-defined modules of secondary structure, appears to be a characteristic of long polypeptide chains, with a specific patterning of coded amino acid residues along the length of sequence. Cooperative hydrogen bond driven secondary structure formation and solvent forces, which contribute favorably to the entropy of folding, by promoting compaction of the polymeric chain, have long been discussed as major determinants of the folding process. First principles design approaches, which use non-coded amino acids, employ an alternative structure directing strategy, by using amino acid residues which exhibit a strong conformational bias for specific regions of the Ramachandran map. This overview of ongoing studies in the authors' laboratory, attempts to explore the use of conformationally restricted amino acid residues in the design of peptides with well-defined secondary structures. Short peptides composed of 20 genetically coded amino acids usually exist in solution as an ensemble of equilibrating conformations. Apolar peptide sequences, which are readily soluble in organic solvents like chloroform and methanol, facilitate formation of structures which are predominately driven by intramolecular hydrogen bond formation. The choice of sequences containing residues with a limited range of conformational choices strongly favors formation of local turn structures, stabilized by short range intramolecular hydrogen bonds. Two residue beta-turns can nucleate either helical or hairpin folding, depending on the precise conformation of the turn segment Restriction of the conformational space available to amino acid residues is easily achieved by introduction of an additional alkyl group at the C alpha carbon atom or by side chain backbone cyclization, as in proline. Studies of synthetic sequences incorporating two prototype residues alpha-aminoisobutyric acid (Aib) and D-proline (DPro) illustrate the utility of the strategy in construction of helices and hairpins. Extensions to the design of conformationally switchable sequences and structurally defined hybrid peptides containing backbone homologated residues are also surveyed.

Propensity to Form Amyloid Fibrils Is Encoded as Excitations in the Free Energy Landscape of Monomeric Proteins

Protein aggregation, linked to many of diseases, is initiated when monomers access rogue conformations that are poised to form amyloid fibrils. We show, using simulations of src SH3 domain, that mechanical force enhances the population of the aggregation-prone (N*) states, which are rarely populated under force free native conditions but are encoded in the spectrum of native fluctuations. The folding phase diagrams of SH3 as a function of denaturant concentration (C]), mechanical force (f), and temperature exhibit an apparent two-state behavior, without revealing the presence of the elusive N* states. Interestingly, the phase boundaries separating the folded and unfolded states at all C] and f fall on a master curve, which can be quantitatively described using an analogy to superconductors in a magnetic field. The free energy profiles as a function of the molecular extension (R), which are accessible in pulling experiments, (R), reveal the presence of a native-like N* with a disordered solvent-exposed amino-terminal beta-strand. The structure of the N* state is identical with that found in Fyn SH3 by NMR dispersion experiments. We show that the timescale for fibril formation can be estimated from the population of the N* state, determined by the free energy gap separating the native structure and the N* state, a finding that can be used to assess fibril forming tendencies of proteins. The structures of the N* state are used to show that oligomer formation and likely route to fibrils occur by a domain-swap mechanism in SH3 domain. (C) 2014 Elsevier Ltd. All rights reserved.

Tautomeric Preference and Conformation Locking in Fenobam, Thiofenobam, and Their Analogues: The Decisive Role of Hydrogen Bond Hierarchy

The crystal and molecular structures of the potential antidepressant drug fenobam and its derivatives are examined in terms of the preferred form among the two possible tautomeric structures. In this study, chemical derivatization has been utilized as a means to ``experimentally simulate'' the tautomeric preference and conformational variability in fenobam. Eight new derivatives of fenobam have been synthesized, and structural features have been characterized by single-crystal X-ray diffraction and NMR spectroscopy. The specific tautomeric preference found in all of these compounds and their known crystal forms have been construed in terms of the stabilizing intramolecular N-H center dot center dot center dot O and N-H center dot center dot center dot S hydrogen bonding. The hierarchy of intramolecular hydrogen bonds evidenced as the preference of the C-H center dot center dot center dot O hydrogen bond over C-H center dot center dot center dot N and that of the C-H center dot center dot center dot N hydrogen bond over C-H center dot center dot center dot S explains the two distinct conformations adopted by fenobam and thiofenobam derivatives. The relative energy values of different molecular conformations have been calculated and compared.

beta gamma-fused turn structures in sugar amino acid (SAA) containing cyclic tetrapeptides with alpha 3 delta architecture

The current manuscript describes conformational analysis of 15-membered cyclic tetrapeptides (CTPs), with alpha 3 delta architecture, containing sugar amino acids (SAA) having variation in the stereocenter at C5 carbon. Conformational analyses of both the series, in protected and deprotected forms, were carried out in DMSO-d(6) using various NMR techniques, supported by restrained MD calculations. It was intriguing to notice that the alpha 3 delta macrocycles got stabilized by both 10-membered beta-turn as well as a seven-membered gamma-turn, fused within the same macrocycle. The presence of fused sub-structures within a 15-membered macrocycle is rare to see. Also, the stereocenter variation at C5 did not affect the fused turn structures and exhibited similar conformations in both the series. The design becomes highly advantageous as fused reverse turn structures are occurring in the cyclic structure with minimalistic size macrocycle and this can be applied to develop suitable pharmacophores in the drug development process. (C) 2014 Elsevier Ltd. All rights reserved.

Succinate esters: odd-even effects in melting points

Dialkyl succinates show a pattern of alternating behavior in their melting points, as the number of C atoms in the alkane side chain increases, unlike in the dialkyl oxalates Joseph et al. (2011). Acta Cryst. B67, 525-534]. Dialkyl succinates with odd numbers of C atoms in the alkyl side chain show higher melting points than the immediately adjacent analogues with even numbers. The crystal structures and their molecular packing have been analyzed for a series of dialkyl succinates with 1 - 4 C atoms in the alkyl side chain. The energy difference (Delta E) between the optimized and observed molecular conformations, density, Kitaigorodskii packing index (KPI) and C-H center dot center dot center dot O interactions are considered to rationalize this behavior. In contrast to the dialkyl oxalates where a larger number of moderately strong C-H center dot center dot center dot O interactions were characteristic of oxalates with elevated melting points, here the molecular packing and the density play a major role in raising the melting point. On moving from oxalate to succinate esters the introduction of the C2 spacer adds two activated H atoms to the asymmetric unit, resulting in the formation of stronger C-H center dot center dot center dot O hydrogen bonds in all succinates. As a result the crystallinity of long-chain alkyl substituted esters improves enormously in the presence of hydrogen bonds from activated donors.

Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

Ligand 5,10,15,20-Tetra(N-methyl-4-pyridyl)porphine (TMPyP4) Prefers the Parallel Propeller-Type Human Telomeric G-Quadruplex DNA over Its Other Polymorphs

The binding of ligand 5,10,15,20-tetra(N-methyl-4-pyridyl)porphine (TMPyP4) with telomeric and genomic G-quadruplex DNA has been extensively studied. However, a comparative study of interactions of TMPyP4 with different conformations of human telomeric G-quadruplex DNA, namely, parallel propeller-type (PP), antiparallel basket-type (AB), and mixed hybrid-type (MH) G-quadruplex DNA, has not been done. We considered all the possible binding sites in each of the G-quadruplex DNA structures and docked TMPyP4 to each one of them. The resultant most potent sites for binding were analyzed from the mean binding free energy of the complexes. Molecular dynamics simulations were then carried out, and analysis of the binding free energy of the TMPyP4-G-quadruplex complex showed that the binding of TMPyP4 with parallel propeller-type G-quadruplex DNA is preferred over the other two G-quadruplex DNA conformations. The results obtained from the change in solvent excluded surface area (SESA) and solvent accessible surface area (SASA) also support the more pronounced binding of the ligand with the parallel propeller-type G-quadruplex DNA.

Directing peptide conformation with centrally positioned pre-organized dipeptide segments: studies of a 12-residue helix and beta-hairpin

Secondary structure formation in oligopeptides can be induced by short nucleating segments with a high propensity to form hydrogen bonded turn conformations. Type I/III turns facilitate helical folding while type II'/I' turns favour hairpin formation. This principle is experimentally verified by studies of two designed dodecapeptides, Boc-Val-Phe-Leu-Phe-Val-Aib-Aib-Val-Phe-Leu-Phe-Val-OMe 1 and Boc-Val-Phe-Leu-Phe-Val- (D) Pro- (L) Pro-Val-Phe-Leu-Phe-Val-OMe 2. The N- and C-terminal flanking pentapeptide sequences in both cases are identical. Peptide 1 adopts a largely alpha-helical conformation in crystals, with a small 3(10) helical segment at the N-terminus. The overall helical fold is maintained in methanol solution as evidenced by NMR studies. Peptide 2 adopts an antiparallel beta-hairpin conformation stabilized by 6 interstrand hydrogen bonds. Key nuclear Overhauser effects (NOEs) provide evidence for the antiparallel beta-hairpin structure. Aromatic proton chemical shifts provide a clear distinction between the conformation of peptides 1 (helical) and 2 (beta-hairpin). The proximity of facing aromatic residues positioned at non-hydrogen bonding positions in the hairpin results in extensively ring current shifted proton resonances in peptide 2.

Exploring the Consequences of a Representative ``Disallowed'' Conformation of Aib on a 3(10)-Helical Fold

The structural effects of a representative disallowed conformation of Aib on the 3(10)-helical fold of an octapeptidomimetic are explored. The 1D (H-1, C-13) & 2D NMR, FT-IR and CD data reveal that the octapeptide 1, adopts a 3(10)-helical conformation in solution, as it does in its crystal structure. The C-terminal methyl carboxylate (CO2Me) of 1 was modified into an 1,3-oxazine (Oxa) functional group in the peptidomimetic 2. This modification results in the stabilization of the backbone of the C-terminal Aib (Aib*-Oxa) of 2, in a conformation (phi, =180, 0) that is natively disallowed to Aib. Consequent to the presence of this natively disallowed conformation, the 3(10)-helical fold is not disrupted in the body of the peptidomimetic 2. But the structural distortions that do occur in 2 are primarily in residues in the immediate vicinity of the natively disallowed conformation, rather than in the whole peptide body. Non-native electronic effects resulting from modifications in backbone functional groups can be at the origin of stabilizing residues in natively disallowed conformations. (c) 2014 Wiley Periodicals, Inc. Biopolymers (Pept Sci) 104: 21-36, 2015.

Hierarchical sampling for metastable conformers determines biomolecular recognition: the case of malectin and diglucosylated N-glycan interactions

Structural information over the entire course of binding interactions based on the analyses of energy landscapes is described, which provides a framework to understand the events involved during biomolecular recognition. Conformational dynamics of malectin's exquisite selectivity for diglucosylated N-glycan (Dig-N-glycan), a highly flexible oligosaccharide comprising of numerous dihedral torsion angles, are described as an example. For this purpose, a novel approach based on hierarchical sampling for acquiring metastable molecular conformations constituting low-energy minima for understanding the structural features involved in a biologic recognition is proposed. For this purpose, four variants of principal component analysis were employed recursively in both Cartesian space and dihedral angles space that are characterized by free energy landscapes to select the most stable conformational substates. Subsequently, k-means clustering algorithm was implemented for geometric separation of the major native state to acquire a final ensemble of metastable conformers. A comparison of malectin complexes was then performed to characterize their conformational properties. Analyses of stereochemical metrics and other concerted binding events revealed surface complementarity, cooperative and bidentate hydrogen bonds, water-mediated hydrogen bonds, carbohydrate-aromatic interactions including CH-pi and stacking interactions involved in this recognition. Additionally, a striking structural transition from loop to beta-strands in malectin CRD upon specific binding to Dig-N-glycan is observed. The interplay of the above-mentioned binding events in malectin and Dig-N-glycan supports an extended conformational selection model as the underlying binding mechanism.

A Gaussian network model study suggests that structural fluctuations are higher for inactive states than active states of protein kinases

We performed Gaussian network model based normal mode analysis of 3-dimensional structures of multiple active and inactive forms of protein kinases. In 14 different kinases, a more number of residues (1095) show higher structural fluctuations in inactive states than those in active states (525), suggesting that, in general, mobility of inactive states is higher than active states. This statistically significant difference is consistent with higher crystallographic B-factors and conformational energies for inactive than active states, suggesting lower stability of inactive forms. Only a small number of inactive conformations with the DFG motif in the ``in'' state were found to have fluctuation magnitudes comparable to the active conformation. Therefore our study reports for the first time, intrinsic higher structural fluctuation for almost all inactive conformations compared to the active forms. Regions with higher fluctuations in the inactive states are often localized to the aC-helix, aG-helix and activation loop which are involved in the regulation and/or in structural transitions between active and inactive states. Further analysis of 476 kinase structures involved in interactions with another domain/protein showed that many of the regions with higher inactive-state fluctuation correspond to contact interfaces. We also performed extensive GNM analysis of (i) insulin receptor kinase bound to another protein and (ii) holo and apo forms of active and inactive conformations followed by multi-factor analysis of variance. We conclude that binding of small molecules or other domains/proteins reduce the extent of fluctuation irrespective of active or inactive forms. Finally, we show that the perceived fluctuations serve as a useful input to predict the functional state of a kinase.

Composition dependent multiple structural transformations of myoglobin in aqueous ethanol solution: A combined experimental and theoretical study

Experimental studies (circular dichroism and ultra-violet (UV) absorption spectra) and large scale atomistic molecular dynamics simulations (accompanied by order parameter analyses) are combined to establish a number of remarkable (and unforeseen) structural transformations of protein myoglobin in aqueous ethanol mixture at various ethanol concentrations. The following results are particularly striking. (1) Two well-defined structural regimes, one at x(EtOH) similar to 0.05 and the other at x(EtOH) similar to 0.25, characterized by formation of distinct partially folded conformations and separated by a unique partially unfolded intermediate state at x(EtOH) similar to 0.15, are identified. (2) Existence of non-monotonic composition dependence of (i) radius of gyration, (ii) long range contact order, (iii) residue specific solvent accessible surface area of tryptophan, and (iv) circular dichroism spectra and UV-absorption peaks are observed. Interestingly at x(EtOH) similar to 0.15, time averaged value of the contact order parameter of the protein reaches a minimum, implying that this conformational state can be identified as a molten globule state. Multiple structural transformations well known in water-ethanol binary mixture appear to have considerably stronger effects on conformation and dynamics of the protein. We compare the present results with studies in water-dimethyl sulfoxide mixture where also distinct structural transformations are observed along with variation of co-solvent composition. (C) 2015 AIP Publishing LLC.