The Genetic Algorithm and the Conformational Search of Polypeptides and Proteins

Abstract

The genetic algorithm is a technique of function optimization derived from the principles of evolutionary theory. We have adapted it to perform conformational search on polypeptides and proteins. The algorithm was first tested on several small polypeptides and the 46 amino acid protein crambin under the AMBER potential energy function. The probable global minimum conformations of the polypeptides were located 90% of the time and a non-native conformation of crambin was located that was 150kcal/mol lower in potential energy than the minimized crystal structure conformation. Next, we used a knowledge-based potential function to predict the structures of melittin, pancreatic polypeptide, and crambin. A 2.31 Å ΔRMS conformation of melittin and a 5.33 Å ΔRMS conformation of pancreatic polypeptide were located by genetic algorithm-based conformational search under the knowledge-based potential function. Although the ΔRMS of pancreatic polypeptide was somewhat high, most of the secondary structure was correct. The secondary structure of crambin was predicted correctly, but the potential failed to promote packing interactions. Finally, we tested the packing aspects of our potential function by attempting to predict the tertiary structure of cytochrome b ₅₆₂ given correct secondary structure as a constraint. The final predicted conformation of cytochrome b ₅₆₂ was an almost completely extended continuous helix which indicated that the knowledge-based potential was useless for tertiary structure prediction. This work serves as a warning against testing potential functions designed for tertiary structure prediction on small proteins.     

The design of a peptide sequence to inhibit HIV replication: a search algorithm combining Monte Carlo and self-consistent mean field techniques

We developed a search algorithm combining Monte Carlo (MC) and self-consistent mean field techniques to evolve a peptide sequence that has good binding capability to the anticodon stem and loop (ASL) of human lysine tRNA species, tRNA^Lys3, with the ultimate purpose of breaking the replication cycle of human immunodeficiency virus-1. The starting point is the 15-amino-acid sequence, RVTHHAFLGAHRTVG, found experimentally by Agris and co-workers to bind selectively to hypermodified tRNA^Lys3. The peptide backbone conformation is determined via atomistic simulation of the peptide-ASL^Lys3 complex and then held fixed throughout the search. The proportion of amino acids of various types (hydrophobic, polar, charged, etc.) is varied to mimic different peptide hydration properties. Three different sets of hydration properties were examined in the search algorithm to see how this affects evolution to the best-binding peptide sequences. Certain amino acids are commonly found at fixed sites for all three hydration states, some necessary for binding affinity and some necessary for binding specificity. Analysis of the binding structure and the various contributions to the binding energy shows that: 1) two hydrophilic residues (asparagine at site 11 and the cysteine at site 12) “recognize” the ASL^Lys3 due to the VDW energy, and thereby contribute to its binding specificity and 2) the positively charged arginines at sites 4 and 13 preferentially attract the negatively charged sugar rings and the phosphate linkages, and thereby contribute to the binding affinity.     

Towards Precision Medicine: Advances in Computational Approaches for the Analysis of Human Variants

Variations and similarities in our individual genomes are part of our history, our heritage, and our identity. Some human genomic variants are associated with common traits such as hair and eye color, while others are associated with susceptibility to disease or response to drug treatment. Identifying the human variations producing clinically relevant phenotypic changes is critical for providing accurate and personalized diagnosis, prognosis, and treatment for diseases. Furthermore, a better understanding of the molecular underpinning of disease can lead to development of new drug targets for precision medicine. Several resources have been designed for collecting and storing human genomic variations in highly structured, easily accessible databases. Unfortunately, a vast amount of information about these genetic variants and their functional and phenotypic associations is currently buried in the literature, only accessible by manual curation or sophisticated text text-mining technology to extract the relevant information. In addition, the low cost of sequencing technologies coupled with increasing computational power has enabled the development of numerous computational methodologies to predict the pathogenicity of human variants. This review provides a detailed comparison of current human variant resources, including HGMD, OMIM, ClinVar, and UniProt/Swiss-Prot, followed by an overview of the computational methods and techniques used to leverage the available data to predict novel deleterious variants. We expect these resources and tools to become the foundation for understanding the molecular details of genomic variants leading to disease, which in turn will enable the promise of precision medicine.     

Methylerythritol phosphate pathway to isoprenoids: Kinetic modeling and in silico enzyme inhibitions in Plasmodium falciparum

The methylerythritol phosphate (MEP) pathway of Plasmodium falciparum (P. falciparum) has become an attractive target for anti-malarial drug discovery. This study describes a kinetic model of this pathway, its use in validating 1-deoxy-d-xylulose 5-phosphate reductoisomerase (DXR) as drug target from the systemic perspective, and additional target identification, using metabolic control analysis and in silico inhibition studies. In addition to DXR, 1-deoxy-d-xylulose 5-phosphate synthase (DXS) can be targeted because it is the first enzyme of the pathway and has the highest flux control coefficient followed by that of DXR. In silico inhibition of both enzymes caused large decrement in the pathway flux. An added advantage of targeting DXS is its influence on vitamin B1 and B6 biosynthesis. Two more potential targets, 2-C-methyl-d-erythritol 2,4-cyclodiphosphate synthase and 1-hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate synthase, were also identified. Their inhibition caused large accumulation of their substrates causing instability of the system.     

