Dynamic Behavior of Trigger Factor on the Ribosome.

Trigger factor (TF) is the only ribosome-associated chaperone in bacteria. It interacts with hydrophobic segments in nascent chain (NCs) as they emerge from the ribosome. TF binds via its N-terminal ribosome-binding domain (RBD) mainly to ribosomal protein uL23 at the tunnel exit on the large ribosomal subunit. Whereas earlier structural data suggested that TF binds as a rigid molecule to the ribosome, recent comparisons of structural data on substrate-bound, ribosome-bound, and TF in solution from different species suggest that this chaperone is a rather flexible molecule. Here, we present two cryo-electron microscopy structures of TF bound to ribosomes translating an mRNA coding for a known TF substrate from Escherichia coli of a different length. The structures reveal distinct degrees of flexibility for the different TF domains, a conformational rearrangement of the RBD upon ribosome binding, and an increase in rigidity within TF when the NC is extended. Molecular dynamics simulations agree with these data and offer a molecular basis for these observations.

How Quantum Coherence Assists Photosynthetic Light Harvesting.

This perspective examines how hundreds of pigment molecules in purple bacteria cooperate through quantum coherence to achieve remarkable light harvesting efficiency. Quantum coherent sharing of excitation, which modifies excited state energy levels and combines transition dipole moments, enables rapid transfer of excitation over large distances. Purple bacteria exploit the resulting excitation transfer to engage many antenna proteins in light harvesting, thereby increasing the rate of photon absorption and energy conversion. We highlight here how quantum coherence comes about and plays a key role in the photosynthetic apparatus of purple bacteria.

Squeezing protein shells: how continuum elastic models, molecular dynamics simulations, and experiments coalesce at the nanoscale.

The current rapid growth in the use of nanosized particles is fueled in part by our increased understanding of their physical properties and ability to manipulate them, which is essential for achieving optimal functionality. Here we report detailed quantitative measurements of the mechanical response of nanosized protein shells (viral capsids) to large-scale physical deformations and compare them with theoretical descriptions from continuum elastic modeling and molecular dynamics (MD). Specifically, we used nanoindentation by atomic force microscopy to investigate the complex elastic behavior of Hepatitis B virus capsids. These capsids are hollow, approximately 30 nm in diameter, and conform to icosahedral (5-3-2) symmetry. First we show that their indentation behavior, which is symmetry-axis-dependent, cannot be reproduced by a simple model based on Föppl-von Kármán thin-shell elasticity with the fivefold vertices acting as prestressed disclinations. However, we can properly describe the measured nonlinear elastic and orientation-dependent force response with a three-dimensional, topographically detailed, finite-element model. Next, we show that coarse-grained MD simulations also yield good agreement with our nanoindentation measurements, even without any fitting of force-field parameters in the MD model. This study demonstrates that the material properties of viral nanoparticles can be correctly described by both modeling approaches. At the same time, we show that even for large deformations, it suffices to approximate the mechanical behavior of nanosized viral shells with a continuum approach, and ignore specific molecular interactions. This experimental validation of continuum elastic theory provides an example of a situation in which rules of macroscopic physics can apply to nanoscale molecular assemblies.

Nanoelectromechanics of methylated DNA in a synthetic nanopore.

Methylation of cytosine is a covalent modification of DNA that can be used to silence genes, orchestrating a myriad of biological processes including cancer. We have discovered that a synthetic nanopore in a membrane comparable in thickness to a protein binding site can be used to detect methylation. We observe a voltage threshold for permeation of methylated DNA through a <2 nm diameter pore, which we attribute to the stretching transition; this can differ by >1 V/20 nm depending on the methylation level, but not the DNA sequence.

Possible pathway for ubiquinone shuttling in Rhodospirillum rubrum revealed by molecular dynamics simulation.

