Theoretical Description of the Primary Proton-Coupled Electron Transfer Reaction in the Cytochrome bc1 Complex.

The cytochrome bc1 complex is a transmembrane enzymatic protein complex that plays a central role in cellular energy production and is present in both photosynthetic and respiratory chain organelles. Its reaction mechanism is initiated by the binding of a quinol molecule to an active site, followed by a series of charge transfer reactions between the quinol and protein subunits. Previous work hypothesized that the primary reaction was a concerted proton-coupled electron transfer (PCET) reaction because of the apparent absence of intermediate states associated with single proton or electron transfer reactions. In the present study, the kinetics of the primary bc1 complex PCET reaction is investigated with a vibronically nonadiabatic PCET theory in conjunction with all-atom molecular dynamics simulations and electronic structure calculations. The computed rate constants and relatively high kinetic isotope effects are consistent with experimental measurements on related biomimetic systems. The analysis implicates a concerted PCET mechanism with significant hydrogen tunneling and nonadiabatic effects in the bc1 complex. Moreover, the employed theoretical framework is shown to serve as a general strategy for describing PCET reactions in bioenergetic systems.

Scalable molecular dynamics on CPU and GPU architectures with NAMD.

NAMDis a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

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.

In Situ Conformational Changes of the Escherichia coli Serine Chemoreceptor in Different Signaling States.

Tsr, the serine chemoreceptor in Escherichia coli, transduces signals from a periplasmic ligand-binding site to its cytoplasmic tip, where it controls the activity of the CheA kinase. To function, Tsr forms trimers of homodimers (TODs), which associate in vivo with the CheA kinase and CheW coupling protein. Together, these proteins assemble into extended hexagonal arrays. Here, we use cryo-electron tomography and molecular dynamics simulation to study Tsr in the context of a near-native array, characterizing its signaling-related conformational changes at both the individual dimer and the trimer level. In particular, we show that individual Tsr dimers within a trimer exhibit asymmetric flexibilities that are a function of the signaling state, highlighting the effect of their different protein interactions at the receptor tips. We further reveal that the dimer compactness of the Tsr trimer changes between signaling states, transitioning at the glycine hinge from a compact conformation in the kinase-OFF state to an expanded conformation in the kinase-ON state. Hence, our results support a crucial role for the glycine hinge: to allow the receptor flexibility necessary to achieve different signaling states while also maintaining structural constraints imposed by the membrane and extended array architecture.IMPORTANCE In Escherichia coli, membrane-bound chemoreceptors, the histidine kinase CheA, and coupling protein CheW form highly ordered chemosensory arrays. In core signaling complexes, chemoreceptor trimers of dimers undergo conformational changes, induced by ligand binding and sensory adaptation, which regulate kinase activation. Here, we characterize by cryo-electron tomography the kinase-ON and kinase-OFF conformations of the E. coli serine receptor in its native array context. We found distinctive structural differences between the members of a receptor trimer, which contact different partners in the signaling unit, and structural differences between the ON and OFF signaling complexes. Our results provide new insights into the signaling mechanism of chemoreceptor arrays and suggest an important functional role for a previously postulated flexible region and glycine hinge in the receptor molecule.

Pro-Nifuroxazide Self-Assembly Leads to Triggerable Nanomedicine for Anti-cancer Therapy.

