A systematic method to calculate anharmonic force constants of crystals is presented. The method employs the direct-method approach, where anharmonic force constants are extracted from the trajectory ...of first-principles molecular dynamics simulations at high temperature. The method is applied to Si where accurate cubic and quartic force constants are obtained. We observe that higher-order correction is crucial to obtain accurate force constants from the trajectory with large atomic displacements. The calculated harmonic and anharmonic force constants are, then, combined with the Boltzmann transport equation (BTE) and non-equilibrium molecular dynamics (NEMD) methods in calculating the thermal conductivity. The BTE approach successfully predicts the lattice thermal conductivity of bulk Si, whereas NEMD shows considerable underestimates. To evaluate the linear extrapolation method employed in NEMD to estimate bulk values, we analyze the size dependence in NEMD based on BTE calculations. We observe strong nonlinearity in the size dependence of NEMD in Si, which can be ascribed to acoustic phonons having long mean-free-paths and carrying considerable heat. Subsequently, we also apply the whole method to a thermoelectric material Mg2Si and demonstrate the reliability of the NEMD method for systems with low thermal conductivities.
The radiation intensity of the outer shell of a carbon onion complex containing an ionized iron atom is calculated according to the scheme proposed by L. D. Landau. A complete description of the ...oscillations arising in the system under consideration is given. Keywords: molecular dynamics, dipole radiation, mathematical modeling.
This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular ...dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.
Cancer cells can evade immune surveillance through binding of its transmembrane receptor CD47 to CD172a on myeloid cells. CD47 is recognized as a promising immune checkpoint for cancer immunotherapy ...inhibiting macrophage phagocytosis. N-terminal post-translated modification (PTM) via glutaminyl cyclase is a landmark event in CD47 function maturation, but the molecular mechanism underlying the mechano-chemical regulation of the modification on CD47/CD172a remains unclear. Here, we performed so-called "ramp-clamp" steered molecular dynamics (SMD) simulations, and found that the N-terminal PTM enhanced interaction of CD172a with CD47 by inducing a dynamics-driven contraction of the binding pocket of the bound CD172a, an additional constraint on CYS15 on CD47 significantly improved the tensile strength of the complex with or without PTM, and a catch bond phenomenon would occur in complex dissociation under tensile force of 25 pN in a PTM-independent manner too. The residues GLN52 and SER66 on CD172a reinforced the H-bonding with their partners on CD47 in responding to PTM, while ARG69 on CD172 with its partner on CD47 might be crucial in the structural stability of the complex. This work might serve as molecular basis for the PTM-induced function improvement of CD47, should be helpful for deeply understanding CD47-relevant immune response and cancer development, and provides a novel insight in developing of new strategies of immunotherapy targeting this molecule interaction.
Display omitted
•Study emphasizes (poly)valine vibrational effects and solvent impact on absorption spectrum and atomic magnetic properties.•Distinctive Magnetic Signatures for carbon atoms in the ...peptide bond exhibit specific 6.2 ± 6.1 ppm magnetic shielding constant.•Peptide chain growth induces a 60 nm red shift in absorption band spectrum.
In this study, we explored how the growth of peptide chains containing valines impacts nuclear magnetic resonance (NMR) and optical absorption spectra. Results stem from equilibrated configurations in water, reflecting solvent influence and valine vibration. Utilizing fully atomistic molecular dynamics (MD), simulations, and density functional theory (DFT) including TD-DFT, we observe a notable shift in optical absorption towards longer wavelengths with increased valine count. Additionally, carbon atom magnetic signatures unveil distinct peptide bond characteristics. This work sheds light on the spectroscopic behavior of valine-rich peptide chains, crucial for understanding structural and functional properties, particularly in biomolecular contexts.
The overexpression of p21-activated kinase 4 (PAK4) is associated with a variety of cancers. In this paper, the binding modes and inhibitory mechanisms of four 7H-pyrrolo2,3-dpyrimidine competitive ...inhibitors of PAK4 were investigated at the molecular level, mainly using molecular dynamics simulations and binding free energy calculations. The results show that the inhibitors had strong interactions with the hinge region, the β-sheets, and the residues with charged side chains around the 4-substituent. The terminal amino group of the inhibitor 5n was different from the other three, which could cause the enhancement of hydrogen bonds or electrostatic interactions formed with the surrounding residues. Thus, inhibitor 5n had the strongest inhibition capacity. The different halogen atoms on the 2-substituents of the inhibitors 5h, 5g, and 5e caused differences in the positions of the 2-benzene rings and affected the interactions of the hinge region. It also affected to some extent the orientations of the 4-imino groups and consequently their affinities for the surrounding charged residues. The combined results lead to the weakest inhibitory capacity of inhibitor 5e.
PDF
