A Retrospective on the Development of Methods for the Analysis of Protein Conformational Ensembles.

Analysing protein conformational ensembles whether from molecular dynamics (MD) simulation or other sources for functionally relevant conformational changes can be very challenging. In the nineteen nineties dimensional reduction methods were developed primarily for analysing MD trajectories to determine dominant motions with the aim of understanding their relationship to function. Coarse-graining methods were also developed so the conformational change between two structures could be described in terms of the relative motion of a small number of quasi-rigid regions rather than in terms of a large number of atoms. When these methods are combined, they can characterize the large-scale motions inherent in a conformational ensemble providing insight into possible functional mechanism. The dimensional reduction methods first applied to protein conformational ensembles were referred to as Quasi-Harmonic Analysis, Principal Component Analysis and Essential Dynamics Analysis. A retrospective on the origin of these methods is presented, the relationships between them explained, and more recent developments reviewed.

Tspan6 stimulates the chemoattractive potential of breast cancer cells for B cells in an EV- and LXR-dependent manner.

The immune microenvironment in breast cancer (BCa) is controlled by a complex network of communication between various cell types. Here, we find that recruitment of B lymphocytes to BCa tissues is controlled via mechanisms associated with cancer cell-derived extracellular vesicles (CCD-EVs). Gene expression profiling identifies the Liver X receptor (LXR)-dependent transcriptional network as a key pathway that controls both CCD-EVs-induced migration of B cells and accumulation of B cells in BCa tissues. The increased accumulation oxysterol ligands for LXR (i.e., 25-hydroxycholesterol and 27-hydroxycholesterol) in CCD-EVs is regulated by the tetraspanin 6 (Tspan6). Tspan6 stimulates the chemoattractive potential of BCa cells for B cells in an EV- and LXR-dependent manner. These results demonstrate that tetraspanins control intercellular trafficking of oxysterols via CCD-EVs. Furthermore, tetraspanin-dependent changes in the oxysterol composition of CCD-EVs and the LXR signaling axis play a key role in specific changes in the tumor immune microenvironment.

Characterization of the Immune Microenvironment in Inflammatory Breast Cancer Using Multiplex Immunofluorescence.

INTRODUCTION: Inflammatory breast cancer (IBC) is an aggressive form of breast cancer with a poorly characterized immune microenvironment. METHODS: We used a five-colour multiplex immunofluorescence panel, including CD68, CD4, CD8, CD20, and FOXP3 for immune microenvironment profiling in 93 treatment-naïve IBC samples. RESULTS: Lower grade tumours were characterized by decreased CD4+ cells but increased accumulation of FOXP3+ cells. Increased CD20+ cells correlated with better response to neoadjuvant chemotherapy and increased CD4+ cells infiltration correlated with better overall survival. Pairwise analysis revealed that both ER+ and triple-negative breast cancer were characterized by co-infiltration of CD20 + cells with CD68+ and CD4+ cells, whereas co-infiltration of CD8+ and CD68+ cells was only observed in HER2+ IBC. Co-infiltration of CD20+, CD8+, CD4+, and FOXP3+ cells, and co-existence of CD68+ with FOXP3+ cells correlated with better therapeutic responses, while resistant tumours were characterized by co-accumulation of CD4+, CD8+, FOXP3+, and CD68+ cells and co-expression of CD68+ and CD20+ cells. In a Cox regression model, response to therapy was the most significant factor associated with improved patient survival. CONCLUSION: Those results reveal a complex unique pattern of distribution of immune cell subtypes in IBC and provide an important basis for detailed characterization of molecular pathways that govern the formation of IBC immune landscape and potential for immunotherapy.

Interactive Flexible-Receptor Molecular Docking in Virtual Reality Using DockIT.

