An official website of the United States government
The pKa values and associated protonation states of ionizable lipids in lipid nanoparticle (LNP) formulations are strongly dependent on their chemical environment. This phenomenon leads to poorly understood structure-function relationships that influence payload delivery, tissue-selective biodistribution, and manufacturing. For example, the charge- and biodistribution of an mRNA-loaded LNP can vary based on the type of ionizable lipid used, the molar ratio of its components, and its cargo. Yet, the spatial resolution of experimental protonation state measurements is currently limited to the apparent charge of an ionizable lipid averaged over all environments/conformations of an LNP — best represented by its apparent pKa value. Such measurements are too coarse to capture the heterogenous charge distributions of ionizable lipids in LNPs, which influence biocorona formation and interactions with the payload. Similar limitations are inherent to classical fixed protonation-state in silico models that cannot capture the environment-dependent protonation states and pKa values determining local pKa. To address this gap in experimental and computational tools available to accurately determine the local charge distributions in LNPs, this work now incorporates a scalable continuous constant pH molecular dynamics (CpHMD) model to simulate the self-assembly processes of five reported distinct LNP formulations. Parameters for ionizable lipids were generated from fitting fixed lambda-state calculations performed with Hamiltonian replica exchange (HREX) to improve conformational sampling during parameterization. Simulated systems were composed of 100 ionizable lipids (50 mol%), cholesterol (40 mol%), distearoylphosphatidylcholine (10 mol%), and mRNA (20 nucleotides) to model the interior of an LNP. Self-assembly was simulated for 100 ns at different pH values to validate the apparent pKa for each system. The theoretically calculated apparent pKa values matched reasonably well with those measured experimentally (mean absolute error = 0.5 pKa units), and all systems exhibited pH-dependent structures. Overall, this work provides a new computational platform technology to (i) predict the pKa values of ionizable lipids in different chemical environments and (ii) enable a structure-based way to model the heterogeneous, environment-dependent charge distributions of ionizable lipids in LNP systems typically encountered during LNP manufacturing and delivery.
more »
« less
pKa Determination of a Histidine Residue in a Short Peptide Using Raman Spectroscopy
Determining the pKa of key functional groups is critical to understanding the pH-dependent behavior of biological proteins and peptide-based biomaterials. Traditionally, 1H NMR spectroscopy has been used to determine the pKa of amino acids; however, for larger molecules and aggregating systems, this method can be practically impossible. Previous studies concluded that the C-D stretches in Raman are a useful alternative for determining the pKa of histidine residues. In this study, we report on the Raman application of the C2-D probe on histidine’s imidazole side chain to determining the pKa of histidine in a short peptide sequence. The pKa of the tripeptide was found via difference Raman spectroscopy to be 6.82, and this value was independently confirmed via 1H NMR spectroscopy on the same peptide. The C2-D probe was also compared to other Raman reporters of the protonation state of histidine and was determined to be more sensitive and reliable than other protonation-dependent signals. The C2-D Raman probe expands the tool box available to chemists interested in directly interrogating the pKa’s of histidine-containing peptide and protein systems.
more »
« less
- Award ID(s):
- 1800080
- PAR ID:
- 10120833
- Date Published:
- Journal Name:
- Molecules
- Volume:
- 24
- Issue:
- 3
- ISSN:
- 1420-3049
- Page Range / eLocation ID:
- 405
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Density functional theory (DFT) benchmark studies of 1H and 13C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured 1H and 13C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for 1H and 13C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for 1H and 0.5 to 2.9 for 13C. Maximum deviations were less than 0.5 and 6.5 ppm for 1H and 13C, respectively.more » « less
-
Block polymers show promise as solid-state battery electrolytes due to the optimization of conductive and mechanical properties enabled via tuning of block chemistry and length. We investigate a polystyrene-block-poly(oligo-oxyethylene methacrylate) (PS-b-POEM) electrolyte doped with various lithium salts to investigate the role of molecular structure on ion transport properties and on local ion dynamics and associations. Anion charge becomes more delocalized with increasing size, reducing the coupling between salt ions while increasing coupling between ion and polymer chain motions and creating a more mobile overall environment. We observe support for this ion-polymer coupling via 1H, 7Li and 19F NMR spectroscopy, from which we obtain ion-specific mobility transition temperatures that differ from the polymer glass transition temperature. We also note faster transport and weaker local energetic interactions with anion size using temperature-dependent NMR diffusometry. 1H NMR spectroscopy further elucidates polymer chain dynamics and enables quantification of the temperature-dependent fraction of the conducting block that is immobile near the PS-POEM domain interface. NMR thus represents a species-specific and timescale-specific platform to quantify phase and interface behavior, and to correlate ion-specific transport with polymer chain dynamics.more » « less
-
Chain entanglements play a crucial role in polymer crystallization, yet their effects on crystallization remain not fully understood. Freeze-drying is one way to potentially preserve disentangled states of long polymer chains. In fact, it is known that freeze-drying (FD) significantly accelerates the crystallization kinetics of semicrystalline polymers. However, the chain-level structure of the FD polymer chains without a long-range order (glass) has been a debatable matter. In this study, we investigate the effect of freeze-drying on single chain-level structures of 13CH3 enriched poly(L-lactic Acid) and 13CH enriched poly(D-lactic acid) racemate by using 1H-1H spin diffusion via 13C detection solid-state NMR spectroscopy. Spatial distributions of PLLA and PDLA glassy chains in the range of a few Å – 30 nm are evaluated via 1H-1H spin diffusion. This analysis provides core-shell morphology of single chains where the outer shell layers include both PDLA and PLLA mixture and the inner core possess a single component.more » « less
-
NMR spectroscopy is the most important technique for understanding the structure of peptides and proteins in solution, providing information at the single-residue and single-atom level. However, written instruction in the interpretation of NMR spectra of peptides and proteins is generally focused on advanced techniques and highly complex spectra, with a lack of simple spectra and guides available for beginning students. In order to address this instructional limitation, we have generated a dataset of 1H NMR spectra of a series of simple peptides that include all canonical amino acids. Peptides examined include Ac-X(S/pS)-NH2, Ac-X(T/pT)-NH2, and Ac-XPPGY-NH2, where X = all encoded amino acids, pS = phosphorylated Ser, and pT = phosphorylated Thr. The characterization of each peptide includes a 1-D spectrum and a TOCSY spectrum, with both the raw and processed available. The spectra can be used for instructional applications including analysis of regions of the spectra (e.g. amide HN, aromatic, Hα, and aliphatic regions); identification of spin systems and residue assignment via TOCSY spectra; analysis of conformational features including amide HN chemical shift dispersion and changes due to hydrogen bonding or post-translational modifications; the 3JαN coupling constant that reports on the φ torsion angle and on order versus disorder at a given residue; conformational preferences at Hα via chemical shift index analysis; understanding of diastereotopic hydrogens; dynamic processes, including hydrogen exchange; and identification of proline cis-trans isomerism. In addition, for a limited number of peptides, NOESY spectra are included to allow sequential resonance assignment and for assignment of trans versus cis proline conformations. Spectra from closely related peptides allow the analysis of the relative effects of single amino acid changes. The paper is written to be directly accessible to students as a tutorial guide. In addition, the data can be used by instructors for problem sets and exams.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/molecules24030405
