Molecules with new topologies derived from hydrogen-bonded dimers of tetraurea calix[4]arenes

Tetraurea calix[4]arenes 2 have been synthesized in which two adjacent aryl urea residues are connected to a loop by an aliphatic chain -O-(CH(2))(n)-O-. The remaining urea residues have a bulky 3,5-di-tert-butylphenyl residue and an omega-alkenyloxyphenyl residue. Since this bulky residue cannot pass through the loop, only one homodimer (22) is formed in apolar solvents, for steric reasons, in which the two alkenyl residues penetrate the two macrocyclic loops. Covalent connection of these alkenyl groups by olefin metathesis followed by hydrogenation creates compounds 3, which consist of molecules with hitherto unknown topology. Their molecular structure was confirmed by (1)H NMR spectroscopy and ESIMS, and for one example by single-crystal X-ray analysis.     

Molecular recognition properties of FN3 monobodies that bind the Src SH3 domain

We have constructed a phage-displayed library based on the human fibronectin tenth type III domain (FN3) scaffold by randomizing residues in its FG and BC loops. Screening against the SH3 domain of human c-Src yielded six different clones. Five of these contained proline-rich sequences in their FG loop that resembled class I (i.e., +xxPxxP) peptide ligands for the Src SH3 domain. The sixth clone lacked the proline-rich sequence and showed particularly high binding specificity to the Src SH3 domain among various SH3 domains tested. Competitive binding, loop replacement, and NMR perturbation experiments were conducted to analyze the recognition properties of selected binders. The strongest binder was able to pull down full-length c-Src from murine fibroblast cell extracts, further demonstrating the potential of this scaffold for use as an antibody mimetic.     

Predicting three-dimensional structures of transmembrane domains of β-barrel membrane proteins

β-Barrel membrane proteins are found in the outer membrane of gram-negative bacteria, mitochondria, and chloroplasts. They are important for pore formation, membrane anchoring, and enzyme activity. These proteins are also often responsible for bacterial virulence. Due to difficulties in experimental structure determination, they are sparsely represented in the protein structure databank. We have developed a computational method for predicting structures of the transmembrane (TM) domains of β-barrel membrane proteins. Based on physical principles, our method can predict structures of the TM domain of β-barrel membrane proteins of novel topology, including those from eukaryotic mitochondria. Our method is based on a model of physical interactions, a discrete conformational state space, an empirical potential function, as well as a model to account for interstrand loop entropy. We are able to construct three-dimensional atomic structure of the TM domains from sequences for a set of 23 nonhomologous proteins (resolution 1.8-3.0 ?). The median rmsd of TM domains containing 75-222 residues between predicted and measured structures is 3.9 ? for main chain atoms. In addition, stability determinants and protein-protein interaction sites can be predicted. Such predictions on eukaryotic mitochondria outer membrane protein Tom40 and VDAC are confirmed by independent mutagenesis and chemical cross-linking studies. These results suggest that our model captures key components of the organization principles of β-barrel membrane protein assembly.     

Proteome-wide classification and identification of mammalian-type GPCRs by binary topology pattern

G protein-coupled receptors (GPCRs), a large eukaryotic protein family, have proved difficult to comprehensively detect and functionally identify by homology searches and domain detection, because they are highly divergent and their sequences share strikingly little similarity. Transmembrane (TM) topology pattern analysis has been used to classify TM proteins, and such patterns are conserved within GPCRs of similar function. Here, we developed a stepwise binary topology pattern (BTP) method for GPCR classification and identification and used it to identify and classify mammalian-type GPCRs in the genomes of 10 different eukaryotic species. A binary topology pattern was obtained for each functional class or group by assigning binary loop threshold lengths of "0" (short loop) or "1" (long loop). The GPCR-classification ability of the BTP method had quite high accuracies for classifying GPCR functions at the class level (Classes A, B, C, Frizzled/Smoothened, Non-GPCR, based on the GPCRDB classification scheme), with many classes being classified with 100% accuracy. Sufficiently high accuracies were also maintained at the functional group level, 0.945 over 15 functional groups. Proteome-wide mammalian-type GPCR searches in 10 eukaryotic genomes (H. sapiens, M. musculus, F. rubripes, C. intestinalis, A. thaliana, D. melanogaster, A. gambiae, C. elegans, P. falciparum, S. cerevisiae) using the BTP method showed much higher classification/identification in non-mammalian genomes than typical BLAST searches, in which a higher number of sequences were classified as Non-GPCR. This stepwise BTP method should prove useful for the identification and functional classification of GPCRs from the genomes of a wide range of species.     

