Sequence composition and environment effects on residue fluctuations in protein structures

The Effect of Pulling and Twisting Forces on Chameleon Sequence Peptides

Chameleon sequences are amino acid sequences found in several distinct configurations in experiment. They challenge our understanding of the link between sequence and structure, and provide insight into structural competition in proteins. Here, we study the energy landscapes for three such sequences, and interrogate how pulling and twisting forces impact the available structural ensembles. Chameleon sequences do not necessarily exhibit multiple structural ensembles on a multifunnel energy landscape when we consider them in isolation. The application of even small forces leads to drastic changes in the energy landscapes. For pulling forces, we observe transitions from helical to extended structures in a very small span of forces. For twisting forces, the picture is much more complex, and highly dependent on the magnitude and handedness of the applied force as well as the reference angle for the twist. Depending on these parameters, more complex and more simplistic energy landscapes are observed alongside more and less diverse structural ensembles. The impact of even small forces is significant, confirming their likely role in folding events. In addition, small forces exerted by the remaining scaffold of a protein may be sufficient to lead to the adoption of a specific structural ensemble by a chameleon sequence.

Mobilization of iron stored in bacterioferritin, a new target for perturbing iron homeostasis and developing antibacterial and antibiofilm molecules

Antibiotic resistance is a global public health threat. The care of chronic infections is complicated by bacterial biofilms. Biofilm embedded cells can be up to 1000-fold more tolerant to antibiotic treatment than planktonic cells. Antibiotic tolerance is a condition which does not involve mutation and enables bacteria to survive in the presence of antibiotics. The antibiotic tolerance of biofilm-cells often renders antibiotics ineffective, even against strains that do not carry resistance-impairing mutations. This review discusses bacterial iron homeostasis and the strategies being developed to target this bacterial vulnerability, with emphasis on a recently proposed approach which aims at targeting the iron storage protein bacterioferritin (Bfr) and its physiological partner, the ferredoxin Bfd. Bfr regulates cytosolic iron concentrations by oxidizing Fe2+ and storing Fe3+ in its internal cavity, and by forming a complex with Bfd to reduce Fe3+ in the internal cavity and release Fe2+ to the cytosol. Blocking the Bfr-Bfd complex in P. aeruginosa cells causes an irreversible accumulation of Fe3+ in BfrB and simultaneous cytosolic iron depletion, which leads to impaired biofilm maintenance and biofilm cell death. Recently discovered small molecule inhibitors of the Bfr-Bfd complex, which bind Bfr at the Bfd binding site, inhibit iron mobilization, and elicit biofilm cell death.

Deimmunization of flagellin adjuvant for clinical application

Flagellin is the cognate ligand for host pattern recognition receptors, toll-like receptor 5 (TLR5) in the cell surface, and NAIP5/NLRC4 inflammasome in the cytosol. TLR5-binding domain is located in D1 domain, where crucial amino acid sequences are conserved among diverse bacteria. The highly conserved C-terminal 35 amino acids of flagellin were proved to be responsible for the inflammasome activation by binding to NAIP5. D2/D3 domains, located in the central region and exposed to the outside surface of flagellar filament, are heterogeneous across bacterial species and highly immunogenic. Taking advantage of TLR5- and NLRC4-stimulating activities, flagellin has been actively developed as a vaccine adjuvant and immunotherapeutic. Because of its immunogenicity, there exist worries concerning diminished efficacy and possible reactogenicity after repeated administration. Deimmunization of flagellin derivatives while preserving the TLR5/NLRC4-mediated immunomodulatory activity should be the most reasonable option for clinical application. This review describes strategies and current achievements in flagellin deimmunization.

Acetylcholinesterase inhibition of Alzheimer's disease: identification of potential phytochemicals and designing more effective derivatives to manage disease condition

