Modeling supramolecular assemblages

The Rac1 regulator ELMO1 controls vascular morphogenesis in zebrafish

RATIONALE

Angiogenesis is regulated by the small GTPase Rac1. The ELMO1/DOCK180 complex forms a guanine nucleotide exchange factor for Rac1, regulating its activation during cell migration in different biological systems.

OBJECTIVE

To investigate the function of ELMO1/DOCK180 in vascular development.

METHODS AND RESULTS

In situ hybridization studies for elmo1 identified a vascular and neuronal expression in zebrafish. Morpholino-based expression silencing of elmo1 severely impaired the formation of the vasculature, including intersomitic vessels, the dorsal longitudinal anastomotic vessel, the parachordal vessel, and the development of the thoracic duct in tg(fli1:EGFP) embryos. Mechanistically, we identified Netrin-1 and its receptor Unc5B as upstream activators of the ELMO1/DOCK180 complex, regulating its functional interaction and leading to Rac1 activation in endothelial cells and vessel formation in zebrafish.

CONCLUSIONS

Our data have identified a novel signaling cascade regulating vasculature formation in zebrafish.

Guidance of vascular development: lessons from the nervous system

The vascular system of vertebrates consists of an organized, branched network of arteries, veins, and capillaries that penetrates all the tissues of the body. One of the most striking features of the vascular system is that its branching pattern is highly stereotyped, with major and secondary branches forming at specific sites and developing highly conserved organ-specific vascular patterns. The factors controlling vascular patterning are not yet completely understood. Recent studies have highlighted the anatomic and structural similarities between blood vessels and nerves. The 2 networks are often aligned, with nerve fibers and blood vessels following parallel routes. Furthermore, both systems require precise control over their guidance and growth. Several molecules with attractive and repulsive properties have been found to modulate the proper guidance of both nerves and blood vessels. These include the Semaphorins, the Slits, and the Netrins and their receptors. In this review, we describe the molecular mechanisms by which blood vessels and axons achieve proper path finding and the molecular cues that are involved in their guidance.

Direct computer simulation of water-mediated force between supported phospholipid membranes

The grand canonical Monte Carlo technique is used to calculate the water-mediated force operating between two supported 1,2-dilauroyl-DL-phosphatidylethanolamine (DLPE) membranes in the short separation range. The intra- and intermolecular interactions in the system are described with a combination of an AMBER-based force field for DLPE and a TIP4P model for water. The long range contributions to the electrostatic interaction energy are treated in the dipole-dipole group-based approximation. The total water-mediated force is analyzed in terms of its hydration component and the component due to the direct interaction between the membranes. The latter is, in addition, partitioned into the electrostatic, van der Waals, and steric repulsion contributions to give an idea of their relative significance in the water-mediated interaction of the membranes.

Modeling 2D and 3D diffusion

Modeling obstructed diffusion is essential to the understanding of diffusion-mediated processes in the crowded cellular environment. Simple Monte Carlo techniques for modeling obstructed random walks are explained and related to Brownian dynamics and more complicated Monte Carlo methods. Random number generation is reviewed in the context of random walk simulations. Programming techniques and event-driven algorithms are discussed as ways to speed simulations.

High performance computing in biology: multimillion atom simulations of nanoscale systems

Computational methods have been used in biology for sequence analysis (bioinformatics), all-atom simulation (molecular dynamics and quantum calculations), and more recently for modeling biological networks (systems biology). Of these three techniques, all-atom simulation is currently the most computationally demanding, in terms of compute load, communication speed, and memory load. Breakthroughs in electrostatic force calculation and dynamic load balancing have enabled molecular dynamics simulations of large biomolecular complexes. Here, we report simulation results for the ribosome, using approximately 2.64 million atoms, the largest all-atom biomolecular simulation published to date. Several other nano-scale systems with different numbers of atoms were studied to measure the performance of the NAMD molecular dynamics simulation program on the Los Alamos National Laboratory Q Machine. We demonstrate that multimillion atom systems represent a 'sweet spot' for the NAMD code on large supercomputers. NAMD displays an unprecedented 85% parallel scaling efficiency for the ribosome system on 1024 CPUs. We also review recent targeted molecular dynamics simulations of the ribosome that prove useful for studying conformational changes of this large biomolecular complex in atomic detail.

The association of tetrameric acetylcholinesterase with ColQ tail: a block normal mode analysis

