Rare Events Sampling Methods for Quantum and Classical Ab Initio Molecular Dynamics.

We provide an approach to sample rare events during classical ab initio molecular dynamics and quantum wavepacket dynamics. For classical AIMD, a set of fictitious degrees of freedom are introduced that may harmonically interact with the electronic and nuclear degrees of freedom to steer the dynamics in a conservative fashion toward energetically forbidden regions. A similar approach when introduced for quantum wavepacket dynamics has the effect of biasing the trajectory of the wavepacket centroid toward the regions of the potential surface that are difficult to sample.

On the Possibility That Bond Strain Is the Mechanism of RING E3 Activation in the E2-Catalyzed Ubiquitination Reaction.

Ubiquitination is a type of post-translational modification wherein the small protein ubiquitin (Ub) is covalently bound to a lysine on a target protein. Ubiquitination can signal for several regulatory pathways including protein degradation. Ubiquitination occurs by a series of reactions catalyzed by three types of enzymes: ubiquitin activating enzymes, E1; ubiquitin conjugating enzymes, E2; and ubiquitin ligases, E3.

A conserved asparagine in a ubiquitin-conjugating enzyme positions the substrate for nucleophilic attack.

The mechanism used by the ubiquitin-conjugating enzyme, Ubc13, to catalyze ubiquitination is probed with three computational techniques: Born-Oppenheimer molecular dynamics, single point quantum mechanics/molecular mechanics energies, and classical molecular dynamics. These simulations support a long-held hypothesis and show that Ubc13-catalyzed ubiquitination uses a stepwise, nucleophilic attack mechanism. Furthermore, they show that the first step-the formation of a tetrahedral, zwitterionic intermediate-is rate limiting.

Far-Infrared Synchrotron Spectroscopy and Quantum Chemical Calculations of the Potentially Important Interstellar Molecule, 2-Chloroethanol.

The high brightness of the Australian synchrotron allowed for detailed spectra to be collected at high resolution (0.00096 cm) in the vicinity of the a/ b/ c-type ν band of 2-chloroethanol, which involves O-H torsional motion about the C-O bond. A rovibrational analysis was performed for both chlorine isotopologues in the ν fundamental (centered at ∼344 cm) which involved the assignment of 7153 lines ( J ≤ 90, K ≤ 41).

Molecular dynamics simulations reveal a new role for a conserved active site asparagine in a ubiquitin-conjugating enzyme.

The role of a highly conserved active site asparagine (N79) in the ubiquitin conjugating enzyme, Ubc13, is probed using molecular dynamics simulations. Both wild type and mutant enzymes (N79A and N79D) are studied. Contrary to a popular hypothesis, we show that it is unlikely that N79 stabilizes a reaction intermediate, but instead preferentially hydrogen bonds to a loop near the active site.

The Disulfide Bonds within BST-2 Enhance Tensile Strength during Viral Tethering.

Human BST-2/tetherin is a host factor that inhibits the release of enveloped viruses, including HIV-1, HIV-2, and SIV, from the cell surface by tethering viruses to the host cell membrane. BST-2 has an α-helical ectodomain that forms disulfide-linked dimers between two monomers forming a coiled coil. The ectodomain contains three cysteine residues that can participate in disulfide bond formation and are critical for viral tethering.

Mechanical dissociation of the M-band titin/obscurin complex is directionally dependent.

Titin and obscurin, two giant muscle proteins, bind to each other in an antiparallel Ig-Ig fashion at the M-band. This interaction must be able to withstand the mechanical strain that the M-band typically experiences and remain intact. The mechanical force on these domains is likely exerted along one of two axes: a longitudinal axis, resulting in a 'shearing' force, or a lateral axis, resulting in a 'peeling' force.

Gauging the flexibility of the active site in soybean lipoxygenase-1 (SLO-1) through an atom-centered density matrix propagation (ADMP) treatment that facilitates the sampling of rare events.

We present a computational methodology to sample rare events in large biological enzymes that may involve electronically polarizing, reactive processes. The approach includes simultaneous dynamical treatment of electronic and nuclear degrees of freedom, where contributions from the electronic portion are computed using hybrid density functional theory and the computational costs are reduced through a hybrid quantum mechanics/molecular mechanics (QM/MM) treatment. Thus, the paper involves a QM/MM dynamical treatment of rare events.

Proton transport pathways in [NiFe]-hydrogenase.

Hydrogenases reversibly catalyze the production of molecular hydrogen. Current interest in these enzymes is focused on understanding the catalysis, since this may prove useful for hydrogen-based fuel cell and photosynthetic hydrogen production cell technologies. A key step in the hydrogenase catalytic cycle and the focus of this work is proton transport (PT) to and from the active site.

Analysis of Hydrogen Tunneling in an Enzyme Active Site using von Neumann Measurements.

We build on our earlier quantum wavepacket study of hydrogen transfer in the biological enzyme, soybean lipoxygenase-1, by using von Neumann quantum measurement theory to gain qualitative insights into the transfer event. We treat the enzyme active site as a measurement device which acts on the tunneling hydrogen nucleus via the potential it exerts at each configuration. A series of changing active site geometries during the tunneling process effects a sequential projection of the initial, reactant state onto the final, product state.

Combining quantum wavepacket ab initio molecular dynamics with QM/MM and QM/QM techniques: Implementation blending ONIOM and empirical valence bond theory.

We discuss hybrid quantum-mechanics/molecular-mechanics (QM/MM) and quantum mechanics/quantum mechanics (QM/QM) generalizations to our recently developed quantum wavepacket ab initio molecular dynamics methodology for simultaneous dynamics of electrons and nuclei. The approach is a synergy between a quantum wavepacket dynamics, ab initio molecular dynamics, and the ONIOM scheme. We utilize this method to include nuclear quantum effects arising from a portion of the system along with a simultaneous description of the electronic structure.

Hydrogen tunneling in an enzyme active site: a quantum wavepacket dynamical perspective.

We study the hydrogen tunneling problem in a model system that represents the active site of the biological enzyme, soybean lipoxygenase-1. Toward this, we utilize quantum wavepacket dynamics performed on potential surfaces obtained by using hybrid density functional theory under the influence of a dynamical active site. The kinetic isotope effect is computed by using the transmission amplitude of the wavepacket, and the experimental value is reproduced.

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects.

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically improved (by several orders of magnitude) through the utilization of wavelet-based techniques combined with the previously introduced time-dependent deterministic sampling procedure measure to achieve stable, picosecond length, quantum-classical dynamics of electrons and nuclei in clusters.

Computational Improvements to Quantum Wave Packet ab Initio Molecular Dynamics Using a Potential-Adapted, Time-Dependent Deterministic Sampling Technique.

In a recent publication, we introduced a computational approach to treat the simultaneous dynamics of electrons and nuclei. The method is based on a synergy between quantum wave packet dynamics and ab initio molecular dynamics. Atom-centered density-matrix propagation or Born-Oppenheimer dynamics can be used to perform ab initio dynamics.