Efficient Characterization of Protein Cavities within Molecular Simulation Trajectories: trj_cavity

The role of hydrophobicity in the cold denaturation of proteins under high pressure: A study on apomyoglobin

An exciting debate arises when microscopic mechanisms involved in the denaturation of proteins at high pressures are explained. In particular, the issue emerges when the hydrophobic effect is invoked, given that hydrophobicity cannot elucidate by itself the volume changes measured during protein unfolding. In this work, we study by the use of molecular dynamics simulations and essential dynamics analysis the relation between the solvation dynamics, volume, and water structure when apomyoglobin is subjected to a hydrostatic pressure regime. Accordingly, the mechanism of cold denaturation of proteins under high-pressure can be related to the disruption of the hydrogen-bond network of water favoring the coexistence of two states, low-density and high-density water, which directly implies in the formation of a molten globule once the threshold of 200 MPa has been overcome.

Microscopic Structure of Compacted Polyelectrolyte Complexes: Insights from Molecular Dynamics Simulations

We utilize atomistic molecular dynamics (MD) simulations to study local structural changes inside a polyelectrolyte complex consisting of poly(styrenesulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) upon densification, in analogy to ultracentrifugation in experiments. In particular, we focus on the water content and on the reinforcement of the PSS-PDADMA network for various external accelerations. We demonstrate that apart from the formation of mesoscopic pores observed experimentally also the microscopic structure and the local relaxation processes likely affect the unique rheological properties of compacted polyelectrolyte complexes, as densification increases both the number of PSS-PDADMA coordinations and the intermixing of PSS and PDADMA. These processes slow down local rearrangements, thus further stabilizing the compacted state. We find that the concept of binary PSS-PDADMA salt bonds-relevant for theoretical models-is not strictly valid in the dense limit.

Mechanics and nanovoid nucleation dynamics: effects of polar functionality in glassy polymer networks

We use molecular simulations and experiments to rationalize the properties of a class of networks based on dicyclopentadiene (DCPD), a polymer with excellent fracture toughness and a high glass transition temperature (Tg), copolymerized with 5-norbornene-2-methanol (NBOH). DCPD is a highly non-polar hydrocarbon, while NBOH contains a hydroxy group, introducing polar functionality and hydrogen bonds (H-bonds). NBOH thus represents a possible route to improve the chemical compatibility of DCPD-based networks with less-hydrophobic materials. We systematically vary the NBOH content (polar chemistry) in DCPD networks, while keeping other network parameters nearly constant, including the molecular weight between cross-links, chain rigidity, and Tg. Using molecular dynamics (MD) simulations, we quantify the thermovolumetric and mechanical properties, including Tg, cohesive energy density, stiffness, and yield strength. We compare these results with experiments on networks of similar composition, finding good agreement. The relation between these properties and polar chemistry are studied by examining a secondary network of physical cross-links, formed by hydrogen bonds between NBOH units. Further, we examine nanovoid formation, an energy dissipation mechanism hypothesized to contribute to the toughness of pDCPD. Using metadynamics to accelerate sampling, we quantify the nanovoid nucleation rate under hydrostatic tension, similar to the stress state in the plastic zone preceding a crack tip. Small additions of NBOH have minimal effect, but the rate drops steeply with larger amounts. Several properties are mapped at nanometer scales, including stiffness and mobility, and associated with void nucleation. Estimates of the length- and time-scale of the plastic zone near a crack tip are used in discussing nanovoid formation as a plausible toughening mechanism in these materials. Overall, the results suggest that pDCPD tolerates the addition of some polar chemistry without degrading its excellent mechanical properties.

Role of Extracellular Loops and Membrane Lipids for Ligand Recognition in the Neuronal Adenosine Receptor Type 2A: An Enhanced Sampling Simulation Study

Human G-protein coupled receptors (GPCRs) are important targets for pharmaceutical intervention against neurological diseases. Here, we use molecular simulation to investigate the key step in ligand recognition governed by the extracellular domains in the neuronal adenosine receptor type 2A (hA_2AR), a target for neuroprotective compounds. The ligand is the high-affinity antagonist (4-(2-(7-amino-2-(furan-2-yl)-[1,2,4]triazolo[1,5-a][1,3,5]triazin-5-ylamino)ethyl)phenol), embedded in a neuronal membrane mimic environment. Free energy calculations, based on well-tempered metadynamics, reproduce the experimentally measured binding affinity. The results are consistent with the available mutagenesis studies. The calculations identify a vestibular binding site, where lipids molecules can actively participate to stabilize ligand binding. Bioinformatic analyses suggest that such vestibular binding site and, in particular, the second extracellular loop, might drive the ligand toward the orthosteric binding pocket, possibly by allosteric modulation. Taken together, these findings point to a fundamental role of the interaction between extracellular loops and membrane lipids for ligands’ molecular recognition and ligand design in hA_2AR.

An Algorithm for Computing Side Chain Conformational Variations of a Protein Tunnel/Channel

In this paper, a novel method to compute side chain conformational variations for a protein molecule tunnel (or channel) is proposed. From the conformational variations, we compute the flexibly deformed shapes of the initial tunnel, and present a way to compute the maximum size of the ligand that can pass through the deformed tunnel. By using the two types of graphs corresponding to amino acids and their side chain rotamers, the suggested algorithm classifies amino acids and rotamers which possibly have collisions. Based on the divide and conquer technique, local side chain conformations are computed first, and then a global conformation is generated by combining them. With the exception of certain cases, experimental results show that the algorithm finds up to 327,680 valid side chain conformations from 12⁸~12³³ conformation candidates within three seconds.

Effect of surfactants on the interaction of phenol with laccase: Molecular docking and molecular dynamics simulation studies

Some surfactants can enhance the removal of phenol by laccase (Lac) in various industrial effluents. Their behavior and function in the biodegradation of phenolic wastewater have been experimentally reported by many researchers, but the underlying molecular mechanism is still unclear. Therefore, the interaction mechanisms of phenol with Lac from Trametes versicolor were investigated in the presence or absence of Triton X-100 (TX100) or rhamnolipid (RL) by molecular docking and molecular dynamics (MD) simulations. The results indicate that phenol contacts with an active site of Lac by hydrogen bonds (HBs) and van der Waals (vdW) interactions in aqueous solution for maintaining its stability. The presence of TX100 or RL results in the significant changes of enzymatic conformations. Meanwhile, the hydrophobic parts of surfactants contact with the outside surface of Lac. These changes lead to the decrease of binding energy between phenol and Lac. The migration behavior of water molecules within hydration shell is also inevitably affected. Therefore, the amphipathic TX100 or RL may influence the phenol degradation ability of Lac by modulating their interactions and water environment. This study offers molecular level of understanding on the function of surfactants in biosystem.

Structural Prediction of the Dimeric Form of the Mammalian Translocator Membrane Protein TSPO: A Key Target for Brain Diagnostics