Energy landscapes of dynamic ensembles of rolling triplet repeat bulge loops: implications for DNA expansion associated with disease states

DNA repeat domains can form ensembles of canonical and noncanonical states, including stable and metastable DNA secondary structures. Such sequence-induced structural diversity creates complex conformational landscapes for DNA processing pathways, including those triplet expansion events that accompany replication, recombination, and/or repair. Here we demonstrate further levels of conformational complexity within repeat domains. Specifically, we show that bulge loop structures within an extended repeat domain can form dynamic ensembles containing a distribution of loop positions, thereby yielding families of positional loop isomers, which we designate as "rollamers". Our fluorescence, absorbance, and calorimetric data are consistent with loop migration/translocation between sites within the repeat domain ("rollamerization"). We demonstrate that such "rollameric" migration of bulge loops within repeat sequences can invade and disrupt previously formed base-paired domains via an isoenthalpic, entropy-driven process. We further demonstrate that destabilizing abasic lesions alter the loop distributions so as to favor "rollamers" with the lesion positioned at the duplex/loop junction, sites where the flexibility of the abasic "universal hinge" relaxes unfavorable interactions and/or facilitates topological accommodation. Another strategic siting of an abasic site induces directed loop migration toward denaturing domains, a phenomenon that merges destabilizing domains. In the aggregate, our data reveal that dynamic ensembles within repeat domains profoundly impact the overall energetics of such DNA constructs as well as the distribution of states by which they denature/renature. These static and dynamic influences within triplet repeat domains expand the conformational space available for selection and targeting by the DNA processing machinery. We propose that such dynamic ensembles and their associated impact on DNA properties influence pathways that lead to DNA expansion.     

A Thermodynamic Perspective on Potential G-Quadruplex Structures as Silencer Elements in the MYC Promoter

Multiple G-tracts within the promoter region of the c-myc oncogene may fold into various G-quadruplexes with the recruitment of different tracts and guanosine residues for the G-core assembly. Thermodynamic profiles for the folding of wild-type and representative truncated as well as mutated sequences were extracted by comprehensive DSC experiments. The unique G-quadruplex involving consecutive G-tracts II–V with formation of two one-nucleotide and one central two-nucleotide propeller loop, previously proposed to be the biologically most relevant species, was found to be the most stable fold in terms of its Gibbs free energy of formation at ambient temperatures. Its stability derives from its short propeller loops but also from the favorable type of loop residues. Whereas quadruplex folds with long propeller loops are significantly disfavored, a snap-back loop structure formed by incorporating a 3’-terminal guanosine into the empty position of a tetrad seems highly competitive based on its thermodynamic stability. However, its destabilization by extending the 3’-terminus questions the significance of such a species under in vivo conditions.     

Molecular dynamics of proteorhodopsin in lipid bilayers by solid-state NMR

