Force-producing ADP state of myosin bound to actin.

Molecular motors produce force when they interact with their cellular tracks. For myosin motors, the primary force-generating state has MgADP tightly bound, whereas myosin is strongly bound to actin. We have generated an 8-Å cryoEM reconstruction of this state for myosin V and used molecular dynamics flexed fitting for model building.

Molecular dissection of human Argonaute proteins by DNA shuffling.

A paramount task in RNA interference research is to decipher the complex biology of cellular effectors, exemplified in humans by four pleiotropic Argonaute proteins (Ago1-Ago4). Here, we exploited DNA family shuffling, a molecular evolution technology, to generate chimeric Ago protein libraries for dissection of intricate phenotypes independently of prior structural knowledge. Through shuffling of human Ago2 and Ago3, we discovered two N-terminal motifs that govern RNA cleavage in concert with the PIWI domain.

Negative protein-protein interaction datasets derived from large-scale two-hybrid experiments.

Negative protein-protein interaction datasets are needed for training and evaluation of interaction prediction methods, as well as validation of high-throughput interaction discovery experiments. In large-scale two-hybrid assays, the direct interaction of a large number of protein pairs is systematically probed. We present a simple method to harness two-hybrid data to obtain negative protein-protein interaction datasets, which we validated using other available experimental data.

Cryo-electron microscopy modeling by the molecular dynamics flexible fitting method.

The increasing power and popularity of cryo-electron microscopy (cryo-EM) in structural biology brought about the development of so-called hybrid methods, which permit the interpretation of cryo-EM density maps beyond their nominal resolution in terms of atomic models. The Cryo-EM Modeling Challenge 2010 is the first community effort to bring together developers of hybrid methods as well as cryo-EM experimentalists. Participating in the challenge, the molecular dynamics flexible fitting (MDFF) method was applied to a number of cryo-EM density maps.

PepSite: prediction of peptide-binding sites from protein surfaces.

Complex biological functions emerge through intricate protein-protein interaction networks. An important class of protein-protein interaction corresponds to peptide-mediated interactions, in which a short peptide stretch from one partner interacts with a large protein surface from the other partner. Protein-peptide interactions are typically of low affinity and involved in regulatory mechanisms, dynamically reshaping protein interaction networks.

Stereochemical errors and their implications for molecular dynamics simulations.

Background: Biological molecules are often asymmetric with respect to stereochemistry, and correct stereochemistry is essential to their function. Molecular dynamics simulations of biomolecules have increasingly become an integral part of biophysical research. However, stereochemical errors in biomolecular structures can have a dramatic impact on the results of simulations.

Structural insights into cognate versus near-cognate discrimination during decoding.

The structural basis of the tRNA selection process is investigated by cryo-electron microscopy of ribosomes programmed with UGA codons and incubated with ternary complex (TC) containing the near-cognate Trp-tRNA(Trp) in the presence of kirromycin. Going through more than 350 000 images and employing image classification procedures, we find ∼8% in which the TC is bound to the ribosome. The reconstructed 3D map provides a means to characterize the arrangement of the near-cognate aa-tRNA with respect to elongation factor Tu (EF-Tu) and the ribosome, as well as the domain movements of the ribosome.

Molecular dynamics of EF-G during translocation.

Elongation factor G (EF-G) plays a crucial role in two stages of mRNA-(tRNA)(2) translocation. First, EF-G•GTP enters the pre-translocational ribosome in its intersubunit-rotated state, with tRNAs in their hybrid (P/E and A/P) positions. Second, a conformational change in EF-G's Domain IV induced by GTP hydrolysis disengages the mRNA-anticodon stem-loops of the tRNAs from the decoding center to advance relative to the small subunit when the ribosome undergoes a backward inter-subunit rotation.

Applications of the molecular dynamics flexible fitting method.

In recent years, cryo-electron microscopy (cryo-EM) has established itself as a key method in structural biology, permitting the structural characterization of large biomolecular complexes in various functional states. The data obtained through single-particle cryo-EM has recently seen a leap in resolution thanks to landmark advances in experimental and computational techniques, resulting in sub-nanometer resolution structures being obtained routinely. The remaining gap between these data and revealing the mechanisms of molecular function can be closed through hybrid modeling tools that incorporate known atomic structures into the cryo-EM data.

The role of L1 stalk-tRNA interaction in the ribosome elongation cycle.

The ribosomal L1 stalk is a mobile structure implicated in directing tRNA movement during translocation through the ribosome. This article investigates three aspects of L1 stalk-tRNA interaction. First, by combining data from cryo electron microscopy, X-ray crystallography, and molecular dynamics simulations through the molecular dynamics flexible fitting method, we obtained atomic models of different tRNAs occupying the hybrid P/E state interacting with the L1 stalk.

Recognition of the regulatory nascent chain TnaC by the ribosome.

Regulatory nascent chains interact with the ribosomal exit tunnel and modulate their own translation. To characterize nascent chain recognition by the ribosome at the atomic level, extensive molecular dynamics simulations of TnaC, the leader peptide of the tryptophanase operon, inside the exit tunnel were performed for an aggregate time of 2.1 mus.