Positron emission tomography (PET) radioligands targeting the human translocator membrane protein (TSPO) are broadly used for the investigations of neuroinflammatory conditions associated with neurological disorders. Structural information on the mammalian protein homodimers-the suggested functional state of the protein-is limited to a solid-state nuclear magnetic resonance (NMR) study and to a model based on the previously-deposited solution NMR structure of the monomeric mouse protein. Computational studies performed here suggest that the NMR-solved structure in the presence of detergents is not prone to dimer formation and is furthermore unstable in its native membrane environment. We, therefore, propose a new model of the functionally-relevant dimeric form of the mouse protein, based on a prokaryotic homologue. The model, fully consistent with solid-state NMR data, is very different from the previous predictions. Hence, it provides, for the first time, structural insights into this pharmaceutically-important target which are fully consistent with experimental data.

Exploring the entry route of palmitic acid and palmitoylcarnitine into myoglobin

Myoglobin, besides its role in oxygen turnover, has gained recognition as a potential regulator of lipid metabolism. Previously, we confirmed the interaction of fatty acids and acylcarnitines with Oxy-Myoglobin, using both molecular dynamic simulations and Isothermal Titration Calorimetry studies. However, those studies were limited to testing only the binding sites derived from homology to fatty acid binding proteins and predictions using automated docking. To explore the entry mechanisms of the lipid ligands into myoglobin, we conducted molecular dynamic simulations of murine Oxy- and Deoxy-Mb structures with palmitate or palmitoylcarnitine starting at different positions near the protein surface. The simulations indicated that both ligands readily (under ∼10-20 ns) enter the Oxy-Mb structure through a dynamic area ("portal region") near heme, known to be the entry point for small molecule gaseous ligands like O₂, CO and NO. The entry is not observed with Deoxy-Mb where lipid ligands move away from protein surface, due to a compaction of the entry portal and the heme-containing crevice in the Mb protein upon O₂ removal. The results suggest quick spontaneous binding of lipids to Mb driven by hydrophobic interactions, strongly enhanced by oxygenation, and consistent with the emergent role of Mb in lipid metabolism.

Interaction of the human erythrocyte Band 3 anion exchanger 1 (AE1, SLC4A1) with lipids and glycophorin A: Molecular organization of the Wright (Wr) blood group antigen

The Band 3 (AE1, SLC4A1) membrane protein is found in red blood cells and in kidney where it functions as an electro-neutral chloride/bicarbonate exchanger. In this study, we have used molecular dynamics simulations to provide the first realistic model of the dimeric membrane domain of human Band 3 in an asymmetric lipid bilayer containing a full complement of phospholipids, including phosphatidylinositol 4,5–bisphosphate (PIP₂) and cholesterol, and its partner membrane protein Glycophorin A (GPA). The simulations show that the annular layer in the inner leaflet surrounding Band 3 was enriched in phosphatidylserine and PIP₂ molecules. Cholesterol was also enriched around Band 3 but also at the dimer interface. The interaction of these lipids with specific sites on Band 3 may play a role in the folding and function of this anion transport membrane protein. GPA associates with Band 3 to form the Wright (Wr) blood group antigen, an interaction that involves an ionic bond between Glu658 in Band 3 and Arg61 in GPA. We were able to recreate this complex by performing simulations to allow the dimeric transmembrane portion of GPA to interact with Band 3 in a model membrane. Large-scale simulations showed that the GPA dimer can bridge Band 3 dimers resulting in the dynamic formation of long strands of alternating Band 3 and GPA dimers.

Human Band 3 (AE1, SLC4A1), an abundant 911 amino acid glycoprotein, catalyzes the exchange of bicarbonate and chloride across the red blood cell membrane, a process necessary for efficient respiration. Malfunction of Band 3 leads to inherited diseases such as Southeast Asian Ovalocytosis, hereditary spherocytosis and distal renal tubular acidosis. Despite much available structural and functional data about Band 3, key questions about the conformational changes associated with transport and the molecular details of its interaction with lipids and other proteins remain unanswered. In this study, we have used computer simulations to investigate the dynamics of Band 3 in lipid bilayers that resemble the red blood cell plasma membrane. Our results suggest that negatively charged phospholipids and cholesterol interact strongly with Band 3 forming an annulus around the protein. Glycophorin A (GPA) interacts with Band 3 to form the Wright (Wr) blood group antigen. We were able to recreate this complex and show that GPA promotes the clustering of Band 3 in red blood cell membranes. Understanding the molecular details of the interaction of Band 3 with GPA has provided new insights into the nature of the Wright blood group antigen.

Partial Intrinsic Disorder Governs the Dengue Capsid Protein Conformational Ensemble

The 11 kDa, positively charged dengue capsid protein (C protein) exists stably as a homodimer and colocalizes with the viral genome within mature viral particles. Its core is composed of four alpha helices encompassing a small hydrophobic patch that may interact with lipids, but approximately 20% of the protein at the N-terminus is intrinsically disordered, making it challenging to elucidate its conformational landscape. Here, we combine small-angle X-ray scattering (SAXS), amide hydrogen-deuterium exchange mass spectrometry (HDXMS), and atomic-resolution molecular dynamics (MD) simulations to probe the dynamics of dengue C proteins. We show that the use of MD force fields (FFs) optimized for intrinsically disordered proteins (IDPs) is necessary to capture their conformational landscape and validate the computationally generated ensembles with reference to SAXS and HDXMS data. Representative ensembles of the C protein dimer are characterized by alternating, clamp-like exposure and occlusion of the internal hydrophobic patch, as well as by residual helical structure at the disordered N-terminus previously identified as a potential source of autoinhibition. Such dynamics are likely to determine the multifunctionality of the C protein during the flavivirus life cycle and hence impact the design of novel antiviral compounds.

The architecture of the OmpC-MlaA complex sheds light on the maintenance of outer membrane lipid asymmetry in Escherichia coli

A distinctive feature of the Gram-negative bacterial cell envelope is the asymmetric outer membrane (OM), where lipopolysaccharides and phospholipids (PLs) reside in the outer and inner leaflets, respectively. This unique lipid asymmetry renders the OM impermeable to external insults, including antibiotics and bile salts. In Escherichia coli, the complex comprising osmoporin OmpC and the OM lipoprotein MlaA is believed to maintain lipid asymmetry by removing mislocalized PLs from the outer leaflet of the OM. How this complex performs this function is unknown. Here, we defined the molecular architecture of the OmpC-MlaA complex to gain insights into its role in PL transport. Using in vivo photo-cross-linking and molecular dynamics simulations, we established that MlaA interacts extensively with OmpC and is located entirely within the lipid bilayer. In addition, MlaA forms a hydrophilic channel, likely enabling PL translocation across the OM. We further showed that flexibility in a hairpin loop adjacent to the channel is critical in modulating MlaA activity. Finally, we demonstrated that OmpC plays a functional role in maintaining OM lipid asymmetry together with MlaA. Our work offers glimpses into how the OmpC-MlaA complex transports PLs across the OM and has important implications for future antibacterial drug development.

