RNA folding using quantum computers.

The 3-dimensional fold of an RNA molecule is largely determined by patterns of intramolecular hydrogen bonds between bases. Predicting the base pairing network from the sequence, also referred to as RNA secondary structure prediction or RNA folding, is a nondeterministic polynomial-time (NP)-complete computational problem. The structure of the molecule is strongly predictive of its functions and biochemical properties, and therefore the ability to accurately predict the structure is a crucial tool for biochemists. Many methods have been proposed to efficiently sample possible secondary structure patterns. Classic approaches employ dynamic programming, and recent studies have explored approaches inspired by evolutionary and machine learning algorithms. This work demonstrates leveraging quantum computing hardware to predict the secondary structure of RNA. A Hamiltonian written in the form of a Binary Quadratic Model (BQM) is derived to drive the system toward maximizing the number of consecutive base pairs while jointly maximizing the average length of the stems. A Quantum Annealer (QA) is compared to a Replica Exchange Monte Carlo (REMC) algorithm programmed with the same objective function, with the QA being shown to be highly competitive at rapidly identifying low energy solutions. The method proposed in this study was compared to three algorithms from literature and, despite its simplicity, was found to be competitive on a test set containing known structures with pseudoknots.

Coexistence of Lipid Phases Stabilizes Interstitial Water in the Outer Layer of Mammalian Skin.

The lipid matrix in the outer layer of mammalian skin, the stratum corneum, has been previously investigated by multiple biophysical techniques aimed at identifying hydrophilic and lipophilic pathways of permeation. Although consensus is developing over the microscopic structure of the lipid matrix, no molecular-resolution model describes the permeability of all chemical species simultaneously. Using molecular dynamics simulations of a model mixture of skin lipids, the self-assembly of the lipid matrix lamellae has been studied. At higher humidity, the resulting lamellar phase is maintained by partitioning excess water into isolated droplets of controlled size and spatial distribution. The droplets may fuse together to form intralamellar water channels, thereby providing a pathway for the permeation of hydrophilic species. These results reconcile competing data on the outer skin's structure and broaden the scope of molecular-based methods to improve the safety of topical products and to advance transdermal drug delivery.

Controlling Association and Separation of Gold Nanoparticles with Computationally Designed Zinc-Coordinating Proteins.

Functionalization of nanoparticles with biopolymers has yielded a wide range of structured and responsive hybrid materials. DNA provides the ability to program length and recognition using complementary oligonucleotide sequences. Nature more often leverages the versatility of proteins, however, where structure, assembly, and recognition are more subtle to engineer. Herein, a protein was computationally designed to present multiple Zn2+ coordination sites and cooperatively self-associate to form an antiparallel helical homodimer. Each subunit was unstructured in the absence of Zn2+ or when the cation was sequestered with a chelating agent. When bound to the surface of gold nanoparticles via cysteine, the protein provided a reversible molecular linkage between particles. Nanoparticle association and changes in interparticle separation were monitored by redshifts in the surface plasmon resonance (SPR) band and by transmission electron microscopy (TEM). Titrations with Zn2+ revealed sigmoidal transitions at submicromolar concentrations. The metal-ion concentration required to trigger association varied with the loading of the proteins on the nanoparticles, the solution ionic strength, and the cation employed. Specifying the number of helical (heptad) repeat units conferred control over protein length and nanoparticle separation. Two different length proteins were designed via extension of the helical structure. TEM and extinction measurements revealed distributions of nanoparticle separations consistent with the expected protein structures. Nanoparticle association, interparticle separation, and SPR properties can be tuned using computationally designed proteins, where protein structure, folding, length, and response to molecular species such as Zn2+ can be engineered.

Self-interrupted synthesis of sterically hindered aliphatic polyamide dendrimers.