In the last decade, the structures of many components of the photosynthetic apparatus of purple bacteria, as well as the mutual organization of these components within the purple membrane, were resolved. One key question that emerged concerned the assembly of the core complex consisting of the reaction center (RC) and the light-harvesting 1 (LH1) complex. In some species, like Rhodobacter sphaeroides, the ring-shaped LH1 complex was found to be open, whereas other species, like Rhodospirillum rubrum, have a closed ring surrounding the reaction center. This poses the question of how the ubiquinone molecule that transports electrons and protons from the RC to the cytochrome bc(1) complex overcomes the apparent barrier of the LH1 ring. In this study, we investigated how, in the case of a closed LH1 ring, the ubiquinone molecule diffuses through the LH1 ring. For this purpose, the LH1 structure of R. rubrum was modeled and the potential of mean force along the diffusion pathway through the LH1 was determined by steered molecular-dynamics simulations. The potential was reconstructed using the fluctuation theorem in combination with the stiff spring approximation. An upper limit for the mean first-passage time for diffusion of ubiquinone through the LH1 ring, based on a worst-case scenario potential, was calculated as approximately 8 x 10(-3) s, which is still in agreement with known turnover rates of RC and RC-LH1 complexes in the range of approximately 1000 Hz.

The electromechanics of DNA in a synthetic nanopore.

We have explored the electromechanical properties of DNA on a nanometer-length scale using an electric field to force single molecules through synthetic nanopores in ultrathin silicon nitride membranes. At low electric fields, E < 200 mV/10 nm, we observed that single-stranded DNA can permeate pores with a diameter >/=1.0 nm, whereas double-stranded DNA only permeates pores with a diameter >/=3 nm. For pores <3.0 nm diameter, we find a threshold for permeation of double-stranded DNA that depends on the electric field and pH. For a 2 nm diameter pore, the electric field threshold is approximately 3.1 V/10 nm at pH = 8.5; the threshold decreases as pH becomes more acidic or the diameter increases. Molecular dynamics indicates that the field threshold originates from a stretching transition in DNA that occurs under the force gradient in a nanopore. Lowering pH destabilizes the double helix, facilitating DNA translocation at lower fields.

Approaches to developing biological H(2)-photoproducing organisms and processes.

The development of efficient biological systems for the direct photoproduction of H(2) gas from water faces several challenges, the more serious of which is the sensitivity of the H(2)-evolving enzymes (hydrogenases) to O(2), an obligatory by-product of photosynthesis. This high sensitivity is common to both FeFe and NiFe hydrogenases, and is caused by O(2) binding to their respective metallocatalytic sites. This overview describes approaches to (i) molecular engineering of algal FeFe-hydrogenase to prevent O(2) access to its catalytic site; (ii) transform a cyanobacterium with an O(2)-tolerant bacterial NiFe hydrogenase or (c) partially inactivate algal O(2)-evolution activity to create physiologically anaerobiosis and induce hydrogenase expression.

Molecular dynamics and experimental investigation of H(2) and O(2) diffusion in [Fe]-hydrogenase.

The [Fe]-hydrogenase enzymes are highly efficient H(2) catalysts found in ecologically and phylogenetically diverse microorganisms, including the photosynthetic green alga, Chlamydomonas reinhardtii. Although these enzymes can occur in several forms, H(2) catalysis takes place at a unique [FeS] prosthetic group or H-cluster, located at the active site. Significant to the function of hydrogenases is how the surrounding protein structure facilitates substrate-product transfer, and protects the active site H-cluster from inactivation. To elucidate the role of protein structure in O(2) inactivation of [Fe]-hydrogenases, experimental and theoretical investigations have been performed. Molecular dynamics was used to comparatively investigate O(2) and H(2) diffusion in CpI ([Fe]-hydrogenase I from Clostridium pasteurianum). Our preliminary results suggest that H(2) diffuses more easily and freely than O(2), which is restricted to a small number of allowed pathways to and from the active site. These O(2) pathways are located in the conserved active site domain, shown experimentally to have an essential role in active site protection.