Interactive docking enables the user to guide and control the docking of two biomolecules into a binding pose. It is of particular use when the binding site is known and is thought to be applicable to structure-based drug design (SBDD) and educating students about biomolecular interactions. For SBDD, it enables expertise and intuition to be brought to bear in the drug design process. In education, it can teach students about the most basic level of biomolecular function. Here, we introduce DockIT for virtual reality (VR) that uses a VR headset and hand-held controllers. Using the method of linear response on explicit solvent molecular dynamics simulations, DockIT can model both global and local conformational changes within the receptor due to forces of interaction with the ligand. It has real-time flexible molecular surface rendering and can show the real-time formation and breaking of hydrogen bonds, both between the ligand and receptor and within the receptor itself as it smoothly changes conformation.

The role of the half-turn in determining structures of Alzheimer's Aß wild-type and mutants.

Half-turns are shown to be the main determinants of many experimental Alzheimer's Aß fibril structures. Fibril structures contain three half-turn types, ßαRß, ßαLß and ß뵧 which each result in a ∼90° bend in a ß-strand. It is shown that only these half-turns enable cross-ß stacking and thus the right-angle fold seen in fibrils is an intrinsic feature of cross-ß. Encoding a strand as a conformational sequence in ß, αR, αL and ε(ßL), pairwise combination rules for consecutive half-turns are used to decode this sequence to give the backbone path. This reveals how structures would be dramatically affected by a deletion. Using a wild-type Aß(42) fibril structure and the pairwise combination rules, the Osaka deletion is predicted to result in exposure of surfaces that are mutually shielding from the solvent. Molecular dynamics simulations on an 11-mer ß-sheet of Alzheimer's Aß(40) of the Dutch (E22Q), Iowa (D23N), Arctic (E22G), and Osaka (E22Δ) mutants, show the crucial role glycine plays in the positioning of ßαRß half-turns. Their "in-phase" positions along the sequence in the wild-type, Dutch mutant and Iowa mutant means that the half-folds all fold to the same side creating the same closed structure. Their out-of-phase positions in Arctic and Osaka mutants creates a flatter structure in the former and an S-shape structure in the latter which, as predicted, exposes surfaces on the inside in the closed wild-type to the outside. This is consistent with the gain of interaction model and indicates how domain swapping might explain the Osaka mutant's unique properties.

Determination of amino acids that favour the αL region using Ramachandran propensity plots. Implications for α-sheet as the possible amyloid intermediate.

In amyloid diseases an insoluble amyloid fibril forms via a soluble oligomeric intermediate. It is this intermediate that mediates toxicity and it has been suggested, somewhat controversially, that it has the α-sheet structure. Nests and α-strands are similar peptide motifs in that alternate residues lie in the αR and γL regions of the Ramachandran plot for nests, or αR and αL regions for α-strands. In nests a concavity is formed by the main chain NH atoms whereas in α-strands the main chain is almost straight. Using "Ramachandran propensity plots" to focus on the αL/γL region, it is shown that glycine favours γL (82% of amino acids are glycine), but disfavours αL (3% are glycine). Most charged and polar amino acids favour αL with asparagine having by far the highest propensity. Thus, glycine favours nests but, contrary to common expectation, should not favour α-sheet. By contrast most charged or polar amino acids should favour α-sheet by their propensity for the αL conformation, which is more discriminating amongst amino acids than the αR conformation. Thus, these results suggest the composition of sequences that favour α-sheet formation and point towards effective prediction of α-sheet from sequence.

DockIT: a tool for interactive molecular docking and molecular complex construction.

SUMMARY: DockIT is a tool that has a unique set of physical and graphical features for interactive molecular docking. It enables the user to bring a ligand and a receptor into a docking pose by controlling relative position and orientation, either with a mouse and keyboard, or with a haptic device. Atomic interactions are modelled using molecular dynamics-based force-fields with the force on the ligand being felt on a haptic device. Real-time calculation and display of intermolecular hydrogen bonds and multipoint collision detection either using maximum force or maximum atomic overlap, mean that together with the ability to monitor selected intermolecular atomic distances, the user can find physically feasible docking poses that satisfy distance constraints derived from experimental methods. With these features and the ability to output and reload docked structures it can be used to accurately build up large multi-component molecular systems in preparation for molecular dynamics simulation. AVAILABILITY AND IMPLEMENTATION: DockIT is available free of charge for non-commercial use at http://www.haptimol.co.uk/downloads.htm. It requires a windows computer with GPU that supports OpenCL 1.2 and OpenGL 4.0. It may be used with a mouse and keyboard, or a haptic device from 3DSystems.

