Studies of folding and misfolding using simplified models

Computer simulations are as vital to our studies of biological systems as experiments. They bridge and rationalize experimental observations, extend the experimental "field of view", which is often limited to a specific time or length scale, and, most importantly, provide novel insights into biological systems, offering hypotheses about yet-to-be uncovered phenomena. These hypotheses spur further experimental discoveries. Simplified molecular models have a special place in the field of computational biology. Branded as less accurate than all-atom protein models, they have offered what all-atom molecular dynamics simulations could not--the resolution of the length and time scales of biological phenomena. Not only have simplified models proven to be accurate in explaining or reproducing several biological phenomena, they have also offered a novel multiscale computational strategy for accessing a broad range of time and length scales upon integration with traditional all-atom simulations. Recent computer simulations of simplified models have shaken or advanced the established understanding of biological phenomena. It was demonstrated that simplified models can be as accurate as traditional molecular dynamics approaches in identifying native conformations of proteins. Their application to protein structure prediction yielded phenomenal accuracy in recapitulating native protein conformations. New studies that utilize the synergy of simplified protein models with all-atom models and experiments yielded novel insights into complex biological processes, such as protein folding, aggregation and the formation of large protein complexes.     

Insights into allosteric control of vinculin function from its large scale conformational dynamics

Vinculin is an important constituent of both cell-cell and cell-matrix junctions, where it plays crucial roles in the regulation of cell adhesion and migration. When activated, it mediates the linkage between cadherins (cell-cell) or integrins (cell-matrix) and the actin cytoskeleton through interactions with various proteins. The activation of vinculin requires structural conversions from an autoinhibited conformation to the "open" conformations in which the occluded binding sites of its different ligands become exposed, while the structural dynamics underlying the vinculin activation remains largely unknown. Here we report the first computational study of large scale conformational dynamics of full-length vinculin. We find that the "holding" and "releasing" motions between vinculin tail and pincer-like structure formed by first three domains of vinculin are the dominant motions near the native state of vinculin, indicating that an inherent flexibility of vinculin has a large influence on its allostery. We also find a cooperative dissociation between the head and tail domains of vinculin with increasing temperature in both thermodynamic and kinetic simulations, implying that vinculin may function as an allosteric switch in response to external signals. We show that the kinetics of vinculin unfolding exhibits specific sequential patterns, suggesting that a sophisticated interplay between domains may synergistically contribute to vinculin activation. We further find that the interaction between vinculin-binding site peptide from talin and vinculin significantly destabilizes the intramolecular head-tail interactions, suggesting a direct role of talin binding in vinculin activation.     

Three-dimensional RNA structure refinement by hydroxyl radical probing

Molecular modeling guided by experimentally derived structural information is an attractive approach for three-dimensional structure determination of complex RNAs that are not amenable to study by high-resolution methods. Hydroxyl radical probing (HRP), which is performed routinely in many laboratories, provides a measure of solvent accessibility at individual nucleotides. HRP measurements have, to date, only been used to evaluate RNA models qualitatively. Here we report the development of a quantitative structure refinement approach using HRP measurements to drive discrete molecular dynamics simulations for RNAs ranging in size from 80 to 230 nucleotides. We first used HRP reactivities to identify RNAs that form extensive helical packing interactions. For these RNAs, we achieved highly significant structure predictions given the inputs of RNA sequence and base pairing. This HRP-directed tertiary structure refinement approach generates robust structural hypotheses that are useful for guiding explorations of structure-function inter-relationships in RNA.     

Structural basis for the sequence-dependent effects of platinum–DNA adducts

The differences in efficacy and molecular mechanisms of platinum based anti-cancer drugs cisplatin (CP) and oxaliplatin (OX) have been hypothesized to be in part due to the differential binding affinity of cellular and damage recognition proteins to CP and OX adducts formed on adjacent guanines in genomic DNA. HMGB1a in particular exhibits higher binding affinity to CP-GG adducts, and the extent of discrimination between CP- and OX-GG adducts is dependent on the bases flanking the adducts. However, the structural basis for this differential binding is not known. Here, we show that the conformational dynamics of CP- and OX-GG adducts are distinct and depend on the sequence context of the adduct. Molecular dynamics simulations of the Pt-GG adducts in the TGGA sequence context revealed that even though the major conformations of CP- and OX-GG adducts were similar, the minor conformations were distinct. Using the pattern of hydrogen bond formation between the Pt–ammines and the adjacent DNA bases, we identified the major and minor conformations sampled by Pt–DNA. We found that the minor conformations sampled exclusively by the CP-GG adduct exhibit structural properties that favor binding by HMGB1a, which may explain its higher binding affinity to CP-GG adducts, while these conformations are not sampled by OX-GG adducts because of the constraints imposed by its cyclohexane ring, which may explain the negligible binding affinity of HMGB1a for OX-GG adducts in the TGGA sequence context. Based on these results, we postulate that the constraints imposed by the cyclohexane ring of OX affect the DNA conformations explored by OX-GG adduct compared to those of CP-GG adduct, which may influence the binding affinities of HMG-domain proteins for Pt-GG adducts, and that these conformations are further influenced by the DNA sequence context of the Pt-GG adduct.     

Protein folding: then and now

Over the past three decades the protein folding field has undergone monumental changes. Originally a purely academic question, how a protein folds has now become vital in understanding diseases and our abilities to rationally manipulate cellular life by engineering protein folding pathways. We review and contrast past and recent developments in the protein folding field. Specifically, we discuss the progress in our understanding of protein folding thermodynamics and kinetics, the properties of evasive intermediates, and unfolded states. We also discuss how some abnormalities in protein folding lead to protein aggregation and human diseases.     

Understanding hierarchical protein evolution from first principles.

