1
|
Zykaj E, Abboud C, Asadi P, Warsame S, Almousa H, Milev MP, Greco BM, López-Sánchez M, Bratkovic D, Kachroo AH, Pérez-Jurado LA, Sacher M. A Humanized Yeast Model for Studying TRAPP Complex Mutations; Proof-of-Concept Using Variants from an Individual with a TRAPPC1-Associated Neurodevelopmental Syndrome. Cells 2024; 13:1457. [PMID: 39273027 PMCID: PMC11394476 DOI: 10.3390/cells13171457] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2024] [Revised: 08/22/2024] [Accepted: 08/25/2024] [Indexed: 09/15/2024] Open
Abstract
Variants in membrane trafficking proteins are known to cause rare disorders with severe symptoms. The highly conserved transport protein particle (TRAPP) complexes are key membrane trafficking regulators that are also involved in autophagy. Pathogenic genetic variants in specific TRAPP subunits are linked to neurological disorders, muscular dystrophies, and skeletal dysplasias. Characterizing these variants and their phenotypes is important for understanding the general and specialized roles of TRAPP subunits as well as for patient diagnosis. Patient-derived cells are not always available, which poses a limitation for the study of these diseases. Therefore, other systems, like the yeast Saccharomyces cerevisiae, can be used to dissect the mechanisms at the intracellular level underlying these disorders. The development of CRISPR/Cas9 technology in yeast has enabled a scar-less editing method that creates an efficient humanized yeast model. In this study, core yeast subunits were humanized by replacing them with their human orthologs, and TRAPPC1, TRAPPC2, TRAPPC2L, TRAPPC6A, and TRAPPC6B were found to successfully replace their yeast counterparts. This system was used for studying the first reported individual with an autosomal recessive disorder caused by biallelic TRAPPC1 variants, a girl with a severe neurodevelopmental disorder and myopathy. We show that the maternal variant (TRAPPC1 p.(Val121Alafs*3)) is non-functional while the paternal variant (TRAPPC1 p.(His22_Lys24del)) is conditional-lethal and affects secretion and non-selective autophagy in yeast. This parallels defects seen in fibroblasts derived from this individual which also showed membrane trafficking defects and altered Golgi morphology, all of which were rescued in the human system by wild-type TRAPPC1. This study suggests that humanized yeast can be an efficient means to study TRAPP subunit variants in the absence of human cells and can assign significance to variants of unknown significance (VUS). This study lays the foundation for characterizing further TRAPP variants through this system, rapidly contributing to disease diagnosis.
Collapse
Affiliation(s)
- Erta Zykaj
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Chelsea Abboud
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Paria Asadi
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Simane Warsame
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Hashem Almousa
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Miroslav P Milev
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Brittany M Greco
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Marcos López-Sánchez
- Department of Medicine and Life Sciences, Universitat Pompeu Fabra, 08003 Barcelona, Spain
- Hospital del Mar, Hospital del Mar Research Institute (IMIM), 08003 Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), ISCIII, 28029 Madrid, Spain
| | - Drago Bratkovic
- Women's and Children's Hospital, Metabolic Clinic, North Adelaide, SA 5006, Australia
| | - Aashiq H Kachroo
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
| | - Luis Alberto Pérez-Jurado
- Department of Medicine and Life Sciences, Universitat Pompeu Fabra, 08003 Barcelona, Spain
- Hospital del Mar, Hospital del Mar Research Institute (IMIM), 08003 Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER), ISCIII, 28029 Madrid, Spain
- Women's and Children's Hospital, Metabolic Clinic, North Adelaide, SA 5006, Australia
| | - Michael Sacher
- Department of Biology, Concordia University, Montreal, QC H4B1R6, Canada
- Department of Anatomy and Cell Biology, McGill University, Montreal, QC H3A 0C7, Canada
| |
Collapse
|
2
|
Chennakesavalu S, Rotskoff GM. Data-Efficient Generation of Protein Conformational Ensembles with Backbone-to-Side-Chain Transformers. J Phys Chem B 2024; 128:2114-2123. [PMID: 38394363 DOI: 10.1021/acs.jpcb.3c08195] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/25/2024]
Abstract
Excitement at the prospect of using data-driven generative models to sample configurational ensembles of biomolecular systems stems from the extraordinary success of these models on a diverse set of high-dimensional sampling tasks. Unlike image generation or even the closely related problem of protein structure prediction, there are currently no data sources with sufficient breadth to parametrize generative models for conformational ensembles. To enable discovery, a fundamentally different approach to building generative models is required: models should be able to propose rare, albeit physical, conformations that may not arise in even the largest data sets. Here we introduce a modular strategy to generate conformations based on "backmapping" from a fixed protein backbone that (1) maintains conformational diversity of the side chains and (2) couples the side-chain fluctuations using global information about the protein conformation. Our model combines simple statistical models of side-chain conformations based on rotamer libraries with the now ubiquitous transformer architecture to sample with atomistic accuracy. Together, these ingredients provide a strategy for rapid data acquisition and hence a crucial ingredient for scalable physical simulation with generative neural networks.
Collapse
Affiliation(s)
| | - Grant M Rotskoff
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
- Institute for Computational and Mathematical Engineering, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
3
|
Kulshreshtha A, Bhatnagar S. Structural effect of the H992D/H418D mutation of angiotensin-converting enzyme in the Indian population: implications for health and disease. J Biomol Struct Dyn 2024:1-18. [PMID: 38411559 DOI: 10.1080/07391102.2024.2321246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 02/14/2024] [Indexed: 02/28/2024]
Abstract
The Non synonymous SNPs (nsSNPs) of the renin-angiotensin-system (RAS) pathway, unique to the Indian population were investigated in view of its importance as an endocrine system. nsSNPs of the RAS pathway genes were mined from the IndiGenome database. Damaging nsSNPs were predicted using SIFT, PredictSNP, SNP and GO, Snap2 and Protein Variation Effect Analyzer. Loss of function was predicted based on protein stability change using I mutant, PremPS and CONSURF. The structural impact of the nsSNPs was predicted using HOPE and Missense3d followed by modeling, refinement, and energy minimization. Molecular Dynamics studies were carried out using Gromacsv2021.1. 23 Indian nsSNPs of the RAS pathway genes were selected for structural analysis and 8 were predicted to be damaging. Further sequence analysis showed that HEMGH zinc binding motif changes to HEMGD in somatic ACE-C domain (sACE-C) H992D and Testis ACE (tACE) H418D resulted in loss of zinc coordination, which is essential for enzymatic activity in this metalloprotease. There was a loss of internal interactions around the zinc coordination residues in the protein structural network. This was also confirmed by Principal Component Analysis, Free Energy Landscape and residue contact maps. Both mutations lead to broadening of the AngI binding cavity. The H992D mutation in sACE-C is likely to be favorable for cardiovascular health, but may lead to renal abnormalities with secondary impact on the heart. H418D in tACE is potentially associated with male infertility.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Akanksha Kulshreshtha
- Computational and Structural Biology Laboratory, Department of Biological Sciences and Engineering, Netaji Subhas University of Technology, Dwarka, New Delhi, India
| | - Sonika Bhatnagar
- Computational and Structural Biology Laboratory, Department of Biological Sciences and Engineering, Netaji Subhas University of Technology, Dwarka, New Delhi, India
| |
Collapse
|
4
|
Luque D, Ortega-Esteban A, Valbuena A, Luis Vilas J, Rodríguez-Huete A, Mateu MG, Castón JR. Equilibrium Dynamics of a Biomolecular Complex Analyzed at Single-amino Acid Resolution by Cryo-electron Microscopy. J Mol Biol 2023; 435:168024. [PMID: 36828271 DOI: 10.1016/j.jmb.2023.168024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2022] [Revised: 02/10/2023] [Accepted: 02/17/2023] [Indexed: 02/24/2023]
Abstract
The biological function of macromolecular complexes depends not only on large-scale transitions between conformations, but also on small-scale conformational fluctuations at equilibrium. Information on the equilibrium dynamics of biomolecular complexes could, in principle, be obtained from local resolution (LR) data in cryo-electron microscopy (cryo-EM) maps. However, this possibility had not been validated by comparing, for a same biomolecular complex, LR data with quantitative information on equilibrium dynamics obtained by an established solution technique. In this study we determined the cryo-EM structure of the minute virus of mice (MVM) capsid as a model biomolecular complex. The LR values obtained correlated with crystallographic B factors and with hydrogen/deuterium exchange (HDX) rates obtained by mass spectrometry (HDX-MS), a gold standard for determining equilibrium dynamics in solution. This result validated a LR-based cryo-EM approach to investigate, with high spatial resolution, the equilibrium dynamics of biomolecular complexes. As an application of this approach, we determined the cryo-EM structure of two mutant MVM capsids and compared their equilibrium dynamics with that of the wild-type MVM capsid. The results supported a previously suggested linkage between mechanical stiffening and impaired equilibrium dynamics of a virus particle. Cryo-EM is emerging as a powerful approach for simultaneously acquiring information on the atomic structure and local equilibrium dynamics of biomolecular complexes.
Collapse
Affiliation(s)
- Daniel Luque
- Spanish National Microbiology Centre, Institute of Health Carlos III, Madrid, Spain
| | - Alvaro Ortega-Esteban
- Department of Structure of Macromolecules, Centro Nacional de Biotecnología (CNB-CSIC), Campus de Cantoblanco, Madrid, Spain
| | - Alejandro Valbuena
- Centro de Biología Molecular "Severo Ochoa" (CSIC-UAM), Universidad Autónoma de Madrid, Madrid, Spain
| | - Jose Luis Vilas
- Biocomputing Unit, Centro Nacional de Biotecnología (CNB-CSIC), Campus de Cantoblanco, Madrid, Spain
| | - Alicia Rodríguez-Huete
- Centro de Biología Molecular "Severo Ochoa" (CSIC-UAM), Universidad Autónoma de Madrid, Madrid, Spain
| | - Mauricio G Mateu
- Centro de Biología Molecular "Severo Ochoa" (CSIC-UAM), Universidad Autónoma de Madrid, Madrid, Spain.
| | - José R Castón
- Department of Structure of Macromolecules, Centro Nacional de Biotecnología (CNB-CSIC), Campus de Cantoblanco, Madrid, Spain.
| |
Collapse
|
5
|
Bolik-Coulon N, Ferrage F. Explicit models of motions to analyze NMR relaxation data in proteins. J Chem Phys 2022; 157:125102. [DOI: 10.1063/5.0095910] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Nuclear Magnetic Resonance (NMR) is a tool of choice to characterize molecular motions. In biological macromolecules, pico- to nano-second motions, in particular, can be probed by nuclear spin relaxation rates which depend on the time fluctuations of the orientations of spin interaction frames. For the past 40 years, relaxation rates have been successfully analyzed using the Model Free (MF) approach which makes no assumption on the nature of motions and reports on the effective amplitude and time-scale of the motions. However, obtaining a mechanistic picture of motions from this type of analysis is difficult at best, unless complemented with molecular dynamics (MD) simulations. In spite of their limited accuracy, such simulations can be used to obtain the information necessary to build explicit models of motions designed to analyze NMR relaxation data. Here, we present how to build such models, suited in particular to describe motions of methyl-bearing protein side-chains and compare them with the MF approach. We show on synthetic data that explicit models of motions are more robust in the presence of rotamer jumps which dominate the relaxation in methyl groups of protein side-chains. We expect this work to motivate the use of explicit models of motion to analyze MD and NMR data.
Collapse
Affiliation(s)
| | - Fabien Ferrage
- Departement de chimie, Ecole Normale Superieure Departement de Chimie, France
| |
Collapse
|
6
|
A single-molecule stochastic theory of protein-ligand binding in the presence of multiple unfolding/folding and ligand binding pathways. Biophys Chem 2022; 285:106803. [DOI: 10.1016/j.bpc.2022.106803] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Revised: 03/17/2022] [Accepted: 03/17/2022] [Indexed: 11/19/2022]
|
7
|
The native state conformational heterogeneity in the energy landscape of protein folding. Biophys Chem 2022; 283:106761. [DOI: 10.1016/j.bpc.2022.106761] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2021] [Revised: 01/10/2022] [Accepted: 01/14/2022] [Indexed: 11/18/2022]
|
8
|
Why are large conformational changes well described by harmonic normal modes? Biophys J 2021; 120:5343-5354. [PMID: 34710378 DOI: 10.1016/j.bpj.2021.10.027] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2021] [Revised: 09/14/2021] [Accepted: 10/20/2021] [Indexed: 12/11/2022] Open
Abstract
Low-frequency normal modes generated by elastic network models tend to correlate strongly with large conformational changes of proteins, despite their reliance on the harmonic approximation, which is only valid in close proximity of the native structure. We consider 12 variants of the torsional network model (TNM), an elastic network model in torsion angle space, that adopt different sets of torsion angles as degrees of freedom and reproduce with similar quality the thermal fluctuations of proteins but present drastic differences in their agreement with conformational changes. We show that these differences are related to the extent of the deviations from the harmonic approximation, assessed through an anharmonic energy function whose harmonic approximation coincides with the TNM. Our results indicate that mode anharmonicity is more strongly related to its collectivity, i.e., the number of atoms displaced by the mode, than to its amplitude; low-frequency modes can remain harmonic even at large amplitudes, provided they are sufficiently collective. Finally, we assess the potential benefits of different strategies to minimize the impact of anharmonicity. The reduction of the number of degrees of freedom or their regularization by a torsional harmonic potential significantly improves the collectivity and harmonicity of normal modes and the agreement with conformational changes. In contrast, the correction of normal mode frequencies to partially account for anharmonicity does not yield substantial benefits. The TNM program is freely available at https://github.com/ugobas/tnm.
Collapse
|
9
|
Autoencoder-based detection of the residues involved in G protein-coupled receptor signaling. Sci Rep 2021; 11:19867. [PMID: 34615896 PMCID: PMC8494915 DOI: 10.1038/s41598-021-99019-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2021] [Accepted: 09/17/2021] [Indexed: 11/08/2022] Open
Abstract
Regulator binding and mutations alter protein dynamics. The transmission of the signal of these alterations to distant sites through protein motion results in changes in protein expression and cell function. The detection of residues involved in signal transmission contributes to an elucidation of the mechanisms underlying processes as vast as cellular function and disease pathogenesis. We developed an autoencoder (AE) based method that detects residues essential for signaling by comparing the fluctuation data, particularly the time fluctuation of the side-chain distances between residues, during molecular dynamics simulations between the ligand-bound and -unbound forms or wild-type and mutant forms of proteins. Here, the AE-based method was applied to the G protein-coupled receptor (GPCR) system, particularly a class A-type GPCR, CXCR4, to detect the essential residues involved in signaling. Among the residues involved in the signaling of the homolog CXCR2, which were extracted from the literature based on the complex structures of the ligand and G protein, our method could detect more than half of the essential residues involved in G protein signaling, including those spanning the fifth and sixth transmembrane helices in the intracellular region, despite the lack of information regarding the interaction with G protein in our CXCR4 models.
Collapse
|
10
|
Chekmarev SF. First-passage times in protein folding: exploring the native-like states vs. overcoming the free energy barrier. Phys Chem Chem Phys 2021; 23:17856-17865. [PMID: 34378547 DOI: 10.1039/d0cp06560a] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Using a model β-hairpin protein as a representative example of simple two-state folders whose kinetics are uncomplicated by the presence of on- and off-pathway intermediates, it is studied how the search for the protein's native state among native-like states affects the folding kinetics. It is revealed that the first-passage time (FPT) distributions are essentially single-exponential not only for the times to overcome the free energy barrier between the unfolded and native-like states but also for the times to find the native state among the native-like ones. The FPT distributions of this type are observed through all studied two-state-like regimes of protein folding, varying from a regime close to two-state folding to a regime close to downhill folding. If the protein explores native-like states for a time much longer than the time to overcome the free energy barrier, which is characteristic of high temperatures, the resulting FPT distribution to reach the native state remains close to exponential but the mean FPT (MFPT) is determined not by the height of the free energy barrier but by the time to explore native-like states. In particular, the mean time to overcome the free energy barrier is in reasonable agreement with the Kramers rate formula and generally far shorter than the overall MFPT to reach the native state. The observed increase of the overall MFPT, as a result of longer exploration of native-like states, may lead to an overestimate of the height of the free energy barrier between the unfolded and folded states when it is calculated from the overall MFPT.
Collapse
|
11
|
Miyashita N, Yonezawa Y. Mutual information analysis of the dynamic correlation between side chains in proteins. J Chem Phys 2021; 155:044107. [PMID: 34340388 DOI: 10.1063/5.0055662] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023] Open
Abstract
Protein dynamics play an essential role in function regulation. In recent years, many experimental and theoretical studies have shown that changes in protein fluctuations in the backbone and side chains fulfill a pivotal role associated with amino acid mutations, chemical modifications, and ligand binding. The dynamic correlations between protein side chains have not been sufficiently studied, and no reliable analysis method has been available so far. Therefore, we developed a method to evaluate the dynamic correlation between protein side chains using mutual information and molecular dynamics simulations. To eliminate the structural superposition errors dealing with conventional analysis methods, and to accurately extract the intrinsic fluctuation properties of the side chains, we employed distance principal component analysis (distPCA). The motion of the side chain was then projected onto the eigenvector space obtained by distPCA, and the mutual information between the projected motions was calculated. The proposed method was then applied to a small protein "eglin c" and the mutants. The results show that even a single mutation significantly changed the dynamic correlations and also suggest that the dynamic change is deeply related to the stability. Those results indicate that our developed method could be useful for analyzing the molecular mechanism of allosteric communication in proteins.
Collapse
Affiliation(s)
- Naoyuki Miyashita
- Department of Computational Systems Biology, Biology-Oriented Science and Technology, Kindai University, 930 Nishimitani, Kinokawa, Wakayama 649-6493, Japan
| | - Yasushige Yonezawa
- High Pressure Protein Research Center, Institute of Advanced Technology, Kindai University, 930 Nishimitani, Kinokawa, Wakayama 649-6493, Japan
| |
Collapse
|
12
|
Dhabal D, Jiang Z, Pallath A, Patel AJ. Characterizing the Interplay between Polymer Solvation and Conformation. J Phys Chem B 2021; 125:5434-5442. [PMID: 33978411 DOI: 10.1021/acs.jpcb.1c02191] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Conformational transitions of flexible molecules, especially those driven by hydrophobic effects, tend to be hindered by desolvation barriers. For such transitions, it is thus important to characterize and understand the interplay between solvation and conformation. Using specialized molecular simulations, here we perform such a characterization for a hydrophobic polymer solvated in water. We find that an external potential, which unfavorably perturbs the polymer hydration waters, can trigger a coil-to-globule or collapse transition, and that the relative stabilities of the collapsed and extended states can be quantified by the strength of the requisite potential. Our results also provide mechanistic insights into the collapse transition, highlighting that the bottleneck to polymer collapse is the formation of a sufficiently large cluster, and the collective dewetting of such a cluster. We also study the collapse of the hydrophobic polymer in octane, a nonpolar solvent, and interestingly, we find that the mechanistic details of the transition are qualitatively similar to that in water.
Collapse
Affiliation(s)
- Debdas Dhabal
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Zhitong Jiang
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Akash Pallath
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Amish J Patel
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
13
|
Prabantu VM, Naveenkumar N, Srinivasan N. Influence of Disease-Causing Mutations on Protein Structural Networks. Front Mol Biosci 2021; 7:620554. [PMID: 33778000 PMCID: PMC7987782 DOI: 10.3389/fmolb.2020.620554] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Accepted: 12/17/2020] [Indexed: 01/18/2023] Open
Abstract
The interactions between residues in a protein tertiary structure can be studied effectively using the approach of protein structure network (PSN). A PSN is a node-edge representation of the structure with nodes representing residues and interactions between residues represented by edges. In this study, we have employed weighted PSNs to understand the influence of disease-causing mutations on proteins of known 3D structures. We have used manually curated information on disease mutations from UniProtKB/Swiss-Prot and their corresponding protein structures of wildtype and disease variant from the protein data bank. The PSNs of the wildtype and disease-causing mutant are compared to analyse variation of global and local dissimilarity in the overall network and at specific sites. We study how a mutation at a given site can affect the structural network at a distant site which may be involved in the function of the protein. We have discussed specific examples of the disease cases where the protein structure undergoes limited structural divergence in their backbone but have large dissimilarity in their all atom networks and vice versa, wherein large conformational alterations are observed while retaining overall network. We analyse the effect of variation of network parameters that characterize alteration of function or stability.
Collapse
Affiliation(s)
| | - Nagarajan Naveenkumar
- Molecular Biophysics Unit, Indian Institute of Science, Bangalore, India.,National Centre for Biological Sciences, TIFR, Bangalore, India.,Bharathidasan University, Tiruchirappalli, India
| | | |
Collapse
|
14
|
Azmi MB, Sultana S, Naeem S, Qureshi SA. In silico investigation on alkaloids of Rauwolfia serpentina as potential inhibitors of 3-hydroxy-3-methyl-glutaryl-CoA reductase. Saudi J Biol Sci 2021; 28:731-737. [PMID: 33424361 PMCID: PMC7783793 DOI: 10.1016/j.sjbs.2020.10.066] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Revised: 10/20/2020] [Accepted: 10/29/2020] [Indexed: 11/26/2022] Open
Abstract
Present work aimed to investigate the in silico activity of the alkaloids of roots of Rauwolfia serpentina as inhibitors of 3-hydroxy-3-methyl-glutaryl-CoA reductase (HMGCR). For this purpose, the three-dimensional (3D) structure of the protein HMGCR (PDB ID: 1HW9) was downloaded from Protein Data Bank (PDB) database, as a target enzyme. The structures of twelve alkaloids from the roots of R. serpentina were selected as ligands and docked with the selected HMGCR enzyme using Molegro Virtual Docker (MVD) software. The software ‘MVD’ computes the binding (atom) energies of selected protein (enzyme) and each ligand at minimum energetic conformation state by using the PLP (Piecewise Linear Potential) scoring mechanism. Docking results of twelve tested alkaloids showed that five alkaloids including compound 1 (ajmalicine), 2 (reserpine), 3 (indobinine), 4 (yohimbine), and 5 (indobine) have displayed the highest MolDock scores and best fit within the prominent active site residues (positioned between 684 and 692 of cis-loop) of HMGCR. According to the lowest MolDock energies obtained through non-covalent interactions of alkaloids with HMGCR, these are characterized to be the potential inhibitors of HMGCR. Therefore, the alkaloids from R. serpentina can effectively suppress the cholesterol biosynthesis pathway through inhibition of HMGCR and can serve as potential lead compounds for the development of new drugs for the treatment of hyperlipidaemia.
Collapse
Affiliation(s)
- Muhammad Bilal Azmi
- Department of Biochemistry, Dow Medical College, Dow University of Health Sciences, Karachi 74200, Pakistan
- Corresponding author.
| | - Saleha Sultana
- Department of Biochemistry, University of Karachi, Karachi 75270, Pakistan
| | - Sadaf Naeem
- Department of Biochemistry, University of Karachi, Karachi 75270, Pakistan
| | | |
Collapse
|
15
|
Novel insights in linking solvent relaxation dynamics and protein conformations utilizing red edge excitation shift approach. Emerg Top Life Sci 2021; 5:89-101. [PMID: 33416893 DOI: 10.1042/etls20200256] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2020] [Revised: 12/14/2020] [Accepted: 12/18/2020] [Indexed: 01/01/2023]
Abstract
Protein hydration dynamics plays an important role in many physiological processes since protein fluctuations, slow solvation, and the dynamics of hydrating water are all intrinsically related. Red edge excitation shift (REES) is a unique and powerful wavelength-selective (i.e. excitation-energy dependent) fluorescence approach that can be used to directly monitor the environment-induced restriction and dynamics around a polar fluorophore in a complex biological system. This review is mainly focused on recent applications of REES and a novel analysis of REES data to monitor the structural dynamics, functionally relevant conformational transitions and to unmask the structural ensembles in proteins. In addition, the novel utility of REES in imaging protein aggregates in a cellular context is discussed. We believe that the enormous potential of REES approach showcased in this review will engage more researchers, particularly from life sciences.
Collapse
|
16
|
Marques BS, Stetz MA, Jorge C, Valentine KG, Wand AJ, Nucci NV. Protein conformational entropy is not slaved to water. Sci Rep 2020; 10:17587. [PMID: 33067552 PMCID: PMC7567893 DOI: 10.1038/s41598-020-74382-5] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2020] [Accepted: 09/23/2020] [Indexed: 12/19/2022] Open
Abstract
Conformational entropy can be an important element of the thermodynamics of protein functions such as the binding of ligands. The observed role for conformational entropy in modulating molecular recognition by proteins is in opposition to an often-invoked theory for the interaction of protein molecules with solvent water. The "solvent slaving" model predicts that protein motion is strongly coupled to various aspects of water such as bulk solvent viscosity and local hydration shell dynamics. Changes in conformational entropy are manifested in alterations of fast internal side chain motion that is detectable by NMR relaxation. We show here that the fast-internal side chain dynamics of several proteins are unaffected by changes to the hydration layer and bulk water. These observations indicate that the participation of conformational entropy in protein function is not dictated by the interaction of protein molecules and solvent water under the range of conditions normally encountered.
Collapse
Affiliation(s)
- Bryan S Marques
- Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA, USA
| | - Matthew A Stetz
- Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA, USA
| | - Christine Jorge
- Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA, USA
| | - Kathleen G Valentine
- Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA, USA
| | - A Joshua Wand
- Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA, USA. .,Department of Biochemistry and Biophysics, Texas A&M University, College Station, TX, 77845-2128, USA.
| | - Nathaniel V Nucci
- Department of Physics and Astronomy and Department of Molecular and Cellular Biosciences, Rowan University, 201 Mullica Hill Road, Glassboro, NJ, 08028, USA.
| |
Collapse
|
17
|
A beginner's guide to molecular dynamics simulations and the identification of cross-correlation networks for enzyme engineering. Methods Enzymol 2020; 643:15-49. [PMID: 32896280 DOI: 10.1016/bs.mie.2020.04.020] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
The functional properties of proteins are decided not only by their relatively rigid overall structures, but even more importantly, by their dynamic properties. In a protein, some regions of structure exhibit highly correlated or anti-correlated motions with others, some are highly dynamic but uncorrelated, while other regions are relatively static. The residues with correlated or anti-correlated motions can form a so-called dynamic cross-correlation network, through which information can be transmitted. Such networks have been shown to be critical to allosteric transitions, and ligand binding, and have also been shown to be able to mediate epistatic interactions between mutations. As a result, they are likely to play a significant role in the development of new enzyme engineering strategies. In this chapter, protocols are provided for the assessment of dynamic cross-correlation networks, and for their application in protein engineering. Transketolase from E. coli is used as a model and the software GROMACS is applied for carrying out MD simulations to generate trajectories containing structural ensembles. The trajectory is then used for a dynamic cross correlation analysis using the R package, Bio3D. A matrix of all atom-wise cross-correlation coefficients is finally obtained, which can be displayed in a graphical representation termed a dynamical cross-correlation matrix.
Collapse
|
18
|
Shi J, Shen Q, Cho JH, Hwang W. Entropy Hotspots for the Binding of Intrinsically Disordered Ligands to a Receptor Domain. Biophys J 2020; 118:2502-2512. [PMID: 32311315 DOI: 10.1016/j.bpj.2020.03.026] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2019] [Revised: 02/28/2020] [Accepted: 03/23/2020] [Indexed: 11/18/2022] Open
Abstract
Proline-rich motifs (PRMs) are widely used for mediating protein-protein interactions with weak binding affinities. Because they are intrinsically disordered when unbound, conformational entropy plays a significant role for the binding. However, residue-level differences of the entropic contribution in the binding of different ligands remain not well understood. We use all-atom molecular dynamics simulation and the maximal information spanning tree formalism to analyze conformational entropy associated with the binding of two PRMs, one from the Abl kinase and the other from the nonstructural protein 1 of the 1918 Spanish influenza A virus, to the N-terminal SH3 (nSH3) domain of the CrkII protein. Side chains of the stably folded nSH3 experience more entropy change upon ligand binding than the backbone, whereas PRMs involve comparable but heterogeneous entropy changes among the backbone and side chains. In nSH3, two conserved nonpolar residues forming contacts with the PRM experience the largest side-chain entropy loss. In contrast, the C-terminal charged residues of PRMs that form polar contacts with nSH3 experience the greatest side-chain entropy loss, although their "fuzzy" nature is attributable to the backbone that remains relatively flexible. Thus, residues that form high-occupancy contacts between nSH3 and PRM do not reciprocally contribute to entropy loss. Furthermore, certain surface residues of nSH3 distal to the interface with PRMs gain entropy, indicating a nonlocal effect of ligand binding. Comparing between the PRMs from cAbl and nonstructural protein 1, the latter involves a larger side-chain entropy loss and forms more contacts with nSH3. Consistent with experiments, this indicates stronger binding of the viral ligand at the expense of losing the flexibility of side chains, whereas the backbone experiences less entropy loss. The entropy "hotspots" as identified in this study will be important for tuning the binding affinity of various ligands to a receptor.
Collapse
Affiliation(s)
- Jie Shi
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas
| | - Qingliang Shen
- Department of Biochemistry and Biophysics, Texas A&M University, College Station, Texas
| | - Jae-Hyun Cho
- Department of Biochemistry and Biophysics, Texas A&M University, College Station, Texas.
| | - Wonmuk Hwang
- Department of Biomedical Engineering, Texas A&M University, College Station, Texas; Department of Materials Science and Engineering, Texas A&M University, College Station, Texas; Department of Physics and Astronomy, Texas A&M University, College Station, Texas; School of Computational Sciences, Korea Institute for Advanced Study, Seoul, South Korea.
| |
Collapse
|
19
|
Bolik-Coulon N, Kadeřávek P, Pelupessy P, Dumez JN, Ferrage F, Cousin SF. Theoretical and computational framework for the analysis of the relaxation properties of arbitrary spin systems. Application to high-resolution relaxometry. JOURNAL OF MAGNETIC RESONANCE (SAN DIEGO, CALIF. : 1997) 2020; 313:106718. [PMID: 32234674 DOI: 10.1016/j.jmr.2020.106718] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2020] [Revised: 03/11/2020] [Accepted: 03/13/2020] [Indexed: 06/11/2023]
Abstract
A wide variety of nuclear magnetic resonance experiments rely on the prediction and analysis of relaxation processes. Recently, innovative approaches have been introduced where the sample travels through a broad range of magnetic fields in the course of the experiment, such as dissolution dynamic nuclear polarization or high-resolution relaxometry. Understanding the relaxation properties of nuclear spin systems over orders of magnitude of magnetic fields is essential to rationalize the results of these experiments. For example, during a high-resolution relaxometry experiment, the absence of control of nuclear spin relaxation pathways during the sample transfers and relaxation delays leads to systematic deviations of polarization decays from an ideal mono-exponential decay with the pure longitudinal relaxation rate. These deviations have to be taken into account to describe quantitatively the dynamics of the system. Here, we present computational tools to (1) calculate analytical expressions of relaxation rates for a broad variety of spin systems and (2) use these analytical expressions to correct the deviations arising in high-resolution relaxometry experiments. These tools lead to a better understanding of nuclear spin relaxation, which is required to improve the sensitivity of many pulse sequences, and to better characterize motions in macromolecules.
Collapse
Affiliation(s)
- Nicolas Bolik-Coulon
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France.
| | - Pavel Kadeřávek
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Philippe Pelupessy
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | | | - Fabien Ferrage
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France.
| | - Samuel F Cousin
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France.
| |
Collapse
|
20
|
Gagnér VA, Lundholm I, Garcia-Bonete MJ, Rodilla H, Friedman R, Zhaunerchyk V, Bourenkov G, Schneider T, Stake J, Katona G. Clustering of atomic displacement parameters in bovine trypsin reveals a distributed lattice of atoms with shared chemical properties. Sci Rep 2019; 9:19281. [PMID: 31848402 PMCID: PMC6917748 DOI: 10.1038/s41598-019-55777-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2019] [Accepted: 11/06/2019] [Indexed: 12/14/2022] Open
Abstract
Low-frequency vibrations are crucial for protein structure and function, but only a few experimental techniques can shine light on them. The main challenge when addressing protein dynamics in the terahertz domain is the ubiquitous water that exhibit strong absorption. In this paper, we observe the protein atoms directly using X-ray crystallography in bovine trypsin at 100 K while irradiating the crystals with 0.5 THz radiation alternating on and off states. We observed that the anisotropy of atomic displacements increased upon terahertz irradiation. Atomic displacement similarities developed between chemically related atoms and between atoms of the catalytic machinery. This pattern likely arises from delocalized polar vibrational modes rather than delocalized elastic deformations or rigid-body displacements. The displacement correlation between these atoms were detected by a hierarchical clustering method, which can assist the analysis of other ultra-high resolution crystal structures. These experimental and analytical tools provide a detailed description of protein dynamics to complement the structural information from static diffraction experiments.
Collapse
Affiliation(s)
- Viktor Ahlberg Gagnér
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
| | - Ida Lundholm
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
| | | | - Helena Rodilla
- Department of Microtechnology and Nanoscience, Chalmers University of Technology, Gothenburg, Sweden
| | - Ran Friedman
- Department of Chemistry and Biomedical Sciences, Linnaeus University, Kalmar, Sweden
| | | | - Gleb Bourenkov
- European Molecular Biology Laboratory Hamburg Outstation, EMBL c/o DESY, Notkestrasse 85, 22603, Hamburg, Germany
| | - Thomas Schneider
- European Molecular Biology Laboratory Hamburg Outstation, EMBL c/o DESY, Notkestrasse 85, 22603, Hamburg, Germany
| | - Jan Stake
- Department of Microtechnology and Nanoscience, Chalmers University of Technology, Gothenburg, Sweden
| | - Gergely Katona
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden.
| |
Collapse
|
21
|
Bastolla U, Dehouck Y. Can Conformational Changes of Proteins Be Represented in Torsion Angle Space? A Study with Rescaled Ridge Regression. J Chem Inf Model 2019; 59:4929-4941. [PMID: 31600071 DOI: 10.1021/acs.jcim.9b00627] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Torsion angles are the natural degrees of freedom of protein structures. The ability to determine torsional variations corresponding to observed changes in Cartesian coordinates is highly valuable, notably to investigate the mechanisms of functional conformational changes or to develop computational models of protein dynamics. This issue is far from trivial in practice since the impact of modifying one torsion angle strongly depends on all other angles, and the compounding effects of small variations in bond lengths and valence angles can completely disrupt a protein fold. We demonstrate that naive strategies, such as directly comparing torsion angles between structures without correcting for variations in bond lengths and valence angles or fitting torsional variations without a proper regularization scheme, fail at producing an adequate representation of conformational changes in internal coordinates. In contrast, rescaled ridge regression, a method recently introduced to regularize multidimensional regressions with correlated explanatory variables, is shown to consistently identify a minimal set of torsion angles variations that closely reproduce changes in Cartesian coordinates. This torsional representation of conformational changes is shown to be robust to the choice of experimental structures. It also provides a better agreement with theoretical models of protein dynamics than the Cartesian representation, regarding notably the predominance of low-frequency normal modes in functional motions and the presence, in predicted equilibrium dynamics, of hints of natural selection for specific functional motions. The software is available at https://github.com/ugobas/tnm .
Collapse
Affiliation(s)
- Ugo Bastolla
- Centro de Biologia Molecular "Severo Ochoa", CSIC-UAM Cantoblanco , 28049 Madrid , Spain
| | - Yves Dehouck
- Centro de Biologia Molecular "Severo Ochoa", CSIC-UAM Cantoblanco , 28049 Madrid , Spain
| |
Collapse
|
22
|
Tsuchiya Y, Taneishi K, Yonezawa Y. Autoencoder-Based Detection of Dynamic Allostery Triggered by Ligand Binding Based on Molecular Dynamics. J Chem Inf Model 2019; 59:4043-4051. [PMID: 31386362 DOI: 10.1021/acs.jcim.9b00426] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Dynamic allostery on proteins, triggered by regulator binding or chemical modifications, transmits information from the binding site to distant regions, dramatically altering protein function. It is accompanied by subtle changes in side-chain conformations of the protein, indicating that the changes in dynamics, and not rigid or large conformational changes, are essential to understand regulation of protein function. Although a lot of experimental and theoretical studies have been dedicated to investigate this issue, the regulation mechanism of protein function is still being debated. Here, we propose an autoencoder-based method that can detect dynamic allostery. The method is based on the comparison of time fluctuations of protein structures, in the form of distance matrices, obtained from molecular dynamics simulations in ligand-bound and -unbound forms. Our method detected that the changes in dynamics by ligand binding in the PDZ2 domain led to the reorganization of correlative fluctuation motions among residue pairs, which revealed a different view of the correlated motions from the PCA and DCCM. In addition, other correlative motions were also found as a result of the dynamic perturbation from the ligand binding, which may lead to dynamic allostery. This autoencoder-based method would be usefully applied to the signal transduction and mutagenesis systems involved in protein functions and severe diseases.
Collapse
Affiliation(s)
- Yuko Tsuchiya
- Artificial Intelligence Research Center , National Institute of Advanced Industrial Science and Technology , 2-4-7 Aomi , Koto-ku , Tokyo 135-0064 , Japan
| | - Kei Taneishi
- Cluster for Science, Technology and Innovation Hub , RIKEN , 6-7-3 Minatojima-minamimachi , Chuo-ku, Kobe , Hyogo 650-0047 , Japan
| | - Yasushige Yonezawa
- High Pressure Protein Research Center, Institute of Advanced Technology , Kindai University , 930 Nishimitani , Kinokawa , Wakayama 649-6493 , Japan
| |
Collapse
|
23
|
Munshi S, Gopi S, Subramanian S, Campos LA, Naganathan AN. Protein plasticity driven by disorder and collapse governs the heterogeneous binding of CytR to DNA. Nucleic Acids Res 2019. [PMID: 29538715 PMCID: PMC5934615 DOI: 10.1093/nar/gky176] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
The amplitude of thermodynamic fluctuations in biological macromolecules determines their conformational behavior, dimensions, nature of phase transitions and effectively their specificity and affinity, thus contributing to fine-tuned molecular recognition. Unique among large-scale conformational changes in proteins are temperature-induced collapse transitions in intrinsically disordered proteins (IDPs). Here, we show that CytR DNA-binding domain, an IDP that folds on binding DNA, undergoes a coil-to-globule transition with temperature in the absence of DNA while exhibiting energetically decoupled local and global structural rearrangements, and maximal thermodynamic fluctuations at the optimal bacterial growth temperature. The collapse is shown to be a continuous transition through a combination of statistical-mechanical modeling and all-atom implicit solvent simulations. Surprisingly, CytR binds single-site cognate DNA with negative cooperativity, described by Hill coefficients less than one, resulting in a graded binding response. We show that heterogeneity arising from varying binding-competent CytR conformations or orientations at the single-molecular level contributes to negative binding cooperativity at the level of bulk measurements due to the conflicting requirements of collapse transition, large fluctuations and folding-upon-binding. Our work reports strong evidence for functionally driven thermodynamic fluctuations in determining the extent of collapse and disorder with implications in protein search efficiency of target DNA sites and regulation.
Collapse
Affiliation(s)
- Sneha Munshi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sandhyaa Subramanian
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Luis A Campos
- National Biotechnology Center, Consejo Superior de Investigaciones Científicas, Darwin 3, Campus de Cantoblanco, 28049 Madrid, Spain
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
24
|
Bolik-Coulon N, Cousin SF, Kadeřávek P, Dumez JN, Ferrage F. Understanding the methyl-TROSY effect over a wide range of magnetic fields. J Chem Phys 2019; 150:224202. [DOI: 10.1063/1.5095757] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nicolas Bolik-Coulon
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Samuel F. Cousin
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Pavel Kadeřávek
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | | | - Fabien Ferrage
- Laboratoire des Biomolécules, LBM, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
25
|
Mishra P, Jha SK. Slow Motion Protein Dance Visualized Using Red-Edge Excitation Shift of a Buried Fluorophore. J Phys Chem B 2019; 123:1256-1264. [PMID: 30640479 DOI: 10.1021/acs.jpcb.8b11151] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
It has been extremely challenging to detect protein structures with a dynamic core, such as dry molten globules, that remain in equilibrium with the tightly packed native (N) state and that are important for a myriad of entropy-driven protein functions. Here, we detect the higher entropy conformations of a human serum protein, using red-edge excitation shift experiments. We covalently introduced a fluorophore inside the protein core and observed that in a subset of native population, the side chains of the polar and buried residues have different spatial arrangements than the mean population and that they solvate the fluorophore on a timescale much slower than the nanosecond timescale of fluorescence. Our results provide direct evidence for the dense fluidity of protein core and show that alternate side-chain packing arrangements exist in the core that might be important for multiple binding functions of this protein.
Collapse
Affiliation(s)
- Prajna Mishra
- Physical and Materials Chemistry Division, Academy of Scientific and Innovative Research (AcSIR) , CSIR-National Chemical Laboratory , Dr. Homi Bhabha Road , Pune 411008 , Maharashtra , India
| | - Santosh Kumar Jha
- Physical and Materials Chemistry Division, Academy of Scientific and Innovative Research (AcSIR) , CSIR-National Chemical Laboratory , Dr. Homi Bhabha Road , Pune 411008 , Maharashtra , India
| |
Collapse
|
26
|
Modulation of allosteric coupling by mutations: from protein dynamics and packing to altered native ensembles and function. Curr Opin Struct Biol 2018; 54:1-9. [PMID: 30268910 PMCID: PMC6420056 DOI: 10.1016/j.sbi.2018.09.004] [Citation(s) in RCA: 64] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2018] [Revised: 08/13/2018] [Accepted: 09/10/2018] [Indexed: 01/12/2023]
Abstract
A large body of work has gone into understanding the effect of mutations on protein structure and function. Conventional treatments have involved quantifying the change in stability, activity and relaxation rates of the mutants with respect to the wild-type protein. However, it is now becoming increasingly apparent that mutational perturbations consistently modulate the packing and dynamics of a significant fraction of protein residues, even those that are located >10–15 Å from the mutated site. Such long-range modulation of protein features can distinctly tune protein stability and the native conformational ensemble contributing to allosteric modulation of function. In this review, I summarize a series of experimental and computational observations that highlight the incredibly pliable nature of proteins and their response to mutational perturbations manifested via the intra-protein interaction network. I highlight how an intimate understanding of mutational effects could pave the way for integrating stability, folding, cooperativity and even allostery within a single physical framework.
Collapse
|
27
|
Cousin SF, Kadeřávek P, Bolik-Coulon N, Gu Y, Charlier C, Carlier L, Bruschweiler-Li L, Marquardsen T, Tyburn JM, Brüschweiler R, Ferrage F. Time-Resolved Protein Side-Chain Motions Unraveled by High-Resolution Relaxometry and Molecular Dynamics Simulations. J Am Chem Soc 2018; 140:13456-13465. [DOI: 10.1021/jacs.8b09107] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Samuel F. Cousin
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Pavel Kadeřávek
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Nicolas Bolik-Coulon
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Yina Gu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Cyril Charlier
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Ludovic Carlier
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Lei Bruschweiler-Li
- Campus Chemical Instrument Center, The Ohio State University, Columbus, Ohio 43210, United States
| | | | - Jean-Max Tyburn
- Bruker BioSpin, 34 rue de l’Industrie BP 10002, 67166 Wissembourg Cedex, France
| | - Rafael Brüschweiler
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Campus Chemical Instrument Center, The Ohio State University, Columbus, Ohio 43210, United States
- Department of Biological Chemistry and Pharmacology, The Ohio State University, Columbus, Ohio 43210, United States
| | - Fabien Ferrage
- Laboratoire des biomolécules, LBM, Département de chimie, École normale supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
28
|
Botlani M, Siddiqui A, Varma S. Machine learning approaches to evaluate correlation patterns in allosteric signaling: A case study of the PDZ2 domain. J Chem Phys 2018; 148:241726. [DOI: 10.1063/1.5022469] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Mohsen Botlani
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, Florida 33620, USA
| | - Ahnaf Siddiqui
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, Florida 33620, USA
| | - Sameer Varma
- Department of Cell Biology, Microbiology and Molecular Biology, University of South Florida, Tampa, Florida 33620, USA
| |
Collapse
|
29
|
Vajda S, Beglov D, Wakefield AE, Egbert M, Whitty A. Cryptic binding sites on proteins: definition, detection, and druggability. Curr Opin Chem Biol 2018; 44:1-8. [PMID: 29800865 PMCID: PMC6088748 DOI: 10.1016/j.cbpa.2018.05.003] [Citation(s) in RCA: 93] [Impact Index Per Article: 15.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2018] [Accepted: 05/03/2018] [Indexed: 12/30/2022]
Abstract
Many proteins in their unbound structures lack surface pockets appropriately sized for drug binding. Hence, a variety of experimental and computational tools have been developed for the identification of cryptic sites that are not evident in the unbound protein but form upon ligand binding, and can provide tractable drug target sites. The goal of this review is to discuss the definition, detection, and druggability of such sites, and their potential value for drug discovery. Novel methods based on molecular dynamics simulations are particularly promising and yield a large number of transient pockets, but it has been shown that only a minority of such sites are generally capable of binding ligands with substantial affinity. Based on recent studies, current methodology can be improved by combining molecular dynamics with fragment docking and machine learning approaches.
Collapse
Affiliation(s)
- Sandor Vajda
- Department of Biomedical Engineering, Boston University, Boston, MA 02215, United States; Department of Chemistry, Boston University, Boston, MA 02215, United States.
| | - Dmitri Beglov
- Department of Biomedical Engineering, Boston University, Boston, MA 02215, United States
| | - Amanda E Wakefield
- Department of Biomedical Engineering, Boston University, Boston, MA 02215, United States; Department of Chemistry, Boston University, Boston, MA 02215, United States
| | - Megan Egbert
- Department of Biomedical Engineering, Boston University, Boston, MA 02215, United States
| | - Adrian Whitty
- Department of Chemistry, Boston University, Boston, MA 02215, United States.
| |
Collapse
|
30
|
Abstract
Molecular recognition by proteins is fundamental to the molecular basis of biology. Dissection of the thermodynamic landscape governing protein-ligand interactions has proven difficult because determination of various entropic contributions is quite challenging. Nuclear magnetic resonance relaxation measurements, theory, and simulations suggest that conformational entropy can be accessed through a dynamical proxy. Here, we review the relationship between measures of fast side-chain motion and the underlying conformational entropy. The dynamical proxy reveals that the contribution of conformational entropy can range from highly favorable to highly unfavorable and demonstrates the potential of this key thermodynamic variable to modulate protein-ligand interactions. The dynamical so-called entropy meter also refines the role of solvent entropy and directly determines the loss in rotational-translational entropy that occurs upon formation of high-affinity complexes. The ability to quantify the roles of entropy through an entropy meter based on measurable dynamical properties promises to highlight its role in protein function.
Collapse
Affiliation(s)
- A Joshua Wand
- Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania 19104-6059, USA; ,
| | - Kim A Sharp
- Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania Perelman School of Medicine, Philadelphia, Pennsylvania 19104-6059, USA; ,
| |
Collapse
|
31
|
Andryushchenko VA, Chekmarev SF. Temperature evolution of Trp-cage folding pathways: An analysis by dividing the probability flux field into stream tubes. J Biol Phys 2017; 43:565-583. [PMID: 28983809 DOI: 10.1007/s10867-017-9470-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2016] [Accepted: 09/01/2017] [Indexed: 11/25/2022] Open
Abstract
Owing to its small size and very fast folding rate, the Trp-cage miniprotein has become a benchmark system to study protein folding. Two folding pathways were found to be characteristic of this protein: pathway I, in which the hydrophobic collapse precedes the formation of α-helix, and pathway II, in which the events occur in the reverse order. At the same time, the relative contribution of these pathways at different temperatures as well as the nature of transition from one pathway to the other remain unclear. To gain insight into this issue, we employ a recently proposed hydrodynamic description of protein folding, in which the process of folding is considered as a motion of a "folding fluid" (Chekmarev et al., Phys. Rev. Lett. 100(1), 018107 2008). Using molecular dynamics simulations, we determine the field of probability fluxes of transitions in a space of collective variables and divide it into stream tubes. Each tube contains a definite fraction of the total folding flow and can be associated with a certain pathway. Specifically, three temperatures were considered, T = 285K, T = 315K, and T = 325K. We have found that as the temperature increases, the contribution of pathway I, which is approximately 90% of the total folding flow at T = 285K, decreases to approximately 10% at T = 325K, i.e., pathway II becomes dominant. At T = 315K, both pathways contribute approximately equally. All these temperatures are found below the calculated melting point, which suggests that the Trp-cage folding mechanism is determined by kinetic factors rather than thermodynamics.
Collapse
Affiliation(s)
- Vladimir A Andryushchenko
- Institute of Thermophysics, SB RAS, 630090, Novosibirsk, Russia
- Department of Physics, Novosibirsk State University, 630090, Novosibirsk, Russia
| | - Sergei F Chekmarev
- Institute of Thermophysics, SB RAS, 630090, Novosibirsk, Russia.
- Department of Physics, Novosibirsk State University, 630090, Novosibirsk, Russia.
| |
Collapse
|
32
|
Rajasekaran N, Sekhar A, Naganathan AN. A Universal Pattern in the Percolation and Dissipation of Protein Structural Perturbations. J Phys Chem Lett 2017; 8:4779-4784. [PMID: 28910120 DOI: 10.1021/acs.jpclett.7b02021] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
Understanding the extent to which information is transmitted through the intramolecular interaction network of proteins upon a perturbation, that is, an allosteric effect, has long remained an unsolved problem. Through an analysis of high-resolution NMR data from the literature on 28 different proteins and 49 structural perturbations, we show that the extent of induced structural changes through mutations and molecular events including protein-protein, protein-peptide, protein-ligand binding, and post-translational modifications exhibit a near-universal exponential functional form. The extent of percolation into the protein structures can be up to 20-25 Å despite no apparent change in the 3D structures. These observations are also consistent with theoretical expectations, elementary graph theoretic analysis of protein structures, detailed molecular dynamics simulations, and experimental double-mutant cycles. Our analysis highlights that most molecular events would contribute to allosteric effects independent of protein structure, topology, or identity and provides a simple avenue to test and potentially model their effects.
Collapse
Affiliation(s)
- Nandakumar Rajasekaran
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras , Chennai 600036, India
| | - Ashok Sekhar
- Departments of Molecular Genetics, Biochemistry, and Chemistry, The University of Toronto , Toronto, Ontario M5S 1A8, Canada
| | - Athi N Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras , Chennai 600036, India
| |
Collapse
|
33
|
Mishra P, Jha SK. An Alternatively Packed Dry Molten Globule-like Intermediate in the Native State Ensemble of a Multidomain Protein. J Phys Chem B 2017; 121:9336-9347. [DOI: 10.1021/acs.jpcb.7b07032] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Prajna Mishra
- Physical and Materials Chemistry
Division, CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune, Maharashtra 411008, India
| | - Santosh Kumar Jha
- Physical and Materials Chemistry
Division, CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune, Maharashtra 411008, India
| |
Collapse
|
34
|
Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea. J Mol Biol 2017; 430:450-464. [PMID: 28782558 DOI: 10.1016/j.jmb.2017.07.022] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Revised: 07/29/2017] [Accepted: 07/31/2017] [Indexed: 11/22/2022]
Abstract
While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the Tm and Cm. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.
Collapse
|
35
|
A self-consistent structural perturbation approach for determining the magnitude and extent of allosteric coupling in proteins. Biochem J 2017; 474:2379-2388. [DOI: 10.1042/bcj20170304] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2017] [Revised: 05/12/2017] [Accepted: 05/17/2017] [Indexed: 11/17/2022]
Abstract
Elucidating the extent of energetic coupling between residues in single-domain proteins, which is a fundamental determinant of allostery, information transfer and folding cooperativity, has remained a grand challenge. While several sequence- and structure-based approaches have been proposed, a self-consistent description that is simultaneously compatible with unfolding thermodynamics is lacking. We recently developed a simple structural perturbation protocol that captures the changes in thermodynamic stabilities induced by point mutations within the protein interior. Here, we show that a fundamental residue-specific component of this perturbation approach, the coupling distance, is uniquely sensitive to the environment of a residue in the protein to a distance of ∼15 Å. With just the protein contact map as an input, we reproduce the extent of percolation of perturbations within the structure as observed in network analysis of intra-protein interactions, molecular dynamics simulations and NMR-observed changes in chemical shifts. Using this rapid protocol that relies on a single structure, we explain the results of statistical coupling analysis (SCA) that requires hundreds of sequences to identify functionally critical sectors, the propagation and dissipation of perturbations within proteins and the higher-order couplings deduced from detailed NMR experiments. Our results thus shed light on the possible mechanistic origins of signaling through the interaction network within proteins, the likely distance dependence of perturbations induced by ligands and post-translational modifications and the origins of folding cooperativity through many-body interactions.
Collapse
|
36
|
Rajasekaran N, Suresh S, Gopi S, Raman K, Naganathan AN. A General Mechanism for the Propagation of Mutational Effects in Proteins. Biochemistry 2016; 56:294-305. [DOI: 10.1021/acs.biochem.6b00798] [Citation(s) in RCA: 42] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Nandakumar Rajasekaran
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | | | - Soundhararajan Gopi
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Karthik Raman
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Athi N. Naganathan
- Department of Biotechnology, Bhupat & Jyoti Mehta School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
37
|
Non-covalent forces tune the electron transfer complex between ferredoxin and sulfite reductase to optimize enzymatic activity. Biochem J 2016; 473:3837-3854. [PMID: 27551107 DOI: 10.1042/bcj20160658] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2016] [Accepted: 08/22/2016] [Indexed: 11/17/2022]
Abstract
Although electrostatic interactions between negatively charged ferredoxin (Fd) and positively charged sulfite reductase (SiR) have been predominantly highlighted to characterize complex formation, the detailed nature of intermolecular forces remains to be fully elucidated. We investigated interprotein forces for the formation of an electron transfer complex between Fd and SiR and their relationship to SiR activity using various approaches over NaCl concentrations between 0 and 400 mM. Fd-dependent SiR activity assays revealed a bell-shaped activity curve with a maximum ∼40-70 mM NaCl and a reverse bell-shaped dependence of interprotein affinity. Meanwhile, intrinsic SiR activity, as measured in a methyl viologen-dependent assay, exhibited saturation above 100 mM NaCl. Thus, two assays suggested that interprotein interaction is crucial in controlling Fd-dependent SiR activity. Calorimetric analyses showed the monotonic decrease in interprotein affinity on increasing NaCl concentrations, distinguished from a reverse bell-shaped interprotein affinity observed from Fd-dependent SiR activity assay. Furthermore, Fd:SiR complex formation and interprotein affinity were thermodynamically adjusted by both enthalpy and entropy through electrostatic and non-electrostatic interactions. A residue-based NMR investigation on the addition of SiR to 15N-labeled Fd at the various NaCl concentrations also demonstrated that a combination of electrostatic and non-electrostatic forces stabilized the complex with similar interfaces and modulated the binding affinity and mode. Our findings elucidate that non-electrostatic forces are also essential for the formation and modulation of the Fd:SiR complex. We suggest that a complex configuration optimized for maximum enzymatic activity near physiological salt conditions is achieved by structural rearrangement through controlled non-covalent interprotein interactions.
Collapse
|
38
|
Bondos SE, Swint-Kruse L, Matthews KS. Flexibility and Disorder in Gene Regulation: LacI/GalR and Hox Proteins. J Biol Chem 2015; 290:24669-77. [PMID: 26342073 DOI: 10.1074/jbc.r115.685032] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
To modulate transcription, a variety of input signals must be sensed by genetic regulatory proteins. In these proteins, flexibility and disorder are emerging as common themes. Prokaryotic regulators generally have short, flexible segments, whereas eukaryotic regulators have extended regions that lack predicted secondary structure (intrinsic disorder). Two examples illustrate the impact of flexibility and disorder on gene regulation: the prokaryotic LacI/GalR family, with detailed information from studies on LacI, and the eukaryotic family of Hox proteins, with specific insights from investigations of Ultrabithorax (Ubx). The widespread importance of structural disorder in gene regulatory proteins may derive from the need for flexibility in signal response and, particularly in eukaryotes, in protein partner selection.
Collapse
Affiliation(s)
- Sarah E Bondos
- From the Department of Molecular and Cellular Medicine, Texas A&M Health Science Center, College Station, Texas 77843
| | - Liskin Swint-Kruse
- the Department of Biochemistry and Molecular Biology, the University of Kansas Medical Center, Kansas City, Kansas 66160, and
| | | |
Collapse
|
39
|
Bowman GR. Accurately modeling nanosecond protein dynamics requires at least microseconds of simulation. J Comput Chem 2015; 37:558-66. [PMID: 26077712 DOI: 10.1002/jcc.23973] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2015] [Revised: 05/03/2015] [Accepted: 05/24/2015] [Indexed: 01/17/2023]
Abstract
Advances in hardware and algorithms have greatly extended the timescales accessible to molecular simulation. This article assesses whether such long timescale simulations improve our ability to calculate order parameters that describe conformational heterogeneity on ps-ns timescales or if force fields are now a limiting factor. Order parameters from experiment are compared with order parameters calculated in three different ways from simulations ranging from 10 ns to 100 μs in length. Importantly, bootstrapping is employed to assess the variability in results for each simulation length. The results of 10-100 ns timescale simulations are highly variable, possibly explaining the variation in levels of agreement between simulation and experiment in published works examining different proteins. Fortunately, microsecond timescale simulations improve both the accuracy and precision of calculated order parameters, at least so long as the full exponential fit or truncated average approximation is used instead of the common long-time limit approximation. The improved precision of these long timescale simulations allows a statistically sound comparison of a number of modern force fields, such as Amber03, Amber99sb-ILDN, and Charmm27. While there is some variation between these force fields, they generally give very similar results for sufficiently long simulations. The fact that so much simulation is required to precisely capture ps-ns timescale processes suggests that extremely extensive simulations are required for slower processes. Advanced sampling techniques could aid greatly, however, such methods will need to maintain accurate kinetics if they are to be of value for calculating dynamical properties like order parameters. © 2015 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Gregory R Bowman
- Department of Biochemistry & Molecular Biophysics, Department of Biomedical Engineering and Center for Biological Systems Engineering, Washington University School of Medicine, St. Louis, Missouri, 63110
| |
Collapse
|