Investigation of sequence features of hinge-bending regions in proteins with domain movements using kernel logistic regression.

BACKGROUND: Hinge-bending movements in proteins comprising two or more domains form a large class of functional movements. Hinge-bending regions demarcate protein domains and collectively control the domain movement. Consequently, the ability to recognise sequence features of hinge-bending regions and to be able to predict them from sequence alone would benefit various areas of protein research. For example, an understanding of how the sequence features of these regions relate to dynamic properties in multi-domain proteins would aid in the rational design of linkers in therapeutic fusion proteins. RESULTS: The DynDom database of protein domain movements comprises sequences annotated to indicate whether the amino acid residue is located within a hinge-bending region or within an intradomain region. Using statistical methods and Kernel Logistic Regression (KLR) models, this data was used to determine sequence features that favour or disfavour hinge-bending regions. This is a difficult classification problem as the number of negative cases (intradomain residues) is much larger than the number of positive cases (hinge residues). The statistical methods and the KLR models both show that cysteine has the lowest propensity for hinge-bending regions and proline has the highest, even though it is the most rigid amino acid. As hinge-bending regions have been previously shown to occur frequently at the terminal regions of the secondary structures, the propensity for proline at these regions is likely due to its tendency to break secondary structures. The KLR models also indicate that isoleucine may act as a domain-capping residue. We have found that a quadratic KLR model outperforms a linear KLR model and that improvement in performance occurs up to very long window lengths (eighty residues) indicating long-range correlations. CONCLUSION: In contrast to the only other approach that focused solely on interdomain hinge-bending regions, the method provides a modest and statistically significant improvement over a random classifier. An explanation of the KLR results is that in the prediction of hinge-bending regions a long-range correlation is at play between a small number amino acids that either favour or disfavour hinge-bending regions. The resulting sequence-based prediction tool, HingeSeek, is available to run through a webserver at hingeseek.cmp.uea.ac.uk.

The CD151-midkine pathway regulates the immune microenvironment in inflammatory breast cancer.

The immune microenvironment in inflammatory breast cancer (IBC) is poorly characterised, and molecular and cellular pathways that control accumulation of various immune cells in IBC tissues remain largely unknown. Here, we discovered a novel pathway linking the expression of the tetraspanin protein CD151 in tumour cells with increased accumulation of macrophages in cancerous tissues. It is notable that elevated expression of CD151 and a higher number of tumour-infiltrating macrophages correlated with better patient responses to chemotherapy. Accordingly, CD151-expressing IBC xenografts were characterised by the increased infiltration of macrophages. In vitro migration experiments demonstrated that CD151 stimulates the chemoattractive potential of IBC cells for monocytes via mechanisms involving midkine (a heparin-binding growth factor), integrin α6ß1, and production of extracellular vesicles (EVs). Profiling of chemokines secreted by IBC cells demonstrated that CD151 increases production of midkine. Purified midkine specifically stimulated migration of monocytes, but not other immune cells. Further experiments demonstrated that the chemoattractive potential of IBC-derived EVs is blocked by anti-midkine antibodies. These results demonstrate for the first time that changes in the expression of a tetraspanin protein by tumour cells can affect the formation of the immune microenvironment by modulating recruitment of effector cells to cancerous tissues. Therefore, a CD151-midkine pathway can be considered as a novel target for controlled changes of the immune landscape in IBC. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Haptic-Assisted Interactive Molecular Docking Incorporating Receptor Flexibility.

