Network approach for capturing ligand-induced subtle global changes in protein structures

Ligand-induced conformational changes in proteins are of immense functional relevance. It is a major challenge to elucidate the network of amino acids that are responsible for the percolation of ligand-induced conformational changes to distal regions in the protein from a global perspective. Functionally important subtle conformational changes (at the level of side-chain noncovalent interactions) upon ligand binding or as a result of environmental variations are also elusive in conventional studies such as those using root-mean-square deviations (r.m.s.d.s). In this article, the network representation of protein structures and their analyses provides an efficient tool to capture these variations (both drastic and subtle) in atomistic detail in a global milieu. A generalized graph theoretical metric, using network parameters such as cliques and/or communities, is used to determine similarities or differences between structures in a rigorous manner. The ligand-induced global rewiring in the protein structures is also quantified in terms of network parameters. Thus, a judicious use of graph theory in the context of protein structures can provide meaningful insights into global structural reorganizations upon perturbation and can also be helpful for rigorous structural comparison. Data sets for the present study include high-resolution crystal structures of serine proteases from the S1A family and are probed to quantify the ligand-induced subtle structural variations.

The structures of bacteriorhodopsin with different retinal-Schiff base orientations--computer modeling and energy minimization studies

Bacteriorhodopsin has been the subject of intense study in order to understand its photochemical function. The recent atomic model proposed by Henderson and coworkers based on electron cryo-microscopic studies has helped in understanding many of the structural and functional aspects of bacteriorhodopsin. However, the accuracy of the positions of the side chains is not very high since the model is based on low-resolution data. In this study, we have minimized the energy of this structure of bacteriorhodopsin and analyzed various types of interactions such as - intrahelical and interhelical hydrogen bonds and retinal environment. In order to understand the photochemical action, it is necessary to obtain information on the structures adopted at the intermediate states. In this direction, we have generated some intermediate structures taking into account certain experimental data, by computer modeling studies. Various isomers of retinal with 13-cis and/or 15-cis conformations and all possible staggered orientations of Lys-216 side chain were generated. The resultant structures were examined for the distance between Lys-216-schiff base nitrogen and the carboxylate oxygen atoms of Asp-96 - a residue which is known to reprotonate the schiff base at later stages of photocycle. Some of the structures were selected on the basis of suitable retinal orientation and the stability of these structures were tested by energy minimization studies. Further, the minimized structures are analyzed for the hydrogen bond interactions and retinal environment and the results are compared with those of the minimized rest state structure. The importance of functional groups in stabilizing the structure of bacteriorhodopsin and in participating dynamically during the photocycle have been discussed.

A novel structural derivative of natural alkaloid ellipticine, MDPSQ, induces necrosis in leukemic cells

DNA intercalating molecules are promising chemotherapeutic agents. In the present study, a novel DNA intercalating compound of pyrimido4',5':4,5]selenolo(2,3-b)quinoline series having 8-methyl-4-(3 diethylaminopropylamino) side chain is studied for its chemotherapeutic properties. Our results showed that 8-methyl-4-(3 diethylaminopropylamino) pyrimido 4',5':4,5] selenolo(2,3-b)quinoline (MDPSQ) induces cytotoxicity in a time- and concentration-dependent manner on leukemic cell lines. Both cell cycle analysis and tritiated thymidine assays revealed that MDPSQ affects DNA replication. Treatment with MDPSQ resulted in both elevated levels of DNA strand breaks and repair proteins, further indicating its cytotoxic effects. Besides, Annexin V/PI staining revealed that MDPSQ induces cell death by triggering necrosis rather than apoptosis.

Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

On the importance of backbone loop and peptide flip in the Walker sequence in F-1-ATPase action

The Walker sequence, GXXXXGKT, present in all the six subunits of F-1-ATPase exists in a folded form, known as phosphate-binding loop (P-loop). Analysis of the Ramachandran angles showed only small RMS deviation between the nucleotide-bound and nucleotide-free forms. This indicated a good overlap of the backbone loops. The catalytic beta-subunits (chains D, E and F) showed significant changes in the Ramachandran angles and the side chain torsion angles, but not the structural alpha-subunits (chains A, B and C). Most striking among these are the changes associated with Val160 and Gly161 corresponding to a flip in the peptide unit between them when a nucleotide is bound (chains D or F compared to nucleotide-free chain E). The conformational analysis further revealed a hitherto unnoticed hydrogen bond between amide-N of the flipped Gly161 and terminal phosphate-O of the nucleotide. This assigns a role for this conserved amino acid, otherwise ignored, of making an unusual direct interaction between the peptide backbone of the enzyme protein and the incoming nucleotide substrate. Significance of this interaction is enhanced, as it is limited only to the catalytic subunits, and also likely to involve a mechanical rotation of bonds of the peptide unit. Hopefully this is part of the overall events that link the chemical hydrolysis of ATP with the mechanical rotation of this molecule, now famous as tiny molecular motor.

