Search

Experimental, theoretical, and computational studies of small proteins suggest that interresidue contacts not present in the folded structure play little or no role in the self-assembly mechanism. Non-native contacts can, however, influence folding kinetics by introducing additional local minima that slow diffusion over the global free-energy barrier between folded and unfolded states. Here, we combine single-molecule fluorescence with all-atom molecular dynamics simulations to discover the structural origin for the slow diffusion that markedly decreases the folding rate for a designed α-helical protein. Our experimental determination of transition path times and our analysis of the simulations point to non-native salt bridges between helices as the source, which provides a quantitative glimpse of how specific intramolecular interactions influence protein folding rates by altering dynamics and not activation free energies., (Copyright © 2015, American Association for the Advancement of Science.)

G protein-coupled receptors (GPCRs) relay diverse extracellular signals into cells by catalyzing nucleotide release from heterotrimeric G proteins, but the mechanism underlying this quintessential molecular signaling event has remained unclear. Here we use atomic-level simulations to elucidate the nucleotide-release mechanism. We find that the G protein α subunit Ras and helical domains-previously observed to separate widely upon receptor binding to expose the nucleotide-binding site-separate spontaneously and frequently even in the absence of a receptor. Domain separation is necessary but not sufficient for rapid nucleotide release. Rather, receptors catalyze nucleotide release by favoring an internal structural rearrangement of the Ras domain that weakens its nucleotide affinity. We use double electron-electron resonance spectroscopy and protein engineering to confirm predictions of our computationally determined mechanism., (Copyright © 2015, American Association for the Advancement of Science.)

The mechanism of ion channel voltage gating-how channels open and close in response to voltage changes-has been debated since Hodgkin and Huxley's seminal discovery that the crux of nerve conduction is ion flow across cellular membranes. Using all-atom molecular dynamics simulations, we show how a voltage-gated potassium channel (KV) switches between activated and deactivated states. On deactivation, pore hydrophobic collapse rapidly halts ion flow. Subsequent voltage-sensing domain (VSD) relaxation, including inward, 15-angstrom S4-helix motion, completes the transition. On activation, outward S4 motion tightens the VSD-pore linker, perturbing linker-S6-helix packing. Fluctuations allow water, then potassium ions, to reenter the pore; linker-S6 repacking stabilizes the open pore. We propose a mechanistic model for the sodium/potassium/calcium voltage-gated ion channel superfamily that reconciles apparently conflicting experimental data.

An outstanding challenge in the field of molecular biology has been to understand the process by which proteins fold into their characteristic three-dimensional structures. Here, we report the results of atomic-level molecular dynamics simulations, over periods ranging between 100 μs and 1 ms, that reveal a set of common principles underlying the folding of 12 structurally diverse proteins. In simulations conducted with a single physics-based energy function, the proteins, representing all three major structural classes, spontaneously and repeatedly fold to their experimentally determined native structures. Early in the folding process, the protein backbone adopts a nativelike topology while certain secondary structure elements and a small number of nonlocal contacts form. In most cases, folding follows a single dominant route in which elements of the native structure appear in an order highly correlated with their propensity to form in the unfolded state.

Molecular dynamics (MD) simulations are widely used to study protein motions at an atomic level of detail, but they have been limited to time scales shorter than those of many biologically critical conformational changes. We examined two fundamental processes in protein dynamics--protein folding and conformational change within the folded state--by means of extremely long all-atom MD simulations conducted on a special-purpose machine. Equilibrium simulations of a WW protein domain captured multiple folding and unfolding events that consistently follow a well-defined folding pathway; separate simulations of the protein's constituent substructures shed light on possible determinants of this pathway. A 1-millisecond simulation of the folded protein BPTI reveals a small number of structurally distinct conformational states whose reversible interconversion is slower than local relaxations within those states by a factor of more than 1000.

Na+/H+ antiporters are central to cellular salt and pH homeostasis. The structure of Escherichia coli NhaA was recently determined, but its mechanisms of transport and pH regulation remain elusive. We performed molecular dynamics simulations of NhaA that, with existing experimental data, enabled us to propose an atomically detailed model of antiporter function. Three conserved aspartates are key to our proposed mechanism: Asp164 (D164) is the Na+-binding site, D163 controls the alternating accessibility of this binding site to the cytoplasm or periplasm, and D133 is crucial for pH regulation. Consistent with experimental stoichiometry, two protons are required to transport a single Na+ ion: D163 protonates to reveal the Na+-binding site to the periplasm, and subsequent protonation of D164 releases Na+. Additional mutagenesis experiments further validated the model.

The urokinase plasminogen activator binds to its cellular receptor with high affinity and initiates signaling cascades that are implicated in pathological processes including tumor growth, metastasis, and inflammation. We report the crystal structure at 1.9 angstroms of the urokinase receptor complexed with the urokinase amino-terminal fragment and an antibody against the receptor. The three domains of urokinase receptor form a concave shape with a central cone-shaped cavity where the urokinase fragment inserts. The structure provides insight into the flexibility of the urokinase receptor that enables its interaction with a wide variety of ligands and a basis for the design of urokinase-urokinase receptor antagonists.