Structural model and excitonic properties of the dimeric RC-LH1-PufX complex from Rhodobacter sphaeroides.

The light-harvesting apparatus of the purple bacterial photosynthetic unit consists of a pool of peripheral light-harvesting complexes that transfer excitation energy to a reaction center (RC) via the surrounding pigment-protein complex LH1. Recent electron microscopy and atomic force microscopy studies have revealed that RC-LH1 units of Rhodobacter sphaeroides form membrane-bending dimeric complexes together with the polypeptide PufX. We present a structural model for these RC-LH1-PufX dimeric complexes constructed using the molecular dynamics flexible fitting method based on an EM density map.

Structural insight into nascent polypeptide chain-mediated translational stalling.

Expression of the Escherichia coli tryptophanase operon depends on ribosome stalling during translation of the upstream TnaC leader peptide, a process for which interactions between the TnaC nascent chain and the ribosomal exit tunnel are critical. We determined a 5.8 angstrom-resolution cryo-electron microscopy and single-particle reconstruction of a ribosome stalled during translation of the tnaC leader gene.

Regulation of the protein-conducting channel by a bound ribosome.

During protein synthesis, it is often necessary for the ribosome to form a complex with a membrane-bound channel, the SecY/Sec61 complex, in order to translocate nascent proteins across a cellular membrane. Structural data on the ribosome-channel complex are currently limited to low-resolution cryo-electron microscopy maps, including one showing a bacterial ribosome bound to a monomeric SecY complex. Using that map along with available atomic-level models of the ribosome and SecY, we have determined, through molecular dynamics flexible fitting (MDFF), an atomic-resolution model of the ribosome-channel complex.

Protein-induced membrane curvature investigated through molecular dynamics flexible fitting.

In the photosynthetic purple bacterium Rhodobacter (Rba.) sphaeroides, light is absorbed by membrane-bound light-harvesting (LH) proteins LH1 and LH2. LH1 directly surrounds the reaction center (RC) and, together with PufX, forms a dimeric (RC-LH1-PufX)2 protein complex.

Molecular dynamics flexible fitting: a practical guide to combine cryo-electron microscopy and X-ray crystallography.

Hybrid computational methods for combining structural data from different sources and resolutions are becoming an essential part of structural biology, especially as the field moves toward the study of large macromolecular assemblies. We have developed the molecular dynamics flexible fitting (MDFF) method for combining high-resolution atomic structures with cryo-electron microscopy (cryo-EM) maps, that results in atomic models representing the conformational state captured by cryo-EM. The method has been applied successfully to the ribosome, a ribonucleoprotein complex responsible for protein synthesis.

Ribosome-induced changes in elongation factor Tu conformation control GTP hydrolysis.

In translation, elongation factor Tu (EF-Tu) molecules deliver aminoacyl-tRNAs to the mRNA-programmed ribosome. The GTPase activity of EF-Tu is triggered by ribosome-induced conformational changes of the factor that play a pivotal role in the selection of the cognate aminoacyl-tRNAs. We present a 6.

Flexible fitting of atomic structures into electron microscopy maps using molecular dynamics.

A novel method to flexibly fit atomic structures into electron microscopy (EM) maps using molecular dynamics simulations is presented. The simulations incorporate the EM data as an external potential added to the molecular dynamics force field, allowing all internal features present in the EM map to be used in the fitting process, while the model remains fully flexible and stereochemically correct. The molecular dynamics flexible fitting (MDFF) method is validated for available crystal structures of protein and RNA in different conformations; measures to assess and monitor the fitting process are introduced.

Dynamics of Recognition between tRNA and elongation factor Tu.

Elongation factor Tu (EF-Tu) binds to all standard aminoacyl transfer RNAs (aa-tRNAs) and transports them to the ribosome while protecting the ester linkage between the tRNA and its cognate amino acid. We use molecular dynamics simulations to investigate the dynamics of the EF-Tu.guanosine 5'-triphosphate.

Accelerating molecular modeling applications with graphics processors.

Molecular mechanics simulations offer a computational approach to study the behavior of biomolecules at atomic detail, but such simulations are limited in size and timescale by the available computing resources. State-of-the-art graphics processing units (GPUs) can perform over 500 billion arithmetic operations per second, a tremendous computational resource that can now be utilized for general purpose computing as a result of recent advances in GPU hardware and software architecture. In this article, an overview of recent advances in programmable GPUs is presented, with an emphasis on their application to molecular mechanics simulations and the programming techniques required to obtain optimal performance in these cases.

ATP-sensitive K+ channels in renal mitochondria.

Isolated kidney mitochondria swell when incubated in hyposmotic solutions containing K+ salts in a manner inhibited by ATP, ADP, 5-hydroxydecanoate, and glibenclamide and stimulated by GTP and diazoxide. These results suggest the existence of ATP-sensitive K+ channels in these mitochondria, similar to those previously described in heart, liver, and brain. Renal mitochondrial ATP-sensitive K+ uptake rates are approximately 140 nmol.