Simulating nanoscale functional motions of biomolecules

We are describing efficient dynamics simulation methods for the characterization of functional motion of biomolecules on the nanometer scale. Multivariate statistical methods are widely used to extract and enhance functional collective motions from molecular dynamics (MD) simulations. A dimension reduction in MD is often realized through a principal component analysis (PCA) or a singular value decomposition (SVD) of the trajectory. Normal mode analysis (NMA) is a related collective coordinate space approach, which involves the decomposition of the motion into vibration modes based on an elastic model. Using the myosin motor protein as an example we describe a hybrid technique termed amplified collective motions (ACM) that enhances sampling of conformational space through a combination of normal modes with atomic level MD. Unfortunately, the forced orthogonalization of modes in collective coordinate space leads to complex dependencies that are not necessarily consistent with the symmetry of biological macromolecules and assemblies. In many biological molecules, such as HIV-1 protease, reflective or rotational symmetries are present that are broken using standard orthogonal basis functions. We present a method to compute the plane of reflective symmetry or the axis of rotational symmetry from the trajectory frames. Moreover, we develop an SVD that best approximates the given trajectory while respecting the symmetry. Finally, we describe a local feature analysis (LFA) to construct a topographic representation of functional dynamics in terms of local features. The LFA representations are low-dimensional, and provide a reduced basis set for collective motions, but unlike global collective modes they are sparsely distributed and spatially localized. This yields a more reliable assignment of essential dynamics modes across different MD time windows.     

Double-labeled rabies virus: live tracking of enveloped virus transport

Here we describe a strategy to fluorescently label the envelope of rabies virus (RV), of the Rhabdoviridae family, in order to track the transport of single enveloped viruses in living cells. Red fluorescent proteins (tm-RFP) were engineered to comprise the N-terminal signal sequence and C-terminal transmembrane spanning and cytoplasmic domain sequences of the RV glycoprotein (G). Two variants of tm-RFP were transported to and anchored in the cell surface membrane, independent of glycosylation. As shown by confocal microscopy, tm-RFP colocalized at the cell surface with the RV matrix and G protein and was incorporated into G gene-deficient virus particles. Recombinant RV expressing the membrane-anchored tm-RFP in addition to G yielded infectious viruses with mosaic envelopes containing both tm-RFP and G. Viable double-labeled virus particles comprising a red fluorescent envelope and a green fluorescent ribonucleoprotein were generated by expressing in addition an enhanced green fluorescent protein-phosphoprotein fusion construct (S. Finke, K. Brzozka, and K. K. Conzelmann, J. Virol. 78:12333-12343, 2004). Individual enveloped virus particles were observed under live cell conditions as extracellular particles and inside endosomal vesicles. Importantly, double-labeled RVs were transported in the retrograde direction over long distances in neurites of in vitro-differentiated NS20Y neuroblastoma cells. This indicates that the typical retrograde axonal transport of RV to the central nervous system involves neuronal transport vesicles in which complete enveloped RV particles are carried as a cargo.     

Functional motions of influenza virus hemagglutinin: a structure-based analytical approach.

Influenza virus hemagglutinin (HA), a homotrimeric integral membrane glycoprotein essential for viral infection, is engaged in two biological functions: recognition of target cells' receptor proteins and fusion of viral and endosomal membranes, both requiring substantial conformational flexibility from the part of the glycoprotein. The different modes of collective motions underlying the functional mobility/adaptability of the protein are determined in the present study using an extension of the Gaussian network model (GNM) to treat concerted anisotropic motions. We determine the molecular mechanisms that may underlie HA function, along with the structural regions or residues whose mutations are expected to impede function. Good agreement between theoretically predicted fluctuations of individual residues and corresponding x-ray crystallographic temperature factors is found, which lends support to the GNM elucidation of the conformational dynamics of HA by focusing upon a subset of dominant modes. The lowest frequency mode indicates a global torsion of the HA trimer about its longitudinal axis, accompanied by a substantial mobility at the viral membrane connection. This mode is proposed to constitute the dominant molecular mechanism for the translocation and aggregation of HAs, and for the opening and dilation of the fusion pore. The second and third collective modes indicate a global bending, allowing for a large lateral surface exposure, which is likely to facilitate the close association of the viral and endosomal membranes before pore opening. The analysis of kinetically hot residues, in contrast, reveals a localization of energy centered around the HA2 residue Asp112, which apparently triggers the solvent exposure of the fusion peptide.     