The relationship between visual working memory and attention: retention of precise colour information in the absence of effects on perceptual selection

We examined the conditions under which a feature value in visual working memory (VWM) recruits visual attention to matching stimuli. Previous work has suggested that VWM supports two qualitatively different states of representation: an active state that interacts with perceptual selection and a passive (or accessory) state that does not. An alternative hypothesis is that VWM supports a single form of representation, with the precision of feature memory controlling whether or not the representation interacts with perceptual selection. The results of three experiments supported the dual-state hypothesis. We established conditions under which participants retained a relatively precise representation of a parcticular colour. If the colour was immediately task relevant, it reliably recruited attention to matching stimuli. However, if the colour was not immediately task relevant, it failed to interact with perceptual selection. Feature maintenance in VWM is not necessarily equivalent with feature-based attentional selection.     

Development of a pharmacophore model for the catecholamine release-inhibitory peptide catestatin: Virtual screening and functional testing identify novel small molecule therapeutics of hypertension

The endogenous catecholamine release-inhibitory peptide catestatin (CST) regulates events leading to hypertension and cardiovascular disease. Earlier we studied the structure of CST by NMR, molecular modeling, and amino acid scanning mutagenesis. That structure has now been exploited for elucidation of interface pharmacophores that mediate binding of CST to its target, with consequent secretory inhibition. Designed pharmacophore models allowed screening of 3D structural domains. Selected compounds were tested on both cultured catecholaminergic cells and an in vivo model of hypertension; in each case, the candidates showed substantial mimicry of native CST actions, with preserved or enhanced potency and specificity. The approach and compounds have thus enabled rational design of novel drug candidates for treatment of hypertension or autonomic dysfunction.     

Foraging success under uncertainty: search tradeoffs and optimal space use

Understanding the structural complexity and the main drivers of animal search behaviour is pivotal to foraging ecology. Yet, the role of uncertainty as a generative mechanism of movement patterns is poorly understood. Novel insights from search theory suggest that organisms should collect and assess new information from the environment by producing complex exploratory strategies. Based on an extension of the first passage time theory, and using simple equations and simulations, we unveil the elementary heuristics behind search behaviour. In particular, we show that normal diffusion is not enough for determining optimal exploratory behaviour but anomalous diffusion is required. Searching organisms go through two critical sequential phases (approach and detection) and experience fundamental search tradeoffs that may limit their encounter rates. Using experimental data, we show that biological search includes elements not fully considered in contemporary physical search theory. In particular, the need to consider search movement as a non‐stationary process that brings the organism from one informational state to another. For example, the transition from remaining in an area to departing from it may occur through an exploratory state where cognitive search is challenged. Therefore, a more comprehensive view of foraging ecology requires including current perspectives about movement under uncertainty.     

Investigating core genetic-and-epigenetic cell cycle networks for stemness and carcinogenic mechanisms,and cancer drug design using big database mining and genome-wide next-generation sequencing data

Recent studies have demonstrated that cell cycle plays a central role in development and carcinogenesis. Thus, the use of big databases and genome-wide high-throughput data to unravel the genetic and epigenetic mechanisms underlying cell cycle progression in stem cells and cancer cells is a matter of considerable interest.

Real genetic-and-epigenetic cell cycle networks (GECNs) of embryonic stem cells (ESCs) and HeLa cancer cells were constructed by applying system modeling, system identification, and big database mining to genome-wide next-generation sequencing data. Real GECNs were then reduced to core GECNs of HeLa cells and ESCs by applying principal genome-wide network projection. In this study, we investigated potential carcinogenic and stemness mechanisms for systems cancer drug design by identifying common core and specific GECNs between HeLa cells and ESCs. Integrating drug database information with the specific GECNs of HeLa cells could lead to identification of multiple drugs for cervical cancer treatment with minimal side-effects on the genes in the common core. We found that dysregulation of miR-29C, miR-34A, miR-98, and miR-215; and methylation of ANKRD1, ARID5B, CDCA2, PIF1, STAMBPL1, TROAP, ZNF165, and HIST1H2AJ in HeLa cells could result in cell proliferation and anti-apoptosis through NFκB, TGF-β, and PI3K pathways. We also identified 3 drugs, methotrexate, quercetin, and mimosine, which repressed the activated cell cycle genes, ARID5B, STK17B, and CCL2, in HeLa cells with minimal side-effects.