Transcription factor STAT3 has been shown to regulate genes that are involved in stem cell self-renewal and thus represents a novel therapeutic target of great biological significance. However, many small-molecule agents with potential effects through STAT3 modulation in cancer therapy lack aqueous solubility and high off-target toxicity, hence impeding efficient bioavailability and activity. This work, for the first time, reports a prodrug-based strategy for selective and safer delivery of STAT3 inhibitors designed toward metastatic and drug-resistant breast cancer. We have synthesized a novel lipase-labile SN-2 phospholipid prodrug from a clinically investigated STAT3 inhibitor, nifuroxazide (Pro-nifuroxazide), which can be regioselectively cleaved by the membrane-abundant enzymes in cancer cells. Pro-nifuroxazide self-assembled to sub 20 nm nanoparticles (NPs), and the cytotoxic ability was screened in ER(+)-MCF-7 and ER(-)-MD-MB231 cells at 48-72 h using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetra-zolium bromide proliferation assay. Results indicated that Pro-nifuroxazide NPs are multifold more effective toward inhibiting cancer cells in a time-dependent manner compared to parent nifuroxazide. A remarkable improvement in the local concentration of drugs to as high as ∼240 fold when assembled into NPs is presumably the reason for this functional improvement. We also introduced molecular dynamics simulations to generate Pro-nifuroxazide nano-assembly, as a model assembly from triggerable anti-cancer drugs, to provide molecular insights correlating physicochemical and anti-cancer properties. In silico properties of Pro-nifuroxazide including size, chemistry of NPs and membrane interactions with individual molecules could be validated by in vitro functional activities in cells of breast cancer origin. The in vivo anti-cancer efficiencies of Pro-nifuroxazide NPs in nude mice xenografts with MCF-7 revealed remarkable growth inhibition of as high as 400%. Histopathological analysis corroborated these findings to show significantly high nuclear fragmentation and retracted cytoplasm. Immunostaining on tumor section demonstrated a significantly lower level of pSTAT-3 by Pro-nifuroxazide NP treatment, establishing the inhibition of STAT-3 phosphorylation. Our strategy for the first time proposes a translatable prodrug agent self-assembled into NPs and demonstrates remarkable enhancement in IC50, induced apoptosis, and reduced cancer cell population through STAT-3 inhibition via reduced phosphorylation.

Fast pressure-jump all-atom simulations and experiments reveal site-specific protein dehydration-folding dynamics.

As theory and experiment have shown, protein dehydration is a major contributor to protein folding. Dehydration upon folding can be characterized directly by all-atom simulations of fast pressure drops, which create desolvated pockets inside the nascent hydrophobic core. Here, we study pressure-drop refolding of three λ-repressor fragment (λ6-85) mutants computationally and experimentally. The three mutants report on tertiary structure formation via different fluorescent helix-helix contact pairs. All-atom simulations of pressure drops capture refolding and unfolding of all three mutants by a similar mechanism, thus validating the nonperturbative nature of the fluorescent contact probes. Analysis of simulated interprobe distances shows that the α-helix 1-3 pair distance displays a slower characteristic time scale than the 1-2 or 3-2 pair distance. To see whether slow packing of α-helices 1 and 3 is reflected in the rate-limiting folding step, fast pressure-drop relaxation experiments captured refolding on a millisecond time scale. These experiments reveal that refolding monitored by 1-3 contact formation indeed is much slower than when monitored by 1-2 or 3-2 contact formation. Unlike the case of the two-state folder [three-α-helix bundle (α3D)], whose drying and core formation proceed in concert, λ6-85 repeatedly dries and rewets different local tertiary contacts before finally forming a solvent-excluded core, explaining the non-two-state behavior observed during refolding in molecular dynamics simulations. This work demonstrates that proteins can explore desolvated pockets and dry globular states numerous times before reaching the native conformation.

Mechanism for the Regulated Control of Bacterial Transcription Termination by a Universal Adaptor Protein.

NusG/Spt5 proteins are the only transcription factors utilized by all cellular organisms. In enterobacteria, NusG antagonizes the transcription termination activity of Rho, a hexameric helicase, during the synthesis of ribosomal and actively translated mRNAs. Paradoxically, NusG helps Rho act on untranslated transcripts, including non-canonical antisense RNAs and those arising from translational stress; how NusG fulfills these disparate functions is unknown. Here, we demonstrate that NusG activates Rho by assisting helicase isomerization from an open-ring, RNA-loading state to a closed-ring, catalytically active translocase. A crystal structure of closed-ring Rho in complex with NusG reveals the physical basis for this activation and further explains how Rho is excluded from translationally competent RNAs. This study demonstrates how a universally conserved transcription factor acts to modulate the activity of a ring-shaped ATPase motor and establishes how the innate sequence bias of a termination factor can be modulated to silence pervasive, aberrant transcription.

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.

