66 research outputs found
Constraint methods for determining pathways and free energy of activated processes
Activated processes from chemical reactions up to conformational transitions
of large biomolecules are hampered by barriers which are overcome only by the
input of some free energy of activation. Hence, the characteristic and
rate-determining barrier regions are not sufficiently sampled by usual
simulation techniques. Constraints on a reaction coordinate r have turned out
to be a suitable means to explore difficult pathways without changing potential
function, energy or temperature. For a dense sequence of values of r, the
corresponding sequence of simulations provides a pathway for the process. As
only one coordinate among thousands is fixed during each simulation, the
pathway essentially reflects the system's internal dynamics. From mean forces
the free energy profile can be calculated to obtain reaction rates and insight
in the reaction mechanism. In the last decade, theoretical tools and computing
capacity have been developed to a degree where simulations give impressive
qualitative insight in the processes at quantitative agreement with
experiments. Here, we give an introduction to reaction pathways and
coordinates, and develop the theory of free energy as the potential of mean
force. We clarify the connection between mean force and constraint force which
is the central quantity evaluated, and discuss the mass metric tensor
correction. Well-behaved coordinates without tensor correction are considered.
We discuss the theoretical background and practical implementation on the
example of the reaction coordinate of targeted molecular dynamics simulation.
Finally, we compare applications of constraint methods and other techniques
developed for the same purpose, and discuss the limits of the approach
The performance of non-polarizable and polarizable force-field parameter sets for ethylene glycol in molecular dynamics simulations of the pure liquid and its aqueous mixtures
Comprehensive and Automated Linear Interaction Energy Based Binding-Affinity Prediction for Multifarious Cytochrome P450 Aromatase Inhibitors
Cytochrome P450 aromatase (CYP19A1) plays a key role in the development of estrogen dependent breast cancer, and aromatase inhibitors have been at the front line of treatment for the past three decades. The development of potent, selective and safer inhibitors is ongoing with in silico screening methods playing a more prominent role in the search for promising lead compounds in bioactivity-relevant chemical space. Here we present a set of comprehensive binding affinity prediction models for CYP19A1 using our automated Linear Interaction Energy (LIE) based workflow on a set of 132 putative and structurally diverse aromatase inhibitors obtained from a typical industrial screening study. We extended the workflow with machine learning methods to automatically cluster training and test compounds in order to maximize the number of explained compounds in one or more predictive LIE models. The method uses protein-ligand interaction profiles obtained from Molecular Dynamics (MD) trajectories to help model search and define the applicability domain of the resolved models. Our method was successful in accounting for 86% of the data set in 3 robust models that show high correlation between calculated and observed values for ligand-binding free energies (RMSE < 2.5 kJ mol(-1)), with good cross-validation statistics
Insights into regioselective metabolism of mefenamic acid by cytochrome P450 BM3 mutants through crystallography, docking, molecular dynamics, and free energy calculations
Biomolecular modeling: Goals, problems, perspectives
Computation based on molecular models is playing an increasingly important role in biology, biological chemistry, and biophysics. Since only a very limited number of properties of biomolecular systems is actually accessible to measurement by experimental means, computer simulation can complement experiment by providing not only averages, but also distributions and time series of any definable quantity, for example, conformational distributions or interactions between parts of systems. Present day biomolecular modeling is limited in its application by four main problems: 1) the force-field problem, 2) the search (sampling) problem, 3) the ensemble (sampling) problem, and 4) the experimental problem. These four problems are discussed and illustrated by practical examples. Perspectives are also outlined for pushing forward the limitations of biomolecular modeling. © 2006 Wiley-VCH Verlag GmbH & Co. KGaA
QM/MM based fitting of atomic polarizabilities for use in condensed-phase biomolecular simulation
Accounting for electronic polarization effects in biomolecular simulation (by using a polarizable force field) can increase the accuracy of simulation results. However, the use of gas-phase estimates of atomic polarizabilities
Semiempirical self-consistent polarization description of bulk water, the liquid-vapor interface, and cubic ice
We have applied an efficient electronic structure approach, the semiempirical self-consistent polarization neglect of diatomic differential overlap (SCP-NDDO) method, previously parametrized to reproduce properties of water clusters by Chang, Schenter, and Garrett J. Chem. Phys. 2008, 128, 164111] and now implemented in the CP2K package, to model ambient liquid water at 300 K (both the bulk and the liquid-vapor interface) and cubic ice at 15 and 250 K The SCP-NDDO potential retains its transferability and good performance across the full range of conditions encountered in the clusters and the bulk phases of water. In particular, we obtain good results for the density, radial distribution functions, enthalpy of vaporization, self-diffusion coefficient, molecular dipole moment distribution, and hydrogen bond populations, in comparison to experimental measurements
- …