Computational Insight Into the Structural Organization of Full-Length Toll-Like Receptor 4 Dimer in a Model Phospholipid Bilayer

Toll-like receptors (TLRs) are a unique category of pattern recognition receptors that recognize distinct pathogenic components, often utilizing the same set of downstream adaptors. Specific molecular features of extracellular, transmembrane (TM), and cytoplasmic domains of TLRs are crucial for coordinating the complex, innate immune signaling pathway. Here, we constructed a full-length structural model of TLR4-a widely studied member of the interleukin-1 receptor/TLR superfamily-using homology modeling, protein-protein docking, and molecular dynamics simulations to understand the differential domain organization of TLR4 in a membrane-aqueous environment. Results showed that each functional domain of the membrane-bound TLR4 displayed several structural transitions that are biophysically essential for plasma membrane integration. Specifically, the extracellular and cytoplasmic domains were partially immersed in the upper and lower leaflets of the membrane bilayer. Meanwhile, TM domains tilted considerably to overcome the hydrophobic mismatch with the bilayer core. Our analysis indicates an alternate dimerization or a potential oligomerization interface of TLR4-TM. Moreover, the helical properties of an isolated TM dimer partly agree with that of the full-length receptor. Furthermore, membrane-absorbed or solvent-exposed surfaces of the toll/interleukin-1 receptor domain are consistent with previous X-ray crystallography and biochemical studies. Collectively, we provided a complete structural model of membrane-bound TLR4 that strengthens our current understanding of the complex mechanism of receptor activation and adaptor recruitment in the innate immune signaling pathway.

Protein Fluctuations and Cavity Changes Relationship

Protein cavities and tunnels are critical for function. Ligand recognition and binding, transport, and enzyme catalysis require cavities rearrangements. Therefore, the flexibility of cavities should be guaranteed by protein vibrational dynamics. Molecular dynamics simulations provide a framework to explore conformational plasticity of protein cavities. Herein, we present a novel procedure to characterize the dynamics of protein cavities in terms of their volume gradient vector. For this purpose, we make use of algorithms for calculation of the cavity volume that result robust for numerical differentiations. Volume gradient vector is expressed in terms of principal component analysis obtained from equilibrated molecular dynamics simulations. We analyze contributions of principal component modes to the volume gradient vector according to their frequency and degree of delocalization. In all our test cases, we find that low frequency modes play a critical role together with minor contributions of high frequency modes. These modes involve concerted motions of significant fractions of the total residues lining the cavities. We make use of variations of the potential energy of a protein in the direction of the volume gradient vector as a measure of flexibility of the cavity. We show that proteins whose collective low frequency fluctuations contribute the most to changes of cavity volume exhibit more flexible cavities.

Potential role of salt-bridges in the hinge-like movement of apicomplexa specific β-hairpin of Plasmodium and Toxoplasma profilins: A molecular dynamics simulation study

Profilin is one of the actin-binding proteins that regulate dynamics of actin polymerization. It plays a key role in cell motility and invasion. It also interacts with several other proteins notably through its poly-L-proline (PLP) binding site. Profilin in apicomplexa is characterized by a unique mini-domain consisting of a large β-hairpin extension and an acidic loop which is relatively longer in Plasmodium species. Profilin is essential for the invasive blood stages of Plasmodium falciparum. In the current study, unbound profilins from Plasmodium falciparum (Pf), Toxoplasma gondii (Tg), and Homo sapiens (Hs) were subjected to molecular dynamics (MD) simulations for a timeframe of 100 ns each to understand the conformational dynamics of these proteins. It was found that the β-hairpin of profilins from Pf and Tg shows a hinge-like movement. This movement in Pf profilin may possibly be driven by the loss of a salt-bridge within profilin. The impact of this conformational change on actin binding was assessed by docking three dimensional (3D) structures of profilin from Pf and Tg with their corresponding actins using ClusPro2.0. The stability of docked Pf profilin-actin complex was assessed through a 50 ns MD simulation. As Hs profilin I does not have the apicomplexa specific mini-domain, MD simulation was performed for this protein and its dynamics was compared to that of profilins from Pf and Tg. Using an immunoinformatics approach, potential epitope regions were predicted for Pf profilin. This has a potential application in the design of vaccines as they mapped to its unique mini-domain.

Geometric Detection Algorithms for Cavities on Protein Surfaces in Molecular Graphics: A Survey

Detecting and analyzing protein cavities provides significant information about active sites for biological processes (e.g., protein-protein or protein-ligand binding) in molecular graphics and modeling. Using the three-dimensional structure of a given protein (i.e., atom types and their locations in 3D) as retrieved from a PDB (Protein Data Bank) file, it is now computationally viable to determine a description of these cavities. Such cavities correspond to pockets, clefts, invaginations, voids, tunnels, channels, and grooves on the surface of a given protein. In this work, we survey the literature on protein cavity computation and classify algorithmic approaches into three categories: evolution-based, energy-based, and geometry-based. Our survey focuses on geometric algorithms, whose taxonomy is extended to include not only sphere-, grid-, and tessellation-based methods, but also surface-based, hybrid geometric, consensus, and time-varying methods. Finally, we detail those techniques that have been customized for GPU (Graphics Processing Unit) computing.

GPU-based detection of protein cavities using Gaussian surfaces

Background

Protein cavities play a key role in biomolecular recognition and function, particularly in protein-ligand interactions, as usual in drug discovery and design. Grid-based cavity detection methods aim at finding cavities as aggregates of grid nodes outside the molecule, under the condition that such cavities are bracketed by nodes on the molecule surface along a set of directions (not necessarily aligned with coordinate axes). Therefore, these methods are sensitive to scanning directions, a problem that we call cavity ground-and-walls ambiguity, i.e., they depend on the position and orientation of the protein in the discretized domain. Also, it is hard to distinguish grid nodes belonging to protein cavities amongst all those outside the protein, a problem that we call cavity ceiling ambiguity.

Results

We solve those two ambiguity problems using two implicit isosurfaces of the protein, the protein surface itself (called inner isosurface) that excludes all its interior nodes from any cavity, and the outer isosurface that excludes most of its exterior nodes from any cavity. Summing up, the cavities are formed from nodes located between these two isosurfaces. It is worth noting that these two surfaces do not need to be evaluated (i.e., sampled), triangulated, and rendered on the screen to find the cavities in between; their defining analytic functions are enough to determine which grid nodes are in the empty space between them.

Conclusion

This article introduces a novel geometric algorithm to detect cavities on the protein surface that takes advantage of the real analytic functions describing two Gaussian surfaces of a given protein.

Allosteric Communication Networks in Proteins Revealed through Pocket Crosstalk Analysis

