Density is an easily adjusted variable in molecular dynamics (MD) simulations. Thus, pressure-jump (P-jump)-induced protein refolding, if it could be made fast enough, would be ideally suited for comparison with MD. Although pressure denaturation perturbs secondary structure less than temperature denaturation, protein refolding after a fast P-jump is not necessarily faster than that after a temperature jump. Recent P-jump refolding experiments on the helix bundle λ-repressor have shown evidence of a <3 μs burst phase, but also of a ~1.5 ms "slow" phase of refolding, attributed to non-native helical structure frustrating microsecond refolding. Here we show that a λ-repressor mutant is nonetheless capable of refolding in a single explicit solvent MD trajectory in about 19 μs, indicating that the burst phase observed in experiments on the same mutant could produce native protein. The simulation reveals that after about 18.5 μs of conformational sampling, the productive structural rearrangement to the native state does not occur in a single swift step but is spread out over a brief series of helix and loop rearrangements that take about 0.9 μs. Our results support the molecular time scale inferred for λ-repressor from near-downhill folding experiments, where transition-state population can be seen experimentally, and also agrees with the transition-state transit time observed in slower folding proteins by single-molecule spectroscopy.

Biological molecules are often asymmetric with respect to stereochemistry, and correct stereochemistry is essential to their function. Molecular dynamics simulations of biomolecules have increasingly become an integral part of biophysical research. However, stereochemical errors in biomolecular structures can have a dramatic impact on the results of simulations.

The proper functioning of biomolecules in living cells requires them to assume particular structures and to undergo conformational changes. Both biomolecular structure and motion can be studied using a wide variety of techniques, but none offers the level of detail as do molecular dynamics (MD) simulations. Integrating two widely used modeling programs, namely NAMD and VMD, we have created a robust, user-friendly software, QwikMD, which enables novices and experts alike to address biomedically relevant questions, where often only molecular dynamics simulations can provide answers. Performing both simple and advanced MD simulations interactively, QwikMD automates as many steps as necessary for preparing, carrying out, and analyzing simulations while checking for common errors and enabling reproducibility. QwikMD meets also the needs of experts in the field, increasing the efficiency and quality of their work by carrying out tedious or repetitive tasks while enabling easy control of every step. Whether carrying out simulations within the live view mode on a small laptop or performing complex and large simulations on supercomputers or Cloud computers, QwikMD uses the same steps and user interface. QwikMD is freely available by download on group and personal computers. It is also available on the cloud at Amazon Web Services.

The knowledge of multiple conformational states is a prerequisite to understand the function of membrane transport proteins. Unfortunately, the determination of detailed atomic structures for all these functionally important conformational states with conventional high-resolution approaches is often difficult and unsuccessful. In some cases, biophysical and biochemical approaches can provide important complementary structural information that can be exploited with the help of advanced computational methods to derive structural models of specific conformational states. In particular, functional and spectroscopic measurements in combination with site-directed mutations constitute one important source of information to obtain these mixed-resolution structural models. A very common problem with this strategy, however, is the difficulty to simultaneously integrate all the information from multiple independent experiments involving different mutations or chemical labels to derive a unique structural model consistent with the data. To resolve this issue, a novel restrained molecular dynamics structural refinement method is developed to simultaneously incorporate multiple experimentally determined constraints (e.g., engineered metal bridges or spin-labels), each treated as an individual molecular fragment with all atomic details. The internal structure of each of the molecular fragments is treated realistically, while there is no interaction between different molecular fragments to avoid unphysical steric clashes. The information from all the molecular fragments is exploited simultaneously to constrain the backbone to refine a three-dimensional model of the conformational state of the protein. The method is illustrated by refining the structure of the voltage-sensing domain (VSD) of the Kv1.2 potassium channel in the resting state and by exploring the distance histograms between spin-labels attached to T4 lysozyme. The resulting VSD structures are in good agreement with the consensus model of the resting state VSD and the spin-spin distance histograms from ESR/DEER experiments on T4 lysozyme are accurately reproduced.