2,2-Bis(azidomethyl)propionic acid was prepared in four steps and 85% yield from the commercially available 2,2-bis(hydroxymethyl)propionic acid and used as the starting building block for the divergent, convergent, and double-stage convergent-divergent iterative methods for the synthesis of dendrimers and dendrons containing ethylenediamine (EDA), piperazine (PPZ), and methyl 2,2-bis(aminomethyl)propionate (COOMe) cores. These cores have the same multiplicity but different conformations. A diversity of synthetic methods were used for the synthesis of dendrimers and dendrons. Regardless of the method used, a self-interruption of the synthesis was observed at generation 4 for the dendrimer with an EDA core and at generation 5 for the one with a PPZ core, whereas for the COOMe core, self-interruption was observed at generation 6 dendron, which is equivalent to generation 5 dendrimer. Molecular modeling and molecular-dynamics simulations demonstrated that the observed self-interruption is determined by the backfolding of the azide groups at the periphery of the dendrimer. The latter conformation inhibits completely the heterogeneous hydrogenation of the azide groups catalyzed by 10% Pd/carbon as well as homogeneous hydrogenation by the Staudinger method. These self-terminated polyamide dendrimers are enzymatically and hydrolytically stable and also exhibit antimicrobial activity. Thus, these nanoscale constructs open avenues for biomedical applications.

Canine CNGA3 Gene Mutations Provide Novel Insights into Human Achromatopsia-Associated Channelopathies and Treatment.

Cyclic nucleotide-gated (CNG) ion channels are key mediators underlying signal transduction in retinal and olfactory receptors. Genetic defects in CNGA3 and CNGB3, encoding two structurally related subunits of cone CNG channels, lead to achromatopsia (ACHM). ACHM is a congenital, autosomal recessive retinal disorder that manifests by cone photoreceptor dysfunction, severely reduced visual acuity, impaired or complete color blindness and photophobia. Here, we report the first canine models for CNGA3-associated channelopathy caused by R424W or V644del mutations in the canine CNGA3 ortholog that accurately mimic the clinical and molecular features of human CNGA3-associated ACHM. These two spontaneous mutations exposed CNGA3 residues essential for the preservation of channel function and biogenesis. The CNGA3-R424W results in complete loss of cone function in vivo and channel activity confirmed by in vitro electrophysiology. Structural modeling and molecular dynamics (MD) simulations revealed R424-E306 salt bridge formation and its disruption with the R424W mutant. Reversal of charges in a CNGA3-R424E-E306R double mutant channel rescued cGMP-activated currents uncovering new insights into channel gating. The CNGA3-V644del affects the C-terminal leucine zipper (CLZ) domain destabilizing intersubunit interactions of the coiled-coil complex in the MD simulations; the in vitro experiments showed incompetent trimeric CNGA3 subunit assembly consistent with abnormal biogenesis of in vivo channels. These newly characterized large animal models not only provide a valuable system for studying cone-specific CNG channel function in health and disease, but also represent prime candidates for proof-of-concept studies of CNGA3 gene replacement therapy for ACHM patients.

Molecular dynamics simulations of cholesterol-rich membranes using a coarse-grained force field for cyclic alkanes.

The architecture of a biological membrane hinges upon the fundamental fact that its properties are determined by more than the sum of its individual components. Studies on model membranes have shown the need to characterize in molecular detail how properties such as thickness, fluidity, and macroscopic bending rigidity are regulated by the interactions between individual molecules in a non-trivial fashion. Simulation-based approaches are invaluable to this purpose but are typically limited to short sampling times and model systems that are often smaller than the required properties. To alleviate both limitations, the use of coarse-grained (CG) models is nowadays an established computational strategy. We here present a new CG force field for cholesterol, which was developed by using measured properties of small molecules, and can be used in combination with our previously developed force field for phospholipids. The new model performs with precision comparable to atomistic force fields in predicting the properties of cholesterol-rich phospholipid bilayers, including area per lipid, bilayer thickness, tail order parameter, increase in bending rigidity, and propensity to form liquid-ordered domains in ternary mixtures. We suggest the use of this model to quantify the impact of cholesterol on macroscopic properties and on microscopic phenomena involving localization and trafficking of lipids and proteins on cellular membranes.

Dehydration of multilamellar fatty acid membranes: towards a computational model of the stratum corneum.