The detection and characterization of binding pockets and allosteric communication in proteins is crucial for studying biological regulation and performing drug design. Nowadays, ever-longer molecular dynamics (MD) simulations are routinely used to investigate the spatiotemporal evolution of proteins. Yet, there is no computational tool that can automatically detect all the pockets and potential allosteric communication networks along these extended MD simulations. Here, we use a novel and fully automated algorithm that examines pocket formation, dynamics, and allosteric communication embedded in microsecond-long MD simulations of three pharmaceutically relevant proteins, namely, PNP, A2A, and Abl kinase. This dynamic analysis uses pocket crosstalk, defined as the temporal exchange of atoms between adjacent pockets, along the MD trajectories as a fingerprint of hidden allosteric communication networks. Importantly, this study indicates that dynamic pocket crosstalk analysis provides new mechanistic understandings on allosteric communication networks, enriching the available experimental data. Thus, our results suggest the prospective use of this unprecedented dynamic analysis to characterize transient binding pockets for structure-based drug design.

POVME 3.0: Software for Mapping Binding Pocket Flexibility

We present a substantial update to the open-source POVME binding pocket analysis software. New capabilities of POVME 3.0 include a flexible chemical coloring scheme for feature identification, postanalysis tools for comparing large ensembles of pockets (e.g., from molecular dynamics simulations), and the introduction of scripts and methods that facilitate binding pocket comparison and analysis. We envision the use of this software for visualization of binding pocket dynamics, selection of representative structures for ensemble docking, and incorporation of molecular dynamics results into ligand design efforts.

Progress in Molecular Dynamics Simulations of Gram-Negative Bacterial Cell Envelopes

Bacteria are protected by complex molecular architectures known as the cell envelope. The cell envelope is composed of regions with distinct chemical compositions and physical properties, namely, membranes and a cell wall. To develop novel antibiotics to combat pathogenic bacteria, molecular level knowledge of the structure, dynamics, and interplay between the chemical components of the cell envelope that surrounds bacterial cells is imperative. In addition, conserved molecular patterns associated with the bacterial envelope are recognized by receptors as part of the mammalian defensive response to infection, and an improved understanding of bacteria-host interactions would facilitate the search for novel immunotherapeutics. This Perspective introduces an emerging area of computational biology: multiscale molecular dynamics simulations of chemically complex models of bacterial lipids and membranes. We discuss progress to date, and identify areas for future development that will enable the study of aspects of the membrane components that are as yet unexplored by computational methods.

Allosteric modulation model of the mu opioid receptor by herkinorin, a potent not alkaloidal agonist

Modulation of opioid receptors is the primary choice for pain management and structural information studies have gained new horizons with the recently available X-ray crystal structures. Herkinorin is one of the most remarkable salvinorin A derivative with high affinity for the mu opioid receptor, moderate selectivity and lack of nitrogen atoms on its structure. Surprisingly, binding models for herkinorin are lacking. In this work, we explore binding models of herkinorin using automated docking, molecular dynamics simulations, free energy calculations and available experimental information. Our herkinorin D-ICM-1 binding model predicted a binding free energy of -11.52 ± 1.14 kcal mol^-1 by alchemical free energy estimations, which is close to the experimental values -10.91 ± 0.2 and -10.80 ± 0.05 kcal mol^-1 and is in agreement with experimental structural information. Specifically, D-ICM-1 molecular dynamics simulations showed a water-mediated interaction between D-ICM-1 and the amino acid H297^6.52, this interaction coincides with the co-crystallized ligands. Another relevant interaction, with N127^2.63, allowed to rationalize herkinorin's selectivity to mu over delta opioid receptors. Our suggested binding model for herkinorin is in agreement with this and additional experimental data. The most remarkable observation derived from our D-ICM-1 model is that herkinorin reaches an allosteric sodium ion binding site near N150^3.35. Key interactions in that region appear relevant for the lack of β-arrestin recruitment by herkinorin. This interaction is key for downstream signaling pathways involved in the development of side effects, such as tolerance. Future SAR studies and medicinal chemistry efforts will benefit from the structural information presented in this work.

Characterizing Structural Stability of Amyloid Motif Fibrils Mediated by Water Molecules

In biological systems, structural confinements of amyloid fibrils can be mediated by the role of water molecules. However, the underlying effect of the dynamic behavior of water molecules on structural stabilities of amyloid fibrils is still unclear. By performing molecular dynamics simulations, we investigate the dynamic features and the effect of interior water molecules on conformations and mechanical characteristics of various amyloid fibrils. We find that a specific mechanism induced by the dynamic properties of interior water molecules can affect diffusion of water molecules inside amyloid fibrils, inducing their different structural stabilities. The conformation of amyloid fibrils induced by interior water molecules show the fibrils' different mechanical features. We elucidate the role of confined and movable interior water molecules in structural stabilities of various amyloid fibrils. Our results offer insights not only in further understanding of mechanical features of amyloids as mediated by water molecules, but also in the fine-tuning of the functional abilities of amyloid fibrils for applications.

Toward Understanding the Molecular Recognition of Albumin by p53-Activating Stapled Peptide ATSP-7041

Reactivation of tumor-suppressing activity of p53 protein by targeting its negative regulator MDM2/MDMX has been pursued as a potential anticancer strategy. A promising dual inhibitor of MDM2/MDMX that has been developed and is currently in clinical trials is the stapled peptide ATSP-7041. The activity of this molecule is reported to be modulated in the presence of serum. Albumin is the most abundant protein in serum and is known to bind reversibly to several molecules. To study this interaction, we develop a protocol combining molecular modeling, docking, and simulations. Exhaustive docking of the peptide with representative simulated structures of human serum albumin led to the identification of probable binding sites on the surface of the protein, including both known canonical and novel binding sites. Sequence differences at putative peptide-binding sites in human and mouse albumin result in differing interaction energies with the peptide and enable us to rationalize the observed differences in vivo. In general, the findings should help in guiding the design of features in such peptides that may affect their distribution and cell permeability, opening a new window in structure-guided design strategies.

Determinants of Macromolecular Specificity from Proteomics-Derived Peptide Substrate Data

BACKGROUND

Recent advances in proteomics methodologies allow for high throughput profiling of proteolytic cleavage events. The resulting substrate peptide distributions provide deep insights in the underlying macromolecular recognition events, as determinants of biomolecular specificity identified by proteomics approaches may be compared to structure-based analysis of corresponding protein-protein interfaces.

METHOD

Here, we present an overview of experimental and computational methodologies and tools applied in the area and provide an outlook beyond the protein class of proteases.

RESULTS AND CONCLUSION

We discuss here future potential, synergies and needs of the emerging overlap disciplines of proteomics and structure-based modelling.

A role for loop G in the β1 strand in GABAA receptor activation

KEY POINTS