Beyond the Gene Chip.

We describe a prospective strategy for reading the encyclopedic information encoded in the genome: using a nanopore in a membrane formed from an MOS-capacitor to sense the charge in DNA. In principle, as DNA permeates the capacitor-membrane through the pore, the electrostatic charge distribution characteristic of the molecule should polarize the capacitor and induce a voltage on the electrodes that can be measured. Silicon nanofabrication and molecular dynamic simulations with atomic detail are technological linchpins in the development of this detector. The sub-nanometer precision available through silicon nanotechnology facilitates the fabrication of the detector, and molecular dynamics provides us with a means to design it and analyze the experimental outcomes.

Modelling DNA loops using continuum and statistical mechanics.

The classical Kirchhoff elastic-rod model applied to DNA is extended to account for sequence-dependent intrinsic twist and curvature, anisotropic bending rigidity, electrostatic force interactions, and overdamped Brownian motion in a solvent. The zero-temperature equilibrium rod model is then applied to study the structural basis of the function of the lac repressor protein in the lac operon of Escherichia coli. The structure of a DNA loop induced by the clamping of two distant DNA operator sites by lac repressor is investigated and the optimal geometries for the loop of length 76 bp are predicted. Further, the mimicked binding of catabolite gene activator protein (CAP) inside the loop provides solutions that might explain the experimentally observed synergy in DNA binding between the two proteins. Finally, a combined Monte Carlo and Brownian dynamics solver for a worm-like chain model is described and a preliminary analysis of DNA loop-formation kinetics is presented.

Multiple probes reveal a native-like intermediate during low-temperature refolding of ubiquitin.

We investigate the refolding of ubiquitin Phe45Trp/Ile61Ala (Ub(*)I61A) in a low-temperature, high-viscosity buffer, where folding is slowed so that apparent two-state and three-state mechanisms are readily distinguishable. Ub(*)I61A forms a compact ensemble rapidly (as judged from stopped-flow, small-angle X-ray scattering) with a secondary structure signature similar to that of the native state (as judged from stopped-flow circular dichroism from 215 nm to 250 nm), but the fluorescence signature still resembles the guanidinium-denatured state. The compact ensemble forms over a range of solvent and temperature conditions. The native fluorescence signature, which requires the tryptophan residue to be packed tightly, is acquired at least 500 times more slowly. Molecular dynamics simulations at 495 K show no contraction of the backbone in ethylene glycol buffer compared to pure aqueous buffer, and no significant effect on the local backbone structure of the unfolded protein. Only at higher simulation temperature does a backbone contraction appear. Thus, it appears unlikely that the aqueous ethylene glycol buffer fundamentally changes the folding mechanism of ubiquitin. We suggest that ubiquitin forms a compact ensemble with native-like secondary structure, but without tight packing, long before the native state.

Forced dissociation of the strand dimer interface between C-cadherin ectodomains.

The force-induced dissociation of the strand dimer interface in C-cadherin has been studied using steered molecular dynamics simulations. The dissociation occurred, without domain unraveling, after the extraction of the conserved trypthophans (Trp2) from their respective hydrophobic pockets. The simulations revealed two stable positions for the Trp2 side chain inside the pocket. The most internal stable position involved a hydrogen bond between the ring Nepsilon of Trp2 and the backbone carbonyl of Glu90. In the second stable position, the aromatic ring is located at the pocket entrance. After extracting the two tryptophans from their pockets, the complex exists in an intermediate bound state that involves a close packing of the tryptophans with residues Asp1 and Asp27 from both domains. Dissociation occurred after this residue association was broken. Simulations carried out with a complex formed between W2A mutants showed that the mutant complex dissociates more easily than the wild type complex does. These results correlate closely with the role of the conserved tryptophans suggested previously by site directed mutagenesis.

Forced detachment of the CD2-CD58 complex.

