Molecular Dynamics and Quantum Chemical Approach for the Estimation of an Intramolecular Hydrogen Bond Strength in Okadaic Acid

Studying the ssDNA loaded adeno-associated virus aggregation using coarse-grained molecular dynamics simulations

The aggregation of adeno-associated viral (AAV) capsids in an aqueous environment was investigated via coarse-grained molecular dynamics (CG-MD) simulations. The primary driving force and mechanism of the aggregation were investigated with or without single-strand DNA (ssDNA) loaded at various process temperatures. Capsid aggregation appeared to involve multiple residue interactions (i.e., hydrophobic, polar and charged residues) leading to complex protein aggregation. In addition, two aggregation mechanisms (i.e., the fivefold face-to-face contact and the edge-to-edge contact) were identified from this study. The ssDNA with its asymmetric structure could be the reason for destabilizing protein subunits and enhancing the interaction between the charged residues, and further result in the non-reversible face-to-face contact. At higher temperature, the capsid structure was found to be unstable with the significant size expansion of the loaded ssDNA which could be attributed to reduced number of intramolecular hydrogen bonds, the increased conformational deviations of protein subunits and the higher residue fluctuations. The CG-MD model was further validated with previous experimental and simulation data, including the full capsid size measurement and the capsid internal pressure. Thus, a good understanding of AAV capsid aggregation, instability and the role of ssDNA were revealed by applying the developed computational model.

From bonds to interactions: comprehensive molecular characterization via polarizable bond-dipole approach

Accurately characterizing molecular interactions stands as a pivotal requirement for ensuring the reliability of molecular dynamics simulations. In line with our bond-dipole-based interaction model proposed by Gao et al. [X.-C. Gao, Q. Hao and C.-S. Wang, J. Chem. Theory Comput., 2017, 13, 2730-2741.], we have implemented an efficient and concise approach to compute electrostatic potential. This methodology capitalizes on the polarizable nature of chemical bond dipoles, resulting in a model of remarkable simplicity. In this study, we have revised the polarizable bond-dipole-based force field (PBFF) through the meticulous curation of quantum chemical data sets. These data sets encompass a comprehensive collection of 40 000 conformations, including those of water, methylamine, methanol, and N-methylacetamide. Additionally, we incorporate 520 hydrogen-bonded dimers into our data sets. In pursuit of enhanced accuracy in molecular dynamics simulations and a more faithful representation of potential energy landscapes, we undertook the re-optimization of the nonbonded parameters within the PBFF framework. Concurrently, we intricately fine-tuned the bonded parameters. The results of our comprehensive evaluation denote that this newly optimized force field method adeptly and efficiently computes structural characteristics, harmonic frequencies, and interaction energies. Overall, this study provides further validation for the applicability of PBFF in molecular dynamics simulations.

Molecular Dynamics Simulation to Uncover the Mechanisms of Protein Instability During Freezing

Freezing is a common process applied in the pharmaceutical industry to store and transport biotherapeutics. Herewith, multi-scale molecular dynamics simulations of Lactate dehydrogenase (LDH) protein in phosphate buffer with/without ice formation performed to uncover the still poorly understood mechanisms and molecular details of protein destabilization upon freezing. Both fast and slow ice growing conditions were simulated at 243 K from one or two-side of the simulation box, respectively. The rate of ice formation at all-atom simulations was crucial to LDH stability, as faster freezing rates resulted in enhanced structural stability maintained by a higher number of intramolecular hydrogen bonds, less flexible protein's residues, lower solvent accessibility and greater structural compactness. Further, protein aggregation investigated by coarse-grained simulations was verified to be initiated by extended protein structures and retained by electrostatic interactions of the salt bridges between charged residues and hydrogen bonds between polar residues of the protein. Lastly, the study of free energy of dissociation through steered molecular dynamics simulation revealed LDH was destabilized by the solvation of the hydrophobic core and the loss of hydrophobic interactions. For the first time, experimentally validated molecular simulations revealed the detailed mechanisms of LDH destabilization upon ice formation and cryoconcentration of solutes.