All-atom virus simulations.

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

Assembly Properties of Hepatitis B Virus Core Protein Mutants Correlate with Their Resistance to Assembly-Directed Antivirals.

The hepatitis B virus (HBV) capsid or core protein (Cp) can self-assemble to form an icosahedral capsid. It is now being pursued as a target for small-molecule antivirals that enhance the rate and extent of its assembly to yield empty and/or aberrant capsids. These small molecules are thus called core protein allosteric modulators (CpAMs). We sought to understand the physical basis of CpAM-resistant mutants and how CpAMs might overcome them. We examined the effects of two closely related CpAMs, HAP12 and HAP13, which differ by a single atom but have drastically different antiviral activities, on the assembly of wild-type Cp and three T109 mutants (T109M, T109I, and T109S) that display a range of resistances. The T109 side chain forms part of the mouth of the CpAM binding pocket. A T109 mutant that has substantial resistance even to a highly active CpAM strongly promotes normal assembly. Conversely, a mutant that weakens assembly is more susceptible to CpAMs. In crystal and cryo-electron microscopy (cryo-EM) structures of T=4 capsids with bound CpAMs, the CpAMs preferentially fit into two of four quasi-equivalent sites. In these static representations of capsid structures, T109 does not interact with the neighboring subunit. However, all-atom molecular dynamics simulations of an intact capsid show that T109 of one of the four classes of CpAM site has a hydrophobic contact with the neighboring subunit at least 40% of the time, providing a physical explanation for the mutation's ability to affect capsid stability, assembly, and sensitivity to CpAMs.IMPORTANCE The HBV core protein and its assembly into capsids have become important targets for development of core protein allosteric modulators (CpAMs) as antivirals. Naturally occurring T109 mutants have been shown to be resistant to some of these CpAMs. We found that mutation of T109 led to changes in capsid stability and recapitulated resistance to a weak CpAM, but much less so than to a strong CpAM. Examination of HBV capsid structures, determined by cryo-EM and crystallography, could not explain how T109 mutations change capsid stability and resistance. However, by mining data from a microsecond-long all-atom molecular dynamics simulation, we found that the capsid was extraordinarily flexible and that T109 can impede entry to the CpAM binding site. In short, HBV capsids are incredibly dynamic and molecular mobility must be considered in discussions of antiviral mechanisms.

All-atom molecular dynamics of the HBV capsid reveals insights into biological function and cryo-EM resolution limits.

The hepatitis B virus capsid represents a promising therapeutic target. Experiments suggest the capsid must be flexible to function; however, capsid structure and dynamics have not been thoroughly characterized in the absence of icosahedral symmetry constraints. Here, all-atom molecular dynamics simulations are leveraged to investigate the capsid without symmetry bias, enabling study of capsid flexibility and its implications for biological function and cryo-EM resolution limits. Simulation results confirm flexibility and reveal a propensity for asymmetric distortion. The capsid's influence on ionic species suggests a mechanism for modulating the display of cellular signals and implicates the capsid's triangular pores as the location of signal exposure. A theoretical image reconstruction performed using simulated conformations indicates how capsid flexibility may limit the resolution of cryo-EM. Overall, the present work provides functional insight beyond what is accessible to experimental methods and raises important considerations regarding asymmetry in structural studies of icosahedral virus capsids.

Molecular Dynamics Simulations of Protein-Drug Complexes: A Computational Protocol for Investigating the Interactions of Small-Molecule Therapeutics with Biological Targets and Biosensors.

MD simulations provide a powerful tool for the investigation of protein-drug complexes. The following chapter uses the aryl acylamidase-acetaminophen system as an example to describe a general protocol for preparing and running simulations of protein-drug complexes, complete with a step-by-step tutorial. The described approach is broadly applicable toward the study of drug interactions in the context of both biological targets and biosensing enzymes.

New insights into influenza A specificity: an evolution of paradigms.

Understanding the molecular origin of influenza receptor specificity is complicated by the paucity of quantitative affinity measurements, and the qualitative and variable nature of glycan array data. Further obstacles arise from the varied impact of viral glycosylation and the relatively narrow spectrum of biologically relevant receptors present on glycan arrays. A survey of receptor conformational properties is presented, leading to the conclusion that conformational entropy plays a key role in defining specificity, as does the newly reported ability of biantennary receptors that terminate in Siaα2-6Gal sequences to form bidentate interactions to two binding sites in a hemagglutinin trimer. Bidentate binding provides a functional explanation for the observation that Siaα2-6 receptors adopt an open-umbrella topology when bound to hemagglutinins from human-infective viruses, and calls for a reassessment of virus avidity and tissue tropism.

Angular measurements of the dynein ring reveal a stepping mechanism dependent on a flexible stalk.

