Calculation of Protein Folding Thermodynamics Using Molecular Dynamics Simulations

Ab initio characterization of protein molecular dynamics with AI2BMD

Biomolecular dynamics simulation is a fundamental technology for life sciences research, and its usefulness depends on its accuracy and efficiency1-3. Classical molecular dynamics simulation is fast but lacks chemical accuracy4,5. Quantum chemistry methods such as density functional theory can reach chemical accuracy but cannot scale to support large biomolecules6. Here we introduce an artificial intelligence-based ab initio biomolecular dynamics system (AI2BMD) that can efficiently simulate full-atom large biomolecules with ab initio accuracy. AI2BMD uses a protein fragmentation scheme and a machine learning force field7 to achieve generalizable ab initio accuracy for energy and force calculations for various proteins comprising more than 10,000 atoms. Compared to density functional theory, it reduces the computational time by several orders of magnitude. With several hundred nanoseconds of dynamics simulations, AI2BMD demonstrated its ability to efficiently explore the conformational space of peptides and proteins, deriving accurate 3J couplings that match nuclear magnetic resonance experiments, and showing protein folding and unfolding processes. Furthermore, AI2BMD enables precise free-energy calculations for protein folding, and the estimated thermodynamic properties are well aligned with experiments. AI2BMD could potentially complement wet-lab experiments, detect the dynamic processes of bioactivities and enable biomedical research that is impossible to conduct at present.

Probing the Molecular Interactions of A22 with Prokaryotic Actin MreB and Eukaryotic Actin: A Computational and Experimental Study

Actin is a major cytoskeletal system that mediates the intricate organization of macromolecules within cells. The bacterial cytoskeletal protein MreB is a prokaryotic actin-like protein governing the cell shape and intracellular organization in many rod-shaped bacteria, including pathogens. MreB stands as a target for antibiotic development, and compounds like A22 and its analogue, MP265, are identified as potent inhibitors of MreB. The bacterial actin MreB shares structural homology with eukaryotic actin despite lacking sequence similarity. It is currently not clear whether small molecules that inhibit MreB can act on eukaryotic actin due to their structural similarity. In this study, we investigate the molecular interactions between A22 and its analogue MP265 with MreB and eukaryotic actin through a molecular dynamics approach. Employing MD simulations and free energy calculations with an all-atom model, we unveil the robust interaction of A22 and MP265 with MreB, and substantial binding affinity is observed for A22 and MP265 with eukaryotic actin. Experimental assays reveal A22's toxicity to eukaryotic cells, including yeast and human glioblastoma cells. Microscopy analysis demonstrates the profound effects of A22 on actin organization in human glioblastoma cells. This integrative computational and experimental study provides new insights into A22's mode of action, highlighting its potential as a versatile tool for probing the dynamics of both prokaryotic and eukaryotic actins.

Energy, water, and protein folding: A molecular dynamics-based quantitative inventory of molecular interactions and forces that make proteins stable

Protein folding energetics can be determined experimentally on a case-by-case basis but it is not understood in sufficient detail to provide deep control in protein design. The fundamentals of protein stability have been outlined by calorimetry, protein engineering, and biophysical modeling, but these approaches still face great difficulty in elucidating the specific contributions of the intervening molecules and physical interactions. Recently, we have shown that the enthalpy and heat capacity changes associated to the protein folding reaction can be calculated within experimental error using molecular dynamics simulations of native protein structures and their corresponding unfolded ensembles. Analyzing in depth molecular dynamics simulations of four model proteins (CI2, barnase, SNase, and apoflavodoxin), we dissect here the energy contributions to ΔH (a key component of protein stability) made by the molecular players (polypeptide and solvent molecules) and physical interactions (electrostatic, van der Waals, and bonded) involved. Although the proteins analyzed differ in length, isoelectric point and fold class, their folding energetics is governed by the same quantitative pattern. Relative to the unfolded ensemble, the native conformations are enthalpically stabilized by comparable contributions from protein-protein and solvent-solvent interactions, and almost equally destabilized by interactions between protein and solvent molecules. The native protein surface seems to interact better with water than the unfolded one, but this is outweighed by the unfolded surface being larger. From the perspective of physical interactions, the native conformations are stabilized by van de Waals and Coulomb interactions and destabilized by conformational strain arising from bonded interactions. Also common to the four proteins, the sign of the heat capacity change is set by interactions between protein and solvent molecules or, from the alternative perspective, by Coulomb interactions.