Alzheimer's disease (AD) is a brain disease characterized by gradual memory loss and cognitive impairments. Acetylcholinesterase (AChE) inhibitors-such as donepezil, memantine, and tacrine-are FDA-approved medications for AD treatment. Due to the lack of their efficacy and higher side effects, many researchers have been searching for effective and safer alternatives. In this study, experimentally proved phytochemicals against brain diseases were screened based on their binding energies to the target site of AChE, pharmacokinetic properties, and drug-likeness. Although some phytochemicals showed higher binding affinities than the control drug (donepezil), they did not show permeability across the blood-brain barrier (BBB). However, berberine, anthocyanin, and diterpene alkaloid can cross the BBB and showed good binding affinities of -10.3, -10.1, and -10.2 kcal/mol, respectively. MD simulation and PCA of the simulation data of the protein and protein-ligand complexes proved that the complexes are stable in the biological environment. A total of 16 derivatives of berberine and 3 derivatives of anthocyanin also showed higher binding energies compared to the binding affinity (-11.5 kcal/mol) of the donepezil. The derivatives were designed by substituting -F, -CF3, -CN, and -NH2, and provided higher docking scores due to increasing of nonbonding interactions. MM/GBSA calculations show that the binding free energies of the best predicted derivatives of diterpene alkaloid, anthocyanin, and berberine (DA22, AC11, and BB40) are -100.4 ± 8.4, -79.3 ± 8.7, and -78.3 ± 10.7 kcal/mol, respectively, with the protein. Overall, this study was successful in finding new, highly effective, and possibly safer inhibitors of AChE.Communicated by Ramaswamy H. Sarma.

Deimmunization of flagellin for repeated administration as a vaccine adjuvant

Flagellin, a protein-based Toll-like receptor agonist, is a versatile adjuvant applicable to wide spectrum of vaccines and immunotherapies. Given reiterated treatments of immunogenic biopharmaceuticals should lead to antibody responses precluding repeated administration, the development of flagellin not inducing specific antibodies would greatly expand the chances of clinical applications. Here we computationally identified immunogenic regions in Vibrio vulnificus flagellin B and deimmunized by simply removing a B cell epitope region. The recombinant deimmunized FlaB (dFlaB) maintains stable TLR5-stimulating activity. Multiple immunization of dFlaB does not induce FlaB-specific B cell responses in mice. Intranasally co-administered dFlaB with influenza vaccine enhanced strong Ag-specific immune responses in both systemic and mucosal compartments devoid of FlaB-specific Ab production. Notably, dFlaB showed better protective immune responses against lethal viral challenge compared with wild type FlaB. The deimmunizing B cell epitope deletion did not compromise stability and adjuvanticity, while suppressing unwanted antibody responses that may negatively affected vaccine antigen-directed immune responses in repeated vaccinations. We explain the underlying mechanism of deimmunization by employing molecular dynamics analysis.

Structural Changes beyond the EF-Hand Contribute to Apparent Calcium Binding Affinities: Insights from Parvalbumins

Members of the parvalbumin (PV) family of calcium (Ca²⁺) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca²⁺ affinities and selectivities against competing ions like Mg²⁺ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca²⁺ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca²⁺ affinity (≈-18 kcal/mol) relative to the wild type (WT) (≈-22 kcal/mol) and an S55D/E59D variant with enhanced affinity (≈-24 kcal/mol). Additionally, we also examined the binding of Mg²⁺ to these isoforms, which is much weaker than Ca²⁺. We used mean spherical approximation (MSA) theory to evaluate ion binding thermodynamics within the proteins' EF-hand domains to account for the impact of ions' finite sizes and the surrounding electrolyte composition. While the MSA scores differentiated Mg²⁺ from Ca²⁺, they did not indicate that Ca²⁺ binding affinities at the binding loop differed between the PV isoforms. Instead, molecular mechanics generalized Born surface area (MM/GBSA) approximation energies, which we used to quantify the thermodynamic cost of structural rearrangement of the proteins upon binding ions, indicated that S55D/E59D αPV favored Ca²⁺ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg²⁺ binding was favored for the S55D/E59D αPV and D94S/G98E αPV variants by -18.32 and -1.65 kcal/mol, respectively. These energies implicate significant contributions to ion binding beyond oxygen coordination at the binding loop, which stemmed from changes in α-helicity, β-sheet character, and hydrogen bonding. Hence, Ca²⁺ affinity and selectivity against Mg²⁺ are emergent properties stemming from both local effects within the proteins' ion binding sites as well as non-local contributions elsewhere. Our findings broaden our understanding of the molecular bases governing αPV ion binding that are likely shared by members of the broad family of CBPs.