The force-induced detachment of the adhesion protein complex CD2-CD58 was studied by steered molecular dynamics simulations. The forced detachment of CD2 and CD58 shows that the system can respond to an external force by two mechanisms, which depend on the loading rate. At the rapid loading rates of 70 and 35 pN/ps (pulling speeds of 1 and 0.5 A/ps) the two proteins unfold before they separate, whereas at slower loading rates of 7 and 3.5 pN/ps (pulling speeds of 0.1 and 0.05 A/ps), the proteins separate before the domains can unfold. When subjected to a constant force of 400 pN, the two proteins separated without significant structural distortion. These findings suggest that protein unfolding is not coupled to the adhesive function of CD2 and CD58. The simulations further confirm that salt bridges primarily determine the tensile strength of the protein-to-protein bond, and that the order of salt bridge rupture depends mainly on the position of the bond, relative to the line of action of the applied force. Salt bridges close to this line break first. The importance of each of the salt bridges for adhesion, determined from the simulations, correlates closely with their role in cell-to-cell adhesion and equilibrium binding determined by site-directed mutagenesis experiments.

The role of intersection topography in bond selectivity of cis-trans photoisomerization.

Ab initio methods are used to characterize the ground and first excited state of the chromophore in the rhodopsin family of proteins: retinal protonated Schiff base. Retinal protonated Schiff base has five double bonds capable of undergoing isomerization. Upon absorption of light, the chromophore isomerizes and the character of the photoproducts (e.g., 13-cis and 11-cis) depends on the environment, protein vs. solution. Our ab initio calculations show that, in the absence of any specific interactions with the environment (e.g., discrete ordered charges in a protein), energetic considerations cannot explain the observed bond selectivity. We instead attribute the origin of bond selectivity to the shape (topography) of the potential energy surfaces in the vicinity of points of true degeneracy (conical intersections) between the ground and first excited electronic states. This provides a molecular example where a competition between two distinct but nearly isoenergetic photochemical reaction pathways is resolved by a topographical difference between two conical intersections.

The mechanism of glycerol conduction in aquaglyceroporins.

BACKGROUND: The E. coli glycerol facilitator, GlpF, selectively conducts glycerol and water, excluding ions and charged solutes. The detailed mechanism of the glycerol conduction and its relationship to the characteristic secondary structure of aquaporins and to the NPA motifs in the center of the channel are unknown. RESULTS: Molecular dynamics simulations of GlpF reveal spontaneous glycerol and water conduction driven, on a nanosecond timescale, by thermal fluctuations. The bidirectional conduction, guided and facilitated by the secondary structure, is characterized by breakage and formation of hydrogen bonds for which water and glycerol compete. The conduction involves only very minor changes in the protein structure, and cooperativity between the GlpF monomers is not evident. The two conserved NPA motifs are strictly linked together by several stable hydrogen bonds and their asparagine side chains form hydrogen bonds with the substrates passing the channel in single file. CONCLUSIONS: A complete conduction of glycerol through the GlpF was deduced from molecular dynamics simulations, and key residues facilitating the conduction were identified. The nonhelical parts of the two half-membrane-spanning segments expose carbonyl groups towards the channel interior, establishing a curve-linear pathway. The conformational stability of the NPA motifs is important in the conduction and critical for selectivity. Water and glycerol compete in a random manner for hydrogen bonding sites in the protein, and their translocations in single file are correlated. The suggested conduction mechanism should apply to the whole family.

Molecular dynamics study of aquaporin-1 water channel in a lipid bilayer.

The aquaporin-1 water channel was modeled in a palmitoyl-oleoyl-phosphatidyl-choline lipid bilayer, by means of molecular dynamics simulations. Interaction of the protein with the membrane and inter-monomer interactions were analyzed. Structural features of the channel important for its biological function, including the Asn-Pro-Ala (NPA) motifs, and the diffusion of water molecules into the channels, were investigated. Simulations revealed the formation of single file water inside the channels for certain relative positions of the NPA motifs.

Simulated refolding of stretched titin immunoglobulin domains.

