Novel Computational Protein Design Algorithms with Applications to Cystic Fibrosis and HIV

Proteins are essential components of cells and are crucial for catalyzing reactions, signaling, recognition, motility, recycling, and structural stability. This diversity of function suggests that nature is only scratching the surface of protein functional space. Protein function is determined by structure, which in turn is determined predominantly by amino acid sequence. Protein design aims to explore protein sequence and conformational space to design novel proteins with new or improved function. The vast number of possible protein sequences makes exploring the space a challenging problem.

Computational structure-based protein design (CSPD) allows for the rational design of proteins. Because of the large search space, CSPD methods must balance search accuracy and modeling simplifications. We have developed algorithms that allow for the accurate and efficient search of protein conformational space. Specifically, we focus on algorithms that maintain provability, account for protein flexibility, and use ensemble-based rankings. We present several novel algorithms for incorporating improved flexibility into CSPD with continuous rotamers. We applied these algorithms to two biomedically important design problems. We designed peptide inhibitors of the cystic fibrosis agonist CAL that were able to restore function of the vital cystic fibrosis protein CFTR. We also designed improved HIV antibodies and nanobodies to combat HIV infections.

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey.

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson-Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1'-C1' distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1'-C1' distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5'-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.

Zone based systems design framework for the realisation of efficient block cipher based message authentication code algorithms

A zone based systems design framework is described and utilised in the implementation of a message authentication code (MAC) algorithm based on symmetric key block ciphers. The resulting block cipher based MAC algorithm may be used to provide assurance of the authenticity and, hence, the integrity of binary data. Using software simulation to benchmark against the de facto cipher block chaining MAC (CBC-MAC) variant used in the TinySec security protocol for wireless sensor networks and the NIST cipher block chaining MAC standard, CMAC; we show that our zone based systems design framework can lead to block cipher based MAC constructs that point to improvements in message processing efficiency, processing throughput and processing latency.

Sequence and Structure Based Protein Folding Studies With Implications

As the expression of the genetic blueprint, proteins are at the heart of all biological systems. The ever increasing set of available protein structures has taught us that diversity is the hallmark of their architecture, a fundamental characteristic that enables them to perform the vast array of functionality upon which all of life depends. This diversity, however, is central to one of the most challenging problems in molecular biology: how does a folding polypeptide chain navigate its way through all of the myriad of possible conformations to find its own particular biologically active form? With few overarching structural principles to draw upon that can be applied to all protein architecture, the search for a solution to the protein folding problem has yet to produce an algorithm that can explain and duplicate this fundamental biological process. In this thesis, we take a two-pronged approach for investigating the protein folding process. Our initial statistical studies of the distributions of hydrophobic and hydrophilic residues within α-helices and β-sheets suggest (i) that hydrophobicity plays a critical role in helix and sheet formation; and (ii) that the nucleation of these motifs may result in largely unidirectional growth. Most tellingly, from an examination of the amino acids found in the smallest β-sheets, we do not find any evidence of a β-nucleating code in the primary protein sequence. Complementing these statistical analyses, we have analyzed the structural environments of several ever-widening aspects of protein topology. Our examination of the gaps between strands in the smallest β-sheets reveals a common organizational principle underlying β-formation involving strands separated by large sequential gaps: with very few exceptions, these large gaps fold into single, compact structural modules, bringing the β-strands that are otherwise far apart in the sequence close together in space. We conclude, therefore, that β-nucleation in the smallest sheets results from the co-location of two strands that are either local in sequence, or local in space following prior folding events. A second study of larger β-sheets both corroborates and extends these findings: virtually all large sequential gaps between pairs of β-strands organize themselves into an hierarchical arrangement, creating a bread-crumb model of go-and-come-back structural organization that ultimately juxtaposes two strands of a parental β-structure that are far apart in the sequence in close spatial proximity. In a final study, we have formalized this go-and-come-back notion into the concept of anti-parallel double-strandedness (DS), and measure this property across protein architecture in general. With over 90% of all residues in a large, non-redundant set of protein structures classified as DS, we conclude that DS is a unifying structural principle that underpins all globular proteins. We postulate, moreover, that this one simple principle, anti-parallel double-strandedness, unites protein structure, protein folding and protein evolution.

Coronavirus Main Proteinase (3CLpro) Structure: Basis for Design of Anti-SARS Drugs

A novel coronavirus has been identified as the causative agent of severe acute respiratory syndrome (SARS). The viral main proteinase (Mpro, also called 3CLpro), which controls the activities of the coronavirus replication complex, is an attractive target for therapy. We determined crystal structures for human coronavirus (strain 229E) Mpro and for an inhibitor complex of porcine coronavirus [transmissible gastroenteritis virus (TGEV)] Mpro, and we constructed a homology model for SARS coronavirus (SARS-CoV) Mpro. The structures reveal a remarkable degree of conservation of the substrate-binding sites, which is further supported by recombinant SARS-CoV Mpro-mediated cleavage of a TGEV Mpro substrate. Molecular modeling suggests that available rhinovirus 3Cpro inhibitors may be modified to make them useful for treating SARS.