Epitopes of Protein Binders Are Related to the Structural Flexibility of a Target Protein Surface

Protein binders including antibodies are known not to bind to random sites of target proteins, and their functional effectiveness mainly depends on the binding region, called the epitope. For the development of protein binders with desired functions, it is thus critical to understand which surface region protein binders prefer (or do not prefer) to bind. The current methods for epitope prediction focus on static indicators such as structural geometry or amino acid propensity, whereas protein binding events are in fact a consequence of dynamic interactions. Here, we demonstrate that the preference for a binding site by protein binders is strongly related to the structural flexibility of a target protein surface. Molecular dynamics simulations on unbound forms of antigen structures revealed that the antigen surface in direct contact with antibodies is less flexible than the rest of the surface. This tendency was shown to be similar in other non-antibody protein binders such as affibody, DARPin, monobody, and repebody. We also found that the relatedness of epitopes to the structural flexibility of a target protein surface is dependent on the secondary structure elements of paratopes. Monobody and repebody, whose binding sites are composed of β-strands, distinctively prefer to bind to a relatively more rigid region of a target protein. These observations enabled us to develop a simple epitope prediction method which shows a comparable performance to the commonly used ones.

Importance of Fluctuating Amino Acid Residues in Folding and Binding of Proteins

BACKGROUND

Conformational flexibility of proteins remains as one of the major events in protein-protein/DNA/ligand/small molecule binding to achieve its biological function in the cell. The availability of high-resolution structures of protein complexes is a valuable resource for researchers to understand the mechanisms behind such interactions and it is found that the flexibility of amino acid residues at binding sites is crucial for many important functions in the cell.

METHODS

In this article, our statistical method (PreFRP) developed based on fluctuating amino acid residues and various amino acid indices related to flexibility/rigidity were used to study the importance of fluctuating amino acid residues in thermonucleases from pathogenic bacteria, cell penetrating peptides and intrinsically disordered proteins responsible for many neural disorders.

RESULTS

The results from our analysis reveal the importance of fluctuating amino acid residues in folding and binding of proteins. The role of moderate and high fluctuating residues in themonucleases, cell penetrating peptide and disordered regions are discussed in detail.

CONCLUSION

Therefore, our analysis will help in understanding the importance of fluctuating amino acid residues in proteins which undergo a conformation change phenomenon.

Challenges in structural approaches to cell modeling

Computational modeling is essential for structural characterization of biomolecular mechanisms across the broad spectrum of scales. Adequate understanding of biomolecular mechanisms inherently involves our ability to model them. Structural modeling of individual biomolecules and their interactions has been rapidly progressing. However, in terms of the broader picture, the focus is shifting toward larger systems, up to the level of a cell. Such modeling involves a more dynamic and realistic representation of the interactomes in vivo, in a crowded cellular environment, as well as membranes and membrane proteins, and other cellular components. Structural modeling of a cell complements computational approaches to cellular mechanisms based on differential equations, graph models, and other techniques to model biological networks, imaging data, etc. Structural modeling along with other computational and experimental approaches will provide a fundamental understanding of life at the molecular level and lead to important applications to biology and medicine. A cross section of diverse approaches presented in this review illustrates the developing shift from the structural modeling of individual molecules to that of cell biology. Studies in several related areas are covered: biological networks; automated construction of three-dimensional cell models using experimental data; modeling of protein complexes; prediction of non-specific and transient protein interactions; thermodynamic and kinetic effects of crowding; cellular membrane modeling; and modeling of chromosomes. The review presents an expert opinion on the current state-of-the-art in these various aspects of structural modeling in cellular biology, and the prospects of future developments in this emerging field.

PreFRP: Prediction and visualization of fluctuation residues in proteins

Aim:

The PreFRP web server extracts sequence and basic information of a protein structure and groups amino acid residues in a protein into three important types such as high, moderate, and weak fluctuating residues.

Materials and Methods:

The server takes a protein data bank file or an amino acid sequence as input and prints the probability of amino acid residues to fluctuate. The server also provides a link to Jmol, a molecular visualization program to visualize the high, moderate, and weak fluctuating residues in three different colors.

Results:

Prediction and visualization of fluctuating amino acid residues in proteins may help to understand the complex three-dimensional structure of proteins and may further help in docking and mutation experiments.

Availability:

The web server is freely accessible through the web page of the author's institution http://www.mpi.edu.in/prefrp/link.html.

Chameleon sequences in neurodegenerative diseases

Chameleon sequences can adopt either alpha helix sheet or a coil conformation. Defining chameleon sequences in PDB (Protein Data Bank) may yield to an insight on defining peptides and proteins responsible in neurodegeneration. In this research, we benefitted from the large PDB and performed a sequence analysis on Chameleons, where we developed an algorithm to extract peptide segments with identical sequences, but different structures. In order to find new chameleon sequences, we extracted a set of 8315 non-redundant protein sequences from the PDB with an identity less than 25%. Our data was classified to "helix to strand (HE)", "helix to coil (HC)" and "strand to coil (CE)" alterations. We also analyzed the occurrence of singlet and doublet amino acids and the solvent accessibility in the chameleon sequences; we then sorted out the proteins with the most number of chameleon sequences and named them Chameleon Flexible Proteins (CFPs) in our dataset. Our data revealed that Gly, Val, Ile, Tyr and Phe, are the major amino acids in Chameleons. We also found that there are proteins such as Insulin Degrading Enzyme IDE and GTP-binding nuclear protein Ran (RAN) with the most number of chameleons (640 and 405 respectively). These proteins have known roles in neurodegenerative diseases. Therefore it can be inferred that other CFP's can serve as key proteins in neurodegeneration, and a study on them can shed light on curing and preventing neurodegenerative diseases.

Protein-protein docking: from interaction to interactome

The protein-protein docking problem is one of the focal points of activity in computational biophysics and structural biology. The three-dimensional structure of a protein-protein complex, generally, is more difficult to determine experimentally than the structure of an individual protein. Adequate computational techniques to model protein interactions are important because of the growing number of known protein structures, particularly in the context of structural genomics. Docking offers tools for fundamental studies of protein interactions and provides a structural basis for drug design. Protein-protein docking is the prediction of the structure of the complex, given the structures of the individual proteins. In the heart of the docking methodology is the notion of steric and physicochemical complementarity at the protein-protein interface. Originally, mostly high-resolution, experimentally determined (primarily by x-ray crystallography) protein structures were considered for docking. However, more recently, the focus has been shifting toward lower-resolution modeled structures. Docking approaches have to deal with the conformational changes between unbound and bound structures, as well as the inaccuracies of the interacting modeled structures, often in a high-throughput mode needed for modeling of large networks of protein interactions. The growing number of docking developers is engaged in the community-wide assessments of predictive methodologies. The development of more powerful and adequate docking approaches is facilitated by rapidly expanding information and data resources, growing computational capabilities, and a deeper understanding of the fundamental principles of protein interactions.

Local packing modulates diversity of iron pathways and cooperative behavior in eukaryotic and prokaryotic ferritins

Ferritin-like molecules show a remarkable combination of the evolutionary conserved activity of iron uptake and release that engage different pores in the conserved ferritin shell. It was hypothesized that pore selection and iron traffic depend on dynamic allostery with no conformational changes in the backbone. In this study, we detect the allosteric networks in Pseudomonas aeruginosa bacterioferritin (BfrB), bacterial ferritin (FtnA), and bullfrog M and L ferritins (Ftns) by a network-weaving algorithm (NWA) that passes threads of an allosteric network through highly correlated residues using hierarchical clustering. The residue-residue correlations are calculated in the packing-on elastic network model that introduces atom packing into the common packing-off model. Applying NWA revealed that each of the molecules has an extended allosteric network mostly buried inside the ferritin shell. The structure of the networks is consistent with experimental observations of iron transport: The allosteric networks in BfrB and FtnA connect the ferroxidase center with the 4-fold pores and B-pores, leaving the 3-fold pores unengaged. In contrast, the allosteric network directly links the 3-fold pores with the 4-fold pores in M and L Ftns. The majority of the network residues are either on the inner surface or buried inside the subunit fold or at the subunit interfaces. We hypothesize that the ferritin structures evolved in a way to limit the influence of functionally unrelated events in the cytoplasm on the allosteric network to maintain stability of the translocation mechanisms. We showed that the residue-residue correlations and the resultant long-range cooperativity depend on the ferritin shell packing, which, in turn, depends on protein sequence composition. Switching from the packing-on to the packing-off model reduces correlations by 35%-38% so that no allosteric network can be found. The influence of the side-chain packing on the allosteric networks explains the diversity in mechanisms of iron traffic suggested by experimental approaches.