Pyrene appended bile acid conjugates: Synthesis and a structure-gelation property study

A wide variety of novel compounds obtained by combining two types of known organogelators, viz., bile acid alkyl amides and pyrene alkanoic acids, were synthesized and screened for their gelation ability. The 3 alpha esters of 1-pyrene butyric acid (PBA) of alkylamides of deoxycholic acid (DCA) turned out to be effective in the gel formation with many organic solvents although the gelation has to be triggered by the addition of a charge transfer (CT) agent 2,4,7-trinitrofluorenone (TNF). The special feature of these molecules is that the organogelation is achieved only after derivatizing the acid moiety of the 1-pyrenealkanoic acids. Additionally, the gelation properties can be fine-tuned by inserting different functional groups at the bile acid side chain. The gels obtained are deep red in colour and optically transparent up to 2% w/v. The SEM studies of the obtained xerogels revealed bundled rod-like morphology without specialized branching.

Revisiting the Mechanism of the Triosephosphate Isomerase Reaction: The Role of the Fully Conserved Glutamic Acid 97 Residue

An analysis of 503 available triosephosphate isomerase sequences revealed nine fully conserved residues. Of these, four residues-K12, H95, E97 and E165-are capable of proton transfer and are all arrayed around the dihydroxyacetone phosphate substrate in the three-dimensional structure. Specific roles have been assigned to the residues K12, H95 and E165, but the nature of the involvement of E97 has not been established. Kinetic and structural characterization is reported for the E97Q and E97D mutants of Plasmodium falciparum triosephosphate isomerase (Pf TIM). A 4000-fold reduction in k(cat) is observed for E97Q, whereas the E97D mutant shows a 100-fold reduction. The control mutant, E165A, which lacks the key catalytic base, shows an approximately 9000-fold drop in activity. The integrity of the overall fold and stability of the dimeric structure have been demonstrated by biophysical studies. Crystal structures of E97Q and E97D mutants have been determined at 2.0 angstrom resolution. In the case of the isosteric replacement of glutamic acid by glutamine in the E97Q mutant a large conformational change for the critical K12 side chain is observed, corresponding to a trans-to-gauche transition about the C gamma-C delta (chi(3)) bond. In the E97D mutant, the K12 side chain maintains the wild-type orientation, but the hydrogen bond between K12 and D97 is lost. The results are interpreted as a direct role for E97 in the catalytic proton transfer cycle. The proposed mechanism eliminates the need to invoke the formation of the energetically unfavourable imidazolate anion at H95, a key feature of the classical mechanism.

Structural basis for the functional and inhibitory mechanisms of beta-hydroxyacyl-acyl carrier protein dehydratase (FabZ) of Plasmodium falciparum

The beta-hydroxyacyl-acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) catalyzes the third and important reaction of the fatty acid elongation cycle. The crystal structure of PfFabZ is available in hexameric (active) and dimeric (inactive) forms. However, PfFabZ has not been crystallized with any bound inhibitors until now. We have designed a new condition to crystallize PfFabZ with its inhibitors bound in the active site, and determined the crystal structures of four of these complexes. This is the first report on any FabZ enzyme with active site inhibitors that interact directly with the catalytic residues. Inhibitor binding not only stabilized the substrate binding loop but also revealed that the substrate binding tunnel has an overall shape of ``U''. In the crystal structures, residue Phe169 located in the middle of the tunnel was found to be in two different conformations, open and closed. Thus, Phe169, merely by changing its side chain conformation, appears to be controlling the length of the tunnel to make it suitable for accommodating longer substrates. The volume of the substrate binding tunnel is determined by the sequence as well as by the conformation of the substrate binding loop region and varies between organisms for accommodating fatty acids of different chain lengths. This report on the crystal structures of the complexes of PfFabZ provides the structural basis of the inhibitory mechanism of the enzyme that could be used to improve the potency of inhibitors against an important component of fatty acid synthesis common to many infectious organisms. (C) 2011 Elsevier Inc. All rights reserved.

Investigating Diproline Segments in Proteins: Occurrences, Conformation and Classification

