Computational strategies to include protein flexibility in Ligand Docking and Virtual Screening

The dynamic character of proteins strongly influences biomolecular recognition mechanisms. With the development of the main models of ligand recognition (lock-and-key, induced fit, conformational selection theories), the role of protein plasticity has become increasingly relevant. In particular, major structural changes concerning large deviations of protein backbones, and slight movements such as side chain rotations are now carefully considered in drug discovery and development. It is of great interest to identify multiple protein conformations as preliminary step in a screening campaign. Protein flexibility has been widely investigated, in terms of both local and global motions, in two diverse biological systems. On one side, Replica Exchange Molecular Dynamics has been exploited as enhanced sampling method to collect multiple conformations of Lactate Dehydrogenase A (LDHA), an emerging anticancer target. The aim of this project was the development of an Ensemble-based Virtual Screening protocol, in order to find novel potent inhibitors. On the other side, a preliminary study concerning the local flexibility of Opioid Receptors has been carried out through ALiBERO approach, an iterative method based on Elastic Network-Normal Mode Analysis and Monte Carlo sampling. Comparison of the Virtual Screening performances by using single or multiple conformations confirmed that the inclusion of protein flexibility in screening protocols has a positive effect on the probability to early recognize novel or known active compounds.

High protein flexibility and reduced hydration water dynamics are key pressure adaptive strategies in prokaryotes

Water and protein dynamics on a nanometer scale were measured by quasi-elastic neutron scattering in the piezophile archaeon Thermococcus barophilus and the closely related pressure-sensitive Thermococcus kodakarensis, at 0.1 and 40 MPa. We show that cells of the pressure sensitive organism exhibit higher intrinsic stability. Both the hydration water dynamics and the fast protein and lipid dynamics are reduced under pressure. In contrast, the proteome of T. barophilus is more pressure sensitive than that of T. kodakarensis. The diffusion coefficient of hydration water is reduced, while the fast protein and lipid dynamics are slightly enhanced with increasing pressure. These findings show that the coupling between hydration water and cellular constituents might not be simply a master-slave relationship. We propose that the high flexibility of the T. barophilus proteome associated with reduced hydration water may be the keys to the molecular adaptation of the cells to high hydrostatic pressure.

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.

Galactokinase promiscuity: a question of flexibility?

Galactokinase catalyses the first committed step of the Leloir pathway, i.e. the ATP-dependent phosphorylation of α-D-galactose at C1-OH. Reduced galactokinase activity results in the inherited metabolic disease type II galactosaemia. However, inhibition of galactokinase is considered a viable approach to treating more severe forms of galactosaemia (types I and III). Considerable progress has been made in the identification of high affinity, selective inhibitors. Although the structure of galactokinase from a variety of species is known, its catalytic mechanism remains uncertain. Although the bulk of evidence suggests that the reaction proceeds via an active site base mechanism, some experimental and theoretical studies contradict this. The enzyme has potential as a biocatalyst in the production of sugar 1-phosphates. This potential is limited by its high specificity. A variety of approaches have been taken to identify galactokinase variants which are more promiscuous. These have broadened galactokinase's specificity to include a wide range of D- and L-sugars. Initial studies suggest that some of these alterations result in increased flexibility at the active site. It is suggested that modulation of protein flexibility is at least as important as structural modifications in determining the success or failure of enzyme engineering.

Flexibility and inhibitor binding in Cdc25 phosphatases