iGNM: a database of protein functional motions based on Gaussian Network Model

MOTIVATION: The knowledge of protein structure is not sufficient for understanding and controlling its function. Function is a dynamic property. Although protein structural information has been rapidly accumulating in databases, little effort has been invested to date toward systematically characterizing protein dynamics. The recent success of analytical methods based on elastic network models, and in particular the Gaussian Network Model (GNM), permits us to perform a high-throughput analysis of the collective dynamics of proteins. RESULTS: We computed the GNM dynamics for 20 058 structures from the Protein Data Bank, and generated information on the equilibrium dynamics at the level of individual residues. The results are stored on a web-based system called iGNM and configured so as to permit the users to visualize or download the results through a standard web browser using a simple search engine. Static and animated images for describing the conformational mobility of proteins over a broad range of normal modes are accessible, along with an online calculation engine available for newly deposited structures. A case study of the dynamics of 20 non-homologous hydrolases is presented to illustrate the utility of the iGNM database for identifying key residues that control the cooperative motions and revealing the connection between collective dynamics and catalytic activity.     

Plant membrane transport.

Recent developments in plant membrane transport, particularly concerning the vacuolar and plasma membranes, have increased our understanding of molecular aspects of primary pumps, carrier systems and ion channels.     

Interconversion of functional motions between mesophilic and thermophilic adenylate kinases

Dynamic properties are functionally important in many proteins, including the enzyme adenylate kinase (AK), for which the open/closed transition limits the rate of catalytic turnover. Here, we compare our previously published coarse-grained (double-well Gō) simulation of mesophilic AK from E. coli (AKmeso) to simulations of thermophilic AK from Aquifex aeolicus (AKthermo). In AKthermo, as with AKmeso, the LID domain prefers to close before the NMP domain in the presence of ligand, but LID rigid-body flexibility in the open (O) ensemble decreases significantly. Backbone foldedness in O and/or transition state (TS) ensembles increases significantly relative to AKmeso in some interdomain backbone hinges and within LID. In contact space, the TS of AKthermo has fewer contacts at the CORE-LID interface but a stronger contact network surrounding the CORE-NMP interface than the TS of AKmeso. A "heated" simulation of AKthermo at 375K slightly increases LID rigid-body flexibility in accordance with the "corresponding states" hypothesis. Furthermore, while computational mutation of 7 prolines in AKthermo to their AKmeso counterparts produces similar small perturbations, mutation of these sites, especially positions 8 and 155, to glycine is required to achieve LID rigid-body flexibility and hinge flexibilities comparable to AKmeso. Mutating the 7 sites to proline in AKmeso reduces some hinges' flexibilities, especially hinge 2, but does not reduce LID rigid-body flexibility, suggesting that these two types of motion are decoupled in AKmeso. In conclusion, our results suggest that hinge flexibility and global functional motions alike are correlated with but not exclusively determined by the hinge residues. This mutational framework can inform the rational design of functionally important flexibility and allostery in other proteins toward engineering novel biochemical pathways.     

Plant hormones: Ins and outs of auxin transport

Regulated transport has long been known to play a key part in action of the plant hormone auxin. Now, at last, a family of auxin efflux carriers has been identified, and the characterisation of one family member has provided strong evidence in support of models that have been proposed to explain gravitropic curvature in roots.     

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.     

Plant cell biology: the ABC of monolignol transport

Lignins are complex aromatic heteropolymers that reinforce the cell walls of terrestrial plants. A new study identifies an ATP-binding cassette ABC transporter that pumps a monolignol lignin precursor across the plasma membrane.     

Plant functional types and climatic change: Introduction

Abstract. Plant functional types are a necessary device for reducing the complex and often uncharted characteristics of species diversity in function and structure when attempting to project the nature and function of species assemblages into future environments. A workshop was held to review the current methods commonly used for defining plant functional types, either globally or for particular biomes, and to compare them with the field experiences of specialists for specific biomes of the world. The methods fall into either an objective and inductive approach or a subjective and deductive approach. When the different methods were tested, it was generally found that the classification for one site or environment was not wholly applicable to a different site or environment. However, the degree of change which is necessary for adjustment between environments may not prove to be a major limitation in the use of functional types.     

