1
|
Huang B, Yin T, Fu S, Liu L, Yang C, Zhou L, Liu X, Zhuang H, Cao Z, Hua Z. Inflammation-oriented montmorillonite adjuvant enhanced oral delivery of anti-TNF-α nanobody against inflammatory bowel disease. Proc Natl Acad Sci U S A 2024; 121:e2320482121. [PMID: 39226349 DOI: 10.1073/pnas.2320482121] [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/03/2023] [Accepted: 07/02/2024] [Indexed: 09/05/2024] Open
Abstract
Oral delivery of proteins faces challenges due to the harsh conditions of the gastrointestinal (GI) tract, including gastric acid and intestinal enzyme degradation. Permeation enhancers are limited in their ability to deliver proteins with high molecular weight and can potentially cause toxicity by opening tight junctions. To overcome these challenges, we propose the use of montmorillonite (MMT) as an adjuvant that possesses both inflammation-oriented abilities and the ability to regulate gut microbiota. This adjuvant can be used as a universal protein oral delivery technology by fusing with advantageous binding amino acid sequences. We demonstrated that anti-TNF-α nanobody (VII) can be intercalated into the MMT interlayer space. The carboxylate groups (-COOH) of aspartic acid (D) and glutamic acid (E) interact with the MMT surface through electrostatic interactions with sodium ions (Na+). The amino groups (NH2) of asparagine (N) and glutamine (Q) are primarily attracted to the MMT layers through hydrogen bonding with oxygen atoms on the surface. This binding mechanism protects VII from degradation and ensures its release in the intestinal tract, as well as retaining biological activity, leading to significantly enhanced therapeutic effects on colitis. Furthermore, VII@MMT increases the abundance of short-chain fatty acids (SCFAs)-producing strains, including Clostridia, Prevotellaceae, Alloprevotella, Oscillospiraceae, Clostridia_vadinBB60_group, and Ruminococcaceae, therefore enhance the production of SCFAs and butyrate, inducing regulatory T cells (Tregs) production to modulate local and systemic immune homeostasis. Overall, the MMT adjuvant provides a promising universal strategy for protein oral delivery by rational designed protein.
Collapse
Affiliation(s)
- Baolian Huang
- School of Biopharmacy, China Pharmaceutical University, Nanjing 211198, People's Republic of China
| | - Te Yin
- The State Key Laboratory of Pharmaceutical Biotechnology, College of Life Sciences, Nanjing University, Nanjing, People's Republic of China
| | - Shuilian Fu
- The State Key Laboratory of Pharmaceutical Biotechnology, College of Life Sciences, Nanjing University, Nanjing, People's Republic of China
| | - Lina Liu
- The State Key Laboratory of Pharmaceutical Biotechnology, College of Life Sciences, Nanjing University, Nanjing, People's Republic of China
| | - Chen Yang
- School of Biopharmacy, China Pharmaceutical University, Nanjing 211198, People's Republic of China
| | - Lulu Zhou
- School of Biopharmacy, China Pharmaceutical University, Nanjing 211198, People's Republic of China
| | - Xing Liu
- The State Key Laboratory of Pharmaceutical Biotechnology, College of Life Sciences, Nanjing University, Nanjing, People's Republic of China
| | - Hongqin Zhuang
- The State Key Laboratory of Pharmaceutical Biotechnology, College of Life Sciences, Nanjing University, Nanjing, People's Republic of China
| | - Zhiting Cao
- School of Biopharmacy, China Pharmaceutical University, Nanjing 211198, People's Republic of China
| | - Zichun Hua
- School of Biopharmacy, China Pharmaceutical University, Nanjing 211198, People's Republic of China
- The State Key Laboratory of Pharmaceutical Biotechnology, College of Life Sciences, Nanjing University, Nanjing, People's Republic of China
- Changzhou High-Tech Research Institute of Nanjing University and Jiangsu TargetPharma Laboratories Inc., Changzhou 213164, People's Republic of China
| |
Collapse
|
2
|
Roth M, Toker Y, Major DT. Monte Carlo-Simulated Annealing and Machine Learning-Based Funneled Approach for Finding the Global Minimum Structure of Molecular Clusters. ACS OMEGA 2024; 9:1298-1309. [PMID: 38222530 PMCID: PMC10785639 DOI: 10.1021/acsomega.3c07600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/01/2023] [Revised: 11/09/2023] [Accepted: 11/15/2023] [Indexed: 01/16/2024]
Abstract
Understanding the physical underpinnings and geometry of molecular clusters is of great importance in many fields, ranging from studying the beginning of the universe to the formation of atmospheric particles. To this end, several approaches have been suggested, yet identifying the most stable cluster geometry (i.e., global potential energy minimum) remains a challenge, especially for highly symmetric clusters. Here, we suggest a new funneled Monte Carlo-based simulated annealing (SA) approach, which includes two key steps: generation of symmetrical clusters and classification of the clusters according to their geometry using machine learning (MCSA-ML). We demonstrate the merits of the MCSA-ML method in comparison to other approaches on several Lennard-Jones (LJ) clusters and four molecular clusters-Ser8(Cl-)2, H+(H2O)6, Ag+(CO2)8, and Bet4Cl-. For the latter of these clusters, the correct structure is unknown, and hence, we compare the experimental and simulated fragmentation patterns, and the fragmentation of the proposed global minimum matches experiments closely. Additionally, based on the fragmentation of the predicted betaine cluster, we were able to identify hitherto unknown neutral fragmentation channels. In comparison to results obtained with other methods, we demonstrated a superior ability of MCSA-ML to predict clusters with high symmetry and similar abilities to predict clusters with asymmetrical structures.
Collapse
Affiliation(s)
- Michal Roth
- Department
of Physics, Bar-Ilan University, Ramat-Gan 5290002, Israel
- Institute
of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Yoni Toker
- Department
of Physics, Bar-Ilan University, Ramat-Gan 5290002, Israel
- Institute
of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Dan T. Major
- Institute
of Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
- Department
of Chemistry, Bar-Ilan University, Ramat-Gan 5290002, Israel
| |
Collapse
|
3
|
Hsieh YC, Delarue M, Orland H, Koehl P. Analyzing the Geometry and Dynamics of Viral Structures: A Review of Computational Approaches Based on Alpha Shape Theory, Normal Mode Analysis, and Poisson-Boltzmann Theories. Viruses 2023; 15:1366. [PMID: 37376665 DOI: 10.3390/v15061366] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Revised: 06/05/2023] [Accepted: 06/09/2023] [Indexed: 06/29/2023] Open
Abstract
The current SARS-CoV-2 pandemic highlights our fragility when we are exposed to emergent viruses either directly or through zoonotic diseases. Fortunately, our knowledge of the biology of those viruses is improving. In particular, we have more and more structural information on virions, i.e., the infective form of a virus that includes its genomic material and surrounding protective capsid, and on their gene products. It is important to have methods that enable the analyses of structural information on such large macromolecular systems. We review some of those methods in this paper. We focus on understanding the geometry of virions and viral structural proteins, their dynamics, and their energetics, with the ambition that this understanding can help design antiviral agents. We discuss those methods in light of the specificities of those structures, mainly that they are huge. We focus on three of our own methods based on the alpha shape theory for computing geometry, normal mode analyses to study dynamics, and modified Poisson-Boltzmann theories to study the organization of ions and co-solvent and solvent molecules around biomacromolecules. The corresponding software has computing times that are compatible with the use of regular desktop computers. We show examples of their applications on some outer shells and structural proteins of the West Nile Virus.
Collapse
Affiliation(s)
- Yin-Chen Hsieh
- Institute for Arctic and Marine Biology, Department of Biosciences, Fisheries, and Economics, UiT The Arctic University of Norway, 9037 Tromso, Norway
| | - Marc Delarue
- Institut Pasteur, Université Paris-Cité and CNRS, UMR 3528, Unité Architecture et Dynamique des Macromolécules Biologiques, 75015 Paris, France
| | - Henri Orland
- Institut de Physique Théorique, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Patrice Koehl
- Department of Computer Science, University of California, Davis, CA 95616, USA
| |
Collapse
|
4
|
Koehl P, Akopyan A, Edelsbrunner H. Computing the Volume, Surface Area, Mean, and Gaussian Curvatures of Molecules and Their Derivatives. J Chem Inf Model 2023; 63:973-985. [PMID: 36638318 PMCID: PMC9930125 DOI: 10.1021/acs.jcim.2c01346] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Geometry is crucial in our efforts to comprehend the structures and dynamics of biomolecules. For example, volume, surface area, and integrated mean and Gaussian curvature of the union of balls representing a molecule are used to quantify its interactions with the water surrounding it in the morphometric implicit solvent models. The Alpha Shape theory provides an accurate and reliable method for computing these geometric measures. In this paper, we derive homogeneous formulas for the expressions of these measures and their derivatives with respect to the atomic coordinates, and we provide algorithms that implement them into a new software package, AlphaMol. The only variables in these formulas are the interatomic distances, making them insensitive to translations and rotations. AlphaMol includes a sequential algorithm and a parallel algorithm. In the parallel version, we partition the atoms of the molecule of interest into 3D rectangular blocks, using a kd-tree algorithm. We then apply the sequential algorithm of AlphaMol to each block, augmented by a buffer zone to account for atoms whose ball representations may partially cover the block. The current parallel version of AlphaMol leads to a 20-fold speed-up compared to an independent serial implementation when using 32 processors. For instance, it takes 31 s to compute the geometric measures and derivatives of each atom in a viral capsid with more than 26 million atoms on 32 Intel processors running at 2.7 GHz. The presence of the buffer zones, however, leads to redundant computations, which ultimately limit the impact of using multiple processors. AlphaMol is available as an OpenSource software.
Collapse
Affiliation(s)
- Patrice Koehl
- Department
of Computer Science, University of California, Davis, California95616, United States,
| | | | | |
Collapse
|
5
|
Spiropyran/Merocyanine Amphiphile in Various Solvents: A Joint Experimental–Theoretical Approach to Photophysical Properties and Self-Assembly. Int J Mol Sci 2022; 23:ijms231911535. [PMID: 36232836 PMCID: PMC9569490 DOI: 10.3390/ijms231911535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Revised: 09/21/2022] [Accepted: 09/24/2022] [Indexed: 11/20/2022] Open
Abstract
This joint experimental-theoretical work focuses on molecular and photophysical properties of the spiropyran-containing amphiphilic molecule in organic and aqueous solutions. Being dissolved in tested organic solvents, the system demonstrates positive photochromism, i.e., upon UV stimulus the colorless spiropyran form is transformed into colorful merocyanine isomer. However, the aqueous solution of the amphiphile possesses a negative photochromism: the orange-red merocyanine form becomes thermodynamically more stable in water, and both UV and vis stimuli lead to the partial or complete photobleaching of the solution. The explanation of this phenomenon is given on the basis of density functional theory calculations and classical modeling including thermodynamic integration. The simulations reveal that stabilization of merocyanine in water proceeds with the energy of ca. 70 kJ mol−1, and that the Helmholtz free energy of hydration of merocyanine form is 100 kJ mol−1 lower as compared to the behavior of SP isomer in water. The explanation of such a difference lies in the molecular properties of the merocyanine: after ring-opening reaction this molecule transforms into a zwitterionic form, as evidenced by the electrostatic potential plotted around the opened form. The presence of three charged groups on the periphery of a flat conjugated backbone stimulates the self-assembly of merocyanine molecules in water, ending up with the formation of elongated associates with stack-like building blocks, as shown in molecular dynamics simulations of the aqueous solution with the concentration above critical micelle concentration. Our quantitative evaluation of the hydrophilicity switching in spiropyran/merocyanine containing surfactants may prompt the search for new systems, including colloidal and polymeric ones, aiming at remote tuning of their morphology, which could give new promising shapes and patterns for the needs of modern nanotechnology.
Collapse
|
6
|
Maity H, Baidya L, Reddy G. Salt-Induced Transitions in the Conformational Ensembles of Intrinsically Disordered Proteins. J Phys Chem B 2022; 126:5959-5971. [PMID: 35944496 DOI: 10.1021/acs.jpcb.2c03476] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Salts modulate the behavior of intrinsically disordered proteins (IDPs) and influence the formation of membraneless organelles through liquid-liquid phase separation (LLPS). In low ionic strength solutions, IDP conformations are perturbed by the screening of electrostatic interactions, independent of the salt identity. In this regime, insight into the IDP behavior can be obtained using the theory for salt-induced transitions in charged polymers. However, salt-specific interactions with the charged and uncharged residues, known as the Hofmeister effect, influence IDP behavior in high ionic strength solutions. There is a lack of reliable theoretical models in high salt concentration regimes to predict the salt effect on IDPs. We propose a simulation methodology using a coarse-grained IDP model and experimentally measured water to salt solution transfer free energies of various chemical groups that allowed us to study the salt-specific transitions induced in the IDPs conformational ensemble. We probed the effect of three different monovalent salts on five IDPs belonging to various polymer classes based on charged residue content. We demonstrate that all of the IDPs of different polymer classes behave as self-avoiding walks (SAWs) at physiological salt concentration. In high salt concentrations, the transitions observed in the IDP conformational ensembles are dependent on the salt used and the IDP sequence and composition. Changing the anion with the cation fixed can result in the IDP transition from a SAW-like behavior to a collapsed globule. An important implication of these results is that a suitable salt can be identified to induce condensation of an IDP through LLPS.
Collapse
Affiliation(s)
- Hiranmay Maity
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka, India 560012
| | - Lipika Baidya
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka, India 560012
| | - Govardhan Reddy
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bengaluru, Karnataka, India 560012
| |
Collapse
|
7
|
Nicy, Chakraborty D, Wales DJ. Energy Landscapes for Base-Flipping in a Model DNA Duplex. J Phys Chem B 2022; 126:3012-3028. [PMID: 35427136 PMCID: PMC9098180 DOI: 10.1021/acs.jpcb.2c00340] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2022] [Revised: 03/24/2022] [Indexed: 12/31/2022]
Abstract
We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5' side of the backbone. We find that bases align toward the 3' side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5' side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.
Collapse
Affiliation(s)
- Nicy
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge, CB2 1EW, U.K.
| | - Debayan Chakraborty
- Department
of Chemistry, The University of Texas at
Austin, Austin, Texas 78712, United States
| | - David J. Wales
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield Road, Cambridge, CB2 1EW, U.K.
| |
Collapse
|
8
|
Li W, Nforneh B, Whitcomb KL, Warncke K. Resolution and characterization of confinement- and temperature-dependent dynamics in solvent phases that surround proteins in frozen aqueous solution by using spin-probe EPR spectroscopy. Methods Enzymol 2022; 666:25-57. [PMID: 35465922 DOI: 10.1016/bs.mie.2022.02.009] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
Spin probe electron paramagnetic resonance spectroscopy is applied to characterize the dynamics of concentric hydration and mesophase solvent domains that surround proteins within the ice boundary in frozen aqueous solutions. The solvent dynamics are tuned by variation of temperature (190-265K) and by the degree of ice boundary confinement, which is modulated by the volume of added cryosolvent (0-~50Å separation distance from protein surface). Goals are to: (1) characterize the protein-coupled solvent dynamics on correlation time scales of ~10-10<τ<10-7s, and spatial scales from protein surface to periphery of the surrounding solution, from the perspective of a free, small-molecule (~7Å diameter) probe, and (2) reveal properties of the solvent-protein coupling that can be correlated with protein functions, that are measureable under the same conditions. Rotational mobility of the nitroxide spin probe, TEMPOL, resolves and tracks two solvent components, the protein-associated domain (PAD; akin to hydration layer) and surrounding mesodomain, through their distinct temperature- and confinement-dependent values of τ and normalized weight. Detailed protocols are described for simulation of two-component nitroxide EPR spectra, which are categorized by line shape regime and guided by a library of template spectra and simulation parameters derived from two model soluble globular proteins. The order-disorder transition in the PAD, which is a universal feature of protein-coupled solvent dynamics, provides a well-defined, tunable property for elucidating mechanism in solvent-protein-function dynamical coupling. The low-temperature mesodomain system and EPR spin probe method are generally applicable to reveal solvent contributions to a broad range of macromolecule-mediated biological processes.
Collapse
Affiliation(s)
- Wei Li
- Department of Physics, Emory University, Atlanta, GA, United States
| | - Benjamen Nforneh
- Department of Physics, Emory University, Atlanta, GA, United States
| | - Katie L Whitcomb
- Department of Physics, Emory University, Atlanta, GA, United States
| | - Kurt Warncke
- Department of Physics, Emory University, Atlanta, GA, United States.
| |
Collapse
|
9
|
Molecular dynamics simulation or structure refinement of proteins: are solvent molecules required? A case study using hen lysozyme. EUROPEAN BIOPHYSICS JOURNAL 2022; 51:265-282. [PMID: 35303138 PMCID: PMC9035012 DOI: 10.1007/s00249-022-01593-1] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Revised: 02/25/2022] [Accepted: 02/26/2022] [Indexed: 12/04/2022]
Abstract
In protein simulation or structure refinement based on values of observable quantities measured in (aqueous) solution, solvent (water) molecules may be explicitly treated, omitted, or represented by a potential of mean-solvation-force term, depending on protein coordinates only, in the force field used. These three approaches are compared for hen egg white lysozyme (HEWL). This 129-residue non-spherical protein contains a variety of secondary-structure elements, and ample experimental data are available: 1630 atom–atom Nuclear Overhauser Enhancement (NOE) upper distance bounds, 213 3 J-couplings and 200 S2 order parameters. These data are used to compare the performance of the three approaches. It is found that a molecular dynamics (MD) simulation in explicit water approximates the experimental data much better than stochastic dynamics (SD) simulation in vacuo without or with a solvent-accessible-surface-area (SASA) implicit-solvation term added to the force field. This is due to the missing energetic and entropic contributions and hydrogen-bonding capacities of the water molecules and the missing dielectric screening effect of this high-permittivity solvent. Omission of explicit water molecules leads to compaction of the protein, an increased internal strain, distortion of exposed loop and turn regions and excessive intra-protein hydrogen bonding. As a consequence, the conformation and dynamics of groups on the surface of the protein, which may play a key role in protein–protein interactions or ligand or substrate binding, may be incorrectly modelled. It is thus recommended to include water molecules explicitly in structure refinement of proteins in aqueous solution based on nuclear magnetic resonance (NMR) or other experimentally measured data.
Collapse
|
10
|
Kanchi PK, Dasmahapatra AK. Enhancing the binding of the β-sheet breaker peptide LPFFD to the amyloid-β fibrils by aromatic modifications: A molecular dynamics simulation study. Comput Biol Chem 2021; 92:107471. [PMID: 33706107 DOI: 10.1016/j.compbiolchem.2021.107471] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2020] [Revised: 02/25/2021] [Accepted: 02/27/2021] [Indexed: 11/25/2022]
Abstract
Alzheimer's is a fatal neurodegenerative disease for which there is no cure at present. The disease is characterized by the presence of plaques in the brains of a patient, which are composed mainly of aggregates of the amyloid-β peptide in the form of β-sheet fibrils. Here, we investigated the possibility of exploiting the superior binding ability of aromatic amino acids to a particular model of the amyloid-β fibrils. which is a difficult target for drug design. The β-sheet breaker peptide LPFFD was modified with aromatic amino acids and its binding to these fibrils was studied. We found that the orientation and the electrostatic complementarity of the modified peptide with respect to the fibrils played a crucial role in determining whether its binding was improved by the aromatic amino acids. The modified LPFFD peptides were able to bind to those fibril residues. which are important in the aggregation of amyloid-β peptides and thus can potentially inhibit the further aggregation of the amyloid-beta peptides by blocking their interactions. We found that the tryptophan modified LPFFD peptides had the best binding affinities. In most cases, the aromatic amino acids in the N-terminus of the modified peptides made more contacts with the fibrils than those in the C-terminus. We also found that increasing the aromatic content did not significantly improve the binding of the LPFFD peptide to the fibrils. Our study can serve as a basis for the design of novel peptide-based drugs for Alzheimer's disease in which aromatic interactions play an important role.
Collapse
Affiliation(s)
- Pavan Krishna Kanchi
- Department of Chemical Engineering, Indian Institute of Technology Guwahati, Guwahati, 781039, Assam, India
| | - Ashok Kumar Dasmahapatra
- Department of Chemical Engineering, Indian Institute of Technology Guwahati, Guwahati, 781039, Assam, India; Center for Nanotechnology, Indian Institute of Technology Guwahati, Guwahati, 781039, Assam, India.
| |
Collapse
|
11
|
PupStruct: Prediction of Pupylated Lysine Residues Using Structural Properties of Amino Acids. Genes (Basel) 2020; 11:genes11121431. [PMID: 33260770 PMCID: PMC7761138 DOI: 10.3390/genes11121431] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2020] [Revised: 11/23/2020] [Accepted: 11/23/2020] [Indexed: 12/23/2022] Open
Abstract
Post-translational modification (PTM) is a critical biological reaction which adds to the diversification of the proteome. With numerous known modifications being studied, pupylation has gained focus in the scientific community due to its significant role in regulating biological processes. The traditional experimental practice to detect pupylation sites proved to be expensive and requires a lot of time and resources. Thus, there have been many computational predictors developed to challenge this issue. However, performance is still limited. In this study, we propose another computational method, named PupStruct, which uses the structural information of amino acids with a radial basis kernel function Support Vector Machine (SVM) to predict pupylated lysine residues. We compared PupStruct with three state-of-the-art predictors from the literature where PupStruct has validated a significant improvement in performance over them with statistical metrics such as sensitivity (0.9234), specificity (0.9359), accuracy (0.9296), precision (0.9349), and Mathew’s correlation coefficient (0.8616) on a benchmark dataset.
Collapse
|
12
|
Singh N, Chaput L, Villoutreix BO. Fast Rescoring Protocols to Improve the Performance of Structure-Based Virtual Screening Performed on Protein-Protein Interfaces. J Chem Inf Model 2020; 60:3910-3934. [PMID: 32786511 DOI: 10.1021/acs.jcim.0c00545] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Protein-protein interactions (PPIs) are attractive targets for drug design because of their essential role in numerous cellular processes and disease pathways. However, in general, PPIs display exposed binding pockets at the interface, and as such, have been largely unexploited for therapeutic interventions with low-molecular weight compounds. Here, we used docking and various rescoring strategies in an attempt to recover PPI inhibitors from a set of active and inactive molecules for 11 targets collected in ChEMBL and PubChem. Our focus is on the screening power of the various developed protocols and on using fast approaches so as to be able to apply such a strategy to the screening of ultralarge libraries in the future. First, we docked compounds into each target using the fast "pscreen" mode of the structure-based virtual screening (VS) package Surflex. Subsequently, the docking poses were postprocessed to derive a set of 3D topological descriptors: (i) shape similarity and (ii) interaction fingerprint similarity with a co-crystallized inhibitor, (iii) solvent-accessible surface area, and (iv) extent of deviation from the geometric center of a reference inhibitor. The derivatized descriptors, together with descriptor-scaled scoring functions, were utilized to investigate possible impacts on VS performance metrics. Moreover, four standalone scoring functions, RF-Score-VS (machine-learning), DLIGAND2 (knowledge-based), Vinardo (empirical), and X-SCORE (empirical), were employed to rescore the PPI compounds. Collectively, the results indicate that the topological scoring algorithms could be valuable both at a global level, with up to 79% increase in areas under the receiver operating characteristic curve for some targets, and in early stages, with up to a 4-fold increase in enrichment factors at 1% of the screened collections. Outstandingly, DLIGAND2 emerged as the best scoring function on this data set, outperforming all rescoring techniques in terms of VS metrics. The described methodology could help in the rational design of small-molecule PPI inhibitors and has direct applications in many therapeutic areas, including cancer, CNS, and infectious diseases such as COVID-19.
Collapse
Affiliation(s)
- Natesh Singh
- Université de Lille, Inserm, Institut Pasteur de Lille, U1177-Drugs and Molecules for Living Systems, F-59000 Lille, France
| | - Ludovic Chaput
- Université de Lille, Inserm, Institut Pasteur de Lille, U1177-Drugs and Molecules for Living Systems, F-59000 Lille, France
| | - Bruno O Villoutreix
- Université de Lille, Inserm, Institut Pasteur de Lille, U1177-Drugs and Molecules for Living Systems, F-59000 Lille, France
| |
Collapse
|
13
|
Vassetti D, Civalleri B, Labat F. Analytical calculation of the solvent-accessible surface area and its nuclear gradients by stereographic projection: A general approach for molecules, polymers, nanotubes, helices, and surfaces. J Comput Chem 2020; 41:1464-1479. [PMID: 32212337 DOI: 10.1002/jcc.26191] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2020] [Revised: 03/07/2020] [Accepted: 03/09/2020] [Indexed: 01/19/2023]
Abstract
In this article, we explore an alternative to the analytical Gauss-Bonnet approach for computing the solvent-accessible surface area (SASA) and its nuclear gradients. These two key quantities are required to evaluate the nonelectrostatic contribution to the solvation energy and its nuclear gradients in implicit solvation models. We extend a previously proposed analytical approach for finite systems based on the stereographic projection technique to infinite periodic systems such as polymers, nanotubes, helices, or surfaces and detail its implementation in the Crystal code. We provide the full derivation of the SASA nuclear gradients, and introduce an iterative perturbation scheme of the atomic coordinates to stabilize the gradients calculation for certain difficult symmetric systems. An excellent agreement of computed SASA with reference analytical values is found for finite systems, while the SASA size-extensivity is verified for infinite periodic systems. In addition, correctness of the analytical gradients is confirmed by the excellent agreement obtained with numerical gradients and by the translational invariance achieved, both for finite and infinite periodic systems. Overall therefore, the stereographic projection approach appears as a general, simple, and efficient technique to compute the key quantities required for the calculation of the nonelectrostatic contribution to the solvation energy and its nuclear gradients in implicit solvation models applicable to both finite and infinite periodic systems.
Collapse
Affiliation(s)
- Dario Vassetti
- Chimie ParisTech, PSL University, CNRS, Institute of Chemistry for Life and Health Sciences, Chemical Theory and Modelling Group, F-75005 Paris, France
| | - Bartolomeo Civalleri
- Department of Chemistry, NIS and INSTM Reference Centre, University of Turin, Via P. Giuria 7, I-10125 Torino, Italy
| | - Frédéric Labat
- Chimie ParisTech, PSL University, CNRS, Institute of Chemistry for Life and Health Sciences, Chemical Theory and Modelling Group, F-75005 Paris, France
| |
Collapse
|
14
|
Novichkova DA, Lushchekina SV, Dym O, Masson P, Silman I, Sussman JL. The four-helix bundle in cholinesterase dimers: Structural and energetic determinants of stability. Chem Biol Interact 2019; 309:108699. [PMID: 31202688 DOI: 10.1016/j.cbi.2019.06.012] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2019] [Revised: 05/30/2019] [Accepted: 06/06/2019] [Indexed: 11/28/2022]
Abstract
The crystal structures of truncated forms of cholinesterases provide good models for assessing the role of non-covalent interactions in dimer assembly in the absence of cross-linking disulfide bonds. These structures identify the four-helix bundle that serves as the interface for formation of acetylcholinesterase and butyrylcholinesterase dimers. Here we performed a theoretical comparison of the structural and energetic factors governing dimerization. This included identification of inter-subunit and intra-subunit hydrogen bonds and hydrophobic interactions, evaluation of solvent-accessible surfaces, and estimation of electrostatic contributions to dimerization. To reveal the contribution to dimerization of individual amino acids within the contact area, free energy perturbation alanine screening was performed. Markov state modelling shows that the loop between the α13 and α14 helices in BChE is unstable, and occupies 4 macro-states. The order of magnitude of mean first passage times between these macrostates is ~10-8 s. Replica exchange molecular dynamics umbrella sampling calculations revealed that the free energy of human BChE dimerization is -15.5 kcal/mol, while that for human AChE is -26.4 kcal/mol. Thus, the C-terminally truncated human butyrylcholinesterase dimer is substantially less stable than that of human acetylcholinesterase. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:CHEMBIOINT:1.
Collapse
Affiliation(s)
- Dana A Novichkova
- N.M. Emanuel Institute of Biochemical Physics of Russian Academy of Sciences, 4 Kosygina St., Moscow, 119334, Russia
| | - Sofya V Lushchekina
- N.M. Emanuel Institute of Biochemical Physics of Russian Academy of Sciences, 4 Kosygina St., Moscow, 119334, Russia.
| | - Orly Dym
- Weizmann Institute of Science, 234 Herzl St., Rehovot, 7610001, Israel
| | - Patrick Masson
- Kazan Federal University, Neuropharmacology Laboratory, 18 Kremlevskaya St., Kazan, 420008, Russia
| | - Israel Silman
- Weizmann Institute of Science, 234 Herzl St., Rehovot, 7610001, Israel
| | - Joel L Sussman
- Weizmann Institute of Science, 234 Herzl St., Rehovot, 7610001, Israel
| |
Collapse
|
15
|
Famiglini V, Silvestri R. Indolylarylsulfones, a fascinating story of highly potent human immunodeficiency virus type 1 non-nucleoside reverse transcriptase inhibitors. Antivir Chem Chemother 2019; 26:2040206617753443. [PMID: 29417826 PMCID: PMC5890576 DOI: 10.1177/2040206617753443] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Indolylarylsulfones are a potent class of human immunodeficiency virus type 1 non-nucleoside reverse transcriptase inhibitors. In this review, the structure activity relationship (SAR) studies to improve the profile of sulfone L-737,126 discovered by Merck AG have been analysed with focus on introduction of the 3′,5′-dimethyl groups at the 3-phenylsulfonyl moiety, the 2-hydroxyethyl tail at the indole-2-carboxamide nitrogen, coupling of the carboxamide nitrogen with one or two glycinamide and alaninamide units, a fluorine atom at position 4 of the indole ring and correlation between configuration of the asymmetric centre and linker length. IAS derivatives look like promising drug candidates for the treatment of AIDS and related infections in combination with other antiretroviral agents.
Collapse
Affiliation(s)
- Valeria Famiglini
- Department of Drug Chemistry and Technologies, Sapienza University of Rome, Laboratory affiliated to Istituto Pasteur Italia - Fondazione Cenci Bolognetti, Roma, Italy
| | - Romano Silvestri
- Department of Drug Chemistry and Technologies, Sapienza University of Rome, Laboratory affiliated to Istituto Pasteur Italia - Fondazione Cenci Bolognetti, Roma, Italy
| |
Collapse
|
16
|
Aytenfisu AH, Simon R, MacKerell AD. Impact of branching on the conformational heterogeneity of the lipopolysaccharide from Klebsiella pneumoniae: Implications for vaccine design. Carbohydr Res 2019; 475:39-47. [PMID: 30818097 DOI: 10.1016/j.carres.2019.02.003] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2019] [Revised: 02/09/2019] [Accepted: 02/13/2019] [Indexed: 02/08/2023]
Abstract
Resistance of Klebsiella pneumoniae (KP) to antibiotics has motivated the development of an efficacious KP human vaccine that would not be subject to antibiotic resistance. Klebsiella lipopolysaccharide (LPS) associated O polysaccharide (OPS) types have provoked broad interest as a vaccine antigen as there are only 4 that predominate worldwide (O1, O2a, O3, O5). Klebsiella O1 and O2 OPS are polygalactans that share a common D-Gal-I structure, for which a variant D-Gal-III was recently discovered. To understand the potential impact of this variability on antigenicity, a detailed molecular picture of the conformational differences associated with the addition of the D-Gal-III (1 → 4)-α-Galp branch is presented using enhanced-sampling molecular dynamics simulations. In D-Gal-I two major conformational states are observed while the presence of the 1 → 4 branch in D-Gal-III resulted in only a single dominant extended state. Stabilization of the more folded states in D-Gal-I is due to a O4-H⋯O2 hydrogen bond in the linear backbone that cannot occur in D-Gal-III as the O4 is in the Galp(1 → 4)Galp glycosidic linkage. The impact of branching in D-Gal-III also significantly decreases the accessibility of the monosaccharides in the linear backbone region of D-Gal-I, while the accessibility of the terminal D-Gal-II region of the OPS is not substantially altered. The present results suggest that a vaccine that targets both the D-Gal-I and D-Gal-III LPS can be developed by using D-Gal-III as the antigen combined with cross-reactivity experiments using the Gal-II polysaccharide to assure that this region of the LPS is the primary epitope of the antigen.
Collapse
Affiliation(s)
- Asaminew H Aytenfisu
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD, 21201, United States
| | - Raphael Simon
- Center for Vaccine Development, Institute for Global Health, School of Medicine, University of Maryland, Baltimore, MD, 21201, United States
| | - Alexander D MacKerell
- University of Maryland Computer-Aided Drug Design Center, Department of Pharmaceutical Sciences, School of Pharmacy, University of Maryland, Baltimore, MD, 21201, United States.
| |
Collapse
|
17
|
Steinbach PJ. Peptide and Protein Structure Prediction with a Simplified Continuum Solvent Model. J Phys Chem B 2018; 122:11355-11362. [PMID: 30230838 DOI: 10.1021/acs.jpcb.8b07264] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
A continuum solvent model based on screened Coulomb potentials has been simplified and parametrized to sample native-like structures in replica-exchange simulations of each of six different peptides and miniproteins. Low-energy, native, and non-native structures were used to iteratively refine 11 parameter values. The centroid of the largest cluster of structures sampled in simulations initiated from an extended conformation represents the predicted structure. The main-chain rms deviation of this prediction from the experimental structure was 0.47 Å for the 12-residue Trp-zip2, 0.86 Å for the 14-residue MBH12, 2.53 Å for the 17-residue U(1-17)T9D, 2.03 Å for the 20-residue BS1, 1.08 Å for the 20-residue Trp-cage, and 3.64 Å for the 35-residue villin headpiece subdomain HP35. The centroid of the sixth largest cluster sampled for HP35 deviated by 0.91 Å. The CHARMM22/CMAP force field was used, with an additional ψ torsion term for residues other than glycine and proline. Six parameters govern the dielectric response of the continuum solvent, and four values of surface tension approximate nonpolar effects. An atom's self-energy and interaction energies are screened independently, each depending on whether the atom is part of a charged group, a neutral hydrogen-bonding main-chain group, or any other neutral group. The parameters inferred result in strong main-chain hydrogen bonds, consistent with the view that protein folding is dominated by the formation of these bonds. (1,2) Conformations of MBH12 and BS1 were excluded from the energy-function refinement, suggesting the parameters, referred to as SCP18, are transferable. An efficient estimate of solvent-accessible surface area is also described.
Collapse
Affiliation(s)
- Peter J Steinbach
- Center for Molecular Modeling, Center for Information Technology , National Institutes of Health , Bethesda , Maryland 20892 , United States
| |
Collapse
|
18
|
Huang H, Simmerling C. Fast Pairwise Approximation of Solvent Accessible Surface Area for Implicit Solvent Simulations of Proteins on CPUs and GPUs. J Chem Theory Comput 2018; 14:5797-5814. [PMID: 30303377 DOI: 10.1021/acs.jctc.8b00413] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We propose a pairwise and readily parallelizable SASA-based nonpolar solvation approach for protein simulations, inspired by our previous pairwise GB polar solvation model development. In this work, we developed a novel function to estimate the atomic and molecular SASAs of proteins, which results in comparable accuracy as the LCPO algorithm in reproducing numerical icosahedral-based SASA values. Implemented in Amber software and tested on consumer GPUs, our pwSASA method reasonably reproduces LCPO simulation results, but accelerates MD simulations up to 30 times compared to the LCPO implementation, which is greatly desirable for protein simulations facing sampling challenges. The value of incorporating the nonpolar term in implicit solvent simulations is explored on a peptide fragment containing the hydrophobic core of HP36 and evaluating thermal stability profiles of four small proteins.
Collapse
Affiliation(s)
- He Huang
- Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States.,Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| | - Carlos Simmerling
- Laufer Center for Physical and Quantitative Biology , Stony Brook University , Stony Brook , New York 11794 , United States.,Department of Chemistry , Stony Brook University , Stony Brook , New York 11794 , United States
| |
Collapse
|
19
|
Binding Direction-Based Two-Dimensional Flattened Contact Area Computing Algorithm for Protein–Protein Interactions. Molecules 2017; 22:molecules22101722. [PMID: 29027921 PMCID: PMC6151622 DOI: 10.3390/molecules22101722] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Revised: 10/06/2017] [Accepted: 10/12/2017] [Indexed: 11/16/2022] Open
Abstract
Interactions between protein molecules are essential for the assembly, function, and regulation of proteins. The contact region between two protein molecules in a protein complex is usually complementary in shape for both molecules and the area of the contact region can be used to estimate the binding strength between two molecules. Although the area is a value calculated from the three-dimensional surface, it cannot represent the three-dimensional shape of the surface. Therefore, we propose an original concept of two-dimensional contact area which provides further information such as the ruggedness of the contact region. We present a novel algorithm for calculating the binding direction between two molecules in a protein complex, and then suggest a method to compute the two-dimensional flattened area of the contact region between two molecules based on the binding direction.
Collapse
|
20
|
Koch M, Saphiannikova M, Santer S, Guskova O. Photoisomers of Azobenzene Star with a Flat Core: Theoretical Insights into Multiple States from DFT and MD Perspective. J Phys Chem B 2017; 121:8854-8867. [PMID: 28832166 DOI: 10.1021/acs.jpcb.7b07350] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This study focuses on comparing physical properties of photoisomers of an azobenzene star with benzene-1,3,5-tricarboxamide core. Three azobenzene arms of the molecule undergo a reversible trans-cis isomerization upon UV-vis light illumination giving rise to multiple states from the planar all-trans one, via two mixed states to the kinked all-cis isomer. Employing density functional theory, we characterize the structural and photophysical properties of each state indicating a role the planar core plays in the coupling between azobenzene chromophores. To characterize the light-triggered switching of solvophilicity/solvophobicity of the star, the difference in solvation free energy is calculated for the transfer of an azobenzene star from its gas phase to implicit or explicit solvents. For the latter case, classical all-atom molecular dynamics simulations of aqueous solutions of azobenzene star are performed employing the polymer consistent force field to shed light on the thermodynamics of explicit hydration as a function of the isomerization state and on the structuring of water around the star. From the analysis of two contributions to the free energy of hydration, the nonpolar van der Waals and the electrostatic terms, it is concluded that isomerization specificity largely determines the polarity of the molecule and the solute-solvent electrostatic interactions. This convertible hydrophilicity/hydrophobicity together with readjustable occupied volume and the surface area accessible to water, affects the self-assembly/disassembly of the azobenzene star with a flat core triggered by light.
Collapse
Affiliation(s)
- Markus Koch
- Institute Theory of Polymers, Leibniz Institute of Polymer Research Dresden , Hohe Str. 6, 01069 Dresden, Germany.,Institute of Theoretical Physics, Technische Universität Dresden , Zellescher Weg 17, 01069 Dresden, Germany
| | - Marina Saphiannikova
- Institute Theory of Polymers, Leibniz Institute of Polymer Research Dresden , Hohe Str. 6, 01069 Dresden, Germany.,Dresden Center for Computational Materials Science (DCMS), Technische Universität Dresden , 01069 Dresden, Germany
| | - Svetlana Santer
- Institute of Physics and Astronomy, University of Potsdam , Karl-Liebknecht-Str. 24/25, 14476 Potsdam, Germany
| | - Olga Guskova
- Institute Theory of Polymers, Leibniz Institute of Polymer Research Dresden , Hohe Str. 6, 01069 Dresden, Germany.,Dresden Center for Computational Materials Science (DCMS), Technische Universität Dresden , 01069 Dresden, Germany
| |
Collapse
|
21
|
Wei S, Brooks CL, Frank AT. A rapid solvent accessible surface area estimator for coarse grained molecular simulations. J Comput Chem 2017; 38:1270-1274. [DOI: 10.1002/jcc.24709] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Revised: 11/08/2016] [Accepted: 11/17/2016] [Indexed: 11/06/2022]
Affiliation(s)
- Shuai Wei
- Department of Chemistry; University of Michigan; Ann Arbor MI 48109
| | - Charles L. Brooks
- Departments of Biophysics and Chemistry; University of Michigan; Ann Arbor MI 48109
| | - Aaron T. Frank
- Departments of Biophysics and Chemistry; University of Michigan; Ann Arbor MI 48109
| |
Collapse
|
22
|
Zhang B, Kilburg D, Eastman P, Pande VS, Gallicchio E. Efficient gaussian density formulation of volume and surface areas of macromolecules on graphical processing units. J Comput Chem 2017; 38:740-752. [PMID: 28160511 DOI: 10.1002/jcc.24745] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2016] [Revised: 01/05/2017] [Accepted: 01/08/2017] [Indexed: 11/07/2022]
Abstract
We present an algorithm to efficiently compute accurate volumes and surface areas of macromolecules on graphical processing unit (GPU) devices using an analytic model which represents atomic volumes by continuous Gaussian densities. The volume of the molecule is expressed by means of the inclusion-exclusion formula, which is based on the summation of overlap integrals among multiple atomic densities. The surface area of the molecule is obtained by differentiation of the molecular volume with respect to atomic radii. The many-body nature of the model makes a port to GPU devices challenging. To our knowledge, this is the first reported full implementation of this model on GPU hardware. To accomplish this, we have used recursive strategies to construct the tree of overlaps and to accumulate volumes and their gradients on the tree data structures so as to minimize memory contention. The algorithm is used in the formulation of a surface area-based non-polar implicit solvent model implemented as an open source plug-in (named GaussVol) for the popular OpenMM library for molecular mechanics modeling. GaussVol is 50 to 100 times faster than our best optimized implementation for the CPUs, achieving speeds in excess of 100 ns/day with 1 fs time-step for protein-sized systems on commodity GPUs. © 2017 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Baofeng Zhang
- Department of Chemistry, Brooklyn College of the City University of New York, 2900 Bedford Avenue, Brooklyn, New York, 11210
| | - Denise Kilburg
- Department of Chemistry, Brooklyn College of the City University of New York, 2900 Bedford Avenue, Brooklyn, New York, 11210.,Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York, 10016
| | - Peter Eastman
- Department of Bioengineering, Stanford University, Stanford, California, 94035
| | - Vijay S Pande
- Department of Chemistry, Stanford University, Stanford, California, 94035
| | - Emilio Gallicchio
- Department of Chemistry, Brooklyn College of the City University of New York, 2900 Bedford Avenue, Brooklyn, New York, 11210.,Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York, 10016
| |
Collapse
|
23
|
Abstract
Computational protein sequence design is the rational design based on computer simulation of new protein molecules to fold to target three-dimensional structures, with the ultimate goal of designing novel functions. It requires a good understanding of the thermodynamic equilibrium properties of the protein of interest. Here, we consider the contribution of the solvent to the stability of the protein. We describe implicit solvent models, focusing on approximations of their nonpolar components using geometric potentials. We consider the surface area (SA) model in which the nonpolar solvation free energy is expressed as a sum of the contributions of all atoms, assumed to be proportional to their accessible surface areas (ASAs). We briefly review existing numerical and analytical approaches that compute the ASA. We describe in more detail the alpha shape theory as it provides a unifying mathematical framework that enables the analytical calculations of the surface area of a macromolecule represented as a union of balls.
Collapse
Affiliation(s)
- Jie Li
- Computational and Systems Biology Group, Genome Institute of Singapore, Agency for Science, Technology and Research, Singapore, Singapore
| | - Patrice Koehl
- Department of Computer Science and Genome Center, University of California, Davis, CA, 95616, USA.
| |
Collapse
|
24
|
Abstract
Molecular replacement is a method for solving the crystallographic phase problem using an atomic model for the target structure. State-of-the-art methods have moved the field significantly from when it was first envisaged as a method for solving cases of high homology and completeness between a model and target structure. Improvements brought about by application of maximum likelihood statistics mean that various errors in the model and pathologies in the data can be accounted for, so that cases hitherto thought to be intractable are standardly solvable. As a result, molecular replacement phasing now accounts for the lion's share of structures deposited in the Protein Data Bank. However, there will always be cases at the fringes of solvability. I discuss here the approaches that will help tackle challenging molecular replacement cases.
Collapse
Affiliation(s)
- Airlie J McCoy
- Department of Haematology, Cambridge Institute for Medical Research, University of Cambridge, Wellcome Trust/MRC Building, Hills Road, Cambridge, CB2 0XY, UK.
| |
Collapse
|
25
|
Abstract
We consider two random sequential packing processes in which spheres of unit radius are randomly attached to the surface of a fixed unit sphere. Independent random spheres are generated and added successively, provided there is no overlap with previous spheres. In model 1, the process stops when a trial sphere intersects one of the previously-accepted spheres. In model 2, random sequential packing, any such overlapping trial sphere is discarded and the next random sphere is tried, until it is impossible to add any further spheres.Previous workers have conjectured convincingly that no exact analytical solution is possible for this type of problem. We use Monte Carlo simulation methods to estimate transition probabilities for the two models. Because some probabilities are extremely small, a simulation using independent repetitions of the model would be inefficient. We designed a branching process of conditionally binomial trials, and performed over 108 trials on a supercomputer.
Collapse
|
26
|
Tavousi P, Behandish M, Ilieş HT, Kazerounian K. Protofold II: Enhanced Model and Implementation for Kinetostatic Protein Folding1. J Nanotechnol Eng Med 2016. [DOI: 10.1115/1.4032759] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
A reliable prediction of three-dimensional (3D) protein structures from sequence data remains a big challenge due to both theoretical and computational difficulties. We have previously shown that our kinetostatic compliance method (KCM) implemented into the Protofold package can overcome some of the key difficulties faced by other de novo structure prediction methods, such as the very small time steps required by the molecular dynamics (MD) approaches or the very large number of samples needed by the Monte Carlo (MC) sampling techniques. In this paper, we improve the free energy formulation used in Protofold by including the typically underrated entropic effects, imparted due to differences in hydrophobicity of the chemical groups, which dominate the folding of most water-soluble proteins. In addition to the model enhancement, we revisit the numerical implementation by redesigning the algorithms and introducing efficient data structures that reduce the expected complexity from quadratic to linear. Moreover, we develop and optimize parallel implementations of the algorithms on both central and graphics processing units (CPU/GPU) achieving speed-ups up to two orders of magnitude on the GPU. Our simulations are consistent with the general behavior observed in the folding process in aqueous solvent, confirming the effectiveness of model improvements. We report on the folding process at multiple levels, namely, the formation of secondary structural elements and tertiary interactions between secondary elements or across larger domains. We also observe significant enhancements in running times that make the folding simulation tractable for large molecules.
Collapse
Affiliation(s)
- Pouya Tavousi
- Kinematics Design Laboratory, Department of Mechanical Engineering, University of Connecticut, Storrs, CT 06269 e-mail:
| | - Morad Behandish
- Computational Design Laboratory, Department of Mechanical Engineering, University of Connecticut, Storrs, CT 06269 e-mail:
| | - Horea T. Ilieş
- Department of Mechanical Engineering, University of Connecticut, Storrs, CT 06269 e-mail:
| | - Kazem Kazerounian
- Department of Mechanical Engineering, University of Connecticut, Storrs, CT 06269 e-mail:
| |
Collapse
|
27
|
Focus on Chirality of HIV-1 Non-Nucleoside Reverse Transcriptase Inhibitors. Molecules 2016; 21:molecules21020221. [PMID: 26891289 PMCID: PMC6273187 DOI: 10.3390/molecules21020221] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2016] [Revised: 02/04/2016] [Accepted: 02/08/2016] [Indexed: 02/02/2023] Open
Abstract
Chiral HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) are of great interest since one enantiomer is often more potent than the corresponding counterpart against the HIV-1 wild type (WT) and the HIV-1 drug resistant mutant strains. This review exemplifies the various studies made to investigate the effect of chirality on the antiretroviral activity of top HIV-1 NNRTI compounds, such as nevirapine (NVP), efavirenz (EFV), alkynyl- and alkenylquinazolinone DuPont compounds (DPC), diarylpyrimidine (DAPY), dihydroalkyloxybenzyloxopyrimidine (DABO), phenethylthiazolylthiourea (PETT), indolylarylsulfone (IAS), arylphosphoindole (API) and trifluoromethylated indole (TFMI) The chiral separation, the enantiosynthesis, along with the biological properties of these HIV-1 NNRTIs, are discussed.
Collapse
|
28
|
Duignan TT, Parsons DF, Ninham BW. A continuum solvent model of ion-ion interactions in water. Phys Chem Chem Phys 2015; 16:22014-27. [PMID: 25205066 DOI: 10.1039/c4cp02822h] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The calculation of ion-ion interactions in water is a problem of long standing importance. Modelling these interactions is a prerequisite to explaining Hofmeister (specific ion) effects. We here generalize our solvation model of ions to calculate the free energy of two ions in water as a function of separation. The same procedure has previously been applied to calculate ion interactions with the air-water interface successfully. The Conductor like Screening Model (COSMO) is used. This treats the ions on a quantum mechanical level and calculates numerically the electrostatic response of the surrounding solvent. Estimates of the change in the cavity formation energy and the change in the ion-water dispersion energy as the ions approach are included separately. The calculated interaction potentials are too attractive and this is a significant issue. However, they do reproduce the affinity of similarly sized ions for each other, which is a crucial property of these potentials. They are also oscillatory, another important property. We normalize the potentials to reduce the over-attraction, and good correlation with experimental values is achieved. We identify the driving contributions to this like-prefers-like behaviour. We then put forward a plausible hypothesis for the over-attraction of the potentials. An agreeable feature of our approach is that it does not rely on salt specific parameters deliberately adjusted to reproduce experimental values.
Collapse
Affiliation(s)
- Timothy T Duignan
- Applied Mathematics Department, Australian National University, Canberra, ACT 0200, Australia.
| | | | | |
Collapse
|
29
|
Faller CE, Guvench O. Sulfation and cation effects on the conformational properties of the glycan backbone of chondroitin sulfate disaccharides. J Phys Chem B 2015; 119:6063-73. [PMID: 25906376 DOI: 10.1021/jp511431q] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Chondroitin sulfate (CS) is one of several glycosaminoglycans that are major components of proteoglycans. A linear polymer consisting of repeats of the disaccharide -4GlcAβ1-3GalNAcβ1-, CS undergoes differential sulfation resulting in five unique sulfation patterns. Because of the dimer repeat, the CS glycosidic "backbone" has two distinct sets of conformational degrees of freedom defined by pairs of dihedral angles: (ϕ1, ψ1) about the β1-3 glycosidic linkage and (ϕ2, ψ2) about the β1-4 glycosidic linkage. Differential sulfation and the possibility of cation binding, combined with the conformational flexibility and biological diversity of CS, complicate experimental efforts to understand CS three-dimensional structures at atomic resolution. Therefore, all-atom explicit-solvent molecular dynamics simulations with Adaptive Biasing Force sampling of the CS backbone were applied to obtain high-resolution, high-precision free energies of CS disaccharides as a function of all possible backbone geometries. All 10 disaccharides (β1-3 vs β1-4 linkage × five different sulfation patterns) were studied; additionally, ion effects were investigated by considering each disaccharide in the presence of either neutralizing sodium or calcium cations. GlcAβ1-3GalNAc disaccharides have a single, broad, thermodynamically important free-energy minimum, whereas GalNAcβ1-4GlcA disaccharides have two such minima. Calcium cations but not sodium cations bind to the disaccharides, and binding is primarily to the GlcA -COO(-) moiety as opposed to sulfate groups. This binding alters the glycan backbone thermodynamics in instances where a calcium cation bound to -COO(-) can act to bridge and stabilize an interaction with an adjacent sulfate group, whereas, in the absence of this cation, the proximity of a sulfate group to -COO(-) results in two like charges being both desolvated and placed adjacent to each other and is found to be destabilizing. In addition to providing information on sulfation and cation effects, the present results can be applied to building models of CS polymers and as a point of comparison in studies of CS polymer backbone dynamics and thermodynamics.
Collapse
Affiliation(s)
- Christina E Faller
- Department of Pharmaceutical Sciences, University of New England College of Pharmacy, 716 Stevens Avenue, Portland, Maine 04103, United States
| | - Olgun Guvench
- Department of Pharmaceutical Sciences, University of New England College of Pharmacy, 716 Stevens Avenue, Portland, Maine 04103, United States
| |
Collapse
|
30
|
Janin J, Wodak SJ, Lensink MF, Velankar S. Assessing Structural Predictions of Protein-Protein Recognition: The CAPRI Experiment. REVIEWS IN COMPUTATIONAL CHEMISTRY 2015. [DOI: 10.1002/9781118889886.ch4] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
31
|
Li L, Wang L, Alexov E. On the energy components governing molecular recognition in the framework of continuum approaches. Front Mol Biosci 2015; 2:5. [PMID: 25988173 PMCID: PMC4429657 DOI: 10.3389/fmolb.2015.00005] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Accepted: 02/04/2015] [Indexed: 01/14/2023] Open
Abstract
Molecular recognition is a process that brings together several biological macromolecules to form a complex and one of the most important characteristics of the process is the binding free energy. Various approaches exist to model the binding free energy, provided the knowledge of the 3D structures of bound and unbound molecules. Among them, continuum approaches are quite appealing due to their computational efficiency while at the same time providing predictions with reasonable accuracy. Here we review recent developments in the field emphasizing on the importance of adopting adequate description of physical processes taking place upon the binding. In particular, we focus on the efforts aiming at capturing some of the atomistic details of the binding phenomena into the continuum framework. When possible, the energy components are reviewed independently of each other. However, it is pointed out that rigorous approaches should consider all energy contributions on the same footage. The two major schemes for utilizing the individual energy components to predict binding affinity are outlined as well.
Collapse
Affiliation(s)
- Lin Li
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University Clemson, SC, USA
| | - Lin Wang
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University Clemson, SC, USA
| | - Emil Alexov
- Computational Biophysics and Bioinformatics, Department of Physics, Clemson University Clemson, SC, USA
| |
Collapse
|
32
|
Mondal S, Ghosh S, De S. Atomistic level molecular dynamics simulation on the solubilization mechanism of aromatic molecules in anionic micelles. RSC Adv 2015. [DOI: 10.1039/c5ra20855f] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Benzene is solubilized in the palisade layer, toluene in the core (or central region), phenol at the surface and pyridine in micelle–water interface. During micellar enhanced ultrafiltration, phenol will have least rejection while benzene is maximum.
Collapse
Affiliation(s)
- Sourav Mondal
- Department of Chemical Engineering
- Indian Institute of Technology Kharagpur
- Kharagpur – 721302
- India
| | - Sudipto Ghosh
- Department of Metallurgical and Materials Engineering
- Indian Institute of Technology Kharagpur
- Kharagpur – 721302
- India
| | - Sirshendu De
- Department of Chemical Engineering
- Indian Institute of Technology Kharagpur
- Kharagpur – 721302
- India
| |
Collapse
|
33
|
Famiglini V, La Regina G, Coluccia A, Pelliccia S, Brancale A, Maga G, Crespan E, Badia R, Riveira-Muñoz E, Esté JA, Ferretti R, Cirilli R, Zamperini C, Botta M, Schols D, Limongelli V, Agostino B, Novellino E, Silvestri R. Indolylarylsulfones carrying a heterocyclic tail as very potent and broad spectrum HIV-1 non-nucleoside reverse transcriptase inhibitors. J Med Chem 2014; 57:9945-57. [PMID: 25418038 DOI: 10.1021/jm5011622] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We synthesized new indolylarylsulfone (IAS) derivatives carrying a heterocyclic tail at the indole-2-carboxamide nitrogen as potential anti-HIV/AIDS agents. Several new IASs yielded EC50 values <1.0 nM against HIV-1 WT and mutant strains in MT-4 cells. The (R)-11 enantiomer proved to be exceptionally potent against the whole viral panel; in the reverse transcriptase (RT) screening assay, it was remarkably superior to NVP and EFV and comparable to ETV. The binding poses were consistent with the one previously described for the IAS non-nucleoside reverse transcriptase inhibitors. Docking studies showed that the methyl group of (R)-11 points toward the cleft created by the K103N mutation, different from the corresponding group of (S)-11. By calculating the solvent-accessible surface, we observed that the exposed area of RT in complex with (S)-11 was larger than the area of the (R)-11 complex. Compounds 6 and 16 and enantiomer (R)-11 represent novel robust lead compounds of the IAS class.
Collapse
Affiliation(s)
- Valeria Famiglini
- Istituto Pasteur - Fondazione Cenci Bolognetti, Dipartimento di Chimica e Tecnologie del Farmaco, Sapienza Università di Roma , Piazzale Aldo Moro 5, I-00185 Roma, Italy
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
34
|
Drechsel NJD, Fennell CJ, Dill KA, Villà-Freixa J. TRIFORCE: Tessellated Semianalytical Solvent Exposed Surface Areas and Derivatives. J Chem Theory Comput 2014; 10:4121-4132. [PMID: 25221446 PMCID: PMC4159216 DOI: 10.1021/ct5002818] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2014] [Indexed: 12/01/2022]
Abstract
We present a new approach to the calculation of solvent-accessible surface areas of molecules with potential application to surface area based methods for determination of solvation free energies. As in traditional analytical and statistical approaches, this new algorithm, called TRIFORCE, reports both component areas and derivatives as a function of the atomic coordinates and radii. Unique to TRIFORCE are the rapid and scalable approaches for the determination of sphere intersection points and numerical estimation of the surface areas, derivatives, and other properties that can be associated with the surface area facets. The algorithm performs a special tessellation and semianalytical integration that uses a precomputed look-up table. This provides a simple way to balance numerical accuracy and memory usage. TRIFORCE calculates derivatives in the same manner, enabling application in force-dependent activities such as molecular geometry minimization. TRIFORCE is available free of charge for academic purposes as both a C++ library, which can be directly interfaced to existing molecular simulation packages, and a web-accessible application.
Collapse
Affiliation(s)
- Nils J. D. Drechsel
- Computational
Biochemistry and Biophysics Laboratory,
Research Unit on Biomedical Informatics, Universitat Pompeu Fabra, C/Doctor Aiguader, 88, 08003 Barcelona, Catalunya, Spain
- Laufer
Center for Physical and Quantitative Biology, Stony Brook University, Stony
Brook, New York 11794-5252, United States
- Department
of Chemistry, Oklahoma State University, Stillwater, Oklahoma 74078, United States
| | - Christopher J. Fennell
- Department
of Chemistry, Oklahoma State University, Stillwater, Oklahoma 74078, United States
| | - Ken A. Dill
- Laufer
Center for Physical and Quantitative Biology and Departments of Physics
and Chemistry, Stony Brook University, Stony Brook, New York 11794-5252, United States
| | - Jordi Villà-Freixa
- Computational
Biochemistry and Biophysics Laboratory,
Research Unit on Biomedical Informatics, Universitat Pompeu Fabra, C/Doctor Aiguader, 88, 08003 Barcelona, Catalunya, Spain
- Escola
Politècnica Superior, Universitat
de Vic—Universitat Central de Catalunya, C/de la Laura, 13, 08500 Vic, Catalunya, Spain
| |
Collapse
|
35
|
Li J, Mach P, Koehl P. Measuring the shapes of macromolecules - and why it matters. Comput Struct Biotechnol J 2013; 8:e201309001. [PMID: 24688748 PMCID: PMC3962087 DOI: 10.5936/csbj.201309001] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2013] [Revised: 11/22/2013] [Accepted: 11/22/2013] [Indexed: 11/22/2022] Open
Abstract
The molecular basis of life rests on the activity of biological macromolecules, mostly nucleic acids and proteins. A perhaps surprising finding that crystallized over the last handful of decades is that geometric reasoning plays a major role in our attempt to understand these activities. In this paper, we address this connection between geometry and biology, focusing on methods for measuring and characterizing the shapes of macromolecules. We briefly review existing numerical and analytical approaches that solve these problems. We cover in more details our own work in this field, focusing on the alpha shape theory as it provides a unifying mathematical framework that enable the analytical calculations of the surface area and volume of a macromolecule represented as a union of balls, the detection of pockets and cavities in the molecule, and the quantification of contacts between the atomic balls. We have shown that each of these quantities can be related to physical properties of the molecule under study and ultimately provides insight on its activity. We conclude with a brief description of new challenges for the alpha shape theory in modern structural biology.
Collapse
Affiliation(s)
- Jie Li
- Genome Center, University of California, Davis, 451 Health Sciences Drive, Davis, CA 95616, United States
| | - Paul Mach
- Graduate Group of Applied Mathematics, University of California, Davis, 1, Shields Ave, Davis, CA, 95616, United States
| | - Patrice Koehl
- Department of Computer Science and Genome Center, University of California, Davis, 1, Shields Ave, Davis, CA, 95616, United States
| |
Collapse
|
36
|
McCoy AJ, Nicholls RA, Schneider TR. SCEDS: protein fragments for molecular replacement in Phaser. ACTA CRYSTALLOGRAPHICA. SECTION D, BIOLOGICAL CRYSTALLOGRAPHY 2013; 69:2216-25. [PMID: 24189233 PMCID: PMC3817695 DOI: 10.1107/s0907444913021811] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/11/2013] [Accepted: 08/05/2013] [Indexed: 11/30/2022]
Abstract
A method is described for generating protein fragments suitable for use as molecular-replacement (MR) template models. The template model for a protein suspected to undergo a conformational change is perturbed along combinations of low-frequency normal modes of the elastic network model. The unperturbed structure is then compared with each perturbed structure in turn and the structurally invariant regions are identified by analysing the difference distance matrix. These fragments are scored with SCEDS, which is a combined measure of the sphericity of the fragments, the continuity of the fragments with respect to the polypeptide chain, the equality in number of atoms in the fragments and the density of C(α) atoms in the triaxial ellipsoid of the fragment extents. The fragment divisions with the highest SCEDS are then used as separate template models for MR. Test cases show that where the protein contains fragments that undergo a change in juxtaposition between template model and target, SCEDS can identify fragments that lead to a lower R factor after ten cycles of all-atom refinement with REFMAC5 than the original template structure. The method has been implemented in the software Phaser.
Collapse
Affiliation(s)
- Airlie J. McCoy
- Cambridge Institute for Medical Research, Department of Haematology, University of Cambridge, Hills Road, Cambridge CB2 0XY, England
| | - Robert A. Nicholls
- MRC Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge Biomedical Campus, Cambridge CB2 0QH, England
| | - Thomas R. Schneider
- European Molecular Biology Laboratory, Hamburg Unit c/o DESY, Notkestrasse 85, 22603 Hamburg, Germany
| |
Collapse
|
37
|
Sega M, Kantorovich SS, Jedlovszky P, Jorge M. The generalized identification of truly interfacial molecules (ITIM) algorithm for nonplanar interfaces. J Chem Phys 2013; 138:044110. [DOI: 10.1063/1.4776196] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
An analytical method for computing atomic contact areas in biomolecules. J Comput Chem 2012; 34:105-20. [DOI: 10.1002/jcc.23111] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2012] [Revised: 08/07/2012] [Indexed: 11/07/2022]
|
39
|
Rawat K, Bohidar HB. Universal Charge Quenching and Stability of Proteins in 1-Methyl-3-alkyl (Hexyl/Octyl) Imidazolium Chloride Ionic Liquid Solutions. J Phys Chem B 2012; 116:11065-74. [DOI: 10.1021/jp3049108] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Kamla Rawat
- Polymer and Biophysics Laboratory, School of Physical
Sciences, Jawaharlal Nehru University,
New Delhi 110067, India
| | - H. B. Bohidar
- Polymer and Biophysics Laboratory, School of Physical
Sciences, Jawaharlal Nehru University,
New Delhi 110067, India
| |
Collapse
|
40
|
Mach P, Koehl P. Geometric measures of large biomolecules: surface, volume, and pockets. J Comput Chem 2011; 32:3023-38. [PMID: 21823134 PMCID: PMC3188685 DOI: 10.1002/jcc.21884] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2011] [Accepted: 06/19/2011] [Indexed: 11/09/2022]
Abstract
Geometry plays a major role in our attempts to understand the activity of large molecules. For example, surface area and volume are used to quantify the interactions between these molecules and the water surrounding them in implicit solvent models. In addition, the detection of pockets serves as a starting point for predictive studies of biomolecule-ligand interactions. The alpha shape theory provides an exact and robust method for computing these geometric measures. Several implementations of this theory are currently available. We show however that these implementations fail on very large macromolecular systems. We show that these difficulties are not theoretical; rather, they are related to the architecture of current computers that rely on the use of cache memory to speed up calculation. By rewriting the algorithms that implement the different steps of the alpha shape theory such that we enforce locality, we show that we can remediate these cache problems; the corresponding code, UnionBall has an apparent O(n) behavior over a large range of values of n (up to tens of millions), where n is the number of atoms. As an example, it takes 136 sec with UnionBall to compute the contribution of each atom to the surface area and volume of a viral capsid with more than five million atoms on a commodity PC. UnionBall includes functions for computing analytically the surface area and volume of the intersection of two, three and four spheres that are fully detailed in an appendix. UnionBall is available as an OpenSource software.
Collapse
Affiliation(s)
- Paul Mach
- Graduate Group in Applied Mathematics, University of California, Davis, CA 95616
| | - Patrice Koehl
- Department of Computer Science and Genome Center, University of California, Davis, CA 95616
| |
Collapse
|
41
|
Kunz APE, Allison JR, Geerke DP, Horta BAC, Hünenberger PH, Riniker S, Schmid N, van Gunsteren WF. New functionalities in the GROMOS biomolecular simulation software. J Comput Chem 2011; 33:340-53. [PMID: 22076815 DOI: 10.1002/jcc.21954] [Citation(s) in RCA: 96] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2011] [Revised: 08/26/2011] [Accepted: 09/03/2011] [Indexed: 01/17/2023]
Abstract
Since the most recent description of the functionalities of the GROMOS software for biomolecular simulation in 2005 many new functions have been implemented. In this article, the new functionalities that involve modified forces in a molecular dynamics (MD) simulation are described: the treatment of electronic polarizability, an implicit surface area and internal volume solvation term to calculate interatomic forces, functions for the GROMOS coarse-grained supramolecular force field, a multiplicative switching function for nonbonded interactions, adiabatic decoupling of a number of degrees of freedom with temperature or force scaling to enhance sampling, and nonequilibrium MD to calculate the dielectric permittivity or viscosity. Examples that illustrate the use of these functionalities are given.
Collapse
Affiliation(s)
- Anna-Pitschna E Kunz
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology ETH, 8093 Zürich, Switzerland
| | | | | | | | | | | | | | | |
Collapse
|
42
|
Cai Q, Ye X, Wang J, Luo R. On-the-fly Numerical Surface Integration for Finite-Difference Poisson-Boltzmann Methods. J Chem Theory Comput 2011; 7:3608-3619. [PMID: 24772042 PMCID: PMC3998210 DOI: 10.1021/ct200389p] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Most implicit solvation models require the definition of a molecular surface as the interface that separates the solute in atomic detail from the solvent approximated as a continuous medium. Commonly used surface definitions include the solvent accessible surface (SAS), the solvent excluded surface (SES), and the van der Waals surface. In this study, we present an efficient numerical algorithm to compute the SES and SAS areas to facilitate the applications of finite-difference Poisson-Boltzmann methods in biomolecular simulations. Different from previous numerical approaches, our algorithm is physics-inspired and intimately coupled to the finite-difference Poisson-Boltzmann methods to fully take advantage of its existing data structures. Our analysis shows that the algorithm can achieve very good agreement with the analytical method in the calculation of the SES and SAS areas. Specifically, in our comprehensive test of 1,555 molecules, the average unsigned relative error is 0.27% in the SES area calculations and 1.05% in the SAS area calculations at the grid spacing of 1/2Å. In addition, a systematic correction analysis can be used to improve the accuracy for the coarse-grid SES area calculations, with the average unsigned relative error in the SES areas reduced to 0.13%. These validation studies indicate that the proposed algorithm can be applied to biomolecules over a broad range of sizes and structures. Finally, the numerical algorithm can also be adapted to evaluate the surface integral of either a vector field or a scalar field defined on the molecular surface for additional solvation energetics and force calculations.
Collapse
Affiliation(s)
- Qin Cai
- Department of Biomedical Engineering, University of California, Irvine, California
- Department of Molecular Biology and Biochemistry, University of California, Irvine, California
| | - Xiang Ye
- Department of Molecular Biology and Biochemistry, University of California, Irvine, California
- Department of Physics, Shanghai Normal University, Shanghai, China
| | - Jun Wang
- Department of Molecular Biology and Biochemistry, University of California, Irvine, California
| | - Ray Luo
- Department of Biomedical Engineering, University of California, Irvine, California
- Department of Molecular Biology and Biochemistry, University of California, Irvine, California
| |
Collapse
|
43
|
Kim DS, Kim CM, Won CI, Kim JK, Ryu J, Cho Y, Lee C, Bhak J. BetaDock: Shape-Priority Docking Method Based on Beta-Complex. J Biomol Struct Dyn 2011; 29:219-42. [DOI: 10.1080/07391102.2011.10507384] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
44
|
Klenin KV, Tristram F, Strunk T, Wenzel W. Derivatives of molecular surface area and volume: Simple and exact analytical formulas. J Comput Chem 2011; 32:2647-53. [DOI: 10.1002/jcc.21844] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2011] [Revised: 04/26/2011] [Accepted: 04/26/2011] [Indexed: 11/06/2022]
|
45
|
Allison JR, Boguslawski K, Fraternali F, van Gunsteren WF. A Refined, Efficient Mean Solvation Force Model that Includes the Interior Volume Contribution. J Phys Chem B 2011; 115:4547-57. [DOI: 10.1021/jp2017117] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jane R. Allison
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology ETH, 8093 Zürich, Switzerland
| | - Katharina Boguslawski
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology ETH, 8093 Zürich, Switzerland
| | - Franca Fraternali
- Randall Division of Cell and Molecular Biophysics, King’s College London, London SE1 1UL, United Kingdom
| | - Wilfred F. van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology ETH, 8093 Zürich, Switzerland
| |
Collapse
|
46
|
Núñez S, Venhorst J, Kruse CG. Assessment of a Novel Scoring Method Based on Solvent Accessible Surface Area Descriptors. J Chem Inf Model 2010; 50:480-6. [DOI: 10.1021/ci9004628] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Sara Núñez
- Research Laboratories, Solvay Pharmaceuticals, C. J. van Houtenlaan 36, 1381 CP Weesp, The Netherlands
| | - Jennifer Venhorst
- Research Laboratories, Solvay Pharmaceuticals, C. J. van Houtenlaan 36, 1381 CP Weesp, The Netherlands
| | - Chris G. Kruse
- Research Laboratories, Solvay Pharmaceuticals, C. J. van Houtenlaan 36, 1381 CP Weesp, The Netherlands
| |
Collapse
|
47
|
Reddy G, Yethiraj A. Solvent effects in polyelectrolyte adsorption: Computer simulations with explicit and implicit solvent. J Chem Phys 2010; 132:074903. [DOI: 10.1063/1.3319782] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
|
48
|
Pal A, Bahadur RP, Ray PS, Chakrabarti P. Accessibility and partner number of protein residues, their relationship and a webserver, ContPlot for their display. BMC Bioinformatics 2009; 10:103. [PMID: 19356223 PMCID: PMC2680847 DOI: 10.1186/1471-2105-10-103] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2008] [Accepted: 04/08/2009] [Indexed: 11/29/2022] Open
Abstract
Background Depending on chemical features residues have preferred locations – interior or exterior – in protein structures, which also determine how many other residues are found around them. The close packing of residues is the hallmark of protein interior and protein-protein interaction sites. Results The average values of accessible surface area (ASA) and partner number (PN, the number of other residues within a distance of 4.5 Å from any atom of a given residue) of different residues have been determined and a webserver, ContPlot has been designed to display these values (relative to the average values) along the protein sequence. This would be useful to visually identify residues that are densely packed, or those involved in protein-protein interactions. The skewness observed in the distribution of PNs is indicative of the hydrophobic or hydrophilic nature of the residue. The variation of ASA with PN can be analytically expressed in terms of a cubic equation. These equations (one for each residue) can be used to estimate the ASA of a polypeptide chain using the PNs of the individual residues in the structure. Conclusion The atom-based PNs (obtained by counting surrounding atoms) are highly correlated to the residue-based PN, indicating that the latter can adequately capture the atomic details of packing. The average values of ASA and PN associated with each residue should be useful in protein structure prediction or fold-recognition algorithm. ContPlot would provide a handy tool to assess the importance of a residue in the protein structure or interaction site.
Collapse
Affiliation(s)
- Arumay Pal
- Department of Biochemistry, Bose Institute, P-1/12 CIT Scheme VIIM, Kolkata, India.
| | | | | | | |
Collapse
|
49
|
Ezkurdia I, Bartoli L, Fariselli P, Casadio R, Valencia A, Tress ML. Progress and challenges in predicting protein-protein interaction sites. Brief Bioinform 2009; 10:233-46. [PMID: 19346321 DOI: 10.1093/bib/bbp021] [Citation(s) in RCA: 120] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The identification of protein-protein interaction sites is an essential intermediate step for mutant design and the prediction of protein networks. In recent years a significant number of methods have been developed to predict these interface residues and here we review the current status of the field. Progress in this area requires a clear view of the methodology applied, the data sets used for training and testing the systems, and the evaluation procedures. We have analysed the impact of a representative set of features and algorithms and highlighted the problems inherent in generating reliable protein data sets and in the posterior analysis of the results. Although it is clear that there have been some improvements in methods for predicting interacting sites, several major bottlenecks remain. Proteins in complexes are still under-represented in the structural databases and in particular many proteins involved in transient complexes are still to be crystallized. We provide suggestions for effective feature selection, and make it clear that community standards for testing, training and performance measures are necessary for progress in the field.
Collapse
Affiliation(s)
- Iakes Ezkurdia
- Centro Nacional de Biotechnolgia, Spanish National Cancer Research Centre (CNIO), Madrid, Spain
| | | | | | | | | | | |
Collapse
|
50
|
Durham E, Dorr B, Woetzel N, Staritzbichler R, Meiler J. Solvent accessible surface area approximations for rapid and accurate protein structure prediction. J Mol Model 2009; 15:1093-108. [PMID: 19234730 PMCID: PMC2712621 DOI: 10.1007/s00894-009-0454-9] [Citation(s) in RCA: 200] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2008] [Accepted: 01/02/2009] [Indexed: 12/01/2022]
Abstract
The burial of hydrophobic amino acids in the protein core is a driving force in protein folding. The extent to which an amino acid interacts with the solvent and the protein core is naturally proportional to the surface area exposed to these environments. However, an accurate calculation of the solvent-accessible surface area (SASA), a geometric measure of this exposure, is numerically demanding as it is not pair-wise decomposable. Furthermore, it depends on a full-atom representation of the molecule. This manuscript introduces a series of four SASA approximations of increasing computational complexity and accuracy as well as knowledge-based environment free energy potentials based on these SASA approximations. Their ability to distinguish correctly from incorrectly folded protein models is assessed to balance speed and accuracy for protein structure prediction. We find the newly developed “Neighbor Vector” algorithm provides the most optimal balance of accurate yet rapid exposure measures.
Collapse
Affiliation(s)
- Elizabeth Durham
- Department of Chemistry, Center for Structural Biology, Vanderbilt University, 465 21st Ave South, Nashville, TN 37232-8725, USA
| | | | | | | | | |
Collapse
|