Optimization of an Elastic Network Augmented Coarse Grained Model to Study CCMV Capsid Deformation

The major protective coat of most viruses is a highly symmetric protein capsid that forms spontaneously from many copies of identical proteins. Structural and mechanical properties of such capsids, as well as their self-assembly process, have been studied experimentally and theoretically, including modeling efforts by computer simulations on various scales. Atomistic models include specific details of local protein binding but are limited in system size and accessible time, while coarse grained (CG) models do get access to longer time and length scales but often lack the specific local interactions. Multi-scale models aim at bridging this gap by systematically connecting different levels of resolution. Here, a CG model for CCMV (Cowpea Chlorotic Mottle Virus), a virus with an icosahedral shell of 180 identical protein monomers, is developed, where parameters are derived from atomistic simulations of capsid protein dimers in aqueous solution. In particular, a new method is introduced to combine the MARTINI CG model with a supportive elastic network based on structural fluctuations of individual monomers. In the parametrization process, both network connectivity and strength are optimized. This elastic-network optimized CG model, which solely relies on atomistic data of small units (dimers), is able to correctly predict inter-protein conformational flexibility and properties of larger capsid fragments of 20 and more subunits. Furthermore, it is shown that this CG model reproduces experimental (Atomic Force Microscopy) indentation measurements of the entire viral capsid. Thus it is shown that one obvious goal for hierarchical modeling, namely predicting mechanical properties of larger protein complexes from models that are carefully parametrized on elastic properties of smaller units, is achievable.     

Crystal structure of a Vip3B family insecticidal protein reveals a new fold and a unique tetrameric assembly

Vegetatively expressed insecticidal proteins (VIPs) produced by Bacillus thuringiensis fall into several classes of which the third, VIP3, is known for their activity against several key Lepidopteran pests of commercial broad acre crops and because their mode of action does not overlap with that of crystalline insecticidal proteins. The details of the VIP3 structure and mode of action have remained obscure for the quarter century that has passed since their discovery. In the present article, we report the first crystal structure of a full‐length VIP3 protein. Crystallization of this target required multiple rounds of construct optimization and screening—over 200 individual sequences were expressed and tested. This protein adopts a novel global fold that combines domains with hitherto unreported topology and containing elements seemingly borrowed from carbohydrate‐binding domains, lectins, or from other insecticidal proteins.     

Conformational Changes and Ligand Recognition of Escherichia colid-Xylose Binding Protein Revealed

ATP binding cassette transport systems account for most import of necessary nutrients in bacteria. The periplasmic binding component (or an equivalent membrane-anchored protein) is critical to recognizing cognate ligand and directing it to the appropriate membrane permease. Here we report the X-ray structures of d-xylose binding protein from Escherichia coli in ligand-free open form, ligand-bound open form, and ligand-bound closed form at 2.15 Å, 2.2 Å, and 2.2 Å resolutions, respectively. The ligand-bound open form is the first such structure to be reported at high resolution; the combination of the three different forms from the same protein furthermore gives unprecedented details concerning the conformational changes involved in binding protein function. As is typical of the structural family, the protein has two similar globular domains, which are connected by a three-stranded hinge region. The open liganded structure shows that xylose binds first to the C-terminal domain, with only very small conformational changes resulting. After a 34° closing motion, additional interactions are formed with the N-terminal domain; changes in this domain are larger and serve to make the structure more ordered near the ligand. An analysis of the interactions suggests why xylose is the preferred ligand. Furthermore, a comparison with the most closely related proteins in the structural family shows that the conformational changes are distinct in each type of binding protein, which may have implications for how the individual proteins act in concert with their respective membrane permeases.     

A model of gradient interpretation based on morphogen binding

We propose a model in which pattern formation is controlled by several concentration gradients of “morphogens” and by allosteric proteins which bind them. In this model, each protein can bind up to two molecules of each morphogen and has an “active state” when one molecule of each morphogen is bound. The concentration of the active state of such a “morphogen binding protein” varies with position in a way that depends on the values given the binding constants. In a contour map of the active state concentration, the contours can have a variety of simple shapes.Simply-shaped regions of cell differentiation can be defined directly by concentration contours of a morphogen binding protein using a threshold-sensing mechanism. More complex shapes may be generated using several proteins and a “winner-take-all” rule according to which each protein specifies some particular sort of cell differentiation and the differentiation of cells in any position is governed by the protein with the highest active state concentration.We present an application of our model to the vertebrate limb skeleton; we use the “winner-take-all” mechanism and thirteen morphogen binding proteins, eleven of which specify cartilage formation. In this model we use one morphogen binding protein to specify the shaft of a typical long bone and one for each epiphysis. Our model is reasonably successful in imitating the in vivo positions and orientations of developing bones and in generating simple, plausible-looking articular surfaces.In addition to the morphogen-binding model we propose a mechanism which could transform morphogen-binding patterns into high-amplitude patterns capable of controlling the activity of structural genes. This “amplifying mechanism” can account for two previously unexplained features of limb skeletal development: the early formation of the diffusely-bounded “scleroblastema” in the limb bud and the center-to-edge gradations in cartilage formation rate which are later seen within individual chondrification foci.A simple modification of the morphogen-binding model provides an explanation for the general anatomical phenomenon of metamerism: The model can account for the formation of inexactly repeating patterns (such as the pattern of the vertebral column) and suggests a mechanism by which such patterns could (1) evolve from exactly repeating patterns, and (2) acquire, in further evolution, a high degree of specialization of the individual repeating units.The most promising approach for testing the morphogen-binding model would appear to involve experiments in which cytoplasm is transferred between cells at various stages of pattern development. Support for the model could also come from the discovery of certain kinds of hereditary limb defects.     