Integrins may undergo large conformational changes during activation, but the dynamic processes and pathways remain poorly understood. We used molecular dynamics to simulate forced unbending of a complete integrin α(v)β₃ ectodomain in both unliganded and liganded forms. Pulling the head of the integrin readily induced changes in the integrin from a bent to an extended conformation. Pulling at a cyclic RGD ligand bound to the integrin head also extended the integrin, suggesting that force can activate integrins. Interactions at the interfaces between the hybrid and β tail domains and between the hybrid and epidermal growth factor 4 domains formed the major energy barrier along the unbending pathway, which could be overcome spontaneously in ~1 µs to yield a partially-extended conformation that tended to rebend. By comparison, a fully-extended conformation was stable. A newly-formed coordination between the α(v) Asp457 and the α-genu metal ion might contribute to the stability of the fully-extended conformation. These results reveal the dynamic processes and pathways of integrin conformational changes with atomic details and provide new insights into the structural mechanisms of integrin activation.

Interactive molecular dynamics, a new modeling tool for rapid investigation of the physical mechanisms of biological processes at the atomic level, is applied to study selectivity and regulation of the membrane channel protein GlpF and the enzyme glycerol kinase. These proteins facilitate the first two steps of Escherichia coli glycerol metabolism. Despite their different function and architecture the proteins are found to employ common mechanisms for substrate selectivity: an induced geometrical fit by structurally homologous binding sites and an induced rapid dipole moment reversal. Competition for hydrogen bonding sites with water in both proteins is critical for substrate motion. In glycerol kinase, it is shown that the proposed domain motion prevents competition with water, in turn regulating the binding of glycerol.

Human immunodeficiency virus type 1 (HIV-1) infection is highly dependent on its capsid. The capsid is a large container, made of ∼1,300 proteins with altogether 4 million atoms. Although the capsid proteins are all identical, they nevertheless arrange themselves into a largely asymmetric structure made of hexamers and pentamers. The large number of degrees of freedom and lack of symmetry pose a challenge to studying the chemical details of the HIV capsid. Simulations of over 64 million atoms for over 1 μs allow us to conduct a comprehensive study of the chemical-physical properties of an empty HIV-1 capsid, including its electrostatics, vibrational and acoustic properties, and the effects of solvent (ions and water) on the capsid. The simulations reveal critical details about the capsid with implications to biological function.

Two structure determination methods, based on the molecular dynamics flexible fitting (MDFF) paradigm, are presented that resolve sub-5 Å cryo-electron microscopy (EM) maps with either single structures or ensembles of such structures. The methods, denoted cascade MDFF and resolution exchange MDFF, sequentially re-refine a search model against a series of maps of progressively higher resolutions, which ends with the original experimental resolution. Application of sequential re-refinement enables MDFF to achieve a radius of convergence of ~25 Å demonstrated with the accurate modeling of β-galactosidase and TRPV1 proteins at 3.2 Å and 3.4 Å resolution, respectively. The MDFF refinements uniquely offer map-model validation and B-factor determination criteria based on the inherent dynamics of the macromolecules studied, captured by means of local root mean square fluctuations. The MDFF tools described are available to researchers through an easy-to-use and cost-effective cloud computing resource on Amazon Web Services.

Oseltamivir (Tamiflu) is currently the frontline antiviral drug employed to fight the flu virus in infected individuals by inhibiting neuraminidase, a flu protein responsible for the release of newly synthesized virions. However, oseltamivir resistance has become a critical problem due to rapid mutation of the flu virus. Unfortunately, how mutations actually confer drug resistance is not well understood. In this study, we employ molecular dynamics (MD) and steered molecular dynamics (SMD) simulations, as well as graphics processing unit (GPU)-accelerated electrostatic mapping, to uncover the mechanism behind point mutation induced oseltamivir-resistance in both H5N1 "avian" and H1N1pdm "swine" flu N1-subtype neuraminidases. The simulations reveal an electrostatic binding funnel that plays a key role in directing oseltamivir into and out of its binding site on N1 neuraminidase. The binding pathway for oseltamivir suggests how mutations disrupt drug binding and how new drugs may circumvent the resistance mechanisms.

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.