Haptic-assisted interactive docking tools immerse the user in an environment where intuition and knowledge can be used to help guide the docking process. Here we present such a tool where the user "holds" a rigid ligand via a haptic device through which they feel interaction forces with a flexible receptor biomolecule. To ensure forces transmitted through the haptic device are smooth and stable, they must be updated at a rate greater than 500 Hz. Due to this time constraint, the majority of haptic docking tools do not attempt to model the conformational changes that would occur when molecules interact during binding. Our haptic-assisted docking tool, "Haptimol FlexiDock", models a receptor's conformational response to forces of interaction with a ligand while maintaining the required haptic refresh rate. In order to model receptor flexibility we use the method of linear response for which we determine the variance-covariance matrix of atomic fluctuations from the trajectory of an explicit-solvent molecular dynamics simulation of the ligand-free receptor molecule. The key to satisfying the time constraint is an eigenvector decomposition of the variance-covariance matrix which enables a good approximation to the conformational response of the receptor to be calculated rapidly. This exploits a feature of protein dynamics whereby most fluctuation occurs within a relatively small subspace. The method is demonstrated on glutamine binding protein in interaction with glutamine and maltose binding protein in interaction with maltose. For both proteins the movement that occurs when the ligand is docked near to its binding site matches the experimentally determined movement well. It is thought that this tool will be particularly useful for structure-based drug design.

Free Energy Profile of Domain Movement in Ligand-Free Citrate Synthase.

Citrate synthase plays a fundamental role in the metabolic cycle of the cell. Its catalytic mechanism is complex involving the binding of two substrates that cause a domain movement. In this paper, we used classical molecular dynamics simulations and umbrella-sampling simulations to determine the potential of mean force along a reaction coordinate for the domain movement in ligand-free citrate synthase from pig ( Sus scrofa). The results show that, at 293 K, the closed-domain conformation has a ∼4 kb T higher energy than the open-domain conformation. In a simple two-state model, this difference means that the enzyme spends 98% of the time in the open-domain conformation ready to receive the substrate, oxaloacetate, rather than the closed-domain conformation where the binding site would be inaccessible to the substrate. Given that experimental evidence indicates that the binding of oxaloacetate induces at least partial closure, this would imply an induced-fit mechanism which we argue is applicable to all enzymes with a functional domain movement for reasons of catalytic efficiency. A barrier of 4 kb T gives an estimation of the mean first passage time in the range 1-10 µs.

Multi-strand ß-sheet of Alzheimer Aß(1-40) folds to ß-strip helix: implication for protofilament formation.

X-ray fiber diffraction experiments on Alzheimer Aß(1-40) fibrils formed in an assembly process thought to simulate a portion of the pathophysiological process in Alzheimer's disease, indicated protofilaments with tilted ß-strands rather than strands oriented perpendicular to the fibril axis as is usually interpreted from cross-ß patterns. The protofilament width and tilt angle determined by these experiments were used to predict a ß-strip helix model - a ß-helix-like structure in which multiple identical polypeptide molecules assemble in-register to form a helical sheet structure such that the outer strands 1 and m join with a register shift t - with m = 11 and t = 22. Starting from untwisted ß-sheets comprising 10, 11, and 12 strands, multiple explicit solvent molecular dynamics (MD) simulations were performed to determine whether the sheets form ß-strip helices matching the dimensions of the experimentally measured protofilament. In the simulations, the predicted 11-strand sheets curled up to form a closed ß-strip helix-like structure with dimensions matching experimental values, whereas the 10- and 12-strand sheets did not form a closed helical structure. The 12-strand structure did, however, show similarity to a cross-ß structure determined by a solid-state NMR experiment. The 11-strand ß-strip helix resembles a trans-membrane ß-barrel which could explain the ability of small oligomers of Aß(1-40) to form toxic ion channels. A further consequence of opposite sides of the 11-strand strip coming together at a register shift of 22 is end-to-end joins between neighboring ß-strip helices, resulting in a protofilament that keeps growing in both directions. Communicated by Ramaswamy H. Sarma.

Methodological improvements for the analysis of domain movements in large biomolecular complexes.

