Mechanism of glycan receptor recognition and specificity switch for avian, swine, and human adapted influenza virus hemagglutinins: a molecular dynamics perspective.

Hemagglutinins (HA's) from duck, swine, and human influenza viruses have previously been shown to prefer avian and human glycan receptor analogues with distinct topological profiles, pentasaccharides LSTa (alpha-2,3 linkage) and LSTc (alpha-2,6 linkage), in comparative molecular dynamics studies. On the basis of detailed analyses of the dynamic motions of the receptor binding domains (RBDs) and interaction energy profiles with individual glycan residues, we have identified approximately 30 residue positions in the RBD that present distinct profiles with the receptor analogues. Glycan binding constrained the conformational space sampling by the HA. Electrostatic steering appeared to play a key role in glycan binding specificity. The complex dynamic behaviors of the major SSE and trimeric interfaces with or without bound glycans suggested that networks of interactions might account for species specificity in these low affinity and high avidity (multivalent) interactions between different HA and glycans. Contact frequency, energetic decomposition, and H-bond analyses revealed species-specific differences in HA-glycan interaction profiles, not readily discernible from crystal structures alone. Interaction energy profiles indicated that mutation events at the set of residues such as 145, 156, 158, and 222 would favor human or avian receptor analogues, often through interactions with distal asialo-residues. These results correlate well with existing experimental evidence, and suggest new opportunities for simulation-based vaccine and drug development

Comprehensive structural model of the mechanochemical cycle of a mitotic motor highlights molecular adaptations in the kinesin family

Kinesins are responsible for a wide variety of microtubule-based, ATP-dependent functions. Their motor domain drives these activities but the molecular adaptations that specify these diverse and essential cellular activities are poorly understood. It has been assumed that the first identified kinesin - the transport motor kinesin-1 – is the mechanistic paradigm for the entire superfamily, but accumulating evidence suggests that this is not the case. To address the deficits in our understanding of the molecular basis of functional divergence within the kinesin superfamily, we studied kinesin-5s, which are essential mitotic motors whose inhibition blocks cell division. Using cryo-electron microscopy and subnanometer resolution structure determination, we have visualised conformations of microtubule-bound human kinesin-5 motor domain at successive steps in its ATPase cycle. Following ATP hydrolysis, nucleotide-dependent conformational changes in the active site are allosterically propagated into rotations of the motor domain and uncurling of the drugbinding loop L5. In addition, the mechanical neck-linker element that is crucial for motor stepping undergoes discrete, ordered displacements. We also observed large reorientations of the motor N-terminus that indicate its importance for kinesin-5 function through control of neck-linker conformation. A kinesin-5 mutant lacking this N-terminus is enzymatically active, and ATP-dependent neck-linker movement and motility is defective although not ablated. All these aspects of kinesin-5 mechanochemistry are distinct from kinesin-1. Our findings directly demonstrate the regulatory role of the kinesin-5 N-terminus in collaboration with the motor’s structured neck-linker, and highlight the multiple adaptations within kinesin motor domains that tune their mechanochemistries according to distinct functional requirements

CRANKITE: a fast polypeptide backbone conformation sampler

Background: CRANKITE is a suite of programs for simulating backbone conformations of polypeptides and proteins. The core of the suite is an efficient Metropolis Monte Carlo sampler of backbone conformations in continuous three-dimensional space in atomic details. Methods: In contrast to other programs relying on local Metropolis moves in the space of dihedral angles, our sampler utilizes local crankshaft rotations of rigid peptide bonds in Cartesian space. Results: The sampler allows fast simulation and analysis of secondary structure formation and conformational changes for proteins of average length

An electrostatic mechanism for Ca(2+)-mediated regulation of gap junction channels.

Gap junction channels mediate intercellular signalling that is crucial in tissue development, homeostasis and pathologic states such as cardiac arrhythmias, cancer and trauma. To explore the mechanism by which Ca(2+) blocks intercellular communication during tissue injury, we determined the X-ray crystal structures of the human Cx26 gap junction channel with and without bound Ca(2+). The two structures were nearly identical, ruling out both a large-scale structural change and a local steric constriction of the pore. Ca(2+) coordination sites reside at the interfaces between adjacent subunits, near the entrance to the extracellular gap, where local, side chain conformational rearrangements enable Ca(2+)chelation. Computational analysis revealed that Ca(2+)-binding generates a positive electrostatic barrier that substantially inhibits permeation of cations such as K(+) into the pore. Our results provide structural evidence for a unique mechanism of channel regulation: ionic conduction block via an electrostatic barrier rather than steric occlusion of the channel pore

Conformational and functional analysis of molecular dynamics trajectories by Self-Organising Maps

<p>Abstract</p> <p>Background</p> <p>Molecular dynamics (MD) simulations are powerful tools to investigate the conformational dynamics of proteins that is often a critical element of their function. Identification of functionally relevant conformations is generally done clustering the large ensemble of structures that are generated. Recently, Self-Organising Maps (SOMs) were reported performing more accurately and providing more consistent results than traditional clustering algorithms in various data mining problems. We present a novel strategy to analyse and compare conformational ensembles of protein domains using a two-level approach that combines SOMs and hierarchical clustering.</p> <p>Results</p> <p>The conformational dynamics of the α-spectrin SH3 protein domain and six single mutants were analysed by MD simulations. The Cα's Cartesian coordinates of conformations sampled in the essential space were used as input data vectors for SOM training, then complete linkage clustering was performed on the SOM prototype vectors. A specific protocol to optimize a SOM for structural ensembles was proposed: the optimal SOM was selected by means of a Taguchi experimental design plan applied to different data sets, and the optimal sampling rate of the MD trajectory was selected. The proposed two-level approach was applied to single trajectories of the SH3 domain independently as well as to groups of them at the same time. The results demonstrated the potential of this approach in the analysis of large ensembles of molecular structures: the possibility of producing a topological mapping of the conformational space in a simple 2D visualisation, as well as of effectively highlighting differences in the conformational dynamics directly related to biological functions.</p> <p>Conclusions</p> <p>The use of a two-level approach combining SOMs and hierarchical clustering for conformational analysis of structural ensembles of proteins was proposed. It can easily be extended to other study cases and to conformational ensembles from other sources.</p

Kinetic distance and kinetic maps from molecular dynamics simulation

Characterizing macromolecular kinetics from molecular dynamics (MD) simulations requires a distance metric that can distinguish slowly-interconverting states. Here we build upon diffusion map theory and define a kinetic distance for irreducible Markov processes that quantifies how slowly molecular conformations interconvert. The kinetic distance can be computed given a model that approximates the eigenvalues and eigenvectors (reaction coordinates) of the MD Markov operator. Here we employ the time-lagged independent component analysis (TICA). The TICA components can be scaled to provide a kinetic map in which the Euclidean distance corresponds to the kinetic distance. As a result, the question of how many TICA dimensions should be kept in a dimensionality reduction approach becomes obsolete, and one parameter less needs to be specified in the kinetic model construction. We demonstrate the approach using TICA and Markov state model (MSM) analyses for illustrative models, protein conformation dynamics in bovine pancreatic trypsin inhibitor and protein-inhibitor association in trypsin and benzamidine

Encounter complexes and dimensionality reduction in protein-protein association

An outstanding challenge has been to understand the mechanism whereby proteins associate. We report here the results of exhaustively sampling the conformational space in protein–protein association using a physics-based energy function. The agreement between experimental intermolecular paramagnetic relaxation enhancement (PRE) data and the PRE profiles calculated from the docked structures shows that the method captures both specific and non-specific encounter complexes. To explore the energy landscape in the vicinity of the native structure, the nonlinear manifold describing the relative orientation of two solid bodies is projected onto a Euclidean space in which the shape of low energy regions is studied by principal component analysis. Results show that the energy surface is canyon-like, with a smooth funnel within a two dimensional subspace capturing over 75% of the total motion. Thus, proteins tend to associate along preferred pathways, similar to sliding of a protein along DNA in the process of protein-DNA recognition

Evolutionary fine-tuning of conformational ensembles in FimH during host-pathogen interactions

Positive selection in the two-domain type 1 pilus adhesin FimH enhances Escherichia coli fitness in urinary tract infection (UTI). We report a comprehensive atomic-level view of FimH in two-state conformational ensembles in solution, composed of one low-affinity tense (T) and multiple high-affinity relaxed (R) conformations. Positively selected residues allosterically modulate the equilibrium between these two conformational states, each of which engages mannose through distinct binding orientations. A FimH variant that only adopts the R state is severely attenuated early in a mouse model of uncomplicated UTI but is proficient at colonizing catheterized bladders in vivo or bladder transitional-like epithelial cells in vitro. Thus, the bladder habitat has barrier(s) to R state–mediated colonization possibly conferred by the terminally differentiated bladder epithelium and/or decoy receptors in urine. Together, our studies reveal the conformational landscape in solution, binding mechanisms, and adhesive strength of an allosteric two-domain adhesin that evolved “moderate” affinity to optimize persistence in the bladder during UTI