Acetylcholinesterase (AChE) rapidly hydrolyzes acetylcholine in the neuromuscular junctions and other cholinergic synapses to terminate the neuronal signal. In physiological conditions, AChE exists as tetramers associated with the proline-rich attachment domain (PRAD) of either collagen-like Q subunit (ColQ) or proline-rich membrane-anchoring protein. Crystallographic studies have revealed that different tetramer forms may be present, and it is not clear whether one or both are relevant under physiological conditions. Recently, the crystal structure of the tryptophan amphiphilic tetramerization (WAT) domain of AChE associated with PRAD ([WAT]₄PRAD), which mimics the interface between ColQ and AChE tetramer, became available. In this study we built a complete tetrameric mouse [AChE_T]₄–ColQ atomic structure model, based on the crystal structure of the [WAT]₄PRAD complex. The structure was optimized using energy minimization. Block normal mode analysis was done to investigate the low-frequency motions of the complex and to correlate the structure model with the two known crystal structures of AChE tetramer. Significant low-frequency motions among the catalytic domains of the four AChE subunits were observed, while the [WAT]₄PRAD part held the complex together. Normal mode involvement analysis revealed that the two lowest frequency modes were primarily involved in the conformational changes leading to the two crystal structures. The first 30 normal modes can account for more than 75% of the conformational changes in both cases. The evidence further supports the idea of a flexible tetramer model for AChE. This model can be used to study the implications of the association of AChE with ColQ.

Acetylcholinesterase (AChE) breaks down acetylcholine in the neuromuscular junction and other cholinergic synapses to terminate neuronal signals. AChE exists as tetramers anchored by structural subunits to the cell membranes in the brain or the basal lamina in the neuromuscular junction. Based on a crystal structure of the tetramerization domain of AChE with a proline-rich attachment domain of the anchoring proteins, a symmetric model of the complex of AChE tetramer with the anchoring protein tail was constructed. Block normal mode analysis revealed the presence of several low-frequency, low-barrier normal modes corresponding to inter-subunit motions. Previous crystal structures of AChE tetramer could be rationalized using these normal modes. These low-frequency modes are due to the presence of a flexible hinge in the structure of AChE. This study paints a picture of a flexible AChE tetramer with different conformational states interconverting easily under physiological conditions, which has important implications on the function of AChE. In particular, AChE is not trapped in the compact tetramer structure, for which access of substrate to two of the active sites is somewhat limited. Rather, the tetramer fluctuates to expose all four of its active sites to ensure rapid removal of acetylcholine.

Analysis of functional motions in Brownian molecular machines with an efficient block normal mode approach: myosin-II and Ca2+ -ATPase

The structural flexibilities of two molecular machines, myosin and Ca(2+)-ATPase, have been analyzed with normal mode analysis and discussed in the context of their energy conversion functions. The normal mode analysis with physical intermolecular interactions was made possible by an improved implementation of the block normal mode (BNM) approach. The BNM results clearly illustrated that the large-scale conformational transitions implicated in the functional cycles of the two motor systems can be largely captured with a small number of low-frequency normal modes. Therefore, the results support the idea that structural flexibility is an essential part of the construction principle of molecular motors through evolution. Such a feature is expected to be more prevalent in motor proteins than in simpler systems (e.g., signal transduction proteins) because in the former, large-scale conformational transitions often have to occur before the chemical events (e.g., ATP hydrolysis in myosin and ATP binding/phosphorylation in Ca(2+)-ATPase). This highlights the importance of Brownian motions associated with the protein domains that are involved in the functional transitions; in this sense, Brownian molecular machines is an appropriate description of molecular motors, although the normal mode results do not address the origin of the ratchet effect. The results also suggest that it might be more appropriate to describe functional transitions in some molecular motors as intrinsic elastic motions modulating local structural changes in the active site, which in turn gets stabilized by the subsequent chemical events, in contrast with the conventional idea of local changes somehow getting amplified into larger-scale motions. In the case of myosin, for example, we favor the idea that Brownian motions associated with the flexible converter propagates to the Switch I/II region, where the salt-bridge formation gets stabilized by ATP hydrolysis, in contrast with the textbook notion that ATP hydrolysis drives the converter motion. Another useful aspect of the BNM results is that selected low-frequency normal modes have been identified to form a set of collective coordinates that can be used to characterize the progress of a significant fraction of large-scale conformational transitions. Therefore, the present normal mode analysis has provided a stepping-stone toward more elaborate microscopic simulations for addressing critical issues in free energy conversions in molecular machines, such as the coupling and the causal relationship between collective motions and essential local changes at the catalytic active site where ATP hydrolysis occurs.

The electrostatic nature of C3d-complement receptor 2 association