The role of the β1 strand in GABAA receptor function is unclear. It lies anti-parallel to the β2 strand, which is known to participate in receptor activation. Molecular dynamics simulation revealed solvent accessible residues within the β1 strand of the GABAA β3 homopentamer that might be amenable to analysis using the substituted Cys accessibility method. Cys substitutions from Asp43 to Thr47 in the GABAA α1 subunit showed that D43C and T47C reduced the apparent potency of GABA. F45C caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Cys43 and Cys47 were accessible to 2-aminoethyl methanethiosulphonate (MTSEA) modification, whereas Cys45 was not. Both GABA and the allosteric agonist propofol reduced MTSEA modification of Cys43 and Cys47. By contrast, modification of Cys64 in the β2 strand loop D was impeded by GABA but unaffected by propofol. These data reveal movement of β1 strand loop G residues during agonist activation of the GABAA receptor.

ABSTRACT

The GABAA receptor α subunit β1 strand runs anti-parallel to the β2 strand, which contains loop D, known to participate in receptor activation and agonist binding. However, a role for the β1 strand has yet to be established. We used molecular dynamics simulation to quantify the solvent accessible surface area (SASA) of β1 strand residues in the GABAA β3 homopentamer structure. Residues in the complementary interface equivalent to those between Asp43 and Thr47 in the α1 subunit have an alternating pattern of high and low SASA consistent with a β strand structure. We investigated the functional role of these β1 strand residues in the α1 subunit by individually replacing them with Cys residues. D43C and T47C substitutions reduced the apparent potency of GABA at α1β2γ2 receptors by 50-fold and eight-fold, respectively, whereas the F45C substitution caused a biphasic GABA concentration-response relationship and increased spontaneous gating. Receptors with D43C or T47C substitutions were sensitive to 2-aminoethyl methanethiosulphonate (MTSEA) modification. However, GABA-evoked currents mediated by α1(F45C)β2γ2 receptors were unaffected by MTSEA, suggesting that this residue is inaccessible. Both GABA and the allosteric agonist propofol reduced MTSEA modification of α1(D43C)β2γ2 and α1(T47C)β2γ2 receptors, indicating movement of the β1 strand even during allosteric activation. This is in contrast to α1(F64C)β2γ2 receptors, where only GABA, but not propofol, reduced MTSEA modification. These findings provide the first functional evidence for movement of the β1 strand during gating of the receptor and identify residues that are critical for maintaining GABAA receptor function.

Voltage Dependence of Conformational Dynamics and Subconducting States of VDAC-1

The voltage-dependent anion channel 1 (VDAC-1) is an important protein of the outer mitochondrial membrane that transports energy metabolites and is involved in apoptosis. The available structures of VDAC proteins show a wide β-stranded barrel pore, with its N-terminal α-helix (N-α) bound to its interior. Electrophysiology experiments revealed that voltage, its polarity, and membrane composition modulate VDAC currents. Experiments with VDAC-1 mutants identified amino acids that regulate the gating process. However, the mechanisms for how these factors regulate VDAC-1, and which changes they trigger in the channel, are still unknown. In this study, molecular dynamics simulations and single-channel experiments of VDAC-1 show agreement for the current-voltage relationships of an "open" channel and they also show several subconducting transient states that are more cation selective in the simulations. We observed voltage-dependent asymmetric distortions of the VDAC-1 barrel and the displacement of particular charged amino acids. We constructed conformational models of the protein voltage response and the pore changes that consistently explain the protein conformations observed at opposite voltage polarities, either in phosphatidylethanolamine or phosphatidylcholine membranes. The submicrosecond VDAC-1 voltage response shows intrinsic structural changes that explain the role of key gating amino acids and support some of the current gating hypotheses. These voltage-dependent protein changes include asymmetric barrel distortion, its interaction with the membrane, and significant displacement of N-α amino acids.

Molecular Properties of Globin Channels and Pores: Role of Cholesterol in Ligand Binding and Movement

Globins contain one or more cavities that control or affect such functions as ligand movement and ligand binding. Here we report that the extended globin family [cytoglobin (Cygb); neuroglobin (Ngb); myoglobin (Mb); hemoglobin (Hb) subunits Hba(α); and Hbb(β)] contain either a transmembrane (TM) helix or pore-lining region as well as internal cavities. Protein motif/domain analyses indicate that Ngb and Hbb each contain 5 cholesterol- binding (CRAC/CARC) domains and 1 caveolin binding motif, whereas the Cygb dimer has 6 cholesterol-binding domains but lacks caveolin-binding motifs. Mb and Hba each exhibit 2 cholesterol-binding domains and also lack caveolin-binding motifs. The Hb αβ-tetramer contains 14 cholesterol-binding domains. Computer algorithms indicate that Cygb and Ngb cavities display multiple partitions and C-terminal pore-lining regions, whereas Mb has three major cavities plus a C-terminal pore-lining region. The Hb tetramer exhibits a large internal cavity but the subunits differ in that they contain a C-terminal TM helix (Hba) and pore-lining region (Hbb). The cavities include 43 of 190 Cygb residues, 38 of 151 of Ngb residues, 55 of 154 Mb residues, and 137 of 688 residues in the Hb tetramer. Each cavity complex includes 6 to 8 residues of the TM helix or pore-lining region and CRAC/CARC domains exist within all cavities. Erythrocyte Hb αβ-tetramers are largely cytosolic but also bind to a membrane anion exchange protein, "band 3," which contains a large internal cavity and 12 TM helices (5 being pore-lining regions). The Hba TM helix may be the erythrocyte membrane "band 3" attachment site. "Band 3" contributes 4 caveolin binding motifs and 10 CRAC/CARC domains. Cholesterol binding may create lipid-disordered phases that alter globin cavities and facilitate ligand movement, permitting ion channel formation and conformational changes that orchestrate anion and ligand (O2, CO2, NO) movement within the large internal cavities and channels of the globins.

Structure and Function of the Su(H)-Hairless Repressor Complex, the Major Antagonist of Notch Signaling in Drosophila melanogaster

Notch is a conserved signaling pathway that specifies cell fates in metazoans. Receptor-ligand interactions induce changes in gene expression, which is regulated by the transcription factor CBF1/Su(H)/Lag-1 (CSL). CSL interacts with coregulators to repress and activate transcription from Notch target genes. While the molecular details of the activator complex are relatively well understood, the structure-function of CSL-mediated repressor complexes is poorly defined. In Drosophila, the antagonist Hairless directly binds Su(H) (the fly CSL ortholog) to repress transcription from Notch targets. Here, we determine the X-ray structure of the Su(H)-Hairless complex bound to DNA. Hairless binding produces a large conformational change in Su(H) by interacting with residues in the hydrophobic core of Su(H), illustrating the structural plasticity of CSL molecules to interact with different binding partners. Based on the structure, we designed mutants in Hairless and Su(H) that affect binding, but do not affect formation of the activator complex. These mutants were validated in vitro by isothermal titration calorimetry and yeast two- and three-hybrid assays. Moreover, these mutants allowed us to solely characterize the repressor function of Su(H) in vivo.

