During the second half of the 20th century, following
structural
biology hallmark works on DNA and proteins, biochemists shifted their
questions from “what does this molecule look like?”
to “how does this process work?”. Prompted by the theoretical
and practical developments in computational chemistry, this led to
the emergence of biomolecular simulations and, along with the 2013
Nobel Prize in Chemistry, to the development of hybrid QM/MM methods.
QM/MM methods are necessary whenever the problem we want to address
involves chemical reactivity and/or a change in the system’s
electronic structure, with archetypal examples being the studies of
an enzyme’s reaction mechanism and a metalloprotein’s
active site. In the last decades QM/MM methods have seen an increasing
adoption driven by their incorporation in widely used biomolecular
simulation software. However, properly setting up a QM/MM simulation
is not an easy task, and several issues need to be properly addressed
to obtain meaningful results. In the present work, we describe both
the theoretical concepts and practical issues that need to be considered
when performing QM/MM simulations. We start with a brief historical
perspective on the development of these methods and describe when
and why QM/MM methods are mandatory. Then we show how to properly
select and analyze the performance of the QM level of theory, the
QM system size, and the position and type of the boundaries. We show
the relevance of performing prior QM model system (or QM cluster)
calculations in a vacuum and how to use the corresponding results
to adequately calibrate those derived from QM/MM. We also discuss
how to prepare the starting structure and how to select an adequate
simulation strategy, including those based on geometry optimizations
as well as free energy methods. In particular, we focus on the determination
of free energy profiles using multiple steered molecular dynamics
(MSMD) combined with Jarzynski’s equation. Finally, we describe
the results for two illustrative and complementary examples: the reaction
performed by chorismate mutase and the study of ligand binding to
hemoglobins. Overall, we provide many practical recommendations (or
shortcuts) together with important conceptualizations that we hope
will encourage more and more researchers to incorporate QM/MM studies
into their research projects.