Cdc25 phosphatases involved in cell cycle checkpoints are now active targets for the development of anti-cancer therapies. Rational drug design would certainly benefit from detailed structural information for Cdc25s. However, only apo- or sulfate-bound crystal structures of the Cdc25 catalytic domain have been described so far. Together with previously available crystalographic data, results from molecular dynamics simulations, bioinformatic analysis, and computer-generated conformational ensembles shown here indicate that the last 30-40 residues in the C-terminus of Cdc25B are partially unfolded or disordered in solution. The effect of C-terminal flexibility upon binding of two potent small molecule inhibitors to Cdc25B is then analyzed by using three structural models with variable levels of flexibility, including an equilibrium distributed ensemble of Cdc25B backbone conformations. The three Cdc25B structural models are used in combination with flexible docking, clustering, and calculation of binding free energies by the linear interaction energy approximation to construct and validate Cdc25B-inhibitor complexes. Two binding sites are identified on top and beside the Cdc25B active site. The diversity of interaction modes found increases with receptor flexibility. Backbone flexibility allows the formation of transient cavities or compact hydrophobic units on the surface of the stable, folded protein core that are unexposed or unavailable for ligand binding in rigid and densely packed crystal structures. The present results may help to speculate on the mechanisms of small molecule complexation to partially unfolded or locally disordered proteins.

Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding.

We propose a general mean field model of ligand-protein interactions to determine the thermodynamic equilibrium of a system at finite temperature. The method is employed in structural assessments of two human immuno-deficiency virus type 1 protease complexes where the gross effects of protein flexibility are incorporated by utilizing a data base of crystal structures. Analysis of the energy spectra for these complexes has revealed that structural and thermo-dynamic aspects of molecular recognition can be rationalized on the basis of the extent of frustration in the binding energy landscape. In particular, the relationship between receptor-specific binding of these ligands to human immunodeficiency virus type 1 protease and a minimal frustration principle is analyzed.

Prediction of Phosphorylation Sites from Protein Sequences

Phosphorylation is amongst the most crucial and well-studied post-translational modifications. It is involved in multiple cellular processes which makes phosphorylation prediction vital for understanding protein functions. However, wet-lab techniques are labour and time intensive. Thus, computational tools are required for efficiency. This project aims to provide a novel way to predict phosphorylation sites from protein sequences by adding flexibility and Sezerman Grouping amino acid similarity measure to previous methods, as discovering new protein sequences happens at a greater rate than determining protein structures. The predictor – NOPAY - relies on Support Vector Machines (SVMs) for classification. The features include amino acid encoding, amino acid grouping, predicted secondary structure, predicted protein disorder, predicted protein flexibility, solvent accessibility, hydrophobicity and volume. As a result, we have managed to improve phosphorylation prediction accuracy for Homo sapiens by 3% and 6.1% for Mus musculus. Sensitivity at 99% specificity was also increased by 6% for Homo sapiens and for Mus musculus by 5% on independent test sets. In this study, we have managed to increase phosphorylation prediction accuracy for Homo sapiens and Mus musculus. When there is enough data, future versions of the software may also be able to predict other organisms.

Simulations numériques de la dynamique des protéines : translation de ligands, flexibilité et dynamique des boucles

La flexibilité est une caractéristique intrinsèque des protéines qui doivent, dès le mo- ment de leur synthèse, passer d’un état de chaîne linéaire à un état de structure tridimen- sionnelle repliée et enzymatiquement active. Certaines protéines restent flexibles une fois repliées et subissent des changements de conformation de grande amplitude lors de leur cycle enzymatique. D’autres contiennent des segments si flexibles que leur structure ne peut être résolue par des méthodes expérimentales. Dans cette thèse, nous présentons notre application de méthodes in silico d’analyse de la flexibilité des protéines : • À l’aide des méthodes de dynamique moléculaire dirigée et d’échantillonnage pa- rapluie, nous avons caractérisé les trajectoires de liaison de l’inhibiteur Z-pro- prolinal à la protéine Prolyl oligopeptidase et identifié la trajectoire la plus pro- bable. Nos simulations ont aussi identifié un mode probable de recrutement des ligands utilisant une boucle flexible de 19 acides aminés à l’interface des deux domaines de la protéine. • En utilisant les méthodes de dynamique moléculaire traditionnelle et dirigée, nous avons examiné la stabilité de la protéine SAV1866 dans sa forme fermée insérée dans une membrane lipidique et étudié un des modes d’ouverture possibles par la séparation de ses domaines liant le nucléotide. • Nous avons adapté auproblème de la prédiction de la structure des longues boucles flexibles la méthode d’activation et de relaxation ART-nouveau précédemment uti- lisée dans l’étude du repliement et de l’agrégation de protéines. Appliqué au replie- ment de boucles de 8 à 20 acides aminés, la méthode démontre une dépendance quadratique du temps d’exécution sur la longueur des boucles, rendant possible l’étude de boucles encore plus longues.