Maturation of HIV-1 particles encompasses a complex morphological transformation of Gag via an orchestrated series of proteolytic cleavage events. A longstanding question concerns the structure of the C-terminal region of CA and the peptide SP1 (CA-SP1), which represents an intermediate during maturation of the HIV-1 virus. By integrating NMR, cryo-EM, and molecular dynamics simulations, we show that in CA-SP1 tubes assembled in vitro, which represent the features of an intermediate assembly state during maturation, the SP1 peptide exists in a dynamic helix-coil equilibrium, and that the addition of the maturation inhibitors Bevirimat and DFH-055 causes stabilization of a helical form of SP1. Moreover, the maturation-arresting SP1 mutation T8I also induces helical structure in SP1 and further global dynamical and conformational changes in CA. Overall, our results show that dynamics of CA and SP1 are critical for orderly HIV-1 maturation and that small molecules can inhibit maturation by perturbing molecular motions.

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

NR3 subtype glutamate receptors have a unique developmental expression profile, but are the least well-characterized members of the NMDA receptor gene family, which have key roles in synaptic plasticity and brain development. Using ligand binding assays, crystallographic analysis, and all atom MD simulations, we investigate mechanisms underlying the binding by NR3A and NR3B of glycine and D-serine, which are candidate neurotransmitters for NMDA receptors containing NR3 subunits. The ligand binding domains of both NR3 subunits adopt a similar extent of domain closure as found in the corresponding NR1 complexes, but have a unique loop 1 structure distinct from that in all other glutamate receptor ion channels. Within their ligand binding pockets, NR3A and NR3B have strikingly different hydrogen bonding networks and solvent structures from those found in NR1, and fail to undergo a conformational rearrangement observed in NR1 upon binding the partial agonist ACPC. MD simulations revealed numerous interdomain contacts, which stabilize the agonist-bound closed-cleft conformation, and a novel twisting motion for the loop 1 helix that is unique in NR3 subunits.

Phosphatidylinositol 4-kinase IIα (PI4KIIα), a membrane-associated PI kinase, plays a central role in cell signalling and trafficking. Its kinase activity critically depends on palmitoylation of its cysteine-rich motif (-CCPCC-) and is modulated by the membrane environment. Lack of atomic structure impairs our understanding of the mechanism regulating kinase activity. Here we present the crystal structure of human PI4KIIα in ADP-bound form. The structure identifies the nucleotide-binding pocket that differs notably from that found in PI3Ks. Two structural insertions, a palmitoylation insertion and an RK-rich insertion, endow PI4KIIα with the 'integral' membrane-binding feature. Molecular dynamics simulations, biochemical and mutagenesis studies reveal that the palmitoylation insertion, containing an amphipathic helix, contributes to the PI-binding pocket and anchors PI4KIIα to the membrane, suggesting that fluctuation of the palmitoylation insertion affects PI4KIIα's activity. We conclude from our results that PI4KIIα's activity is regulated indirectly through changes in the membrane environment.

Helicases play key roles in genome maintenance, yet it remains elusive how these enzymes change conformations and how transitions between different conformational states regulate nucleic acid reshaping. Here, we developed a computational technique combining structural bioinformatics approaches and atomic-level free-energy simulations to characterize how the Escherichia coli DNA repair enzyme UvrD changes its conformation at the fork junction to switch its function from unwinding to rezipping DNA. The lowest free-energy path shows that UvrD opens the interface between two domains, allowing the bound ssDNA to escape. The simulation results predict a key metastable 'tilted' state during ssDNA strand switching. By simulating FRET distributions with fluorophores attached to UvrD, we show that the new state is supported quantitatively by single-molecule measurements. The present study deciphers key elements for the 'hyper-helicase' behavior of a mutant and provides an effective framework to characterize directly structure-function relationships in molecular machines.

A molecular model of the swine influenza A/H1N1 ( also called H1N1pdm) type-I neuraminidase was built using the pathogenic avian H5N1 type-I neuraminidase as a basis, due to the higher sequence identity between A/H1N1 and H5N1 (91.47%) compared to Spanish H1N1 (88.37%) neuraminidase. All-atom molecular dynamics (MD) simulations of all three neuraminidases were performed, either as apo-structures or with commercial antiviral drugs Tamiflu or Relenza separately bound; the simulations allowed for the identification of both conserved and unique drug-protein interactions across all three proteins. Specifically, conserved networks of hydrogen bonds stabilizing the drugs in the sialic acid binding site of the simulated neuraminidases are analyzed, providing insight into how disruption due to mutations may lead to increased drug resistance. In addition, a possible mechanism through which the residue 294 mutation acquires drug resistance is proposed by mapping the mutation site onto an electrostatic pathway which may play a role in controlling drug access to the binding pocket of neuraminidase, establishing a starting point for further investigations of neuraminidase drug resistance.