Conformational flexibility of the leucine binding protein examined by protein domain coarse-grained molecular dynamics

Periplasmic binding proteins are the initial receptors for the transport of various substrates over the inner membrane of gram-negative bacteria. The binding proteins are composed of two domains, and the substrate is entrapped between these domains. For several of the binding proteins it has been established that a closed-up conformation exists even without substrate present, suggesting a highly flexible apo-structure which would compete with the ligand-bound protein for the transporter interaction. For the leucine binding protein (LBP), structures of both open and closed conformations are known, but no closed-up structure without substrate has been reported. Here we present molecular dynamics simulations exploring the conformational flexibility of LBP. Coarse grained models based on the MARTINI force field are used to access the microsecond timescale. We show that a standard MARTINI model cannot maintain the structural stability of the protein whereas the ELNEDIN extension to MARTINI enables simulations showing a stable protein structure and nanosecond dynamics comparable to atomistic simulations, but does not allow the simulation of conformational flexibility. A modification to the MARTINI-ELNEDIN setup, referred to as domELNEDIN, is therefore presented. The domELNEDIN setup allows the protein domains to move independently and thus allows for the simulation of conformational changes. Microsecond domELNEDIN simulations starting from either the open or the closed conformations consistently show that also for LBP, the apo-structure is flexible and can exist in a closed form.

Comparative Structural Analysis of Human DEAD-Box RNA Helicases

DEAD-box RNA helicases play various, often critical, roles in all processes where RNAs are involved. Members of this family of proteins are linked to human disease, including cancer and viral infections. DEAD-box proteins contain two conserved domains that both contribute to RNA and ATP binding. Despite recent advances the molecular details of how these enzymes convert chemical energy into RNA remodeling is unknown. We present crystal structures of the isolated DEAD-domains of human DDX2A/eIF4A1, DDX2B/eIF4A2, DDX5, DDX10/DBP4, DDX18/myc-regulated DEAD-box protein, DDX20, DDX47, DDX52/ROK1, and DDX53/CAGE, and of the helicase domains of DDX25 and DDX41. Together with prior knowledge this enables a family-wide comparative structural analysis. We propose a general mechanism for opening of the RNA binding site. This analysis also provides insights into the diversity of DExD/H- proteins, with implications for understanding the functions of individual family members.     

Fluorescence studies of rat cellular retinol binding protein II produced in Escherichia coli: an analysis of four tryptophan substitution mutants.

Rat intestinal cellular retinol binding protein II (CRBP II) is an abundant 134-residue protein that binds all-trans-retinol which contains 4 tryptophans in positions 9, 89, 107, and 110. Our ability to express CRBP II in Escherichia coli and to construct individual tryptophan substitution mutants by site-directed mutagenesis has provided a useful model system for studying the fluorescence of a multi-tryptophan protein. Each of the four mutant proteins binds all-trans-retinol with high affinity, although their affinities are less than that of the wild-type protein. Steady-state and time-resolved fluorescence analyses of these proteins indicate that W107 is at the hydrophobic binding site, W110 is in a polar environment, and the remaining two tryptophans are in a hydrophobic environment. Time-resolved fluorescence study indicates that excited-state energy transfer occurs from the hydrophobic tryptophans to W110. The Stern-Volmer analysis with acrylamide of these proteins reveals that static quenching occurs in the W9F mutant protein while others do not. The fluorescence of rat intestinal fatty acid binding protein (I-FABP), a related protein of known X-ray structure, was also studied for comparison. The results of these findings, coupled with those derived from NMR studies and molecular graphics, suggest that CRBP II undergoes minor structural changes in all of the mutant proteins. Since these effects may be cumulative on the protein structure and function, any conclusions derived from higher mutants in this family of proteins must be treated with caution.     

Contribution of the tyrosines to the structure and function of the human U1A N-terminal RNA binding domain.