From local structure to a global framework: recognition of protein folds

Protein folding has been a major area of research for many years. Nonetheless, the mechanisms leading to the formation of an active biological fold are still not fully apprehended. The huge amount of available sequence and structural information provides hints to identify the putative fold for a given sequence. Indeed, protein structures prefer a limited number of local backbone conformations, some being characterized by preferences for certain amino acids. These preferences largely depend on the local structural environment. The prediction of local backbone conformations has become an important factor to correctly identifying the global protein fold. Here, we review the developments in the field of local structure prediction and especially their implication in protein fold recognition.

Stability mechanisms of a thermophilic laccase probed by molecular dynamics

Laccases are highly stable, industrially important enzymes capable of oxidizing a large range of substrates. Causes for their stability are, as for other proteins, poorly understood. In this work, multiple-seed molecular dynamics (MD) was applied to a Trametes versicolor laccase in response to variable ionic strengths, temperatures, and glycosylation status. Near-physiological conditions provided excellent agreement with the crystal structure (average RMSD ∼0.92 Å) and residual agreement with experimental B-factors. The persistence of backbone hydrogen bonds was identified as a key descriptor of structural response to environment, whereas solvent-accessibility, radius of gyration, and fluctuations were only locally relevant. Backbone hydrogen bonds decreased systematically with temperature in all simulations (∼9 per 50 K), probing structural changes associated with enthalpy-entropy compensation. Approaching T_opt (∼350 K) from 300 K, this change correlated with a beginning “unzipping” of critical β-sheets. 0 M ionic strength triggered partial denucleation of the C-terminal (known experimentally to be sensitive) at 400 K, suggesting a general salt stabilization effect. In contrast, F⁻ (but not Cl⁻) specifically impaired secondary structure by formation of strong hydrogen bonds with backbone NH, providing a mechanism for experimentally observed small anion destabilization, potentially remedied by site-directed mutagenesis at critical intrusion sites. N-glycosylation was found to support structural integrity by increasing persistent backbone hydrogen bonds by ∼4 across simulations, mainly via prevention of F⁻ intrusion. Hydrogen-bond loss in distinct loop regions and ends of critical β-sheets suggest potential strategies for laboratory optimization of these industrially important enzymes.

Low-resolution structural modeling of protein interactome

Structural characterization of protein-protein interactions across the broad spectrum of scales is key to our understanding of life at the molecular level. Low-resolution approach to protein interactions is needed for modeling large interaction networks, given the significant level of uncertainties in large biomolecular systems and the high-throughput nature of the task. Since only a fraction of protein structures in interactome are determined experimentally, protein docking approaches are increasingly focusing on modeled proteins. Current rapid advancement of template-based modeling of protein-protein complexes is following a long standing trend in structure prediction of individual proteins. Protein-protein templates are already available for almost all interactions of structurally characterized proteins, and about one third of such templates are likely correct.

Structure fluctuations and conformational changes in protein binding