Steered molecular dynamics (SMD) is used to investigate forced unfolding and spontaneous refolding of immunoglobulin I27, a domain of the muscle protein titin. Previous SMD simulations revealed the events leading to stretch-induced unfolding of I27, the rupture of hydrogen bonds bridging beta-strands A and B, and those bridging beta-strands A' and G, the latter rupture occurring at an extension of approximately 15 A and preceding the complete unfolding. Simulations are now used to study the refolding of partially unfolded I27 domains. The results reveal that stretched domains with ruptured interstrand hydrogen bonds shrink along the extension direction. Two types of refolding patterns are recognized: for separated beta-strands A' and G, in most simulations five of the six hydrogen bonds between A' and G stably reformed in 2 ns, whereas for separated beta-strands A and B hydrogen bonds seldom reformed in eight 2-ns simulations. The mechanical stability of the partially refolded intermediates has been tested by re-stretching.

Phylogenetic analysis of metabolic pathways.

The information provided by completely sequenced genomes can yield insights into the multi-level organization of organisms and their evolution. At the lowest level of molecular organization individual enzymes are formed, often through assembly of multiple polypeptides. At a higher level, sets of enzymes group into metabolic networks. Much has been learned about the relationship of species from phylogenetic trees comparing individual enzymes. In this article we extend conventional phylogenetic analysis of individual enzymes in different organisms to the organisms' metabolic networks. For this purpose we suggest a method that combines sequence information with information about the underlying reaction networks. A distance between pathways is defined as incorporating distances between substrates and distances between corresponding enzymes. The new analysis is applied to electron-transfer and amino acid biosynthesis networks yielding a more comprehensive understanding of similarities and differences between organisms.

Steered molecular dynamics investigations of protein function.

Molecular recognition and mechanical properties of proteins govern molecular processes in the cell that can cause disease and can be targeted for drug design. Single molecule measurement techniques have greatly advanced knowledge but cannot resolve enough detail to be interpreted in terms of protein structure. We seek to complement the observations through so-called Steered Molecular Dynamics (SMD) simulations that link directly to experiments and provide atomic-level descriptions of the underlying events. Such a research program has been initiated in our group and has involved, for example, studies of elastic properties of immunoglobulin and fibronectin domains as well as the binding of biotin and avidin. In this article we explain the SMD method and suggest how it can be applied to the function of three systems that are the focus of modern molecular biology research: force transduction by the muscle protein titin and extracellular matrix protein fibronectin, recognition of antibody-antigene pairs, and ion selective conductivity of the K+ channel.

Comparison of the early stages of forced unfolding for fibronectin type III modules.

The structural changes accompanying stretch-induced early unfolding events were investigated for the four type III fibronectin (FN-III) modules, FN-III(7), FN-III(8), FN-III(9), and FN-III(10) by using steered molecular dynamics. Simulations revealed that two main energy barriers, I and II, have to be overcome to initiate unraveling of FN-III's tertiary structure. In crossing the first barrier, the two opposing beta-sheets of FN-III are rotated against each other such that the beta-strands of both beta-sheets align parallel to the force vector (aligned state). All further events in the unfolding pathway proceed from this intermediate state. A second energy barrier has to be overcome to break the first major cluster of hydrogen bonds between adjacent beta-strands. Simulations revealed that the height of barrier I varied significantly among the four modules studied, being largest for FN-III(7) and lowest for FN-III(10), whereas the height of barrier II showed little variation. Key residues affecting the mechanical stability of FN-III modules were identified. These results suggest that FN-III modules can be prestretched into an intermediate state with only minor changes to their tertiary structures. FN-III(10), for example, extends 12 A from the native "twisted" to the intermediate aligned state, and an additional 10 A from the aligned state to further unfolding where the first beta-strand is peeled away. The implications of the existence of intermediate states regarding the elasticity of fibrillar fibers and the stretch-induced exposure of cryptic sites are discussed.