Environmental factors such as temperature, hydration, and lipid bilayer properties are tightly coupled to the dynamics of membrane proteins. So far, site-resolved data visualizing the protein's response to alterations in these factors are rare, and conclusions had to be drawn from dynamic data averaged over the whole protein structure. In the current study, high-resolution solid-state NMR at high magnetic field was used to investigate their effects on the molecular dynamics of green proteorhodopsin, a bacterial light-driven proton pump. Through-space and through-bond correlation experiments were employed to identify and characterize highly mobile and motionally restricted regions of proteorhodopsin. Our data show that hydration water plays an essential role for enhancing molecular dynamics of residues in tails and interhelical loops, while it is found less important for residues in transmembrane domains or rigid, structured loop segments. In contrast, switching the lipids from the gel to their liquid crystalline phase enhances molecular fluctuations mainly in transmembrane helices on a time scale of 10(-6) s, but has little effect on loop and tail residues. Increased mobility is especially observed in helices C, F, and G, but also in the EF loop. Fluctuations in those regions are relevant to structural dynamics during the photocycle of proteorhodopsin. Our data are important for the functional understanding of proteorhodopsin, but also offer an important contribution to the general understanding of site-resolved effects of water and lipid bilayers onto the dynamic properties of membrane proteins.     

The challenge of annotating protein sequences: The tale of eight domains of unknown function in Pfam

The Pfam database is an important tool in genome annotation, since it provides a collection of curated protein families. However, a subset of these families, known as domains of unknown function (DUFs), remains poorly characterized. We have related sequences from DUF404, DUF407, DUF482, DUF608, DUF810, DUF853, DUF976 and DUF1111 to homologs in PDB, within the midnight zone (9–20%) of sequence identity. These relationships were extended to provide functional annotation by sequence analysis and model building. Also described are examples of residue plasticity within enzyme active sites, and change of function within homologous sequences of a DUF.     

Loop-length-dependent folding of G-quadruplexes

Guanine-rich DNA sequences can form a large number of structurally diverse quadruplexes. These vary in terms of strand polarity, loop composition, and conformation. We have derived guidelines for understanding the influence of loop length on the structure adopted by intramolecular quadruplex-forming sequences, using a combination of experimental (using CD and UV melting data) and molecular modeling and simulation techniques. We find that a parallel-stranded intramolecular quadruplex structure is the only possible fold when three single residue loops are present. When single thymine loops are present in combination with longer length loops, or when all loops are longer than two residues, both parallel- and antiparallel-folded structures are able to form. Multiple conformations of each structure are likely to coexist in solution, as they were calculated to have very similar free energies.     

Insight into G-DNA structural polymorphism and folding from sequence and loop connectivity through free energy analysis

The lengths of G-tracts and their connecting loop sequences determine G-quadruplex folding and stability. Complete understanding of the sequence-structure relationships remains elusive. Here, single-loop G-quadruplexes were investigated using explicit solvent molecular dynamics (MD) simulations to characterize the effect of loop length, loop sequence, and G-tract length on the folding topologies and stability of G-quadruplexes. Eight loop types, including different variants of lateral, diagonal, and propeller loops, and six different loop sequences [d0 (i.e., no intervening residues in the loop), dT, dT(2), dT(3), dTTA, and dT(4)] were considered through MD simulation and free energy analysis. In most cases the free energetic estimates agree well with the experimental observations. The work also provides new insight into G-quadruplex folding and stability. This includes reporting the observed instability of the left propeller loop, which extends the rules for G-quadruplex folding. We also suggest a plausible explanation why human telomere sequences predominantly form hybrid-I and hybrid-II type structures in K(+) solution. Overall, our calculation results indicate that short loops generally are less stable than longer loops, and we hypothesize that the extreme stability of sequences with very short loops could possibly derive from the formation of parallel multimers. The results suggest that free energy differences, estimated from MD and free energy analysis with current force fields and simulation protocols, are able to complement experiment and to help dissect and explain loop sequence, loop length, and G-tract length and orientation influences on G-quadruplex structure.     

End-to-end vs interior loop formation kinetics in unfolded polypeptide chains