The force-generating mechanism of dynein differs from the force-generating mechanisms of other cytoskeletal motors. To examine the structural dynamics of dynein's stepping mechanism in real time, we used polarized total internal reflection fluorescence microscopy with nanometer accuracy localization to track the orientation and position of single motors. By measuring the polarized emission of individual quantum nanorods coupled to the dynein ring, we determined the angular position of the ring and found that it rotates relative to the microtubule (MT) while walking. Surprisingly, the observed rotations were small, averaging only 8.3°, and were only weakly correlated with steps. Measurements at two independent labeling positions on opposite sides of the ring showed similar small rotations. Our results are inconsistent with a classic power-stroke mechanism, and instead support a flexible stalk model in which interhead strain rotates the rings through bending and hinging of the stalk. Mechanical compliances of the stalk and hinge determined based on a 3.3-µs molecular dynamics simulation account for the degree of ring rotation observed experimentally. Together, these observations demonstrate that the stepping mechanism of dynein is fundamentally different from the stepping mechanisms of other well-studied MT motors, because it is characterized by constant small-scale fluctuations of a large but flexible structure fully consistent with the variable stepping pattern observed as dynein moves along the MT.

Computational Methodologies for Real-Space Structural Refinement of Large Macromolecular Complexes.

The rise of the computer as a powerful tool for model building and refinement has revolutionized the field of structure determination for large biomolecular systems. Despite the wide availability of robust experimental methods capable of resolving structural details across a range of spatiotemporal resolutions, computational hybrid methods have the unique ability to integrate the diverse data from multimodal techniques such as X-ray crystallography and electron microscopy into consistent, fully atomistic structures. Here, commonly employed strategies for computational real-space structural refinement are reviewed, and their specific applications are illustrated for several large macromolecular complexes: ribosome, virus capsids, chemosensory array, and photosynthetic chromatophore. The increasingly important role of computational methods in large-scale structural refinement, along with current and future challenges, is discussed.

All-Atom Molecular Dynamics of Virus Capsids as Drug Targets.

Virus capsids are protein shells that package the viral genome. Although their morphology and biological functions can vary markedly, capsids often play critical roles in regulating viral infection pathways. A detailed knowledge of virus capsids, including their dynamic structure, interactions with cellular factors, and the specific roles that they play in the replication cycle, is imperative for the development of antiviral therapeutics. The following Perspective introduces an emerging area of computational biology that focuses on the dynamics of virus capsids and capsid-protein assemblies, with particular emphasis on the effects of small-molecule drug binding on capsid structure, stability, and allosteric pathways. When performed at chemical detail, molecular dynamics simulations can reveal subtle changes in virus capsids induced by drug molecules a fraction of their size. Here, the current challenges of performing all-atom capsid-drug simulations are discussed, along with an outlook on the applicability of virus capsid simulations to reveal novel drug targets.

Calculating binding free energies for protein-carbohydrate complexes.

A variety of computational techniques may be applied to compute theoretical binding free energies for protein-carbohydrate complexes. Elucidation of the intermolecular interactions, as well as the thermodynamic effects, that contribute to the relative strength of receptor binding can shed light on biomolecular recognition, and the resulting initiation or inhibition of a biological process. Three types of free energy methods are discussed here, including MM-PB/GBSA, thermodynamic integration, and a non-equilibrium alternative utilizing SMD. Throughout this chapter, the well-known concanavalin A lectin is employed as a model system to demonstrate the application of these methods to the special case of carbohydrate binding.

Effect of microfibril twisting on theoretical powder diffraction patterns of cellulose Iß

Previous studies of calculated diffraction patterns for cellulose crystallites suggest that distortions that arise once models have been subjected to MD simulation are the result of both microfibril twisting and changes in unit cell dimensions induced by the empirical force field; to date, it has not been possible to separate the individual contributions of these effects. To provide a better understanding of how twisting manifests in diffraction data, the present study demonstrates a method for generating twisted and linear cellulose structures that can be compared without the bias of dimensional changes, allowing assessment of the impact of twisting alone. Analysis of unit cell dimensions, microfibril volume, hydrogen bond patterns, glycosidic torsion angles, and hydroxymethyl group orientations confirmed that the twisted and linear structures collected with this method were internally consistent, and theoretical powder diffraction patterns for the two were shown to be effectively indistinguishable. These results indicate that differences between calculated patterns for the crystal coordinates and twisted structures from MD simulation can result entirely from changes in unit cell dimensions, and not from microfibril twisting alone. Although powder diffraction patterns for models in the 81-chain size regime were shown to be unaffected by twisting, suggesting that a modest degree of twist is not inconsistent with experimental data, it may be that other diffraction techniques are capable of detecting this structural difference. Until such time as definitive experimental evidence comes to light, the results of this study suggest that both twisted and linear microfibrils may represent an appropriate model for cellulose Iß.

Unraveling cellulose microfibrils: a twisted tale.

Molecular dynamics (MD) simulations of cellulose microfibrils are pertinent to the paper, textile, and biofuels industries for their unique capacity to characterize dynamic behavior and atomic-level interactions with solvent molecules and cellulase enzymes. While high-resolution crystallographic data have established a solid basis for computational analysis of cellulose, previous work has demonstrated a tendency for modeled microfibrils to diverge from the linear experimental structure and adopt a twisted conformation. Here, we investigate the dependence of this twisting behavior on computational approximations and establish the theoretical basis for its occurrence. We examine the role of solvent, the effect of nonbonded force field parameters [partial charges and van der Waals (vdW) contributions], and the use of explicitly modeled oxygen lone pairs in both the solute and solvent. Findings suggest that microfibril twisting is favored by vdW interactions, and counteracted by both intrachain hydrogen bonds and solvent effects at the microfibril surface.