RNA binding domains (RBDs) are members of a large family of proteins that share minimal sequence conservation but adopt an alpha beta sandwich global fold. Defining the contributions of specific amino acids to RBD structure and RNA binding is critical to understanding the functions of these proteins. In these experiments with the human U1A N-terminal RNA binding domain (RBD1), the contributions from each of its four tyrosines to protein structure, stability, and RNA binding were measured. Each tyrosine was substituted with phenylalanine and one other selected residue, and the resulting proteins were characterized by chemical denaturation to measure their unfolding free energy, by binding free energies to the wild-type RNA hairpin, and by 19F NMR to probe for structural changes. Features of the protein identified in these experiments include a possible tyrosine/lysine contact in an alpha-helix, which may be an example of an energetically favorable aromatic/amino side chain interaction. One long loop of the protein, which shows unusual 15N backbone and tyrosine side-chain dynamics, is implicated in protein:protein association. The diverse interactions of the four tyrosine residues in the organization of RBD1 illustrate how each member of this family of proteins will have unique molecular details that contribute to function.     

Hydrophilicity Matching – A Potential Prerequisite for the Formation of Protein-Protein Complexes in the Cell

A binding event between two proteins typically consists of a diffusional search of binding partners for one another, followed by a specific recognition of the compatible binding sites resulting in the formation of the complex. However, it is unclear how binding partners find each other in the context of the crowded, constantly fluctuating, and interaction-rich cellular environment. Here we examine the non-specific component of protein-protein interactions, which refers to those physicochemical properties of the binding partners that are independent of the exact details of their binding sites, but which can affect their localization or diffusional search for one another. We show that, for a large set of high-resolution experimental 3D structures of binary, transient protein complexes taken from the DOCKGROUND database, the binding partners display a surprising, statistically significant similarity in terms of their total hydration free energies normalized by a size-dependent variable. We hypothesize that colocalization of binding partners, even within individual cellular compartments such as the cytoplasm, may be influenced by their relative hydrophilicity, potentially in response to local hydrophilic gradients.     

Structural basis for antagonism by suramin of heparin binding to vaccinia complement protein

Suramin is a competitive inhibitor of heparin binding to many proteins, including viral envelope proteins, protein tyrosine phosphatases, and fibroblast growth factors (FGFs). It has been clinically evaluated as a potential therapeutic in treatment of cancers caused by unregulated angiogenesis, triggered by FGFs. Although it has shown clinical promise in treatment of several cancers, suramin has many undesirable side effects. There is currently no experimental structure that reveals the molecular interactions responsible for suramin inhibition of heparin binding, which could be of potential use in structure-assisted design of improved analogues of suramin. We report the structure of suramin, in complex with the heparin-binding site of vaccinia virus complement control protein (VCP), which interacts with heparin in a geometrically similar manner to many FGFs. The larger than anticipated flexibility of suramin manifested in this structure, and other details of VCP-suramin interactions, might provide useful structural information for interpreting interactions of suramin with many proteins.     

Functional roles of S100 proteins, calcium-binding proteins of the EF-hand type.

A multigenic family of Ca2+-binding proteins of the EF-hand type known as S100 comprises 19 members that are differentially expressed in a large number of cell types. Members of this protein family have been implicated in the Ca2+-dependent (and, in some cases, Zn2+- or Cu2+-dependent) regulation of a variety of intracellular activities such as protein phosphorylation, enzyme activities, cell proliferation (including neoplastic transformation) and differentiation, the dynamics of cytoskeleton constituents, the structural organization of membranes, intracellular Ca2+ homeostasis, inflammation, and in protection from oxidative cell damage. Some S100 members are released or secreted into the extracellular space and exert trophic or toxic effects depending on their concentration, act as chemoattractants for leukocytes, modulate cell proliferation, or regulate macrophage activation. Structural data suggest that many S100 members exist within cells as dimers in which the two monomers are related by a two-fold axis of rotation and that Ca2+ binding induces in individual monomers the exposure of a binding surface with which S100 dimers are believed to interact with their target proteins. Thus, any S100 dimer is suggested to expose two binding surfaces on opposite sides, which renders homodimeric S100 proteins ideal for crossbridging two homologous or heterologous target proteins. Although in some cases different S100 proteins share their target proteins, in most cases a high degree of target specificity has been described, suggesting that individual S100 members might be implicated in the regulation of specific activities. On the other hand, the relatively large number of target proteins identified for a single S100 protein might depend on the specific role played by the individual regions that in an S100 molecule contribute to the formation of the binding surface. The pleiotropic roles played by S100 members, the identification of S100 target proteins, the analysis of functional correlates of S100-target protein interactions, and the elucidation of the three-dimensional structure of some S100 members have greatly increased the interest in S100 proteins and our knowledge of S100 protein biology in the last few years. S100 proteins probably are an example of calcium-modulated, regulatory proteins that intervene in the fine tuning of a relatively large number of specific intracellular and (in the case of some members) extracellular activities. Systems, including knock-out animal models, should be now used with the aim of defining the correspondence between the in vitro regulatory role(s) attributed to individual members of this protein family and the in vivo function(s) of each S100 protein.     