Domain movements play a prominent role in the function of many biomolecules such as the ribosome and F0F1-ATP synthase. As more structures of large biomolecules in different functional states become available as experimental techniques for structure determination advance, there is a need to develop methods to understand the conformational changes that occur. DynDom and DynDom3D were developed to analyse two structures of a biomolecule for domain movements. They both used an original method for domain recognition based on clustering of "rotation vectors". Here we introduce significant improvements in both the methodology and implementation of a tool for the analysis of domain movements in large multimeric biomolecules. The main improvement is in the recognition of domains by using all six degrees of freedom required to describe the movement of a rigid body. This is achieved by way of Chasles' theorem in which a rigid-body movement can be described as a screw movement about a unique axis. Thus clustering now includes, in addition to rotation vector data, screw-axis location data and axial climb data. This improves both the sensitivity of domain recognition and performance. A further improvement is the recognition and annotation of interdomain bending regions, something not done for multimeric biomolecules in DynDom3D. This is significant as it is these regions that collectively control the domain movement. The new stand-alone, platform-independent implementation, DynDom6D, can analyse biomolecules comprising protein, DNA and RNA, and employs an alignment method to automatically achieve the required equivalence of atoms in the two structures.

Morphing and docking visualisation of biomolecular structures using Multi-Dimensional Scaling.

Protein structures are often solved at atomic resolution in two states defining a functional movement but intervening conformations are usually unknown. Morphing methods generate intervening conformations between two known structures. When viewed as an animation using molecular graphics, a smooth, direct morph enables the eye to track changes in structure that might be otherwise missed. We present a morphing method that aims to linearly interpolate interatomic distances and which uses SMACOF (Scaling by MAjorisation of COmplicated Function) and multigrid techniques with a cut-off distance based weighting that optimizes the MolProbity score of intervening structures. The all-atom morphs are smooth, move directly between the two structures, and are shown, in general, to pass closer to a set of known intermediates than those generated using other methods. The techniques are also used for docking by putting the unbound structures in a "near-approach pose" and then morphing to the bound complex. The resulting GPU-accelerated tools are available on a webserver, Morphit_Pro, at http://morphit-pro.cmp.uea.ac.uk/ and more than 5000 domains movements available at the DynDom website can now be viewed as morphs http://morphit-pro.cmp.uea.ac.uk/dyndom/.

High quality rendering of protein dynamics in space filling mode.

Producing high quality depictions of molecular structures has been an area of academic interest for years, with visualisation tools such as UCSF Chimera, Yasara and PyMol providing a huge number of different rendering modes and lighting effects. However, no visualisation program supports per-pixel lighting effects with shadows whilst rendering a molecular trajectory in space filling mode. In this paper, a new approach to rendering high quality visualisations of molecular trajectories is presented. To enhance depth, ambient occlusion is included within the render. Shadows are also included to help the user perceive relative motions of parts of the protein as they move based on their trajectories. Our approach requires a regular grid to be constructed every time the molecular structure deforms allowing per-pixel lighting effects and ambient occlusion to be rendered every frame, at interactive refresh rates. Two different regular grids are investigated, a fixed grid and a memory efficient compact grid. The algorithms used allow trajectories of proteins comprising of up to 300,000 atoms in size to be rendered at ninety frames per second on a desktop computer using the GPU for general purpose computations. Regular grid construction was found to only take up a small proportion of the total time to render a frame. It was found that despite being slower to construct, the memory efficient compact grid outperformed the theoretically faster fixed grid when the protein being rendered is large, owing to its more efficient memory access patterns. The techniques described could be implemented in other molecular rendering software.

Geometrical principles of homomeric ß-barrels and ß-helices: Application to modeling amyloid protofilaments.

Examples of homomeric ß-helices and ß-barrels have recently emerged. Here we generalize the theory for the shear number in ß-barrels to encompass ß-helices and homomeric structures. We introduce the concept of the "ß-strip," the set of parallel or antiparallel neighboring strands, from which the whole helix can be generated giving it n-fold rotational symmetry. In this context, the shear number is interpreted as the sum around the helix of the fixed register shift between neighboring identical ß-strips. Using this approach, we have derived relationships between helical width, pitch, angle between strand direction and helical axis, mass per length, register shift, and number of strands. The validity and unifying power of the method is demonstrated with known structures including α-hemolysin, T4 phage spike, cylindrin, and the HET-s(218-289) prion. From reported dimensions measured by X-ray fiber diffraction on amyloid fibrils, the relationships can be used to predict the register shift and the number of strands within amyloid protofilaments. This was used to construct models of transthyretin and Alzheimer ß(40) amyloid protofilaments that comprise a single strip of in-register ß-strands folded into a "ß-strip helix." Results suggest both stabilization of an individual ß-strip helix and growth by addition of further ß-strip helices can involve the same pair of sequence segments associating with ß-sheet hydrogen bonding at the same register shift. This process would be aided by a repeat sequence. Hence, understanding how the register shift (as the distance between repeat sequences) relates to helical dimensions will be useful for nanotube design.