Microneedle-based drug delivery systems: Microfabrication, drug delivery, and safety

Many promising therapeutic agents are limited by their inability to reach the systemic circulation, due to the excellent barrier properties of biological membranes, such as the stratum corneum (SC) of the skin or the sclera/cornea of the eye and others. The outermost layer of the skin, the SC, is the principal barrier to topically-applied medications. The intact SC thus provides the main barrier to exogenous substances, including drugs. Only drugs with very specific physicochemical properties (molecular weight <500 Da, adequate lipophilicity, and low melting point) can be successfully administered transdermally. Transdermal delivery of hydrophilic drugs and macromolecular agents of interest, including peptides, DNA, and small interfering RNA is problematic. Therefore, facilitation of drug penetration through the SC may involve by-pass or reversible disruption of SC molecular architecture. Microneedles (MNs), when used to puncture skin, will by-pass the SC and create transient aqueous transport pathways of micron dimensions and enhance the transdermal permeability. These micropores are orders of magnitude larger than molecular dimensions, and, therefore, should readily permit the transport of hydrophilic macromolecules. Various strategies have been employed by many research groups and pharmaceutical companies worldwide, for the fabrication of MNs. This review details various types of MNs, fabrication methods and, importantly, investigations of clinical safety of MN.

Structure-Based Mutational Analysis of eIF4E in Relation to sbm1 Resistance to Pea Seed-Borne Mosaic Virus in Pea

Background: Pea encodes eukaryotic translation initiation factor eIF4E (eIF4E(S)), which supports the multiplication of Pea seed-borne mosaic virus (PSbMV). In common with hosts for other potyviruses, some pea lines contain a recessive allele (sbm1) encoding a mutant eIF4E (eIF4E(R)) that fails to interact functionally with the PSbMV avirulence protein, VPg, giving genetic resistance to infection.

Non-genetic mechanisms communicating antibiotic resistance:rethinking strategies for antimicrobial drug design

Introduction: Infections by multidrug-resistant bacteria are of great concern worldwide. In many cases, resistance is not due to the presence of specific antibiotic-modifying enzymes, but rather associated with a general impermeability of the bacterial cell envelope. The molecular bases of this intrinsic resistance are not completely understood. Moreover, horizontal gene transfers cannot solely explain the spread of intrinsic resistance among bacterial strains. Areas covered: This review focuses on the increased intrinsic antibiotic resistance mediated by small molecules. These small molecules can either be secreted from bacterial cells of the same or different species (e.g., indole, polyamines, ammonia, and the Pseudomonas quinolone signal) or be present in the bacterial cell milieu, whether in the environment, such as indole acetic acid and other plant hormones, or in human tissues and body fluids, such as polyamines. These molecules are metabolic byproducts that act as infochemicals and modulate bacterial responses toward antibiotics leading to increasing or decreasing resistance levels. Expert opinion: The non-genetic mechanisms of antibiotic response modulation and communication discussed in this review should reorient our thinking of the mechanisms of intrinsic resistance to antibiotics and its spread across bacterial cell populations. The identification of chemical signals mediating increased intrinsic antibiotic resistance will expose novel critical targets for the development of new antimicrobial strategies.

GPCR structure, function, drug discovery and crystallography: report from Academia-Industry International Conference (UK Royal Society) Chicheley Hall, 1-2 September 2014

G-protein coupled receptors (GPCRs) are the targets of over half of all prescribed drugs today. The UniProt database has records for about 800 proteins classified as GPCRs, but drugs have only been developed against 50 of these. Thus, there is huge potential in terms of the number of targets for new therapies to be designed. Several breakthroughs in GPCRs biased pharmacology, structural biology, modelling and scoring have resulted in a resurgence of interest in GPCRs as drug targets. Therefore, an international conference, sponsored by the Royal Society, with world-renowned researchers from industry and academia was recently held to discuss recent progress and highlight key areas of future research needed to accelerate GPCR drug discovery. Several key points emerged. Firstly, structures for all three major classes of GPCRs have now been solved and there is increasing coverage across the GPCR phylogenetic tree. This is likely to be substantially enhanced with data from x-ray free electron sources as they move beyond proof of concept. Secondly, the concept of biased signalling or functional selectivity is likely to be prevalent in many GPCRs, and this presents exciting new opportunities for selectivity and the control of side effects, especially when combined with increasing data regarding allosteric modulation. Thirdly, there will almost certainly be some GPCRs that will remain difficult targets because they exhibit complex ligand dependencies and have many metastable states rendering them difficult to resolve by crystallographic methods. Subtle effects within the packing of the transmembrane helices are likely to mask and contribute to this aspect, which may play a role in species dependent behaviour. This is particularly important because it has ramifications for how we interpret pre-clinical data. In summary, collaborative efforts between industry and academia have delivered significant progress in terms of structure and understanding of GPCRs and will be essential for resolving problems associated with the more difficult targets in the future.