Identifying sequence regions undergoing conformational change via predicted continuum secondary structure

Motivation: Conformational flexibility is essential to the function of many proteins, e.g. catalytic activity. To assist efforts in determining and exploring the functional properties of a protein, it is desirable to automatically identify regions that are prone to undergo conformational changes. It was recently shown that a probabilistic predictor of continuum secondary structure is more accurate than categorical predictors for structurally ambivalent sequence regions, suggesting that such models are suited to characterize protein flexibility. Results: We develop a computational method for identifying regions that are prone to conformational change directly from the amino acid sequence. The method uses the entropy of the probabilistic output of an 8-class continuum secondary structure predictor. Results for 171 unique amino acid sequences with well-characterized variable structure (identified in the 'Macromolecular movements database') indicate that the method is highly sensitive at identifying flexible protein regions, but false positives remain a problem. The method can be used to explore conformational flexibility of proteins (including hypothetical or synthetic ones) whose structure is yet to be determined experimentally.

The 2.1 angstrom crystal structure of the far-red fluorescent protein HcRed: Inherent conformational flexibility of the chromophore

We have determined the crystal structure of HcRed, a far-red fluorescent protein isolated from Heteractis crispa, to 2.1 resolution. HcRed was observed to form a dimer, in contrast to the monomeric form of green fluorescent protein (GFP) or the tetrameric forms of the GFP-like proteins (eqFP611, Rtms5 and DsRed). Unlike the well-defined chromophore conformation observed in GFP and the GFP-like proteins, the HcRed chromophore was observed to be considerably mobile. Within the HcRed structure, the cyclic tripeptide chromophore, Glu64-Tyr65-Gly66, was observed to adopt both a cis coplanar and a tran. non-coplanar conformation. As a result of these two con formations, the hydroxyphenyl moiety of the chromophore makes distinct interactions within the interior of the b-can. These data together with a quantum chemical model of the chromophore, suggest the cis coplanar conformation to be consistent with the fluorescent properties of HcRed, and the trans non-coplanar conformation to be consistent with non-fluorescent properties of hcCP, the chromoprotein parent of HcRed. Moreover, within the GFP-like family, it appears that where conformational freedom is permissible then flexibility in the chromophore conformation is possible. 2005 Elsevier Ltd. All rights reserved.

Kalata B8, a novel antiviral circular protein, exhibits conformational flexibility in the cystine knot motif

The cyclotides are a family of circular proteins with a range of biological activities and potential pharmaceutical and agricultural applications. The biosynthetic mechanism of cyclization is unknown and the discovery of novel sequences may assist in achieving this goal. In the present study, we have isolated a new cyclotide from Oldenlandia affinis, kalata B8, which appears to be a hybrid of the two major subfamilies (Mobius and bracelet) of currently known cyclotides. We have determined the three-dimensional structure of kalata B8 and observed broadening of resonances directly involved in the cystine knot motif, suggesting flexibility in this region despite it being the core structural element of the cyclotides. The cystine knot motif is widespread throughout Nature and inherently stable, making this apparent flexibility a surprising result. Further-more, there appears to be isomerization of the peptide backbone at an Asp-Gly sequence in the region involved in the cyclization process. Interestingly, such isomerization has been previously characterized in related cyclic knottins from Momordica cochinchinensis that have no sequence similarity to kalata B8 apart from the six conserved cysteine residues and may result from a common mechanism of cyclization. Kalata B8 also provides insight into the structure-activity relationships of cyclotides as it displays anti-HIV activity but lacks haemolytic activity. The 'uncoupling' of these two activities has not previously been observed for the cyclotides and may be related to the unusual hydrophilic nature of the peptide.