The conformational search for favorable intramolecular interactions during protein folding is limited by intrachain diffusion processes. Recent studies on the dynamics of loop formation in unfolded polypeptide chains have focused on loops involving residues near the chain ends. During protein folding, however, most contacts are formed between residues in the interior of the chain. We compared the kinetics of end-to-end loop formation (type I loops) to the formation of end-to-interior (type II loops) and interior-to-interior loops (type III loops) using triplet-triplet energy transfer from xanthone to naphthylalanine. The results show that formation of type II and type III loops is slower compared to type I loops of the same size and amino acid sequence. The rate constant for type II loop formation decreases with increasing overall chain dimensions up to a limiting value, at which loop formation is about 2.5-fold slower for type II loops compared to type I loops. Comparing type II loops of different loop size and amino acid sequence shows that the ratio of loop dimension over total chain dimension determines the rate constant for loop formation. Formation of type III loops is 1.7-fold slower than formation of type II loops, indicating that local chain motions are strongly coupled to motions of other chain segments which leads to faster dynamics toward the chain ends. Our results show that differences in the kinetics of formation of type I, type II, and type III loops are mainly caused by differences in internal flexibility at the different positions in the polypeptide chain. Interactions of the polypeptide chain with the solvent contribute to the kinetics of loop formation, which are strongly viscosity-dependent. However, the observed differences in the kinetics of formation of type I, type II, and type III loops are not due to the increased number of peptide-solvent interactions in type II and type III loops compared to type I loops as indicated by identical viscosity dependencies for the kinetics of formation of the different types of loops.     

Formation of (3+1) G-quadruplexes with a long loop by human telomeric DNA spanning five or more repeats

Structural studies of human telomeric repeats represent an active field of research with potential applications toward the development of specific telomeric quadruplex-targeting drugs for anticancer treatment. To date, high-definition structures were limited to DNA sequences containing up to four GGGTTA repeats. Here we investigate the formation of G-quadruplexes in sequences spanning five to seven human telomeric repeats using NMR, UV, and CD spectroscopy. A (3+1) G-quadruplex with a long propeller loop was isolated from a five-repeat sequence utilizing a guanine-to-inosine substitution. A simple approach of selective site-specific labeling of guanine residues was devised to rigorously determine the folding topology of the oligonucleotide. The same scaffold could be extrapolated to six- and seven-repeat sequences. Our results suggest that long human telomeric sequences consisting of five or more GGGTTA repeats could adopt (3+1) G-quadruplex structures harboring one or more repeat(s) within a single loop. We report on the formation of a Watson-Crick duplex within the long propeller loop upon addition of the complementary strand, demonstrating that the long loop could serve as a new recognition motif.     

Site selection in tandem arrays of metal-binding domains

The metal binding properties of peptides corresponding to metal-binding sites spanning regions that normally function as linkers in tandem arrays of metal-binding domain-containing proteins were examined. For a peptide with two His residues from one TFIIIA-like zinc finger domain, a canonical TFIIIA-like linker, and two Cys residues from an adjacent zinc domain, the dissociation constant for the 1:1 peptide to cobalt(II) was found to be 15 +/- 10 microM, compared with 60 nM for the corresponding zinc finger domains themselves. Peptides overlapping two sets of metal-binding domains from human TRAF (tumor necrosis factor receptor-associated factor) proteins were examined. In one case, the affinity of the presumed metal-binding domain and that for the linker region were comparable, while in the second case, the affinity of the linker peptide was higher than that for the corresponding presumed metal-binding domain peptide. These studies revealed that cobalt(II) affinities in the micromolar range can occur even for peptides that do not correspond to natural zinc-binding domains and that the degree of distinction between authentic metal-binding domains and the corresponding linker-spanning peptides may be modest, at least for single domain peptide models.     

Analysis of sequence repeats of proteins in the PDB