The covalent linkage between the side-chain and the backbone nitrogen atom of proline leads to the formation of the five-membered pyrrolidine ring and hence restriction of the backbone torsional angle phi to values of -60 degrees +/- 30 degrees for the L-proline. Diproline segments constitute a chain fragment with considerably reduced conformational choices. In the current study, the conformational states for the diproline segment ((L)Pro-(L)Pro) found in proteins has been investigated with an emphasis on the cis and trans states for the Pro-Pro peptide bond. The occurrence of diproline segments in turns and other secondary structures has been studied and compared to that of Xaa-Pro-Yaa segments in proteins which gives us a better understanding on the restriction imposed on other residues by the diproline segment and the single proline residue. The study indicates that P(II)-P(II) and P(II)-alpha are the most favorable conformational states for the diproline segment. The analysis on Xaa-Pro-Yaa sequences reveals that the XaaPro peptide bond exists preferably as the trans conformer rather than the cis conformer. The present study may lead to a better understanding of the behavior of proline occurring in diproline segments which can facilitate various designed diproline-based synthetic templates for biological and structural studies. (C) 2011 Wiley Periodicals, Inc. Biopolymers 97: 54-64, 2012.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Network properties of protein-decoy structures

Convergence of the vast sequence space of proteins into a highly restricted fold/conformational space suggests a simple yet unique underlying mechanism of protein folding that has been the subject of much debate in the last several decades. One of the major challenges related to the understanding of protein folding or in silico protein structure prediction is the discrimination of non-native structures/decoys from the native structure. Applications of knowledge-based potentials to attain this goal have been extensively reported in the literature. Also, scoring functions based on accessible surface area and amino acid neighbourhood considerations were used in discriminating the decoys from native structures. In this article, we have explored the potential of protein structure network (PSN) parameters to validate the native proteins against a large number of decoy structures generated by diverse methods. We are guided by two principles: (a) the PSNs capture the local properties from a global perspective and (b) inclusion of non-covalent interactions, at all-atom level, including the side-chain atoms, in the network construction accommodates the sequence dependent features. Several network parameters such as the size of the largest cluster, community size, clustering coefficient are evaluated and scored on the basis of the rank of the native structures and the Z-scores. The network analysis of decoy structures highlights the importance of the global properties contributing to the uniqueness of native structures. The analysis also exhibits that the network parameters can be used as metrics to identify the native structures and filter out non-native structures/decoys in a large number of data-sets; thus also has a potential to be used in the protein `structure prediction' problem.

Structural characterization of folded pentapeptides containing centrally positioned beta(R)Val, gamma(R)Val and gamma(S)Val residues

A cylindrical pore of similar to 7.5 angstrom diameter containing a one-dimensional water wire, within the confines of a hydrophobic channel lined with the valine side chain, has been observed in crystals of the peptide Boc-D-Pro-Aib-Val-Aib-Val-OMe (1) (Raghavender et al., 2009, 2010). The synthesis and structural characterization in crystals of three backbone homologated analogues Boc-D-Pro-Aib-beta(3)(R) Val-Aib-Val-OMe (2), Boc-D-Pro-Aib-gamma(4)(R)Val-Aib-Val-OMe (3), Boc-D-Pro-Aib-gamma(4)(S)Val-Aib-Val-OMe (4) are described. Crystal structures of peptides 2, 3 and 4 reveal close-packed arrangements in which no pore was formed. In peptides 2 and 3 the N-terminus D-Pro-Aib segment adopted conformations closely related to Type II' beta-turns, while residues 2-4 form one turn of an alpha beta right-handed C-11 helix in 2 and an alpha gamma C-12 helix in 3. In peptide 4, a continuous left-handed helical structure was observed with the D-Pro-Aib segment forming a Type III' beta-turn, followed by one turn of a left-handed alpha gamma C-12 helix. (C) 2012 Elsevier Ltd. All rights reserved.

A touch of history and a peep into the future of the lipid-quinone known as coenzyme Q and ubiquinone

Coenzyme Q (ubiquinone), a fully substituted benzoquinone with polyprenyl side chain, participates in many cellular redox activities. Paradoxically it was discovered only in 1957, albeit being ubiquitous. It required a person, F. L. Crane, a place, Enzyme Institute, Madison, USA, and a time when D. E. Green was directing vigorous research on mitochondria. Located at the transition of 2-electron flavoproteins and 1-electron cytochrome carriers, it facilitates electron transfer through the elegant Q-cycle in mitochondria to reduce O-2 to H2O, and to H2O2, now a significant signal-transducing agent, as a minor activity in shunt pathway (animals) and alternative oxidase (plants). The ability to form Q-radical by losing an electron and a proton was ingeniously used by Mitchell to explain the formation of the proton gradient, considered the core of energy transduction, and also in acidification in vacuoles. Known to be a mobile membrane constituent (microsomes, plasma membrane and Golgi apparatus), allowing it to reach multiple sites, coenzyme Q is expected to have other activities. Coenzyme Q protects circulating lipoproteins being a better lipid antioxidant than even vitamin E. Binding to proteins such as QPS, QPN, QPC and uncoupling protein in mitochondria, QA and QB in the reaction centre in R. sphaeroides, and disulfide bond-forming protein in E. coli (possibly also in Golgi), coenzyme Q acquires selective functions. A characteristic of orally dosed coenzyme Q is its exclusive absorption into the liver, but not the other tissues. This enrichment of Q is accompanied by significant decrease of blood pressure and of serum cholesterol. Inhibition of formation of mevalonate, the common precursor in the branched isoprene pathway, by the minor product, coenzyme Q, decreases the major product, cholesterol. Relaxation of contracted arterial smooth muscle by a side-chain truncated product of coenzyme Q explains its effect of decreasing blood pressure. Extensive clinical studies carried out on oral supplements of coenzyine Q, initially by K. Folkers and Y. Yamamura and followed many others, revealed a large number of beneficial effects, significantly in cardiovascular diseases. Such a variety of effects by this lipid quinone cannot depend on redox activity alone. The fat-soluble vitamins (A, D, E and K) that bear structural relationship with coenzyme Q are known to be active in their polar forms. A vignette of modified forms of coenzyme Q taking active role in its multiple effects is emerging.

Designed three-stranded beta-sheet in an alpha/beta hybrid peptide

The incorporation of beta-amino acid residues into the antiparallel beta-strand segments of a multi-stranded beta-sheet peptide is demonstrated for a 19-residue peptide, Boc-LV(beta)FV(D)PGL(beta)FVVL(D)PGLVL(beta)FVV-OMe (BBH19). Two centrally positioned (D)Pro-Gly segments facilitate formation of a stable three-stranded beta-sheet, in which beta-phenylalanine ((beta)Phe) residues occur at facing positions 3, 8 and 17. Structure determination in methanol solution is accomplished by using NMR-derived restraints obtained from NOEs, temperature dependence of amide NH chemical shifts, rates of H/D exchange of amide protons and vicinal coupling constants. The data are consistent with a conformationally well-defined three-stranded beta-sheet structure in solution. Cross-strand interactions between (beta)Phe3/(beta)Phe17 and (beta)Phe3/Val15 residues define orientations of these side-chains. The observation of close contact distances between the side-chains on the N- and C-terminal strands of the three-stranded beta-sheet provides strong support for the designed structure. Evidence is presented for multiple side-chain conformations from an analysis of NOE data. An unusual observation of the disappearance of the Gly NH resonances upon prolonged storage in methanol is rationalised on the basis of a slow aggregation step, resulting in stacking of three-stranded beta-sheet structures, which in turn influences the conformational interconversion between type I' and type II' beta-turns at the two (D)Pro-Gly segments. Experimental evidence for these processes is presented. The decapeptide fragment Boc-LV(beta)FV(D)PGL(beta)FVV-OMe (BBH10), which has been previously characterized as a type I' beta-turn nucleated hairpin, is shown to favour a type II' beta-turn conformation in solution, supporting the occurrence of conformational interconversion at the turn segments in these hairpin and sheet structures.

Chemical unfolding of chicken villin headpiece in aqueous dimethyl sulfoxide solution: cosolvent concentration dependence, pathway, and microscopic mechanism

Unfolding of a protein often proceeds through partial unfolded intermediate states (PUIS). PUIS have been detected in several experimental and simulation studies. However, complete analyses of transitions between different PUIS and the unfolding trajectory are sparse. To understand such dynamical processes, we study chemical unfolding of a small protein, chicken villin head piece (HP-36), in aqueous dimethyl sulfoxide (DMSO) solution. We carry out molecular dynamics simulations at various solution compositions under ambient conditions. In each concentration, the initial step of unfolding involves separation of two adjacent native contacts, between phenyl alanine residues (11-18 and 7-18). This first step induces, under appropriate conditions, subsequent separation among other hydrophobic contacts, signifying a high degree of cooperativity in the unfolding process. The observed sequence of structural changes in HP-36 on increasing DMSO concentration and the observed sequence of PUIS, are in approximate agreement with earlier simulation results (in pure water) and experimental observations on unfolding of HP-36. Peculiar to water-DMSO mixture, an intervening structural transformation (around 15% of DMSO) in the binary mixture solvent retards the progression of unfolding as composition is increased. This is reflected in a remarkable nonmonotonic composition dependence of RMSD, radius of gyration and the fraction of native contacts. At 30% mole fraction of DMSO, we find the extended randomly coiled structure of the unfolded protein. The molecular mechanism of DMSO induced unfolding process is attributed to the initial preferential solvation of the hydrophobic side chain atoms through the methyl groups of DMSO, followed by the hydrogen bonding of the oxygen atom of DMSO to the exposed backbone NH groups of HP-36.