Comparing apples with oranges: on the functional equivalence of food products for comparative LCAs

Purpose

In this article, we present an innovative way of deriving comparable functional systems for comparative life cycle assessments (LCAs) of food products. We define the functional unit as the contribution of one or more foods to the nutrient composition of a weekly diet and, after a product substitution, employ a product system expansion approach to search for an alternative set of products which provides an equivalent nutritional composition.

Methods

Replacement is regarded within the context of a weekly diet. The comparable diet is a solution to a linear problem which finds the diet that is most similar to the starting one, subject to nutritional and/or other constraints that guarantee a minimum dietary quality. The formulation gives priority to selecting food products according to popularity.

Results

We illustrate our method with two examples. We show that a baseline diet containing 3.6 servings of apples a week is equivalent to a similar diet in which the apples are replaced with 3.6 servings of oranges and servings of strawberry and kiwi are removed. These changes are necessary mainly because of differences in the vitamin C content between apples and oranges. The second example is a replacement of all meat in a weekly diet by a soy-based meat substitute. In this case, additional fish products need to be consumed to make up for a lack of selenium and essential amino acids.

Conclusions

We present an innovative and objective way to overcome the challenge of comparing two or more food products in a comparative LCA. Our approach is systematic and finds the alternative diet that best meets the nutritional criteria as well as reflecting the food preferences of the population. The method selects products according to the role they play in the dietary pattern. Moreover, the method is flexible enough to allow for different selection criteria and other nutritional and non-nutritional constraints.     

Predicting the functional motions of p97 using symmetric normal modes

p97 is a protein complex of the AAA+ family. Although functions of p97 are well understood, the mechanism by which p97 performs its unfolding activities remains unclear. In this work, we present a novel way of applying normal mode analysis to study this six‐fold symmetric molecular machine. By selecting normal modes that are axial symmetric and give the largest movements at D1 or D2 pore residues, we are able to predict the functional motions of p97, which are then validated by experimentally observed conformational changes. Our results shed light and provide new understandings on several key steps of the p97 functional process that were previously unclear or controversial, and thus are able to reconcile multiple previous findings. Specifically, our results reveal that (i) a venous valve‐like mechanism is used at D2 pore to ensure a one‐way exit‐only traffic of substrates; (ii) D1 pore remains shut during the functional process; (iii) the “swing‐up” motion of the N domain is closely coupled with the vertical motion of the D1 pore along the pore axis; (iv) because of the shut D1 pore and the one‐way traffic at D2 pore, it is highly likely that substrates enter the chamber through the gaps at the D1/D2 interface. The limited chamber volume inside p97 suggests that a substrate may be pulling out from D2 while at the same time being pulling in at the interface; (v) lastly, p97 uses a series of actions that alternate between twisting and pulling to remove the substrate. Proteins 2016; 84:1823–1835. © 2016 Wiley Periodicals, Inc.     

Modeling conformational ensembles of slow functional motions in Pin1-WW

Protein-protein interactions are often mediated by flexible loops that experience conformational dynamics on the microsecond to millisecond time scales. NMR relaxation studies can map these dynamics. However, defining the network of inter-converting conformers that underlie the relaxation data remains generally challenging. Here, we combine NMR relaxation experiments with simulation to visualize networks of inter-converting conformers. We demonstrate our approach with the apo Pin1-WW domain, for which NMR has revealed conformational dynamics of a flexible loop in the millisecond range. We sample and cluster the free energy landscape using Markov State Models (MSM) with major and minor exchange states with high correlation with the NMR relaxation data and low NOE violations. These MSM are hierarchical ensembles of slowly interconverting, metastable macrostates and rapidly interconverting microstates. We found a low population state that consists primarily of holo-like conformations and is a "hub" visited by most pathways between macrostates. These results suggest that conformational equilibria between holo-like and alternative conformers pre-exist in the intrinsic dynamics of apo Pin1-WW. Analysis using MutInf, a mutual information method for quantifying correlated motions, reveals that WW dynamics not only play a role in substrate recognition, but also may help couple the substrate binding site on the WW domain to the one on the catalytic domain. Our work represents an important step towards building networks of inter-converting conformational states and is generally applicable.