force fields
Recently Published Documents





2022 ◽  
Vol 72 ◽  
pp. 187-193
Xibing He ◽  
Brandon Walker ◽  
Viet H. Man ◽  
Pengyu Ren ◽  
Junmei Wang

2022 ◽  
Vol 156 (2) ◽  
pp. 024801
Charly Empereur-mot ◽  
Riccardo Capelli ◽  
Mattia Perrone ◽  
Cristina Caruso ◽  
Giovanni Doni ◽  

2022 ◽  
Petr Jurecka ◽  
Marie Zgarbova ◽  
Filip Cerny ◽  
Jan Salomon

When DNA interacts with a protein, its structure often undergoes significant conformational adaptation. Perhaps the most common is the transition from canonical B-DNA towards the A-DNA form, which is not a two-state, but rather a continuous transition. The A- and B- forms differ mainly in sugar pucker P (north/south) and glycosidic torsion χ (high-anti/anti). The combination of A-like P and B-like χ (and vice versa) represents the nature of the intermediate states lying between the pure A- and B- forms. In this work, we study how the A/B equilibrium and in particular the A/B intermediate states, which are known to be over-represented at protein-DNA interfaces, are modeled by current AMBER force fields. Eight protein-DNA complexes and their naked (unbound) DNAs were simulated with OL15 and bsc1 force fields as well as an experimental combination OL15χOL3. We found that while the geometries of the A-like intermediate states in the molecular dynamics (MD) simulations agree well with the native X-ray geometries found in the protein-DNA complexes, their populations (stabilities) are significantly underestimated. Different force fields predict different propensities for A-like states growing in the order OL15 < bsc1 < OL15χOL3, but the overall populations of the A-like form are too low in all of them. Interestingly, the force fields seem to predict the correct sequence-dependent A-form propensity, as they predict larger populations of the A-like form in naked (unbound) DNA in those steps that acquire A-like conformations in protein-DNA complexes. The instability of A-like geometries in current force fields may significantly alter the geometry of the simulated protein-DNA complex, destabilize the binding motif, and reduce the binding energy, suggesting that refinement is needed to improve description of protein-DNA interactions in AMBER force fields.

Polymers ◽  
2022 ◽  
Vol 14 (2) ◽  
pp. 252
Natalia Lukasheva ◽  
Dmitry Tolmachev ◽  
Hector Martinez-Seara ◽  
Mikko Karttunen

Electrostatic interactions have a determining role in the conformational and dynamic behavior of polyelectrolyte molecules. In this study, anionic polyelectrolyte molecules, poly(glutamic acid) (PGA) and poly(aspartic acid) (PASA), in a water solution with the most commonly used K+ or Na+ counterions, were investigated using atomistic molecular dynamics (MD) simulations. We performed a comparison of seven popular force fields, namely AMBER99SB-ILDN, AMBER14SB, AMBER-FB15, CHARMM22*, CHARMM27, CHARMM36m and OPLS-AA/L, both with their native parameters and using two common corrections for overbinding of ions, the non-bonded fix (NBFIX), and electronic continuum corrections (ECC). These corrections were originally introduced to correct for the often-reported problem concerning the overbinding of ions to the charged groups of polyelectrolytes. In this work, a comparison of the simulation results with existing experimental data revealed several differences between the investigated force fields. The data from these simulations and comparisons with previous experimental data were then used to determine the limitations and strengths of these force fields in the context of the structural and dynamic properties of anionic polyamino acids. Physical properties, such as molecular sizes, local structure, and dynamics, were studied using two types of common counterions, namely potassium and sodium. The results show that, in some cases, both the macroion size and dynamics depend strongly on the models (parameters) for the counterions due to strong overbinding of the ions and charged side chain groups. The local structures and dynamics are more sensitive to dihedral angle parameterization, resulting in a preference for defined monomer conformations and the type of correction used. We also provide recommendations based on the results.

2022 ◽  
Shuai Jiang ◽  
Yi-Rong Liu ◽  
Teng Huang ◽  
Ya-Juan Feng ◽  
Chun-Yu Wang ◽  

Abstract Atmospheric aerosol nucleation contributes to around half of cloud condensation nuclei globally. Despite the importance for climate, detailed nucleation mechanisms are still poorly understood. Understanding aerosol nucleation dynamics is hindered by non-reactivity of force fields and high computational costs due to rare event nature of aerosol nucleation. Developing reactive force fields for nucleation systems are even more challenging than covalently bonded materials because of wide size range and high dimensional characteristics of non-covalent hydrogen bonding bridging clusters. Here we proposes a system transferable framework to train an accurate reactive force field (FF) based on deep neural network (DNN) and further bridges the DNN-FF based molecular dynamics (MD) with cluster kinetics model based on Poisson distributions of reactive events to overcome high computational costs from direct MD. We found that previously reported acid-base formation rates tend to be underestimated several times, emphasizing acid-base nucleation observed in multiple environments should be revisited.

Frank Jensen

We review different models for introducing electrical polarization in force fields, with special focus on methods where polarization is modelled at the atomic charge level. While electric polarization has been...

Sign in / Sign up

Export Citation Format

Share Document