Internal repeats in protein sequences play a significant role in the evolution of protein structure and function. Applications of different bioinformatics tools help in the identification and characterization of these repeats. In the present study, we analyzed sequence repeats in a non-redundant set of proteins available in the Protein Data Bank (PDB). We used RADAR for detecting internal repeats in a protein, PDBeFOLD for assessing structural similarity, PDBsum for finding functional involvement and Pfam for domain assignment of the repeats in a protein. Through the analysis of sequence repeats, we found that identity of the sequence repeats falls in the range of 20–40% and, the superimposed structures of the most of the sequence repeats maintain similar overall folding. Analysis sequence repeats at the functional level reveals that most of the sequence repeats are involved in the function of the protein through functionally involved residues in the repeat regions. We also found that sequence repeats in single and two domain proteins often contained conserved sequence motifs for the function of the domain.     

Structural basis of telomeric RNA quadruplex--acridine ligand recognition

Human telomeric DNA is now known to be transcribed into noncoding RNA sequences, termed TERRA. These sequences, which are believed to play roles in the regulation of telomere function, can form higher-order quadruplex structures and may themselves be the target of therapeutic intervention. The crystal structure of a TERRA quadruplex-acridine small-molecule complex at a resolution of 2.60 ?, is reported here and contrasts remarkably with the structure of the analogous DNA quadruplex complex. The bimolecular RNA complex has a parallel-stranded topology with propeller-like UUA loops. These loops are held in particular conformations by multiple hydrogen bonds involving the O2' hydroxyl groups of the ribonucleotide sugars and play an active role in binding the acridine molecules to the RNA quadruplex. By contrast, the analogous DNA quadruplex complex has simpler 1:1 acridine binding, with no loop involvement. There are significant loop conformational changes in the RNA quadruplex compared to the native TERRA quadruplex (Collie, G. W.; Haider, S. M.; Neidle, S.; Parkinson, G. N. Nucleic Acids Res. 2010, 38, 5569 - 5580), which have implications for the future design of small molecules targeting TERRA quadruplexes, and RNA quadruplexes more generally.     

Topology variation and loop structural homology in crystal and simulated structures of a bimolecular DNA quadruplex

The topology of DNA quadruplexes depends on the nature and number of the nucleotides linking G-quartet motifs. To assess the effects of a three-nucleotide TTT linker, the crystal structure of the DNA sequence d(G(4)T(3)G(4)) has been determined at 1.5 A resolution, together with that of the brominated analogue d(G(4)(Br)UTTG(4)) at 2.4 A resolution. Both sequences form bimolecular intermolecular G-quadruplexes with lateral loops. d(G(4)(Br)UTTG(4)) crystallized in the monoclinic space group P2(1) with three quadruplex molecules in the asymmetric unit, two associating together as a head-to-head stacked dimer, and the third as a single head-to-tail dimer. The head-to-head dimers have two lateral loops on the same G-quadruplex face and form an eight-G-quartet stack, with a linear array of seven K(+) ions between the quartets. d(G(4)T(3)G(4)) crystallized in the orthorhombic space group C222 and has a structure very similar to the head-to-tail dimer in the P2(1) unit cell. The sequence studied here is able to form several different folds; however, all four quadruplexes in the two structures have lateral loops, in contrast to the diagonal loops reported for the analogous quadruplex with T(4) loops. A total of seven independent T(3) loops were observed in the two structures. These can be classified into two discrete conformational classes, suggesting that these represent preferred loop conformations that are independent of crystal-packing forces.     

Logical Models of Molecular Shapes and Their Families

A new discrete mathematical model of molecular shape is proposed, making use of the partition property of a representation of molecular shape. According to its geometrical and topological structure, a molecular surface can be partitioned into unbounded two-dimensional subsets (domains) and some common subsets of closures of two or more domains. The sets of these domains as a base of a finite topology, containing the Boolean n-cube as a lower Boolean sub-lattice of this topology, defines the domain of the proposed logical model. A logical function can be obtained that reflects the properties of the topological domains as well as the interrelations on the set of domains. Based on classical or quantum-chemical representations of molecular shape, these models allow one the implementation of methods of logical diagnostics in chemistry, and the definition of a metric on the set of molecular shape equivalence classes. The families of molecular shapes can be considered as sets of logical models. The proposed model is unified in the sense that the structures of differentiable and non-differentiable surfaces can be represented in the same mathematical framework. These logical models will also work for interpenetrations of the above types of surfaces.     