The transcription factor CSL regulates gene expression in response to Notch pathway signaling. The X-ray structure of the complex between the fruit fly version of CSL, Su(H), and its antagonist, Hairless, reveals a novel binding mode and unanticipated structural plasticity.

Notch signaling is a form of cell-to-cell communication, in which extracellular receptor-ligand interactions ultimately result in changes in gene expression. The Notch pathway is highly conserved from the model organism Drosophila melanogaster to humans. When mutations occur within Notch pathway components, this often leads to human disease, such as certain types of cancers and birth defects. Transcription of Notch target genes is regulated by the transcription factor CSL (for CBF1/RBP-J in mammals, Su(H) in Drosophila, and Lag-1 in Caenorhabditis elegans). CSL functions as both a transcriptional activator and repressor by forming complexes with coactivator and corepressor proteins, respectively. Here we determine the high-resolution X-ray structure of Su(H) (the fly CSL ortholog) in complex with the corepressor Hairless, which is the major antagonist of Notch signaling in Drosophila. The structure unexpectedly reveals that Hairless binding results in a dramatic conformational change in Su(H). In parallel, we designed mutations in Su(H) and Hairless based on our structure and showed that these mutants are defective in complex formation in vitro and display functional deficiencies in in vivo assays. Taken together, our work provides significant molecular insights into how CSL functions as a transcriptional repressor in the Notch pathway.

Identification of an activation site in Bak and mitochondrial Bax triggered by antibodies

During apoptosis, Bak and Bax are activated by BH3-only proteins binding to the α2–α5 hydrophobic groove; Bax is also activated via a rear pocket. Here we report that antibodies can directly activate Bak and mitochondrial Bax by binding to the α1–α2 loop. A monoclonal antibody (clone 7D10) binds close to α1 in non-activated Bak to induce conformational change, oligomerization, and cytochrome c release. Anti-FLAG antibodies also activate Bak containing a FLAG epitope close to α1. An antibody (clone 3C10) to the Bax α1–α2 loop activates mitochondrial Bax, but blocks translocation of cytosolic Bax. Tethers within Bak show that 7D10 binding directly extricates α1; a structural model of the 7D10 Fab bound to Bak reveals the formation of a cavity under α1. Our identification of the α1–α2 loop as an activation site in Bak paves the way to develop intrabodies or small molecules that directly and selectively regulate these proteins.

During apoptosis, Bak and Bax are activated by BH3-only proteins binding to a specific hydrophobic groove. Here, the authors show that antibodies can also activate Bak and mitochondrial Bax by binding to the α1-α2 loop, thus identifying a potential clinical target.

Protein Binding Pocket Dynamics

The dynamics of protein binding pockets are crucial for their interaction specificity. Structural flexibility allows proteins to adapt to their individual molecular binding partners and facilitates the binding process. This implies the necessity to consider protein internal motion in determining and predicting binding properties and in designing new binders. Although accounting for protein dynamics presents a challenge for computational approaches, it expands the structural and physicochemical space for compound design and thus offers the prospect of improved binding specificity and selectivity. A cavity on the surface or in the interior of a protein that possesses suitable properties for binding a ligand is usually referred to as a binding pocket. The set of amino acid residues around a binding pocket determines its physicochemical characteristics and, together with its shape and location in a protein, defines its functionality. Residues outside the binding site can also have a long-range effect on the properties of the binding pocket. Cavities with similar functionalities are often conserved across protein families. For example, enzyme active sites are usually concave surfaces that present amino acid residues in a suitable configuration for binding low molecular weight compounds. Macromolecular binding pockets, on the other hand, are located on the protein surface and are often shallower. The mobility of proteins allows the opening, closing, and adaptation of binding pockets to regulate binding processes and specific protein functionalities. For example, channels and tunnels can exist permanently or transiently to transport compounds to and from a binding site. The influence of protein flexibility on binding pockets can vary from small changes to an already existent pocket to the formation of a completely new pocket. Here, we review recent developments in computational methods to detect and define binding pockets and to study pocket dynamics. We introduce five different classes of protein pocket dynamics: (1) appearance/disappearance of a subpocket in an existing pocket; (2) appearance/disappearance of an adjacent pocket on the protein surface in the direct vicinity of an already existing pocket; (3) pocket breathing, which may be caused by side-chain fluctuations or backbone or interdomain vibrational motion; (4) opening/closing of a channel or tunnel, connecting a pocket inside the protein with solvent, including lid motion; and (5) the appearance/disappearance of an allosteric pocket at a site on a protein distinct from an already existing pocket with binding of a ligand to the allosteric binding site affecting the original pocket. We suggest that the class of pocket dynamics, as well as the type and extent of protein motion affecting the binding pocket, should be factors considered in choosing the most appropriate computational approach to study a given binding pocket. Furthermore, we examine the relationship between pocket dynamics classes and induced fit, conformational selection, and gating models of ligand binding on binding kinetics and thermodynamics. We discuss the implications of protein binding pocket dynamics for drug design and conclude with potential future directions for computational analysis of protein binding pocket dynamics.

Nanovoid formation and mechanics: a comparison of poly(dicyclopentadiene) and epoxy networks from molecular dynamics simulations

Protective equipment in civilian and military applications requires the use of polymer materials that are both stiff and tough over a wide range of strain rates. However, typical structural materials, like tightly cross-linked epoxies, are very brittle. Recent experiments demonstrated that cross-linked poly(dicyclopentadiene) (pDCPD) networks can circumvent this trade-off by providing structural properties such as a high glass transition temperature and glassy modulus, while simultaneously exhibiting excellent toughness and high-rate impact resistance. The greater performance of pDCPD was attributed to more facile plastic deformation and nano-scale void formation, but the chemical and structural mechanisms underlying this response were not clear. Here, we use atomistic molecular dynamics to compare the molecular- and chain-level properties of pDCPD and epoxy networks undergoing high strain rate deformation. We quantify the tensile modulus and yield strength of the networks as well as the prevalence and characteristics of nanovoids that form during deformation. Networks of similar molecular weight between cross-links are compared. Two key molecular-level properties are identified - monomer flexibility and polar chemistry - that influence the behavior of the networks. Increasing monomer flexibility reduces the modulus and yield strength, while strong non-covalent interactions (e.g., hydrogen bonds) that accompany polar moieties provide higher modulus and yield strength. The lack of strong non-covalent interactions in pDCPD was found to account for its lower modulus and yield strength compared to the epoxies. We examine the molecular-level properties of nanovoids, such as shape, alignment, and local stress distribution, as well as the local chemical environment, finding that nanovoid formation and growth are increased by the monomer rigidity but decreased by polar chemistry. As a result, the pDCPD network, which has a stiff chain backbone with nonpolar alkane chemistry, exhibits more and larger nanovoids that grow more readily during deformation, which could account for the higher toughness and more ductile behavior observed in pDCPD.