The level of hydration controls the cohesion between apposed lamellae of saturated free fatty acids found in the lipid matrix of stratum corneum, the outermost layer of mammalian skin. This multilamellar lipid matrix is highly impermeable to water and ions, so that the local hydration shell of its fatty acids may not always be in equilibrium with the acidity and relative humidity, which significantly change over a course of days during skin growth. The homeostasis of the stratum corneum at each moment of its growth likely requires a balance between two factors, which affect in opposite ways the diffusion of hydrophilic species through the stratum corneum: (i) an increase in water order as the lipid lamellae come in closer contact, and (ii) a decrease in water order as the fraction of charged fatty acids is lowered by pH. Herein molecular dynamics simulations are employed to estimate the impact of both effects on water molecules confined between lamellae of fatty acids. Under conditions where membrane undulations are energetically favorable, the charged fatty acids are able to sequester cations around points of contact between lamellae that are fully dehydrated, while essentially maintaining a multilamellar structure for the entire system. This observation suggests that the undulations of the fatty acid lamellae control the diffusion of hydrophilic species through the water phase by altering the positional and rotational order of water molecules in the embedded/occluded "droplets."

Engineering complementary hydrophobic interactions to control ß-hairpin peptide self-assembly, network branching, and hydrogel properties.

The MAX1 ß-hairpin peptide (VKVKVKVK-V(D)PPT-KVKVKVKV-NH2) has been shown to form nanofibrils having a cross-section of two folded peptides forming a hydrophobic, valine-rich core, and the polymerized fibril exhibits primarily ß-sheet hydrogen bonding.1-7 These nanofibrils form hydrogel networks through fibril entanglements as well as fibril branching.8 Fibrillar branching in MAX1 hydrogel networks provide the ability to flow under applied shear stress and immediately reform a hydrogel solid on cessation of shear. New ß-hairpins were designed to limit branching during nanofibril growth because of steric specificity in the assembled fibril hydrophobic core. The nonturn valines of MAX1 were substituted by 2-naphthylalanine (Nal) and alanine (A) residues, with much larger and smaller side chain volumes, respectively, to obtain LNK1 (Nal)K(Nal)KAKAK-V(D)PPT-KAKAK(Nal)K(Nal)-NH2. LNK1 was targeted to self-associate with a specific "lock and key" complementary packing in the hydrophobic core in order to accommodate the Nal and Ala residue side chains. The experimentally observable manifestation of reduced fibrillar branching in the LNK1 peptide is the lack of solid hydrogel formation after shear in stark contrast to the MAX1 branched fibril system. Molecular dynamics simulations provide a molecular picture of interpeptide interactions within the assembly that is consistent with the branching propensity of MAX1 vs LNK1 and in agreement with experimental observations.

Magnesium-dependent RNA binding to the PA endonuclease domain of the avian influenza polymerase.

Influenza A viruses are highly pathogenic and pose an unpredictable public health danger to humans. An attractive target for developing new antiviral drugs is the PA N-terminal domain (PAN) of influenza polymerase, which is responsible for the endonuclease activity and essential for viral replication. Recently, the crystal structures of the holo form of PAN as well as PAN bound to different inhibitors have been reported, but the potency and selectivity of these inhibitors still need to be improved. New drug design can be guided by a better understanding of the endonuclease activity of PAN. However, this requires the structure of PAN in complex with the host mRNA, which has not been determined yet. In particular, divalent metal ions are known to be essential for RNA cleavage, but it is not clear whether there is either one or two Mg ions in the PAN active site. In the present work, we have modeled the complex of the PAN endonuclease domain with the host mRNA in the presence of either one or two Mg(2+) by using all-atom molecular dynamics. These simulations identify crucial interactions between the enzyme and the nucleic acid. Moreover, they validate a previous hypothesis that a second metal ion binds in the presence of the RNA substrate and therefore support a two-metal ion mechanism, in which K134 decreases the pKa of the nucleophilic water. Nevertheless, at low Mg concentrations an alternative, one-metal ion mechanism is possible, with K137 as the catalytic lysine and H41 as the general base, rationalizing previous unexpected mutagenesis results. The RNA-enzyme interactions determined here could likely be used to design more specific endonuclease inhibitors to fight influenza viral infections.

Origins of the helical wrapping of phenyleneethynylene polymers about single-walled carbon nanotubes.