Key issues in the computational simulation of GPCR function: representation of loop domains

Some key concerns raised by molecular modeling and computational simulation of functional mechanisms for membrane proteins are discussed and illustrated for members of the family of G protein coupled receptors (GPCRs). Of particular importance are issues related to the modeling and computational treatment of loop regions. These are demonstrated here with results from different levels of computational simulations applied to the structures of rhodopsin and a model of the 5-HT2A serotonin receptor, 5-HT2AR. First, comparative Molecular Dynamics (MD) simulations are reported for rhodopsin in vacuum and embedded in an explicit representation of the membrane and water environment. It is shown that in spite of a partial accounting of solvent screening effects by neutralization of charged side chains, vacuum MD simulations can lead to severe distortions of the loop structures. The primary source of the distortion appears to be formation of artifactual H-bonds, as has been repeatedly observed in vacuum simulations. To address such shortcomings, a recently proposed approach that has been developed for calculating the structure of segments that connect elements of secondary structure with known coordinates, is applied to 5-HT2AR to obtain an initial representation of the loops connecting the transmembrane (TM) helices. The approach consists of a simulated annealing combined with biased scaled collective variables Monte Carlo technique, and is applied to loops connecting the TM segments on both the extra-cellular and the cytoplasmic sides of the receptor. Although this initial calculation treats the loops as independent structural entities, the final structure exhibits a number of interloop interactions that may have functional significance. Finally, it is shown here that in the case where a given loop from two different GPCRs (here rhodopsin and 5-HT2AR) has approximately the same length and some degree of sequence identity, the fold adopted by the loops can be similar. Thus, in such special cases homology modeling might be used to obtain initial structures of these loops. Notably, however, all other loops in these two receptors appear to be very different in sequence and structure, so that their conformations can be found reliably only by ab initio, energy based methods and not by homology modeling.     

The interaction of highly helical structural mutants with the NOP receptor discloses the role of the address domain of nociceptin/orphanin FQ

Nociceptin is a heptadecapeptide whose sequence is similar to that of Dynorphin A, sharing a message domain characterized by two glycines and two aromatic residues, and a highly basic C-terminal address domain but, in spite of these similarities, displays no opioid activity. Establishing the relative importance of the message and address domains of nociceptin has so far been hampered by its extreme conformational flexibility. Here we show that mutants of this peptide, designed to increase the helical content in the address domain, can be employed to explain the mode of interaction with the NOP receptor. Nociceptin analogues in which Ala residues are substituted with aminoisobutyric acid (Aib) show a substantial increment of activity in their interaction with the NOP receptor. The increment of biological activity was attributed to the well-documented ability of Aib to induce helicity. Here we have verified this working hypothesis by a conformational investigation extended to new analogues in which the role of Aib is taken up by Leu. The NMR conformational analysis confirms that all Ala/Aib peptides as well as [Leu(7,11)]-N/OFQ-amide and [Leu(11,15)]-N/OFQ-amide mutants (N/OFQ=nociceptin/orphanin FQ) have comparable helix content in helix-promoting media. We show that the helical address domain of nociceptin can place key basic residues at an optimal distance from complementary acidic groups of the EL(2) loop of the receptor. Our structural data are used to rationalize pharmacological data which show that although [Leu(11,15)]-N/OFQ-amide has an activity comparable to those of Ala/Aib peptides, [Leu(7,11)]-N/OFQ-amide is less active than N/OFQ-amide. We hypothesize that bulky residues cannot be hosted in or near the hinge region (Thr(5)-Gly(6)-Ala(7)) without severe steric clash with the receptor. This hypothesis is also consistent with previous data on this hinge region obtained by systematic substitution of Thr, Gly, and Ala with Pro.