Differential single-stranded DNA binding properties of the paralogous SsbA and SsbB proteins from Streptococcus pneumoniae

The naturally transformable Gram-positive bacterium Streptococcus pneumoniae has two single-stranded DNA-binding (SSB) proteins, designated SsbA and SsbB. The SsbA protein is similar in size to the well characterized SSB protein from Escherichia coli (SsbEc). The SsbB protein, in contrast, is a smaller protein that is specifically induced during natural transformation and has no counterpart in E. coli. In this report, the single-stranded DNA (ssDNA) binding properties of the SsbA and SsbB proteins were examined and compared with those of the SsbEc protein. The ssDNA binding characteristics of the SsbA protein were similar to those of the SsbEc protein in every ssDNA binding assay used in this study. The SsbB protein differed from the SsbA and SsbEc proteins, however, both in its binding to short homopolymeric dT(n) oligomers (as judged by polyacrylamide gel-shift assays) and in its binding to the longer naturally occurring X and M13 ssDNAs (as judged by agarose gel-shift assays and electron microscopic analysis). The results indicate that an individual SsbB protein binds to ssDNA with an affinity that is similar or higher than that of the SsbA and SsbEc proteins. However, the manner in which multiple SsbB proteins assemble onto a ssDNA molecule differs from that observed with the SsbA and SsbEc proteins. These results represent the first analysis of paralogous SSB proteins from any bacterial species and provide a foundation for further investigations into the biological roles of these proteins.     

Assembling the pieces

Transition metals function as cofactors in specific proteins, catalyzing electron exchange reactions, binding substrates and stabilizing protein structure. Studies of human diseases and of model organisms have defined many of the molecular details of metal uptake, trafficking, and excretion. The current challenge is to integrate these details into a systematic view of metal content, speciation, localization and use within organisms and ecosystems.     

Calmodulin-binding proteins in the model organism Dictyostelium: a complete & critical review

Calmodulin is an essential protein in the model organism Dictyostelium discoideum. As in other organisms, this small, calcium-regulated protein mediates a diversity of cellular events including chemotaxis, spore germination, and fertilization. Calmodulin works in a calcium-dependent or -independent manner by binding to and regulating the activity of target proteins called calmodulin-binding proteins. Profiling suggests that Dictyostelium has 60 or more calmodulin-binding proteins with specific subcellular localizations. In spite of the central importance of calmodulin, the study of these target proteins is still in its infancy. Here we critically review the history and state of the art of research into all of the identified and presumptive calmodulin-binding proteins of Dictyostelium detailing what is known about each one with suggestions for future research. Two individual calmodulin-binding proteins, the classic enzyme calcineurin A (CNA; protein phosphatase 2B) and the nuclear protein nucleomorphin (NumA), which is a regulator of nuclear number, have been particularly well studied. Research on the role of calmodulin in the function and regulation of the various myosins of Dictyostelium, especially during motility and chemotaxis, suggests that this is an area in which future active study would be particularly valuable. A general, hypothetical model for the role of calmodulin in myosin regulation is proposed.     

Protein stability, flexibility and function

Proteins rely on flexibility to respond to environmental changes, ligand binding and chemical modifications. Potentially, a perturbation that changes the flexibility of a protein may interfere with its function. Millions of mutations have been performed on thousands of proteins in quests for a delineation of the molecular details of their function. Several of these mutations interfered with the binding of a specific ligand with a concomitant effect on the stability of the protein scaffold. It has been ambiguous and not straightforward to recognize if any relationships exist between the stability of a protein and the affinity for its ligand. In this review, we present examples of proteins where changes in stability results in changes in affinity and of proteins where stability and affinity are uncorrelated. We discuss the possibility for a relationship between stability and binding. From the data presented is it clear that there are specific sites (flexibility hotspots) in proteins that are important for both binding and stability. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

The third leg: Molecular dynamics simulations of lipid binding proteins

Molecular dynamics computer simulations can provide a third leg which balances the contributions of both structural biology and binding studies performed on the lipid binding protein family. In this context, these calculations help to establish a dialogue between all three communities, by relating experimental observables with details of structure. Working towards this connection is important, since experience has shown the difficulty of inferring thermodynamic properties from a single static conformation. The challenge is exemplified by ongoing attempts to interpret the impact of mutagenesis on structure and function (i.e. binding). A detailed atomic-level understanding of this system could be achieved with the support of all three legs, paving the way towards rational design of proteins with novel specificities. This paper provides an outline of the connections possible between experiment and theory concerning lipid binding proteins.