Structure fluctuations and conformational changes accompany all biological processes involving macromolecules. The paper presents a classification of protein residues based on the normalized equilibrium fluctuations of the residue centers of mass in proteins and a statistical analysis of conformation changes in the side-chains upon binding. Normal mode analysis and an elastic network model were applied to a set of protein complexes to calculate the residue fluctuations and develop the residue classification. Comparison with a classification based on normalized B-factors suggests that the B-factors may underestimate protein flexibility in solvent. Our classification shows that protein loops and disordered fragments are enriched with highly fluctuating residues and depleted with weakly fluctuating residues. Strategies for engineering thermostable proteins are discussed. To calculate the dihedral angles distribution functions, the configuration space was divided into cells by a cubic grid. The effect of protein association on the distribution functions depends on the amino acid type and a grid step in the dihedral angles space. The changes in the dihedral angles increase from the near-backbone dihedral angle to the most distant one, for most residues. On average, one fifth of the interface residues change the rotamer state upon binding, whereas the rest of the interface residues undergo local readjustments within the same rotamer.

The structure of the BfrB-Bfd complex reveals protein-protein interactions enabling iron release from bacterioferritin

Ferritin-like molecules are unique to cellular iron homeostasis because they can store iron at concentrations much higher than those dictated by the solubility of Fe(3+). Very little is known about the protein interactions that deliver iron for storage or promote the mobilization of stored iron from ferritin-like molecules. Here, we report the X-ray crystal structure of Pseudomonas aeruginosa bacterioferritin (Pa-BfrB) in complex with bacterioferritin-associated ferredoxin (Pa-Bfd) at 2.0 Å resolution. As the first example of a ferritin-like molecule in complex with a cognate partner, the structure provides unprecedented insight into the complementary interface that enables the [2Fe-2S] cluster of Pa-Bfd to promote heme-mediated electron transfer through the BfrB protein dielectric (~18 Å), a process that is necessary to reduce the core ferric mineral and facilitate mobilization of Fe(2+). The Pa-BfrB-Bfd complex also revealed the first structure of a Bfd, thus providing a first view to what appears to be a versatile metal binding domain ubiquitous to the large Fer2_BFD family of proteins and enzymes with diverse functions. Residues at the Pa-BfrB-Bfd interface are highly conserved in Bfr and Bfd sequences from a number of pathogenic bacteria, suggesting that the specific recognition between Pa-BfrB and Pa-Bfd is of widespread significance to the understanding of bacterial iron homeostasis.

Analysis of conformational variation in macromolecular structural models

Experimental conditions or the presence of interacting components can lead to variations in the structural models of macromolecules. However, the role of these factors in conformational selection is often omitted by in silico methods to extract dynamic information from protein structural models. Structures of small peptides, considered building blocks for larger macromolecular structural models, can substantially differ in the context of a larger protein. This limitation is more evident in the case of modeling large multi-subunit macromolecular complexes using structures of the individual protein components. Here we report an analysis of variations in structural models of proteins with high sequence similarity. These models were analyzed for sequence features of the protein, the role of scaffolding segments including interacting proteins or affinity tags and the chemical components in the experimental conditions. Conformational features in these structural models could be rationalized by conformational selection events, perhaps induced by experimental conditions. This analysis was performed on a non-redundant dataset of protein structures from different SCOP classes. The sequence-conformation correlations that we note here suggest additional features that could be incorporated by in silico methods to extract dynamic information from protein structural models.

Side-chain conformational changes upon Protein-Protein Association

Conformational changes upon protein-protein association are the key element of the binding mechanism. The study presents a systematic large-scale analysis of such conformational changes in the side chains. The results indicate that short and long side chains have different propensities for the conformational changes. Long side chains with three or more dihedral angles are often subject to large conformational transition. Shorter residues with one or two dihedral angles typically undergo local conformational changes not leading to a conformational transition. A relationship between the local readjustments and the equilibrium fluctuations of a side chain around its unbound conformation is suggested. Most of the side chains undergo larger changes in the dihedral angle most distant from the backbone. The frequencies of the core-to-surface interface transitions of six nonpolar residues and Tyr are larger than the frequencies of the opposite surface-to-core transitions. The binding increases both polar and nonpolar interface areas. However, the increase of the nonpolar area is larger for all considered classes of protein complexes, suggesting that the protein association perturbs the unbound interfaces to increase the hydrophobic contribution to the binding free energy. To test modeling approaches to side-chain flexibility in protein docking, conformational changes in the X-ray set were compared with those in the docking decoy sets. The results lead to a better understanding of the conformational changes in proteins and suggest directions for efficient conformational sampling in docking protocols.