Selective inhibitors of trypanosomal uridylyl transferase RET1 establish druggability of RNA post-transcriptional modifications

Non-coding RNAs are crucial regulators for a vast array of cellular processes and have been implicated in human disease. These biological processes represent a hitherto untapped resource in our fight against disease. In this work we identify small molecule inhibitors of a non-coding RNA uridylylation pathway. The TUTase family of enzymes is important for modulating non-coding RNA pathways in both human cancer and pathogen systems. We demonstrate that this new class of drug target can be accessed with traditional drug discovery techniques. Using the Trypanosoma brucei TUTase, RET1, we identify TUTase inhibitors and lay the groundwork for the use of this new target class as a therapeutic opportunity for the under-served disease area of African Trypanosomiasis. In a broader sense this work demonstrates the therapeutic potential for targeting RNA post-transcriptional modifications with small molecules in human disease.

Molecular principle of the cyclin-dependent kinase selectivity of 4-(thiazol-5-yl)-2-(phenylamino) pyrimidine-5-carbonitrile derivatives revealed by molecular modeling studies

Due to the high sequence identity of the binding pockets of cyclin-dependent kinases (CDKs), designing highly selective inhibitors towards a specific CDK member remains a big challenge.

Analysis of water channels by molecular dynamics simulation of heterotetrameric sarcosine oxidase

A precise 100-ns molecular dynamics simulation in aquo was performed for the heterotetrameric sarcosine oxidase bound with a substrate analogue, dimethylglycine. The spatial region including the protein was divided into small rectangular cells. The average number of the water molecules locating within each cell was calculated based on the simulation trajectory. The clusters of the cells filled with water molecules were used to determine the water channels. The narrowness of the channels, the average hydropathy indices of the residues of the channels, and the number of migration events of water molecules through the channels were consistent with the selective transport hypothesis whereby tunnel T3 is the pathway for the exit of the iminium intermediate of the enzyme reaction.

Multiple Binding Poses in the Hydrophobic Cavity of Bee Odorant Binding Protein AmelOBP14

In the first step of olfaction, odorants are bound and solubilized by small globular odorant binding proteins (OBPs) which shuttle them to the membrane of a sensory neuron. Low ligand affinity and selectivity at this step enable the recognition of a wide range of chemicals. Honey bee Apis mellifera’s OBP14 (AmelOBP14) binds different plant odorants in a largely hydrophobic cavity. In long molecular dynamics simulations in the presence and absence of ligand eugenol, we observe a highly dynamic C-terminal region which forms one side of the ligand-binding cavity, and the ligand drifts away from its crystallized orientation. Hamiltonian replica exchange simulations, allowing exchanges of conformations sampled by the real ligand with those sampled by a noninteracting dummy molecule and several intermediates, suggest an alternative, quite different ligand pose which is adopted immediately and which is stable in long simulations. Thermodynamic integration yields binding free energies which are in reasonable agreement with experimental data.

Energetics of Endotoxin Recognition in the Toll-Like Receptor 4 Innate Immune Response

Bacterial outer membrane lipopolysaccharide (LPS) potently stimulates the mammalian innate immune system, and can lead to sepsis, the primary cause of death from infections. LPS is sensed by Toll-like receptor 4 (TLR4) in complex with its lipid-binding coreceptor MD-2, but subtle structural variations in LPS can profoundly modulate the response. To better understand the mechanism of LPS-induced stimulation and bacterial evasion, we have calculated the binding affinity to MD-2 of agonistic and antagonistic LPS variants including lipid A, lipid IVa, and synthetic antagonist Eritoran, and provide evidence that the coreceptor is a molecular switch that undergoes ligand-induced conformational changes to appropriately activate or inhibit the receptor complex. The plasticity of the coreceptor binding cavity is shown to be essential for distinguishing between ligands, whilst similar calculations for a model bacterial LPS bilayer reveal the "membrane-like" nature of the protein cavity. The ability to predict the activity of LPS variants should facilitate the rational design of TLR4 therapeutics.

The Structural Basis for Activation and Inhibition of ZAP-70 Kinase Domain

ZAP–70 (Zeta-chain-associated protein kinase 70) is a tyrosine kinase that interacts directly with the activated T-cell receptor to transduce downstream signals, and is hence a major player in the regulation of the adaptive immune response. Dysfunction of ZAP–70 causes selective T cell deficiency that in turn results in persistent infections. ZAP–70 is activated by a variety of signals including phosphorylation of the kinase domain (KD), and binding of its regulatory tandem Src homology 2 (SH2) domains to the T cell receptor. The present study investigates molecular mechanisms of activation and inhibition of ZAP–70 via atomically detailed molecular dynamics simulation approaches. We report microsecond timescale simulations of five distinct states of the ZAP–70 KD, comprising apo, inhibited and three phosphorylated variants. Extensive analysis of local flexibility and correlated motions reveal crucial transitions between the states, thus elucidating crucial steps in the activation mechanism of the ZAP–70 KD. Furthermore, we rationalize previously observed staurosporine-bound crystal structures, suggesting that whilst the KD superficially resembles an “active-like” conformation, the inhibitor modulates the underlying protein dynamics and restricts it in a compact, rigid state inaccessible to ligands or cofactors. Finally, our analysis reveals a novel, potentially druggable pocket in close proximity to the activation loop of the kinase, and we subsequently use its structure in fragment-based virtual screening to develop a pharmacophore model. The pocket is distinct from classical type I or type II kinase pockets, and its discovery offers promise in future design of specific kinase inhibitors, whilst mutations in residues associated with this pocket are implicated in immunodeficiency in humans.

ZAP–70 is a key protein kinase in the adaptive immune system. It is essential for development and function of T cells and natural killer cells, and associated mutations can lead to conditions such as severe combined immunodeficiency (SCID). Here, simulations of the ZAP–70 kinase domain are used to study its dynamics in response to different mechanistic signals. We identify crucial motions over microsecond timescales, which help to rationalize in atomic detail previous structural and experimental data regarding its biological regulation. We subsequently propose a scheme for the phosphorylation-dependent activation cascade of ZAP–70, and for its ligand-dependent inhibition. Finally, we characterize a novel cryptic pocket adjacent to the active site and activation loop, which is distinct from classical type I or type II kinase sites. The pocket is in close proximity to several residues whose mutations cause SCID in humans, and its identification offers promise in future drug design efforts.

Cheminformatics Research at the Unilever Centre for Molecular Science Informatics Cambridge

The Centre for Molecular Informatics, formerly Unilever Centre for Molecular Science Informatics (UCMSI), at the University of Cambridge is a world-leading driving force in the field of cheminformatics. Since its opening in 2000 more than 300 scientific articles have fundamentally changed the field of molecular informatics. The Centre has been a key player in promoting open chemical data and semantic access. Though mainly focussing on basic research, close collaborations with industrial partners ensured real world feedback and access to high quality molecular data. A variety of tools and standard protocols have been developed and are ubiquitous in the daily practice of cheminformatics. Here, we present a retrospective of cheminformatics research performed at the UCMSI, thereby highlighting historical and recent trends in the field as well as indicating future directions.

