The adaptive-partitioning (AP) schemes for combined quantum-mechanical/molecular-mechanical (QM/MM) calculations allow on-the-fly reclassifications of atoms and molecules as QM or MM in dynamics simulations. The permuted-AP (PAP) scheme (J. Phys. Chem. B 2007, 111, 2231) introduces a thin layer of buffer zone between the QM subsystem (also called active zone) and the MM subsystem (also known as the environmental zone) to provide a continuous and smooth transition and expresses the potential energy in a many-body expansion manner. The PAP scheme has been successfully applied to study small molecules solvated in bulk solvent. Here, we propose two modifications to the original PAP scheme to treat solvent molecules entering and leaving protein binding sites. First, the center of the active zone is placed at a pseudoatom in the binding site, whose position is not affected by the movements of ligand or residues in the binding site. Second, the extra forces due to the smoothing functions are deleted. The modified PAP scheme no longer describes a Hamiltonian system, but it satisfies the conservation of momentum. As a proof-of-concept experiment, the modified PAP scheme is applied to the simulations under the canonical ensemble for two binding sites of the Escherichia coli CLC chloride ion transport protein, in particular, the intracellular binding site Sint discovered by crystallography and one putative additional binding site Sadd suggested by molecular modeling. The exchange of water molecules between the binding sites and bulk solvent is monitored. For comparison, simulations are also carried out using the same model system and setup with only one exception: the extra forces due to the smoothing functions are retained. The simulations are benchmarked against conventional QM/MM simulations with large QM subsystems. The results demonstrate that the active zone centered at the pseudo atom is a reasonable and convenient representation of the binding site. Moreover, the transient extra forces are non-negligible and cause the QM water molecules to move out of the active zone. The modified PAP scheme, where the extra forces are excluded, avoids the artifact, providing a realistic description of the exchange of water molecules between the protein binding sites and bulk solvent.

ClC transport proteins show a distinct "broken-helix" architecture, in which certain α-helices are oriented with their N-terminal ends pointed toward the binding sites where the chloride ions are held extensively by the backbone amide nitrogen atoms from the helices. To understand the effectiveness of such binding structures, we carried out natural bond orbital analysis and energy decomposition analysis employing truncated active-site model systems for the bound chloride ions along the translocation pore of the EcClC proteins. Our results indicated that the chloride ions are stabilized in such a binding environment by electrostatic, polarization, and charge-transfer interactions with the backbone and a few side chains. Up to ~25% of the formal charges of the chloride ions were found smeared out to the surroundings primarily via charge transfer from the chloride's lone pair n(Cl) orbitals to the protein's antibonding σ*(N-H) or σ*(O-H) orbitals; those σ* orbitals are localized at the polar N-H and O-H bonds in the chloride's first solvation shells formed by the backbone amide groups and the side chains of residues Ser107, Arg147, Glu148, and Tyr445. Polarizations by the chloride ions were dominated by the redistribution of charge densities among the π orbitals and lone pair orbitals of the protein atoms, in particular the atoms of the backbone peptide links and of the side chains of Arg147, Glu148, and Tyr445. The substantial amounts of electron density involved in charge transfer and in polarization were consistent with the large energetic contributions by the two processes revealed by the energy decomposition analysis. The significant polarization and charge-transfer effects may have impacts on the mechanisms and dynamics of the chloride transport by the ClC proteins.

Complex double-stranded DNA viruses utilize a terminase enzyme to package their genomes into a preassembled procapsid shell. DNA packaging triggers a major conformational change in the proteins assembled into the shell and most often subsequent addition of a decoration protein that is required to stabilize the structure. In bacteriophage λ, DNA packaging drives a procapsid expansion transition to afford a larger but fragile shell. The gpD decoration protein adds to the expanded shell as trimeric spikes at each of the 140 three-fold axes. The spikes provide mechanical strength to the shell such that it can withstand the tremendous internal forces generated by the packaged DNA in addition to environmental insults. Hydrophobic, electrostatic, and aromatic-proline noncovalent interactions have been proposed to mediate gpD trimer spike assembly at the expanded shell surface. Here, we directly examine each of these interactions and demonstrate that hydrophobic interactions play the dominant role. In the course of this study, we unexpectedly found that Trp308 in the λ major capsid protein (gpE) plays a critical role in shell assembly. The gpE-W308A mutation affords a soluble, natively folded protein that does not further assemble into a procapsid shell, despite the fact that it retains binding interactions with the scaffolding protein, the shell assembly chaparone protein. The data support a model in which the λ procapsid shell assembles via cooperative interaction of monomeric capsid proteins, as observed in the herpesviruses and phages such as P22. The significance of the results with respect to capsid assembly, maturation, and stability is discussed.

Like many dsDNA viruses, bacteriophage λ replicates its genome as a concatemer consisting of multiple copies of covalently linked dsDNA genomes. To encapsidate a single genome within a nascent procapsid, λ must: 1) find its own dsDNA amongst the multitude of host nucleic acids; 2) identify the genomic start site; 3) cut the DNA; 4) bring the excised DNA to a procapsid; 5) translocate DNA into the capsid; 6) cut DNA again at a packaging termination site, 7) disengage from the newly filled capsid; and 8) bring the remainder of the genomic concatemer to fill another empty procapsid. These disparate genome processing tasks are carried out by a single virus-encoded enzyme complex called terminase. While it has been shown that λ terminase initially forms a tetrameric complex to cut DNA, it is not clear whether the same configuration translocates DNA. Here, we describe biophysical and initial structural characterization of a λ terminase translocation complex. Analytical ultracentrifugation (AUC) and small angle X-ray scattering (SAXS) indicate that between 4 and 5 protomeric subunits assemble a cone-shaped terminase complex with a maximum dimension of ~230 and radius of gyration of ~72 Å. Two-dimensional classification of cryoEM images of λ terminase are consistent with these dimensions and show that particles assume a preferred orientation in ice. The orientations appear to be end-on, as terminase rings resemble a starfish with approximate pentameric symmetry. While ~5-fold symmetry is apparent, one of the five arms appears partially displaced with weaker more diffuse density in some classes, suggesting flexibility and/or partial occupancy. Charge detection mass spectrometry (CDMS) is consistent with a pentameric complex, with evidence that one motor subunit is weakly bound. Kinetic analysis indicates that the complex hydrolyzes ATP at a rate comparable to the rates of other phage packaging motors. Together with previously published data, these results suggest that λ terminase assembles conformationally and stoichiometrically distinct complexes to carry out different genome processing tasks. We propose a symmetry resolution pathway to explain how terminase transitions between these structurally and functionally distinct states.

No abstract