Virtual Environment for Studying the Docking Interactions of Rigid Biomolecules with Haptics.

Haptic technology facilitates user interaction with the virtual world via the sense of touch. In molecular docking, haptics enables the user to sense the interaction forces during the docking process. Here we describe a haptics-assisted interactive software tool, called Haptimol_RD, for the study of docking interactions. By utilizing GPU-accelerated proximity querying methods very large systems can now be studied. Methods for force scaling, multipoint collision response and haptic navigation are described that address force stability issues that are particular to the interactive docking of large systems. Thus, Haptimol_RD expands, for the first time, the use of interactive biomolecular haptics to the study of protein-protein interactions. Unlike existing approaches, Haptimol_RD is designed to run on relatively inexpensive consumer-level hardware and is freely available to the community.

Determination of locked interfaces in biomolecular complexes using Haptimol_RD.

Interactive haptics-assisted docking provides a virtual environment for the study of molecular complex formation. It enables the user to interact with the virtual molecules, experience the interaction forces via their sense of touch, and gain insights about the docking process itself. Here we use a recently developed haptics software tool, Haptimol_RD, for the rigid docking of protein subunits to form complexes. Dimers, both homo and hetero, are loaded into the software with their subunits separated in space for the purpose of assessing whether they can be brought back into the correct docking pose via rigid-body movements. Four dimers were classified into two types: two with an interwinding subunit interface and two with a non-interwinding subunit interface. It was found that the two with an interwinding interface could not be docked whereas the two with the non-interwinding interface could be. For the two that could be docked a "sucking" effect could be felt on the haptic device when the correct binding pose was approached which is associated with a minimum in the interaction energy. It is clear that for those that could not be docked, the conformation of one or both of the subunits must change upon docking. This leads to the steric-based concept of a locked or non-locked interface. Non-locked interfaces have shapes that allow the subunits to come together or come apart without the necessity of intra-subunit conformational change, whereas locked interfaces require a conformational change within one or both subunits for them to be able to come apart.

The DynDom3D Webserver for the Analysis of Domain Movements in Multimeric Proteins.

DynDom3D is a program for the analysis of domain movements in multimeric proteins. Its inputs are two structure files that indicate a possible domain movement, but the onus has been on the user to process the files so that there is the necessary one-to-one equivalence between atoms in the two atom lists. This is often a prohibitive task to carry out manually, which has limited the application of DynDom3D. Here we report on a webserver with a preprocessor that automatically creates an equivalence between atoms using sequence alignment methods. The processed structure files are passed to DynDom3D and the results are presented on a webpage that includes molecular graphics for easy visualization.

Monte Carlo Sampling with Linear Inverse Kinematics for Simulation of Protein Flexible Regions.

A Monte Carlo linear inverse-kinematics method for the simulation of protein chains with fixed ends is introduced. It includes backbone bond-angle bending and simultaneous loop and ring closure to allow full proline ring flexibility. An obstacle to linear null-space methods is the eventual drift of the end group. Maintenance of the end group at its initial position by occasional reset is performed in a way that is consistent with the overall methodology and minimally disruptive to the current conformation. The implementation permitted multiple rigid regions within the chain, enabling the simulation of domain movements where domains are rigid bodies connected by flexible interdomain regions. The method was tested on polyalanine, polyglycine, loop 6 of triosephosphate isomerase, and glutamine binding protein. Simulations of glutamine binding protein, where only 11 of the 226 residues at the interdomain bending regions were flexible, accurately reproduced the experimentally determined domain movement.