The association of complement component C3d with B or T cell complement receptor 2 (CR2 or CD21) is a link between innate and adaptive immunity. It has been recognized in experimental studies that the C3d-CR2 association is pH- and ionic strength-dependent. This led us to perform electrostatic calculations to obtain a theoretical understanding of the mechanism of C3d-CR2 association. We used the crystallographic structures of human free C3d, free CR2 (short consensus repeat (SCR)1-2), and the C3d-CR2(SCR1-2) complex, and continuum solvent representation, to obtain a detailed atomic-level picture of the components of the two molecules that contribute to association. Based on the calculation of electrostatic potentials for the free and bound species and apparent pK(a) values for each ionizable residue, we show that C3d-CR2(SCR1-2) recognition is electrostatic in nature and involves not only the association interface, but also the whole molecules. Our results are in qualitative agreement with experimental data that measured the ionic strength and pH dependence of C3d-CR2 association. Also, our results for the native molecules and a number of theoretical mutants of C3d explain experimental mutagenesis studies of amino acid replacements away from the association interface that modulate binding of iC3b with full-length CR2. Finally, we discuss the packing of the two SCR domains. Overall, our data provide global and site-specific explanations of the physical causes that underlie the ionic strength dependence of C3d-CR2 association in a unified model that accounts for all experimental data, some of which were previously thought to be contradictory.

A computational model for the electrostatic sequestration of PI(4,5)P2 by membrane-adsorbed basic peptides

The multivalent acidic phospholipid phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) plays a key role in many biological processes. Recent studies show that unstructured clusters of basic residues from a number of peripheral proteins can laterally sequester PI(4,5)P2 in membranes. Specifically, experiments suggest that the basic effector domain of the myristoylated alanine-rich C kinase substrate (MARCKS), or a peptide corresponding to this domain, MARCKS(151-175), sequesters several PI(4,5)P2 and that this sequestration is due to nonspecific electrostatic interactions. Here, we use the finite difference Poisson-Boltzmann method to test this hypothesis by calculating the electrostatic free energy of lateral sequestration of PI(4,5)P2 by membrane-adsorbed basic peptides: Lys-7, Lys-13, and FA-MARCKS(151-175), a peptide based on MARCKS(151-175). In agreement with experiments, we find that the electrostatic free energy becomes more favorable when: 1), Lys-13 and FA-MARCKS(151-175) sequester several PI(4,5)P2; 2), the linear charge density of the basic peptide increases; 3), the mol percent monovalent acidic lipid in the membrane decreases; and 4), the ionic strength of the solution decreases. In addition, the electrostatic sequestration free energy is in excess of the entropic penalty associated with localizing PI(4,5)P2. Our calculations, thus, provide a structural and quantitative description of the observed interaction of PI(4,5)P2 with membrane-adsorbed basic sequences.

Conformational sampling for the impatient

Several new methods for sampling conformations of biomolecules have appeared recently. A brief review thereof is presented, with particular emphasis on applications that have been published, and suitability for different kinds of systems. Four methods (namely: RESPA, replica-exchange molecular dynamics, CONCOORD and Gaussian network method) are readily applicable for biomolecular systems.

The implementation of slab geometry for membrane-channel molecular dynamics simulations

Slab geometric boundary conditions are applied in the molecular dynamics simulation of a simple membrane-channel system. The results of the simulation were compared to those of an analogous system using normal three-dimensional periodic boundary conditions. Analysis of the dynamics and electrostatics of the system show that slab geometric periodicity eliminates the artificial bulk water orientational polarization that is present while using normal three-dimensional periodicity. Furthermore, even though the water occupancy and volume of our simple channel is the same when using either method, the electrostatic properties are considerably different when using slab geometry. In particular, the orientational polarization of water is seen to be different in the interior of the channel. This gives rise to a markedly different electric field within the channel. We discuss the implications of slab geometry for the future simulation of this type of system and for the study of channel transport properties.

In silico protein recombination: enhancing template and sequence alignment selection for comparative protein modelling

Comparative modelling of proteins is a predictive technique to build an atomic model for a given amino acid sequence, on the basis of the structures of other proteins (templates) that have been determined experimentally. Critical problems arise in this procedure: selecting the correct templates, aligning the query sequence with them and building the non-conserved surface loops. In this work, we apply a genetic algorithm, with crossover and mutation, as a new tool to overcome the first two. In silico protein recombination proves to be an effective way to exploit the variability of templates and sequence alignments to produce populations of optimized models by artificial selection. Despite some limitations, the procedure is shown to be robust to alignment errors, while simplifying the task of selecting templates, making it a good candidate for automatic building of reliable protein models.

Atomic-level observation of macromolecular crowding effects: escape of a protein from the GroEL cage

Experimental work has demonstrated that the efficient operation of the GroEL-GroES chaperonin machinery is sensitive to the presence of macromolecular crowding agents. Here, I describe atomically detailed computer simulations that provide a microscopic view of how crowding effects are exerted. Simulations were performed to compute the free energy required to extract the protein rhodanese from the central cavity of GroEL into solutions containing a range of crowder concentrations. The computed energetics allow the total yield of folded protein to be predicted; the calculated yields show a nonlinear dependence on the concentration of crowding agent identical to that observed experimentally. The close correspondence between simulation and experiment prompts the use of the former in a truly predictive setting: simulations are used to suggest that more effective crowding agents might be designed by exploiting an "agoraphobic effect."