Stress and developmental responses in Arabidopsis thaliana: regulation through the transcription factor-interacting protein RCD1

Plants constantly face adverse environmental conditions, such as drought or extreme temperatures that threaten their survival. They demonstrate astonishing metabolic flexibility in overcoming these challenges and one of the key responses to stresses is changes in gene expression leading to alterations in cellular functions. This is brought about by an intricate network of transcription factors and associated regulatory proteins. Protein-protein interactions and post-translational modifications are important steps in this control system along with carefully regulated degradation of signaling proteins. This work concentrates on the RADICAL-INDUCED CELL DEATH1 (RCD1) protein which is an important regulator of abiotic stress-related and developmental responses in Arabidopsis thaliana. Plants lacking this protein function display pleiotropic phenotypes including sensitivity to apoplastic reactive oxygen species (ROS) and salt, ultraviolet B (UV-B) and paraquat tolerance, early flowering and senescence. Additionally, the mutant plants overproduce nitric oxide, have alterations in their responses to several plant hormones and perturbations in gene expression profiles. The RCD1 gene is transcriptionally unresponsive to environmental signals and the regulation of the protein function is likely to happen post-translationally. RCD1 belongs to a small protein family and, together with its closest homolog SRO1, contains three distinguishable domains: In the N-terminus, there is a WWE domain followed by a poly(ADP-ribose) polymerase-like domain which, despite sequence conservation, does not seem to be functional. The C-terminus of RCD1 contains a novel domain called RST. It is present in RCD1-like proteins throughout the plant kingdom and is able to mediate physical interactions with multiple transcription factors. In conclusion, RCD1 is a key point of signal integration that links ROS-mediated cues to transcriptional regulation by yet unidentified means, which are likely to include post-translational mechanisms. The identification of RCD1-interacting transcription factors, most of whose functions are still unknown, opens new avenues for studies on plant stress as well as developmental responses.

Random network behaviour of protein structures

Geometric and structural constraints greatly restrict the selection of folds adapted by protein backbones, and yet, folded proteins show an astounding diversity in functionality. For structure to have any bearing on function, it is thus imperative that, apart from the protein backbone, other tunable degrees of freedom be accountable. Here, we focus on side-chain interactions, which non-covalently link amino acids in folded proteins to form a network structure. At a coarse-grained level, we show that the network conforms remarkably well to realizations of random graphs and displays associated percolation behavior. Thus, within the rigid framework of the protein backbone that restricts the structure space, the side-chain interactions exhibit an element of randomness, which account for the functional flexibility and diversity shown by proteins. However, at a finer level, the network exhibits deviations from these random graphs which, as we demonstrate for a few specific examples, reflect the intrinsic uniqueness in the structure and stability, and perhaps specificity in the functioning of biological proteins.

Structure of Mycobacterium tuberculosis RuvA, a protein involved in recombination

The process of recombinational repair is crucial for maintaining genomic integrity and generating biological diversity. In association with RuvB and RuvC, RuvA plays a central role in processing and resolving Holliday junctions, which are a critical intermediate in homologous recombination. Here, the cloning, purification and structure determination of the RuvA protein from Mycobacterium tuberculosis (MtRuvA) are reported. Analysis of the structure and comparison with other known RuvA proteins reveal an octameric state with conserved subunit-subunit interaction surfaces, indicating the requirement of octamer formation for biological activity. A detailed analysis of plasticity in the RuvA molecules has led to insights into the invariant and variable regions, thus providing a framework for understanding regional flexibility in various aspects of RuvA function.