Free-energy simulations reveal molecular mechanism for functional switch of a DNA helicase.

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.

Molecular mechanism of extreme mechanostability in a pathogen adhesin.

High resilience to mechanical stress is key when pathogens adhere to their target and initiate infection. Using atomic force microscopy-based single-molecule force spectroscopy, we explored the mechanical stability of the prototypical staphylococcal adhesin SdrG, which targets a short peptide from human fibrinogen ß. Steered molecular dynamics simulations revealed, and single-molecule force spectroscopy experiments confirmed, the mechanism by which this complex withstands forces of over 2 nanonewtons, a regime previously associated with the strength of a covalent bond. The target peptide, confined in a screwlike manner in the binding pocket of SdrG, distributes forces mainly toward the peptide backbone through an intricate hydrogen bond network. Thus, these adhesins can attach to their target with exceptionally resilient mechanostability, virtually independent of peptide side chains.

NAMD goes quantum: an integrative suite for hybrid simulations.

Hybrid methods that combine quantum mechanics (QM) and molecular mechanics (MM) can be applied to studies of reaction mechanisms in locations ranging from active sites of small enzymes to multiple sites in large bioenergetic complexes. By combining the widely used molecular dynamics and visualization programs NAMD and VMD with the quantum chemistry packages ORCA and MOPAC, we created an integrated, comprehensive, customizable, and easy-to-use suite (http://www.ks.uiuc.edu/Research/qmmm). Through the QwikMD interface, setup, execution, visualization, and analysis are streamlined for all levels of expertise.

PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations.

Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as molecular recognition, binding strength, and mechanical and structural stability. Extracting, evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. We have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to investigate biomolecular interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid analysis and data visualization without any additional programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical analysis representation is interactively combined with full mapping of the results on the molecular system through the synergistic connection between PyContact and VMD. We showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, we examine the protein-protein interaction network of the mechanically ultrastable cohesin-dockering complex.

Quenching protein dynamics interferes with HIV capsid maturation.

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.

Constant-pH Molecular Dynamics Simulations for Large Biomolecular Systems.

An increasingly important endeavor is to develop computational strategies that enable molecular dynamics (MD) simulations of biomolecular systems with spontaneous changes in protonation states under conditions of constant pH. The present work describes our efforts to implement the powerful constant-pH MD simulation method, based on a hybrid nonequilibrium MD/Monte Carlo (neMD/MC) technique within the highly scalable program NAMD. The constant-pH hybrid neMD/MC method has several appealing features; it samples the correct semigrand canonical ensemble rigorously, the computational cost increases linearly with the number of titratable sites, and it is applicable to explicit solvent simulations. The present implementation of the constant-pH hybrid neMD/MC in NAMD is designed to handle a wide range of biomolecular systems with no constraints on the choice of force field. Furthermore, the sampling efficiency can be adaptively improved on-the-fly by adjusting algorithmic parameters during the simulation. Illustrative examples emphasizing medium- and large-scale applications on next-generation supercomputing architectures are provided.

CryoEM structure of MxB reveals a novel oligomerization interface critical for HIV restriction.

Human dynamin-like, interferon-induced myxovirus resistance 2 (Mx2 or MxB) is a potent HIV-1 inhibitor. Antiviral activity requires both the amino-terminal region of MxB and protein oligomerization, each of which has eluded structural determination due to difficulties in protein preparation. We report that maltose binding protein-fused, full-length wild-type MxB purifies as oligomers and further self-assembles into helical arrays in physiological salt. Guanosine triphosphate (GTP), but not guanosine diphosphate, binding results in array disassembly, whereas subsequent GTP hydrolysis allows its reformation. Using cryo-electron microscopy (cryoEM), we determined the MxB assembly structure at 4.6 Å resolution, representing the first near-atomic resolution structure in the mammalian dynamin superfamily. The structure revealed previously described and novel MxB assembly interfaces. Mutational analyses demonstrated a critical role for one of the novel interfaces in HIV-1 restriction.

Detection of methylation on dsDNA using nanopores in a MoS2 membrane.

Methylation at the 5-carbon position of the cytosine nucleotide base in DNA has been shown to be a reliable diagnostic biomarker for carcinogenesis. Early detection of methylation and intervention could drastically increase the effectiveness of therapy and reduce the cancer mortality rate. Current methods for detecting methylation involve bisulfite genomic sequencing, which are cumbersome and demand a large sample size of bodily fluids to yield accurate results. Hence, more efficient and cost effective methods are desired. Based on our previous work, we present a novel nanopore-based assay using a nanopore in a MoS2 membrane, and the methyl-binding protein (MBP), MBD1x, to detect methylation on dsDNA. We show that the dsDNA translocation was effectively slowed down using an asymmetric concentration of buffer and explore the possibility of profiling the position of methylcytosines on the DNA strands as they translocate through the 2D membrane. Our findings advance us one step closer towards the possible use of nanopore sensing technology in medical applications such as cancer detection.

Physical properties of the HIV-1 capsid from all-atom molecular dynamics simulations.

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.

Skeletal Dysplasia Mutations Effect on Human Filamins' Structure and Mechanosensing.

Cells' ability to sense mechanical cues in their environment is crucial for fundamental cellular processes, leading defects in mechanosensing to be linked to many diseases. The actin cross-linking protein Filamin has an important role in the conversion of mechanical forces into biochemical signals. Here, we reveal how mutations in Filamin genes known to cause Larsen syndrome and Frontometaphyseal dysplasia can affect the structure and therefore function of Filamin domains 16 and 17. Employing X-ray crystallography, the structure of these domains was first solved for the human Filamin B. The interaction seen between domains 16 and 17 is broken by shear force as revealed by steered molecular dynamics simulations. The effects of skeletal dysplasia associated mutations of the structure and mechanosensing properties of Filamin were studied by combining various experimental and theoretical techniques. The results showed that Larsen syndrome associated mutations destabilize or even unfold domain 17. Interestingly, those Filamin functions that are mediated via domain 17 interactions with other proteins are not necessarily affected as strongly interacting peptide binding to mutated domain 17 induces at least partial domain folding. Mutation associated to Frontometaphyseal dysplasia, in turn, transforms 16-17 fragment from compact to an elongated form destroying the force-regulated domain pair.

Permeability of a Fluid Lipid Bilayer to Short-Chain Alcohols from First Principles.

Computational prediction of membrane permeability to small molecules requires accurate description of both the thermodynamics and kinetics underlying translocation across the lipid bilayer. In this contribution, well-converged, microsecond-long free-energy calculations are combined with a recently developed subdiffusive kinetics framework to describe the membrane permeation of a homologous series of short-tail alcohols, from methanol to 1-butanol, with unprecedented fidelity to the underlying molecular models. While the free-energy profiles exhibit barriers for passage through the center of the bilayer in all cases, the height of these barriers decreases with the length of the aliphatic chain of the alcohol, in quantitative agreement with experimentally determined differential solvation free energies in water and oil. A unique aspect of the subdiffusive model employed herein, which was developed in a previous article, is the determination of a position-dependent fractional order which quantifies the degree to which the motion of the alcohol deviates from classical diffusion along the thickness of the membrane. In the aqueous medium far from the bilayer, this quantity approaches 1.0, the asymptotic limit for purely classical diffusion, whereas it dips below 0.75 near the center of the membrane irrespective of the permeant. Remarkably, the fractional diffusivity near the center of membrane, where its influence on the permeability is the greatest, is similar among the four permeants despite the large difference in molecular weight and lipophilicity between methanol and 1-butanol. The relative permeabilities, which are estimated from the free-energy and fractional diffusivity profiles, are therefore determined predominantly by differences in the former rather than the latter. The predicted relative permeabilities are highly correlated with existing experimental results, albeit they do not agree quantitatively with them. On the other hand, quite unexpectedly, the reported experimental values for the short-tail alcohols are nearly three orders of magnitude lower than the available experimental measurement for water. Plausible explanations for this apparent disagreement between theory and experiment are considered in detail.