The impact of molecular dynamics (MD) simulations in molecular biology and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other ...biomolecules in full atomic detail and at very fine temporal resolution. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of experimental structural data, have increased the appeal of biomolecular simulation to experimentalists—a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomolecules, in uncovering the structural basis for disease, and in the design and optimization of small molecules, peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further experimental work.
Hollingsworth and Dror review modern molecular dynamics (MD) simulations, with an emphasis on how such simulations complement wet-lab experiments. MD simulations capture biomolecular motion in atomic detail and have come into widespread use because of recent technological and scientific advances.
Molecular dynamics simulations provide a vehicle for capturing the structures, motions, and interactions of biological macromolecules in full atomic detail. The accuracy of such simulations, however, ...is critically dependent on the force field--the mathematical model used to approximate the atomic-level forces acting on the simulated molecular system. Here we present a systematic and extensive evaluation of eight different protein force fields based on comparisons of experimental data with molecular dynamics simulations that reach a previously inaccessible timescale. First, through extensive comparisons with experimental NMR data, we examined the force fields' abilities to describe the structure and fluctuations of folded proteins. Second, we quantified potential biases towards different secondary structure types by comparing experimental and simulation data for small peptides that preferentially populate either helical or sheet-like structures. Third, we tested the force fields' abilities to fold two small proteins--one α-helical, the other with β-sheet structure. The results suggest that force fields have improved over time, and that the most recent versions, while not perfect, provide an accurate description of many structural and dynamical properties of proteins.
The mutation and overexpression of the epidermal growth factor receptor (EGFR) are associated with the development of a variety of cancers, making this prototypical dimerization-activated receptor ...tyrosine kinase a prominent target of cancer drugs. Using long-timescale molecular dynamics simulations, we find that the N lobe dimerization interface of the wild-type EGFR kinase domain is intrinsically disordered and that it becomes ordered only upon dimerization. Our simulations suggest, moreover, that some cancer-linked mutations distal to the dimerization interface, particularly the widespread L834R mutation (also referred to as L858R), facilitate EGFR dimerization by suppressing this local disorder. Corroborating these findings, our biophysical experiments and kinase enzymatic assays indicate that the L834R mutation causes abnormally high activity primarily by promoting EGFR dimerization rather than by allowing activation without dimerization. We also find that phosphorylation of EGFR kinase domain at Tyr845 may suppress the intrinsic disorder, suggesting a molecular mechanism for autonomous EGFR signaling.
Display omitted
► EGFR kinase is intrinsically disordered in the dimerization interface ► Cancer mutations stabilize the interface, promoting EGFR dimerization ► High activity of EGFR cancer mutants is due primarily to enhanced dimerization ► Phosphorylation similarly promotes dimerization and allows autonomous signaling
A prevalent EGFR mutation linked to lung cancer was previously thought to promote activation of the monomeric form of the receptor. A combination of biochemical and long-timescale molecular dynamics experiments now provides evidence that this mutation promotes ligand-independent dimerization by stabilizing receptors' dimerization interface.
GPCR Dynamics: Structures in Motion Latorraca, Naomi R; Venkatakrishnan, A J; Dror, Ron O
Chemical reviews,
2017-Jan-11, Letnik:
117, Številka:
1
Journal Article
Recenzirano
Odprti dostop
The function of G protein-coupled receptors (GPCRs)-which represent the largest class of both human membrane proteins and drug targets-depends critically on their ability to change shape, ...transitioning among distinct conformations. Determining the structural dynamics of GPCRs is thus essential both for understanding the physiology of these receptors and for the rational design of GPCR-targeted drugs. Here we review what is currently known about the flexibility and dynamics of GPCRs, as determined through crystallography, spectroscopy, and computer simulations. We first provide an overview of the types of motion exhibited by a GPCR and then discuss GPCR dynamics in the context of ligand binding, activation, allosteric modulation, and biased signaling. Finally, we discuss the implications of GPCR conformational plasticity for drug design.
How Fast-Folding Proteins Fold Lindorff-Larsen, Kresten; Piana, Stefano; Dror, Ron O. ...
Science (American Association for the Advancement of Science),
10/2011, Letnik:
334, Številka:
6055
Journal Article
Recenzirano
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 n ative like 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.
Geometric deep learning of RNA structure Townshend, Raphael J L; Eismann, Stephan; Watkins, Andrew M ...
Science,
08/2021, Letnik:
373, Številka:
6558
Journal Article
Recenzirano
Odprti dostop
RNA molecules adopt three-dimensional structures that are critical to their function and of interest in drug discovery. Few RNA structures are known, however, and predicting them computationally has ...proven challenging. We introduce a machine learning approach that enables identification of accurate structural models without assumptions about their defining characteristics, despite being trained with only 18 known RNA structures. The resulting scoring function, the Atomic Rotationally Equivariant Scorer (ARES), substantially outperforms previous methods and consistently produces the best results in community-wide blind RNA structure prediction challenges. By learning effectively even from a small amount of data, our approach overcomes a major limitation of standard deep neural networks. Because it uses only atomic coordinates as inputs and incorporates no RNA-specific information, this approach is applicable to diverse problems in structural biology, chemistry, materials science, and beyond.
Mechanism of Voltage Gating in Potassium Channels Jensen, Morten Ø.; Jogini, Vishwanath; Borhani, David W. ...
Science (American Association for the Advancement of Science),
04/2012, Letnik:
336, Številka:
6078
Journal Article
Recenzirano
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/caldum voltage-gated ion channel superfamily that reconciles apparently conflicting experimental data.
G-protein-coupled receptors (GPCRs) can modulate diverse signaling pathways, often in a ligand-specific manner. The full range of functionally relevant GPCR conformations is poorly understood. Here, ...we use NMR spectroscopy to characterize the conformational dynamics of the transmembrane core of the β2-adrenergic receptor (β2AR), a prototypical GPCR. We labeled β2AR with 13CH3ε-methionine and obtained HSQC spectra of unliganded receptor as well as receptor bound to an inverse agonist, an agonist, and a G-protein-mimetic nanobody. These studies provide evidence for conformational states not observed in crystal structures, as well as substantial conformational heterogeneity in agonist- and inverse-agonist-bound preparations. They also show that for β2AR, unlike rhodopsin, an agonist alone does not stabilize a fully active conformation, suggesting that the conformational link between the agonist-binding pocket and the G-protein-coupling surface is not rigid. The observed heterogeneity may be important for β2AR’s ability to engage multiple signaling and regulatory proteins.
Display omitted
► NMR using 13CH3-ε-Met reveals dynamics of β2 adrenergic receptor (β2AR) ► NMR and computational approaches show unanticipated conformational states ► Conformational heterogeneity is observed in both unliganded and antagonist-bound β2AR ► Agonist alone does not fully stabilize the active conformation of the β2AR
NMR studies reveal that β2 adrenergic receptor experiences conformational heterogeneity until both agonist and a G protein mimetic bind, shifting the receptor into its active conformation.
Molecular dynamics (MD) simulations have become a powerful and popular method for the study of protein allostery, the widespread phenomenon in which a stimulus at one site on a protein influences the ...properties of another site on the protein. By capturing the motions of a protein's constituent atoms, simulations can enable the discovery of allosteric binding sites and the determination of the mechanistic basis for allostery. These results can provide a foundation for applications including rational drug design and protein engineering. Here, we provide an introduction to the investigation of protein allostery using molecular dynamics simulation. We emphasize the importance of designing simulations that include appropriate perturbations to the molecular system, such as the addition or removal of ligands or the application of mechanical force. We also demonstrate how the bidirectional nature of allostery-the fact that the two sites involved influence one another in a symmetrical manner-can facilitate such investigations. Through a series of case studies, we illustrate how these concepts have been used to reveal the structural basis for allostery in several proteins and protein complexes of biological and pharmaceutical interest.