Preferred WMSA catalytic mechanism of the nucleotidyl transfer reaction in human DNA polymerase κ elucidates error-free bypass of a bulky DNA lesion

Effects of the Y432S Cancer-Associated Variant on the Reaction Mechanism of Human DNA Polymerase κ

Human DNA polymerases are vital for genetic information management. Their function involves catalyzing the synthesis of DNA strands with unparalleled accuracy, which ensures the fidelity and stability of the human genomic blueprint. Several disease-associated mutations and their functional impact on DNA polymerases have been reported. One particular polymerase, human DNA polymerase kappa (Pol κ), has been reported to be susceptible to several cancer-associated mutations. The Y432S mutation in Pol κ, associated with various cancers, is of interest due to its impact on polymerization activity and markedly reduced thermal stability. Here, we have used computational simulations to investigate the functional consequences of the Y432S using classical molecular dynamics (MD) and coupled quantum mechanics/molecular mechanics (QM/MM) methods. Our findings suggest that Y432S induces structural alterations in domains responsible for nucleotide addition and ternary complex stabilization while retaining structural features consistent with possible catalysis in the active site. Calculations of the minimum energy path associated with the reaction mechanism of the wild type (WT) and Y432S Pol κ indicate that, while both enzymes are catalytically competent (in terms of energetics and the active site's geometries), the cancer mutation results in an endoergic reaction and an increase in the catalytic barrier. Interactions with a third magnesium ion and environmental effects on nonbonded interactions, particularly involving key residues, contribute to the kinetic and thermodynamic distinctions between the WT and mutant during the catalytic reaction. The energetics and electronic findings suggest that active site residues favor the catalytic reaction with dCTP3- over dCTP4-.

Probing the general base for DNA polymerization in telomerase: a molecular dynamics investigation

Telomerase is an RNA-dependent DNA polymerase that plays a role in the maintenance of the 3' end of the eukaryotic chromosome, known as a telomere, by catalyzing the DNA polymerization reaction in cancer and embryonic stem cells. The detailed molecular details of the DNA polymerization by telomerase, especially the general base for deprotonating the terminal 3'-hydroxyl, which triggers the chemical reaction, remain elusive. We conducted a computational investigation using hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to probe the detailed mechanism of the reaction. Our simulations started with the telomerase:RNA:DNA:dNTP ternary complex, and by using enhanced sampling QM/MM MD simulations, we probed the general base involved directly in the polymerization. We report the participation of an aspartate (Asp344) coordinated to Mg and an active site water molecule, jointly acting as a base during nucleic acid addition. The Asp344 residue remains transiently protonated during the course of the reaction, and later it deprotonates by transferring its proton to the water at the end of the reaction.

Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian

Molecular dynamics (MD) simulations employing ab initio quantum mechanical and molecular mechanical (ai-QM/MM) potentials are considered to be the state of the art, but the high computational cost associated with the ai-QM calculations remains a theoretical challenge for their routine application. Here, we present a modified protocol of the multiple time step (MTS) method for accelerating ai-QM/MM MD simulations of condensed-phase reactions. Within a previous MTS protocol [Nam J. Chem. Theory Comput. 2014, 10, 4175], reference forces are evaluated using a low-level (semiempirical QM/MM) Hamiltonian and employed at inner time steps to propagate the nuclear motions. Correction forces, which arise from the force differences between high-level (ai-QM/MM) and low-level Hamiltonians, are applied at outer time steps, where the MTS algorithm allows the time-reversible integration of the correction forces. To increase the outer step size, which is bound by the highest-frequency component in the correction forces, the semiempirical QM Hamiltonian is recalibrated in this work to minimize the magnitude of the correction forces. The remaining high-frequency modes, which are mainly bond stretches involving hydrogen atoms, are then removed from the correction forces. When combined with a Langevin or SIN(R) thermostat, the modified MTS-QM/MM scheme remains robust with an up to 8 (with Langevin) or 10 fs (with SIN(R)) outer time step (with 1 fs inner time steps) for the chorismate mutase system. This leads to an over 5-fold speedup over standard ai-QM/MM simulations, without sacrificing the accuracy in the predicted free energy profile of the reaction.

Progress on Exploring the Luminescent Properties of Organic Molecular Aggregates by Multiscale Modeling

Luminescent molecular aggregates have attracted worldwide attention because of their potential applications in many fields. The luminescent properties of organic aggregates are complicated and highly morphology-dependent, unraveling the intrinsic mechanism behind is urgent. This review summarizes recent works on investigating the structure-property relationships of organic molecular aggregates at different environments, including crystal, cocrystal, amorphous aggregate, and doped systems by multiscale modeling protocol. We aim to explore the influence of intermolecular non-covalent interactions on molecular packing and their photophysical properties and then pave the effective way to design, synthesize, and develop advanced organic luminescent materials.

5-Formylcytosine mediated DNA-protein cross-links block DNA replication and induce mutations in human cells

5-Formylcytosine (5fC) is an epigenetic DNA modification introduced via TET protein-mediated oxidation of 5-methyl-dC. We recently reported that 5fC form reversible DNA–protein conjugates (DPCs) with histone proteins in living cells (Ji et al. (2017) Angew. Chem. Int. Ed., 56:14130–14134). We now examined the effects of 5fC mediated DPCs on DNA replication. Synthetic DNA duplexes containing site-specific DPCs between 5fC and lysine-containing proteins and peptides were subjected to primer extension experiments in the presence of human translesion synthesis DNA polymerases η and κ. We found that DPCs containing histones H2A or H4 completely inhibited DNA replication, but the replication block was removed when the proteins were subjected to proteolytic digestion. Cross-links to 11-mer or 31-mer peptides were bypassed by both polymerases in an error-prone manner, inducing targeted C→T transitions and –1 deletions. Similar types of mutations were observed when plasmids containing 5fC-peptide cross-links were replicated in human embryonic kidney (HEK) 293T cells. Molecular simulations of the 11-mer peptide-dC cross-links bound to human polymerases η and κ revealed that the peptide fits well on the DNA major groove side, and the modified dC forms a stable mismatch with incoming dATP via wobble base pairing in the polymerase active site.

Two symmetric arginine residues play distinct roles in Thermus thermophilus Argonaute DNA guide strand-mediated DNA target cleavage

Bacterium Thermus thermophilus Argonaute (Ago; TtAgo) is a prokaryotic Ago (pAgo) that acts as the host defense against the uptake and propagation of foreign DNA by catalyzing the DNA cleavage reaction. The TtAgo active site consists of a plugged-in glutamate finger with two arginine residues (R545 and R486) located symmetrically around it. An interesting challenge is to understand how they can collaboratively facilitate enzymatic catalysis. In Kluyveromyces polysporus Ago, a eukaryotic Ago, the evolutionarily symmetrical residues are arginine and histidine, both of which function to stabilize the plugged-in catalytic tetrad conformation. Surprisingly, our simulation results indicated that, in TtAgo, only R545 is involved in the cleavage reaction by serving as a critical structural anchor to stabilize the catalytic tetrad Asp-Glu-Asp-Asp that is completed by the insertion of the glutamate finger, whereas R486 is not involved in target cleavage. The TtAgo-mediated target DNA cleavage occurs in a substrate-assisted mechanism, in which the pro-Rp (Rp, a tetrahedral phosphorus center with "R-type" chirality) oxygen of scissile phosphate acts as a general base to activate the nucleophilic water. Our unexpected theoretical findings on distinct roles played by R545 and R486 in TtAgo catalysis have been validated by single-point site-mutagenesis experiments, wherein the target cleavage is abolished for all mutants of R545. In sharp contrast, the cleavage activity is maintained for all mutants of R486. Our work provides mechanistic insights on the catalytic specificity of Ago proteins and could facilitate the design of new gene-editing tools in the long term.

Computational Simulations of DNA Polymerases: Detailed Insights on Structure/Function/Mechanism from Native Proteins to Cancer Variants

Genetic information is vital in the cell cycle of DNA-based organisms. DNA polymerases (DNA Pols) are crucial players in transactions dealing with these processes. Therefore, the detailed understanding of the structure, function, and mechanism of these proteins has been the focus of significant effort. Computational simulations have been applied to investigate various facets of DNA polymerase structure and function. These simulations have provided significant insights over the years. This perspective presents the results of various computational studies that have been employed to research different aspects of DNA polymerases including detailed reaction mechanism investigation, mutagenicity of different metal cations, possible factors for fidelity synthesis, and discovery/functional characterization of cancer-related mutations on DNA polymerases.

A new paradigm of DNA synthesis: three-metal-ion catalysis

Enzyme catalysis has been studied for over a century. How it actually occurs has not been visualized until recently. By combining in crystallo reaction and X-ray diffraction analysis of reaction intermediates, we have obtained unprecedented atomic details of the DNA synthesis process. Contrary to the established theory that enzyme-substrate complexes and transition states have identical atomic composition and catalysis occurs by the two-metal-ion mechanism, we have discovered that an additional divalent cation has to be captured en route to product formation. Unlike the canonical two metal ions, which are coordinated by DNA polymerases, this third metal ion is free of enzyme coordination. Its location between the α- and β-phosphates of dNTP suggests that the third metal ion may drive the phosphoryltransfer from the leaving group opposite to the 3'-OH nucleophile. Experimental data indicate that binding of the third metal ion may be the rate-limiting step in DNA synthesis and the free energy associated with the metal-ion binding can overcome the activation barrier to the DNA synthesis reaction.

A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases

The enzymatic polymerization of DNA and RNA is the basis for genetic inheritance for all living organisms. It is catalyzed by the DNA/RNA polymerase (Pol) superfamily. Here, bioinformatics analysis reveals that the incoming nucleotide substrate always forms an H-bond between its 3'-OH and β-phosphate moieties upon formation of the Michaelis complex. This previously unrecognized H-bond implies a novel self-activated mechanism (SAM), which synergistically connects the in situ nucleophile formation with subsequent nucleotide addition and, importantly, nucleic acid translocation. Thus, SAM allows an elegant and efficient closed-loop sequence of chemical and physical steps for Pol catalysis. This is markedly different from previous mechanistic hypotheses. Our proposed mechanism is corroborated via ab initio QM/MM simulations on a specific Pol, the human DNA polymerase-η, an enzyme involved in repairing damaged DNA. The structural conservation of DNA and RNA Pols supports the possible extension of SAM to Pol enzymes from the three domains of life.

Analysis of proton wires in the enzyme active site suggests a mechanism of c-di-GMP hydrolysis by the EAL domain phosphodiesterases

We report for the first time a hydrolysis mechanism of the cyclic dimeric guanosine monophosphate (c-di-GMP) by the EAL domain phosphodiesterases as revealed by molecular simulations. A model system for the enzyme-substrate complex was prepared on the base of the crystal structure of the EAL domain from the BlrP1 protein complexed with c-di-GMP. The nucleophilic hydroxide generated from the bridging water molecule appeared in a favorable position for attack on the phosphorus atom of c-di-GMP. The most difficult task was to find a pathway for a proton transfer to the O3' atom of c-di-GMP to promote the O3'P bond cleavage. We show that the hydrogen bond network extended over the chain of water molecules in the enzyme active site and the Glu359 and Asp303 side chains provides the relevant proton wires. The suggested mechanism is consistent with the structural, mutagenesis, and kinetic experimental studies on the EAL domain phosphodiesterases. Proteins 2016; 84:1670-1680. © 2016 Wiley Periodicals, Inc.

Born-Oppenheimer Ab Initio QM/MM Molecular Dynamics Simulations of Enzyme Reactions

There are two key requirements for reliably simulating enzyme reactions: one is a reasonably accurate potential energy surface to describe the bond-forming/breaking process as well as to adequately model the heterogeneous enzyme environment; the other is to perform extensive sampling since an enzyme system consists of at least thousands of atoms and its energy landscape is very complex. One attractive approach to meet both daunting tasks is Born-Oppenheimer ab initio QM/MM molecular dynamics (aiQM/MM-MD) simulation with umbrella sampling. In this chapter, we describe our recently developed pseudobond Q-Chem-Amber interface, which employs a combined electrostatic-mechanical embedding scheme with periodic boundary condition and the particle mesh Ewald method for long-range electrostatics interactions. In our implementation, Q-Chem and the sander module of Amber are combined at the source code level without using system calls, and all necessary data communications between QM and MM calculations are achieved via computer memory. We demonstrate the applicability of this pseudobond Q-Chem-Amber interface by presenting two examples, one reaction in aqueous solution and one enzyme reaction. Finally, we describe our established aiQM/MM-MD enzyme simulation protocol, which has been successfully applied to study more than a dozen enzymes.

Structure and mechanism of error-free replication past the major benzo[a]pyrene adduct by human DNA polymerase κ

Benzo[a]pyrene (BP) is a well-known and frequently encountered carcinogen which generates a bulky DNA adduct (+)-trans-10S-BP-N(2)-dG (BP-dG) in cells. DNA polymerase kappa (polκ) is the only known Y-family polymerase that bypasses BP-dG accurately and thus protects cells from genotoxic BP. Here, we report the structures of human polκ in complex with DNA containing either a normal guanine (G) base or a BP-dG adduct at the active site and a correct deoxycytidine. The structures and supporting biochemical data reveal a unique mechanism for accurate replication by translesion synthesis past the major bulky adduct. The active site of polκ opens at the minor groove side of the DNA substrate to accommodate the bulky BP-dG that is attached there. More importantly, polκ stabilizes the lesion DNA substrate in the same active conformation as for regular B-form DNA substrates and the bulky BPDE ring in a 5' end pointing conformation. The BP-dG adducted DNA substrate maintains a Watson-Crick (BP-dG:dC) base pair within the active site, governing correct nucleotide insertion opposite the bulky adduct. In addition, polκ's unique N-clasp domain supports the open conformation of the enzyme and the extended conformation of the single-stranded template to allow bypass of the bulky lesion. This work illustrates the first molecular mechanism for how a bulky major adduct is replicated accurately without strand misalignment and mis-insertion.

Cooperative motion of a key positively charged residue and metal ions for DNA replication catalyzed by human DNA Polymerase-η

Trans-lesion synthesis polymerases, like DNA Polymerase-η (Pol-η), are essential for cell survival. Pol-η bypasses ultraviolet-induced DNA damages via a two-metal-ion mechanism that assures DNA strand elongation, with formation of the leaving group pyrophosphate (PPi). Recent structural and kinetics studies have shown that Pol-η function depends on the highly flexible and conserved Arg61 and, intriguingly, on a transient third ion resolved at the catalytic site, as lately observed in other nucleic acid-processing metalloenzymes. How these conserved structural features facilitate DNA replication, however, is still poorly understood. Through extended molecular dynamics and free energy simulations, we unravel a highly cooperative and dynamic mechanism for DNA elongation and repair, which is here described by an equilibrium ensemble of structures that connect the reactants to the products in Pol-η catalysis. We reveal that specific conformations of Arg61 help facilitate the recruitment of the incoming base and favor the proper formation of a pre-reactive complex in Pol-η for efficient DNA editing. Also, we show that a third transient metal ion, which acts concertedly with Arg61, serves as an exit shuttle for the leaving PPi. Finally, we discuss how this effective and cooperative mechanism for DNA repair may be shared by other DNA-repairing polymerases.

Applications of quantum mechanical/molecular mechanical methods to the chemical insertion step of DNA and RNA polymerization

We review theoretical attempts to model the chemical insertion reactions of nucleoside triphosphates catalyzed by the nucleic acid polymerases using combined quantum mechanical/molecular mechanical methodology. Due to an existing excellent database of high-resolution X-ray crystal structures, the DNA polymerase β system serves as a useful template for discussion and comparison. The convergence of structures of high-quality complexes and continued developments of theoretical techniques suggest a bright future for understanding the global features of nucleic acid polymerization.

Structural and dynamic characterization of polymerase κ's minor groove lesion processing reveals how adduct topology impacts fidelity

DNA lesion bypass polymerases process different lesions with varying fidelities, but the structural, dynamic, and mechanistic origins of this phenomenon remain poorly understood. Human DNA polymerase κ (Polκ), a member of the Y family of lesion bypass polymerases, is specialized to bypass bulky DNA minor groove lesions in a predominantly error-free manner, by housing them in its unique gap. We have investigated the role of the unique Polκ gap and N-clasp structural features in the fidelity of minor groove lesion processing with extensive molecular modeling and molecular dynamics simulations to pinpoint their functioning in lesion bypass. Here we consider the N²-dG covalent adduct derived from the carcinogenic aromatic amine, 2-acetylaminofluorene (dG-N²-AAF), that is produced via the combustion of kerosene and diesel fuel. Our simulations reveal how the spacious gap directionally accommodates the lesion aromatic ring system as it transits through the stages of incorporation of the predominant correct partner dCTP opposite the damaged guanine, with preservation of local active site organization for nucleotidyl transfer. Furthermore, flexibility in Polκ’s N-clasp facilitates the significant misincorporation of dTTP opposite dG-N²-AAF via wobble pairing. Notably, we show that N-clasp flexibility depends on lesion topology, being markedly reduced in the case of the benzo[a]pyrene-derived major adduct to N²-dG, whose bypass by Polκ is nearly error-free. Thus, our studies reveal how Polκ’s unique structural and dynamic properties can regulate its bypass fidelity of polycyclic aromatic lesions and how the fidelity is impacted by lesion structures.

An overview of Y-Family DNA polymerases and a case study of human DNA polymerase η

Y-Family DNA polymerases specialize in translesion synthesis, bypassing damaged bases that would otherwise block the normal progression of replication forks. Y-Family polymerases have unique structural features that allow them to bind damaged DNA and use a modified template base to direct nucleotide incorporation. Each Y-Family polymerase is unique and has different preferences for lesions to bypass and for dNTPs to incorporate. Y-Family polymerases are also characterized by a low catalytic efficiency, a low processivity, and a low fidelity on normal DNA. Recruitment of these specialized polymerases to replication forks is therefore regulated. The catalytic center of the Y-Family polymerases is highly conserved and homologous to that of high-fidelity and high-processivity DNA replicases. In this review, structural differences between Y-Family and A- and B-Family polymerases are compared and correlated with their functional differences. A time-resolved X-ray crystallographic study of the DNA synthesis reaction catalyzed by the Y-Family DNA polymerase human polymerase η revealed transient elements that led to the nucleotidyl-transfer reaction.

Variants of mouse DNA polymerase κ reveal a mechanism of efficient and accurate translesion synthesis past a benzo[a]pyrene dG adduct

DNA polymerase κ (Polκ) is the only known Y-family DNA polymerase that bypasses the 10S (+)-trans-anti-benzo[a]pyrene diol epoxide (BPDE)-N(2)-deoxyguanine adducts efficiently and accurately. The unique features of Polκ, a large structure gap between the catalytic core and little finger domain and a 90-residue addition at the N terminus known as the N-clasp, may give rise to its special translesion capability. We designed and constructed two mouse Polκ variants, which have reduced gap size on both sides [Polκ Gap Mutant (PGM) 1] or one side flanking the template base (PGM2). These Polκ variants are nearly as efficient as WT in normal DNA synthesis, albeit with reduced accuracy. However, PGM1 is strongly blocked by the 10S (+)-trans-anti-BPDE-N(2)-dG lesion. Steady-state kinetic measurements reveal a significant reduction in efficiency of dCTP incorporation opposite the lesion by PGM1 and a moderate reduction by PGM2. Consistently, Polκ-deficient cells stably complemented with PGM1 GFP-Polκ remained hypersensitive to BPDE treatment, and complementation with WT or PGM2 GFP-Polκ restored BPDE resistance. Furthermore, deletion of the first 51 residues of the N-clasp in mouse Polκ (mPolκ(52-516)) leads to reduced polymerization activity, and the mutant PGM2(52-516) but not PGM1(52-516) can partially compensate the N-terminal deletion and restore the catalytic activity on normal DNA. However, neither WT nor PGM2 mPolκ(52-516) retains BPDE bypass activity. We conclude that the structural gap physically accommodates the bulky aromatic adduct and the N-clasp is essential for the structural integrity and flexibility of Polκ during translesion synthesis.

Molecular modeling of nucleic acid structure: energy and sampling

An overview of computer simulation techniques as applied to nucleic acid systems is presented. This unit discusses methods used to treat energy and to sample representative configurations. Emphasis is placed on molecular mechanics and empirical force fields.

Quantum mechanical modeling: a tool for the understanding of enzyme reactions

Most enzyme reactions involve formation and cleavage of covalent bonds, while electrostatic effects, as well as dynamics of the active site and surrounding protein regions, may also be crucial. Accordingly, special computational methods are needed to provide an adequate description, which combine quantum mechanics for the reactive region with molecular mechanics and molecular dynamics describing the environment and dynamic effects, respectively. In this review we intend to give an overview to non-specialists on various enzyme models as well as established computational methods and describe applications to some specific cases. For the treatment of various enzyme mechanisms, special approaches are often needed to obtain results, which adequately refer to experimental data. As a result of the spectacular progress in the last two decades, most enzyme reactions can be quite precisely treated by various computational methods.

DNA cytosine methylation: structural and thermodynamic characterization of the epigenetic marking mechanism

DNA cytosine methyltransferases regulate the expression of the genome through the precise epigenetic marking of certain cytosines with a methyl group, and aberrant methylation is a hallmark of human diseases including cancer. Targeting these enzymes for drug design is currently a high priority. We have utilized ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to investigate extensively the reaction mechanism of the representative DNA methyltransferase HhaI (M.HhaI) from prokaryotes, whose overall mechanism is shared with the mammalian enzymes. We obtain for the first time full free energy profiles for the complete reaction, together with reaction dynamics in atomistic detail. Our results show an energetically preferred mechanism in which nucleophilic attack of cytosine C5 on the S-adenosyl-L-methionine (AdoMet) methyl group is concerted with formation of the Michael adduct between a conserved Cys in the active site with cytosine C6. Spontaneous and reversible proton transfer between a conserved Glu in the active site and cytosine N3 at the transition state was observed in our simulations, revealing the chemical participation of this Glu residue in the catalytic mechanism. Subsequently, the β-elimination of the C5 proton utilizes as base an OH(-) derived from a conserved crystal water that is part of a proton wire water channel, and this syn β-elimination reaction is the rate-limiting step. Design of novel cytosine methylation inhibitors would be advanced by our structural and thermodynamic characterization of the reaction mechanism.

Sirtuin Deacetylation Mechanism and Catalytic Role of the Dynamic Cofactor Binding Loop

Sirtuins constitute a novel family of protein deacetylases and play critical roles in epigenetics, cell death, and metabolism. In spite of numerous experimental studies, the key and most complicated stage of its NAD+-dependent catalytic mechanism remains to be elusive. Herein by employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, a state-of-the-art computational approach to study enzyme reactions, we have characterized the complete deacetylation mechanism for a sirtuin enzyme, determined its multistep free-energy reaction profile, and elucidated essential catalytic roles of the conserved dynamic cofactor binding loop. These new detailed mechanistic insights would facilitate the design of novel mechanism-based sirtuin modulators.

Revelation of a catalytic calcium-binding site elucidates unusual metal dependence of a human apyrase

Human soluble calcium-activated nucleotidase 1 (hSCAN-1) represents a new family of apyrase enzymes that catalyze the hydrolysis of nucleotide di- and triphosphates, thereby modulating extracellular purinergic and pyrimidinergic signaling. Among well-characterized phosphoryl transfer enzymes, hSCAN-1 is unique not only in its unusual calcium-dependent activation, but also in its novel phosphate-binding motif. Its catalytic site does not utilize backbone amide groups to bind phosphate, as in the common P-loop, but contains a large cluster of acidic ionizable side chains. By employing a state-of-the-art computational approach, we have revealed a previously uncharacterized catalytic calcium-binding site in hSCAN-1, which elucidates the unusual calcium-dependence of its apyrase activity. In a high-order coordination shell, the newly identified calcium ion organizes the active site residues to mediate nucleotide binding, to orient the nucleophilic water, and to facilitate the phosphoryl transfer reaction. From ab initio QM/MM molecular dynamics simulations with umbrella sampling, we have characterized a reverse protonation catalytic mechanism for hSCAN-1 and determined its free energy reaction profile. Our results are consistent with available experimental studies and provide new detailed insight into the structure-function relationship of this novel calcium-activated phosphoryl transfer enzyme.

Molecular modeling of nucleic acid structure: energy and sampling

An overview of computer simulation techniques as applied to nucleic acid systems is presented. This unit expands an accompanying overview unit (UNIT 7.5) by discussing methods used to treat the energy and sample representative configurations. Emphasis is placed on molecular mechanics and empirical force fields.