We propose a model that explains the hierarchical organization of proteins in fold families. The model, which is based on the evolutionary selection of proteins by their native state stability, reproduces patterns of amino acids conserved across protein families. Due to its dynamic nature, the model sheds light on the evolutionary time-scales. By studying the relaxation of the correlation function between consecutive mutations at a given position in proteins, we observe separation of the evolutionary time-scales: at short time intervals families of proteins with similar sequences and structures are formed, while at long time intervals the families of structurally similar proteins that have low sequence similarity are formed. We discuss the evolutionary implications of our model. We provide a "profile" solution to our model and find agreement between predicted patterns of conserved amino acids and those actually observed in nature.     

Molecular dynamics simulation of amyloid beta dimer formation

Recent experiments with amyloid beta (Abeta) peptide indicate that formation of toxic oligomers may be an important contribution to the onset of Alzheimer's disease. The toxicity of Abeta oligomers depends on their structure, which is governed by assembly dynamics. Due to limitations of current experimental techniques, a detailed knowledge of oligomer structure at the atomic level is missing. We introduce a molecular dynamics approach to study Abeta dimer formation. 1), We use discrete molecular dynamics simulations of a coarse-grained model to identify a variety of dimer conformations; and 2), we employ all-atom molecular mechanics simulations to estimate thermodynamic stability of all dimer conformations. Our simulations of a coarse-grained Abeta peptide model predicts 10 different planar beta-strand dimer conformations. We then estimate the free energies of all dimer conformations in all-atom molecular mechanics simulations with explicit water. We compare the free energies of Abeta(1-42) and Abeta(1-40) dimers. We find that 1), dimer conformations have higher free energies compared to their corresponding monomeric states; and 2), the free-energy difference between the Abeta(1-42) and the corresponding Abeta(1-40) dimer conformation is not significant. Our results suggest that Abeta oligomerization is not accompanied by the formation of thermodynamically stable planar beta-strand dimers.     

Mechanism for the alpha-helix to beta-hairpin transition

The aggregation of alpha-helix-rich proteins into beta-sheet-rich amyloid fibrils is associated with fatal diseases, such as Alzheimer's disease and prion disease. During an aggregation process, protein secondary structure elements-alpha-helices-undergo conformational changes to beta-sheets. The fact that proteins with different sequences and structures undergo a similar transition on aggregation suggests that the sequence nonspecific hydrogen bond interaction among protein backbones is an important factor. We perform molecular dynamics simulations of a polyalanine model, which is an alpha-helix in its native state and observe a metastable beta-hairpin intermediate. Although a beta-hairpin has larger potential energy than an alpha-helix, the entropy of a beta-hairpin is larger because of fewer constraints imposed by the hydrogen bonds. In the vicinity of the transition temperature, we observe the interconversion of the alpha-helix and beta-sheet states via a random coil state. We also study the effect of the environment by varying the relative strength of side-chain interactions for a designed peptide-an alpha-helix in its native state. For a certain range of side-chain interaction strengths, we find that the intermediate beta-hairpin state is destabilized and even disappears, suggesting an important role of the environment in the aggregation propensity of a peptide.     

Multiscale modeling of nucleosome dynamics

Nucleosomes form the fundamental building blocks of chromatin. Subtle modifications of the constituent histone tails mediate chromatin stability and regulate gene expression. For this reason, it is important to understand structural dynamics of nucleosomes at atomic levels. We report a novel multiscale model of the fundamental chromatin unit, a nucleosome, using a simplified model for rapid discrete molecular dynamics simulations and an all-atom model for detailed structural investigation. Using a simplified structural model, we perform equilibrium simulations of a single nucleosome at various temperatures. We further reconstruct all-atom nucleosome structures from simulation trajectories. We find that histone tails bind to nucleosomal DNA via strong salt-bridge interactions over a wide range of temperatures, suggesting a mechanism of chromatin structural organization whereby histone tails regulate inter- and intranucleosomal assemblies via binding with nucleosomal DNA. We identify specific regions of the histone core H2A/H2B-H4/H3-H3/H4-H2B/H2A, termed “cold sites”, which retain a significant fraction of contacts with adjoining residues throughout the simulation, indicating their functional role in nucleosome organization. Cold sites are clustered around H3-H3, H2A-H4 and H4-H2A interhistone interfaces, indicating the necessity of these contacts for nucleosome stability. Essential dynamics analysis of simulation trajectories shows that bending across the H3-H3 is a prominent mode of intranucleosomal dynamics. We postulate that effects of salts on mononucleosomes can be modeled in discrete molecular dynamics by modulating histone-DNA interaction potentials. Local fluctuations in nucleosomal DNA vary significantly along the DNA sequence, suggesting that only a fraction of histone-DNA contacts make strong interactions dominating mononucleosomal dynamics. Our findings suggest that histone tails have a direct functional role in stabilizing higher-order chromatin structure, mediated by salt-bridge interactions with adjacent DNA.     

Stabilization of μ-opioid receptor facilitates its cellular translocation and signaling

The G protein-coupled μ-opioid receptor (μ-OR) mediates the majority of analgesia effects for morphine and other pain relievers. Despite extensive studies of its structure and activation mechanisms, the inherently low maturation efficiency of μ-OR represents a major hurdle to understanding its function. Here we computationally designed μ-OR mutants with altered stability to probe the relationship between cell-surface targeting, signal transduction, and agonist efficacy. The stabilizing mutation T315Y enhanced μ-OR trafficking to the plasma membrane and significantly promoted the morphine-mediated inhibition of downstream signaling. In contrast, the destabilizing mutation R165Y led to intracellular retention of μ-OR and reduced the response to morphine stimulation. These findings suggest that μ-OR stability is an important factor in regulating receptor signaling and provide a viable avenue to improve the efficacy of analgesics.