Modelling ligand selectivity of serine proteases using integrative proteochemometric approaches improves model performance and allows the multi-target dependent interpretation of features

Serine proteases, implicated in important physiological functions, have a high intra-family similarity, which leads to unwanted off-target effects of inhibitors with insufficient selectivity. However, the availability of sequence and structure data has now made it possible to develop approaches to design pharmacological agents that can discriminate successfully between their related binding sites. In this study, we have quantified the relationship between 12,625 distinct protease inhibitors and their bioactivity against 67 targets of the serine protease family (20,213 data points) in an integrative manner, using proteochemometric modelling (PCM). The benchmarking of 21 different target descriptors motivated the usage of specific binding pocket amino acid descriptors, which helped in the identification of active site residues and selective compound chemotypes affecting compound affinity and selectivity. PCM models performed better than alternative approaches (models trained using exclusively compound descriptors on all available data, QSAR) employed for comparison with R(2)/RMSE values of 0.64 ± 0.23/0.66 ± 0.20 vs. 0.35 ± 0.27/1.05 ± 0.27 log units, respectively. Moreover, the interpretation of the PCM model singled out various chemical substructures responsible for bioactivity and selectivity towards particular proteases (thrombin, trypsin and coagulation factor 10) in agreement with the literature. For instance, absence of a tertiary sulphonamide was identified to be responsible for decreased selective activity (by on average 0.27 ± 0.65 pChEMBL units) on FA10. Among the binding pocket residues, the amino acids (arginine, leucine and tyrosine) at positions 35, 39, 60, 93, 140 and 207 were observed as key contributing residues for selective affinity on these three targets.

Structural glycobiology of human α1-acid glycoprotein and its implications for pharmacokinetics and inflammation

Human α1-acid glycoprotein (AGP) is an abundant human plasma glycoprotein that may be N-glycosylated at five positions. AGP plays important roles on pharmacokinetics and can rise up to 5-fold in inflammatory events. In such events, the glycan chains attached to Asn54, Asn75 and Asn85 may become fucosylated, originating a sialyl-Lewis X epitope. This epitope, in turn, can bind selectin proteins. Such interplay is important for immunomodulation. While the X-ray structure of unglycosylated AGP has been reported, the absence of the glycan chains hampered the further insights into its structural biology and, ultimately, into its biological function. Thus, the current work intends to contribute in the characterization of the structural glycobiology and function of AGP by building a structural model of its fully glycosylated form, taking into account the different glycoforms that are found in vivo. The obtained data points to the absence of a major influence of glycosylation on AGP's secondary structure, in agreement with crystallography observations. However, the glycan chains seem able to interfere with the protein dynamics, mainly at the AGP-ligand-binding site, indicating a possible role in its complexation to drugs and other bioactive compounds. By examining the influence of fucosylation on AGP structure and binding to selectins, it is proposed that the latter may bind to glycan chains linked to Asn54 and Asn75, and that this binding may involve other glycans, such as the one attached to Asn15. These results point to an increased participation of carbohydrates on the observed AGP roles in pharmacokinetics and inflammation.

Room temperature crystal structure of the fast switching M159T mutant of the fluorescent protein dronpa

The fluorescent protein Dronpa undergoes reversible photoswitching reactions between the bright "on" and dark "off" states via photoisomerization and proton transfer reactions. We report the room temperature crystal structure of the fast switching Met159Thr mutant of Dronpa at 2.0-Å resolution in the bright on state. Structural differences with the wild type include shifted backbone positions of strand β8 containing Thr159 as well as an altered A-C dimer interface involving strands β7, β8, β10, and β11. The Met159Thr mutation increases the cavity volume for the p-hydroxybenzylidene-imidazolinone chromophore as a result of both the side chain difference and the backbone positional differences.

Understanding the thermostability and activity of Bacillus subtilis lipase mutants: insights from molecular dynamics simulations

Improving the thermostability of industrial enzymes is an important protein engineering challenge. Point mutations, induced to increase thermostability, affect the structure and dynamics of the target protein in several ways and thus can also affect its activity. There appears to be no general rules for improving the thermostabilty of enzymes without adversely affecting their enzymatic activity. We report MD simulations, of wild type Bacillus subtilis lipase (WT) and its six progressively thermostable mutants (2M, 3M, 4M, 6M, 9M, and 12M), performed at different temperatures, to address this issue. Less thermostable mutants (LTMs), 2M to 6M, show WT-like dynamics at all simulation temperatures. However, the two more thermostable mutants (MTMs) show the required flexibility at appropriate temperature ranges and maintain conformational stability at high temperature. They show a deep and rugged free-energy landscape, confining them within a near-native conformational space by conserving noncovalent interactions, and thus protecting them from possible aggregation. In contrast, the LTMs having marginally higher thermostabilities than WT show greater probabilities of accessing non-native conformations, which, due to aggregation, have reduced possibilities of reverting to their respective native states under refolding conditions. Our analysis indicates the possibility of nonadditive effects of point mutations on the conformational stability of LTMs.

Epock: rapid analysis of protein pocket dynamics

Summary: The volume of an internal protein pocket is fundamental to ligand accessibility. Few programs that compute such volumes manage dynamic data from molecular dynamics (MD) simulations. Limited performance often prohibits analysis of large datasets. We present Epock, an efficient command-line tool that calculates pocket volumes from MD trajectories. A plugin for the VMD program provides a graphical user interface to facilitate input creation, run Epock and analyse the results.

Availability and implementation: Epock C++ source code, Python analysis scripts, VMD Tcl plugin, documentation and installation instructions are freely available at http://epock.bitbucket.org.

Contact:benoist.laurent@gmail.com or baaden@smplinux.de

Supplementary information:Supplementary data are available at Bioinformatics online.

POVME 2.0: An Enhanced Tool for Determining Pocket Shape and Volume Characteristics

Analysis of macromolecular/small-molecule binding pockets can provide important insights into molecular recognition and receptor dynamics. Since its release in 2011, the POVME (POcket Volume MEasurer) algorithm has been widely adopted as a simple-to-use tool for measuring and characterizing pocket volumes and shapes. We here present POVME 2.0, which is an order of magnitude faster, has improved accuracy, includes a graphical user interface, and can produce volumetric density maps for improved pocket analysis. To demonstrate the utility of the algorithm, we use it to analyze the binding pocket of RNA editing ligase 1 from the unicellular parasite Trypanosoma brucei, the etiological agent of African sleeping sickness. The POVME analysis characterizes the full dynamics of a potentially druggable transient binding pocket and so may guide future antitrypanosomal drug-discovery efforts. We are hopeful that this new version will be a useful tool for the computational- and medicinal-chemist community.