The highly charged, conjugated polymer poly[p-{2,5-bis(3-propoxysulfonicacidsodiumsalt)}phenylene]ethynylene (PPES) has been shown to wrap single-wall carbon nanotubes (SWNTs), adopting a robust helical superstructure. Surprisingly, PPES adopts a helical rather than a linear conformation when adhered to SWNTs. The complexes formed by PPES and related polymers upon helical wrapping of a SWNT are investigated using atomistic molecular dynamics (MD) simulations in the presence and absence of aqueous solvent. In simulations of the PPES/SWNT system in an aqueous environment, PPES spontaneously takes on a helical conformation. A potential of mean force, ΔA(ξ), is calculated as a function of ξ, the component of the end-to-end vector of the polymer chain projected on the SWNT axis; ξ is a monotonic function of the polymer's helical pitch. ΔA(ξ) provides a means to quantify the relative free energies of helical conformations of the polymer when wrapped about the SWNT. The aqueous system possesses a global minimum in ΔA(ξ) at the experimentally observed value of the helical pitch. The presence of this minimum is associated with preferred side chain conformations, where the side chains adopt conformations that provide van der Waals contact between the tubes and the aliphatic components of the side chains, while exposing the anionic sulfonates for aqueous solvation. The simulations provide a free energy estimate of a 0.2 kcal/mol/monomer preference for the helical over the linear conformation of the PPES/SWNT system in an aqueous environment.

Computational design of a protein crystal.

Protein crystals have catalytic and materials applications and are central to efforts in structural biology and therapeutic development. Designing predetermined crystal structures can be subtle given the complexity of proteins and the noncovalent interactions that govern crystallization. De novo protein design provides an approach to engineer highly complex nanoscale molecular structures, and often the positions of atoms can be programmed with sub-Å precision. Herein, a computational approach is presented for the design of proteins that self-assemble in three dimensions to yield macroscopic crystals. A three-helix coiled-coil protein is designed de novo to form a polar, layered, three-dimensional crystal having the P6 space group, which has a "honeycomb-like" structure and hexameric channels that span the crystal. The approach involves: (i) creating an ensemble of crystalline structures consistent with the targeted symmetry; (ii) characterizing this ensemble to identify "designable" structures from minima in the sequence-structure energy landscape and designing sequences for these structures; (iii) experimentally characterizing candidate proteins. A 2.1 Å resolution X-ray crystal structure of one such designed protein exhibits sub-Å agreement [backbone root mean square deviation (rmsd)] with the computational model of the crystal. This approach to crystal design has potential applications to the de novo design of nanostructured materials and to the modification of natural proteins to facilitate X-ray crystallographic analysis.

Theoretical and computational protein design.

From exponentially large numbers of possible sequences, protein design seeks to identify the properties of those that fold to predetermined structures and have targeted structural and functional properties. The interactions that confer structure and function involve intermolecular forces and large numbers of interacting amino acids. As a result, the identification of sequences can be subtle and complex. Sophisticated methods for characterizing sequences consistent with a particular structure have been developed, assisting the design of novel proteins. Developments in such computational protein design are discussed, along with recent accomplishments, ranging from the redesign of existing proteins to the design of new functionalities and nonbiological applications.

Electrostatic polarization is crucial for reproducing pKa shifts of carboxylic residues in Turkey ovomucoid third domain.

We have computed pKa shifts for carboxylic residues of the serine protease inhibitor turkey ovomucoid third domain (residues Asp7, Glu10, Glu19, Asp27, and Glu43). Both polarizable and nonpolarizable empirical force fields were employed. Hydration was represented by the surface generalized Born and Poisson-Boltzmann continuum model. The calculations were carried out in the most physically straightforward fashion, by directly comparing energies of the protonated and deprotonated protein forms, without any additional parameter fitting or adjustment. Our studies have demonstrated that (i) the Poisson-Boltzmann solvation model is more than adequate in reproducing pKa shifts, most likely due to its intrinsically many-body formalism; (ii) explicit treatment of electrostatic polarization included in our polarizable force field (PFF) calculations appears to be crucial in reproducing the acidity constant shifts. The average error of the PFF results was found to be as low as 0.58 pKa units, with the best fixed-charges average deviation being 3.28 units. Therefore, the pKa shifts phenomena and the governing electrostatics are clearly many-body controlled in their intrinsic nature; (iii) our results confirm previously reported conclusions that pKa shifts for protein residues are controlled by the immediate environment of the residues in question, as opposed to long-range interactions in proteins. We are confident that our confirmation of the importance of explicit inclusion of polarization in empirical force fields for protein studies will be useful far beyond the immediate goal of accurate calculation of acidity constants.