Nuclear pore complexes (NPCs) form gateways for material transfer across the nuclear envelope of eukaryotic cells. Disordered proteins, rich in phenylalanine-glycine repeat motifs (FG-nups), form the central transport channel. Understanding how nups are arranged in the interior of the NPC may explain how NPC functions as a selectivity filter for transport of large molecules and a sieve-like filter for diffusion of small molecules (<9 nm or 40 kDa). We employed molecular dynamics to model the structures formed by various assemblies of one kind of nup, namely the 609-aa-long FG domain of Nsp1 (Nsp1-FG). The simulations started from different initial conformations and geometrical arrangements of Nsp1-FGs. In all cases Nsp1-FGs collectively formed brush-like structures with bristles made of bundles of 2-27 nups, however, the bundles being cross-linked through single nups leaving one bundle and joining a nearby one. The degree of cross-linking varies with different initial nup conformations and arrangements. Structural analysis reveals that FG-repeats of the nups not only involve formation of bundle structures, but are abundantly present in cross-linking regions where the epitopes of FG-repeats are highly accessible. Large molecules that are assisted by transport factors (TFs) are selectively transported through NPC apparently by binding to FG-nups through populated FG-binding pockets on the TF surface. Therefore, our finding suggests that TFs bind concertedly to multiple FGs in cross-linking regions and break-up the bundles to create wide pores for themselves and their cargoes to pass. In addition, the cross-linking between Nsp1-FG bundles, arising from simulations, is found to set a molecular size limit of <9 nm (40 kDa) for passive diffusion of molecules. Our simulations suggest that the NPC central channel, near the periphery where tethering of nups is dominant, features brush-like moderately cross-linked bundles, but in the central region, where tethering loses its effect, features a sieve-like structure of bundles and frequent cross-links.

The protein-conducting channel, or translocon, is an evolutionarily conserved complex that allows nascent proteins to cross a cellular membrane or integrate into it. The crystal structure of an archaeal translocon, the SecY complex, revealed that two elements contribute to sealing the channel: a small "plug" domain blocking the periplasmic region of the channel, and a pore ring composed of six hydrophobic residues acting as a constriction point at the channel's center. To determine the independent functions of these two elements, we have performed molecular dynamics simulations of the native channel as well as of two recently structurally resolved mutants in which portions of their plugs were deleted. We find that in the mutants, the instability in the plug region leads to a concomitant increase in flexibility of the pore ring. The instability is quantified by the rate of water permeation in each system as well as by the force required for oligopeptide translocation. Through a novel simulation in which the interactions between the plug and water were independently controlled, we find that the role of the plug in stabilizing the pore ring is significantly more important than its role as a purely steric barrier.

DNA methylation plays an essential role in transcriptional control of organismal development in epigenetics, from turning off a specific gene to inactivation of entire chromosomes. While the biological function of DNA methylation is becoming increasingly clear, the mechanism of methylation-induced gene regulation is still poorly understood. Through single-molecule force experiments and simulation we investigated the effects of methylation on strand separation of DNA, a crucial step in gene expression. Molecular force assay and single-molecule force spectroscopy revealed a strong methylation dependence of strand separation. Methylation is observed to either inhibit or facilitate strand separation, depending on methylation level and sequence context. Molecular dynamics simulations provided a detailed view of methylation effects on strand separation, suggesting the underlying physical mechanism. According to our study, methylation in epigenetics may regulate gene expression not only through mechanisms already known but also through changing mechanical properties of DNA.

Sequence-specific DNA-binding proteins must quickly and reliably localize specific target sites on DNA. This search process has been well characterized for monomeric proteins, but it remains poorly understood for systems that require assembly into dimers or oligomers at the target site. We present a single-molecule study of the target-search mechanism of protelomerase TelK, a recombinase-like protein that is only active as a dimer. We show that TelK undergoes 1D diffusion on non-target DNA as a monomer, and it immobilizes upon dimerization even in the absence of a DNA target site. We further show that dimeric TelK condenses non-target DNA, forming a tightly bound nucleoprotein complex. Together with theoretical calculations and molecular dynamics simulations, we present a novel target-search model for TelK, which may be generalizable to other dimer and oligomer-active proteins.