Ever since DNA was found with a double-stranded helical structure, people were wondering about how the double strands could be separated apart so that genetic information stored inside the helix could be delivered from generation to generation. A class of protein enzymes known to achieve this function are DNA helicases (see Sept 2006 highlight). Recently, however, researchers have found a novel mechanism that can also serve to split the double strands of DNA apart, namely, by dragging DNA from water to non-polar solvent. This mechanism was demonstrated by means of atomic force microscope-based single molecule force spectroscopy as well as all-atom molecular dynamics simulation, using NAMD, reported in a recent publication. An intriguing idea arising then is that the mechanism may be employed by some DNA helicases, that could achieve the splitting of DNA strands by altering the local-environment from hydrophilic to hydrophobic. More information can be found on our research website.

Since Leeuwenhoek introduced it to biology 300 years ago, the light microscope has brought about multiple discoveries, many achieved through improving over time the instrument's resolution. However, in 1873 Abbe recognized that the resolution has a limit, given by the wave length of light. This limit was considered absolute, until in 1992 Hell suggested a microscope that breaks the limit postulated by Abbe. This is achieved by sending coherent light through two opposing objectives, the resulting interference pattern squeezing the radiation into spots significantly smaller than the light's wavelength. This improvement has already permitted biologists to see a new level of detail in living cells. However, the pattern of light in the Hell microscope is rather complex and certain quantitative measurements require a computational analysis to take advantage of the full benefits of the instrument. Such analysis has been accomplished and validated in a recent study, the validation involving measurements on known systems. The developed numerical algorithms harness the computational power of modern processors, in particular they resort to expoiting the computational power of graphics processors (see also the Oct 2007 highlight). The new methodology combined with the new microscopes opens the avenue to unprecedented measurements in living cells. More information can be found here.

Globins are oxygen-storing proteins, vital to life. In our blood, hemoglobins carry oxygen from our lungs to every cell in our body. In our muscles, myoglobins keep reserves of oxygen to make sure it is available when needed. In some plants, leghemoglobins capture oxygen molecules that would otherwise be harmful to the production of ammonium necessary for the plant's survival. All these globins possess an iron-containig "heme", that grabs on to oxygen for a short time, and share the same protein architecture, despite large variations in their sequences. Since the heme group is buried inside a globin, scientists wondered how oxygen makes its way inside the protein to reach it. Exploring the motion and energetics of globins using the program NAMD researchers learned to gather data that permitted them to visualize, utilizing the VMD software, all the pathways taken by oxygen migrating inside whale myoglobin (see the Aug 2006 highlight and related publication). However, when the researchers turned their attention to the rest of the globin family to compute their oxygen pathways, they found, on their computational spelunking trip, something surprising. Given the conserved architecture of all globins, they expected to see similar oxygen pathways throughout the globin family, but they saw the opposite! Aside from a conserved pocket right at the heme binding site, the distribution of oxygen pathways showed very little similarity from one globin to the next. This result is described in a recent report, which shows that oxygen-pathways are not conserved by evolution, and that their location is not determined by a protein's overall architecture, but rather by its local amino acids. The researchers also learned which amino acids are found more often than others lining oxygen pathways, recognizing that bulky side groups are not hindering, but favoring oxygen passage. More information can be found here.

Modern computers include a massively parallel graphics processing unit (GPU) designed to perform geometry and lighting calculations at blistering speeds. State of the art GPUs can perform 0.5 teraFLOPS with their hundred cores. The tremendous computational power of GPUs was untapped by scientific computations because it could only be accessed with difficulty until now. As reported in the Journal of Computational Chemistry, recent advances allowing GPUs to be used for general purpose computing have boosted the performance of a number of molecular modeling applications, including NAMD simulations and VMD electrostatic potential calculations. The accelerated versions of these applications run five to one hundred times faster than on the best CPU-based hardware, allowing a single desktop computer equipped with a GPU to provide processing power equivalent to an entire, large computing cluster. More information on GPU acceleration of molecular modeling applications is provided here.

The cells of higher organisms store their genetic material, the genome, in the so-called nucleus where they organize transcription of DNA into messenger-RNA, the blueprint for proteins. The messenger-RNA leaves the cell to be decoded by ribosomes that synthesize the respective proteins. Transcription factors, also proteins, control in the nucleus which parts of the cells' genomes are transcribed. Naturally, the access to the nucleus as well as exit from it must be restricted to transcription factors and related biomolecules. This is achieved by the nuclear pores, wide channels lined with brushes of polymers. The polymers are disordered proteins and prevent passage for most cellular proteins, except for so-called transport factors which bus transcription factors, messenger RNA, and certain larger biomolecule into and out of the nucleus. How transport factors are permitted to pass the nuclear pores, despite many studies, has been largely unknown. Molecular dynamics simulations, based on relevant crystallographic structures, using NAMD provided a comprehensive picture on the passage mechanism as reported recently. The simulations, analyzed with VMD, revealed that transport factors are dotted rather regularly on their surface with spots that bind to the brushes of nuclear pore proteins. While any protein may accidentally exhibit such a binding spot or two, only transport factors offer a regular pattern of such spots on their surface that apparently is their passport permitting them movement into and out of the nucleus, i.e., helping them to glide through the pores' protein brushes. More on simulations of transport factors can be found here.

Because oxygen gas is very reactive, it is frequently employed by the cell as a reagent by proteins called enzymes, which build the organic compounds that the cell needs. One such enzyme belongs to the copper amine oxidase family. These proteins transform amine-containing compounds into molecules needed by the cell, by reacting the compounds with oxygen. Researchers have long been interested in finding out how the various reagents reach the buried copper active site before the final oxidation reaction can occur. While copper amine oxidases exhibit an obvious channel for capturing the amine compounds to be modified, it had been unclear until now how oxygen molecules make their way through the enzyme. With the help of computer simulations using NAMD, researchers have identified in a recent publication, the routes taken by oxygen inside various copper amine oxidases from different species. In order to accomplish this, they analyzed simulations of the motions of four copper amine oxidases, using the VMD analysis and visualization software, which can predict the probability of finding oxygen molecules anywhere inside the simulated proteins. This analysis found numerous oxygen conduction routes inside each copper amine oxidase, i.e., oxygen can enter the protein through many routes, as it would in a sponge. More information on finding O2 migration pathways in proteins can be found here.

The five human senses are based on amazingly sensitive molecular processes: smell and taste are based on molecular recognition, hearing and touch on molecular mechanics, vision on molecular electronic excitation. Some animals have additional sensory capabilities; for example, some possess a magnetic sense used for orientation by means of the geomagnetic field. The magnetic sense has long been poorly understood since the underlying molecular process could not be identified, but recently some progress has been made. Surprisingly, animal vision has been implicated and evidence has been accumulated that animals can see the geomagnetic field. A long-hidden receptor in the eye, a protein aptly called cryptochrome, is likely involved. Unfortunately, cryptochrome exists only in minute amounts in animal eyes, e.g., those of migratory birds, so that only behavioral measurements on animals can be taken, but not physical measurements directly on cryptochromes. Fortunately, cryptochromes exist also in plants where they control hypocotyl growth inhibition in seedlings. Experimentalists have observed that cryptochrome-dependent responses in Arabidopsis thaliana seedlings are magnetic-field-dependent. Researchers have now also computationally demonstrated that light activation of plant cryptochrome is magnetic-field-dependent. A recent report showed that light excitation leads to cascading electron transfer in which electron spins are influenced by weak magnetic fields; the spin dynamics was found to influence the activation of cryptochrome. Arabidopsis thaliana cryptochrome can be produced in quantities large enough for physical measurements so that the door is now wide open for cracking the secret behind the long-mysterious magnetic sense of animals. More on our cryptochrome web site.

The nucleus is responsible for storing the genome of eukaryotic cells, isolating it from the cellular cytoplasm. Partitioning the genetic material is very important in protecting it from cellular processes or foreign molecules. However, the nucleus also needs to provide access for the rest of the cell to the information stored in the genome. Numerous nuclear pores in the nuclear envelope offer communication pathways between the nucleoplasm and cytoplasm. The pathways are restricted to so-called transport receptors, proteins that taxi molecules into and out of the nucleus. If a molecule wishes to enter or leave the nucleus, it associates with a transport receptor. The complex passes through the pore and then dissociates. The question is why transport receptors can pass the nuclear pores while other proteins cannot. The answer lies in the role of FG-repeat proteins lining the pores and filling much of their free volume. These proteins are disordered peptides, consisting of repeating phenylalanine-glycine (FG) residues separated by a sequence of hydrophilic linker residues. Only proteins that interact favorably with the FG-repeat regions can pass through, while other proteins are excluded. A recent report used molecular dynamics via NAMD to examine the way in which the transport factor NTF2 interacts with the FG-repeats. The study described binding spots for FG-repeat peptides on the surface of NTF2, confirming known binding spots discovered previously via experimental means, and suggesting the existence of further binding spots. The new binding spots may play a role in steering NTF2, upon import or export, along a particular path through the nuclear pore. See also a previous highlight from January 2006, "Gateway to the Nucleus", as well as our webpage on the nuclear pore complex.

Viruses are the cause of many human diseases, from the common cold to AIDS, and medicine is continuously searching for better ways to battle viruses through vaccination or medication. Detailed knowledge of the life cycles of viruses should be useful in the treatment of viral diseases. A key focus of investigations is the virus capsid, a protein coat that protects the viral genome, but also triggers release of the genome and other viral factors upon contact with the body's cells. X-ray crystallography has resolved the average structures of many types of virus capsids, providing the basis for detailed investigations, for example by means of molecular dynamics methods, of capsid dynamical properties, e.g., in assembly and disassembly. Unfortunately, due to their large size most virus capsids are beyond the reach of molecular dynamics simulations, with one notable exception (see the March 2006 highlight "Simulating an Entire Life Form"). This earlier simulation allowed researchers to develop and test a method for coarse-grained molecular dynamics simulations that glosses over atomic detail and, thereby, permits microsecond descriptions of entire viral particles. As reported recently (see also journal cover) such simulations, employing the program NAMD, were applied to the empty capsids of several viruses. These simulations revealed a variety of behaviors, from rapid collapse to high stability, depending on the strength of interactions between the proteins from which capsids are built. The new method offers unprecedented views of capsid dynamics that may assist in battling viral diseases. More information on the simulations can be found on our virus web page.

Most forms of life need to detect and respond to changes in their environment for survival and optimal growth. For this purpose organisms rely on receptors that are based on sensory proteins. In plants, several sensory proteins detect the ambient light for optimal exposure of their photosynthetic apparatus. One class of plant light sensors, the phototropins, influence photosynthesis and induce the transition between root and stem growth when seedlings emerge out of the ground. Induction is activated through several protein domains, two of which actually absorb light and for their sensitivity to light, oxygen, and voltage, are called LOV1 and LOV2 domains. Understanding the LOV domains' involvement in activation is important for studying the signaling mechanisms of other types of sensory proteins. Strangely, light absorbed by LOV domains is observed to lead to a distinct, but only very minute, structural change that does not explain how activation might come about. NAMD-based molecular dynamics simulations of the LOV domain have now revealed, as reported in a recent publication, that photoactivated LOV domains exhibit altered patterns of motion that can induce a signal for plant cells. More information may be found on our biological photoreceptors website.

Viruses, the cause of many diseases, are the smallest natural organisms known. They are extremely primitive and parasitic such that biologists refer to them as "particles", rather than organisms. Viruses contain in a protein shell, the capsid, their own building plan, the genome, in the form of DNA or RNA. Viruses hijack a biological cell and make it produce from one virus many new ones. Viruses have evolved elaborate mechanisms to infect host cells, to to produce and assemble their own components, and to leave the host cell when it bursts from viral overcrowding. Because of their simplicity and small size, computational biologists selected a virus for their first attempt to reverse engineer in a computer program, NAMD, an entire life form, choosing one of the tiniest viruses for this purpose, the satellite tobacco mosaic virus. As described in a recent report, the researchers simulated the virus in a small drop of salt water, altogether involving over a million atoms. This provided an unprecedented view into the dynamics of the virus for a very brief time, revealing nevertheless the key physical properties of the viral particle as well as providing crucial information on its assembly. It may take still a long time to simulate a dog wagging its tail in the computer, but a big first step has been taken to simulate living organisms. Naturally, this step will assist modern medicine (more on our satellite tobacco mosaic virus web page).

Eukaryotic cells envelop their genetic material in the cell nucleus whose boundary contains numerous pores. Only small molecules can pass through these nuclear pores unhindered. For all larger ones, passage is highly selective and controlled. The control involves import and export proteins (transport receptors) that load and release cargo on the proper side of the nucleus upon interaction with signaling proteins. Researchers are presently solving the structure of the nuclear pore and its transport receptors with increasing resolution, and the first atomic level investigation into the mechanism of nuclear pore selectivity has recently been reported [paper]. The study inspected the interaction between the transport receptor importin-β with key nuclear pore proteins that appear disordered near the center of the pore and contain characteristic phenylalanine-glycine sequence repeats. Molecular dynamics simulations using NAMD and analyzed using VMD revealed a key insight into the selectivity mechanism. The simulations showed that the key sequences of the repeat proteins interact strongly with certain spots on the surface of importin-β. The study confirmed spots that had previously been identified experimentally and, moreover, found numerous binding spots not yet seen in experiment. Further experiments and simulations promise an understanding of the selectivity of entry and exit from the nucleus, a key element of the cell's genetic control. For more information see our nuclear pore complex webpage.

An important means for generating genetic diversity to provide raw material for evolution and maintain genomic stability is sexual reproduction. At the molecular level, the genes of two individuals are mixed through a process called homologous recombination. This process is found also in many simple life forms, even bacteria. At the beginning of recombination, two DNA duplexes, e.g., from mother and father, are aligned next to one another as the result of homology search, i.e., like strands are brought together with like strands. The four single DNA strands, two in each duplex, cross reciprocally two of the strands between the duplexes. The result is a joint molecule that contains DNA crossovers, named Holliday junctions. The Holliday junction is highly polymorphic in moving along two DNA duplexes, exchanging their DNA. Researchers are now investigating the physical mechanism of Holliday junction migration. The polymorphic, dynamic character of this migration makes observations difficult and the researchers resorted to molecular dynamics simulations using NAMD. The results, reported recently, resolved the dynamics of maternal-paternal DNA exchange through Holliday junction transitions in unprecedented detail providing an atomic level view of sexual reproduction. Check a brief review on our website.

The nucleus of the cell is centrally important to an organism. It serves to store and organize genetic information, the atomic blueprint for the organism, while separating and protecting this very important information from the host of other cellular components. While the nucleus requires this protective isolation, it also needs to communicate with the rest of the cell, exchanging proteins and RNA, for a variety of nuclear and cytoplasmic processes which act in concert. The nuclear pore complex (NPC), perhaps the largest protein complex in the cell, is responsible for the protected exchange of components between the nucleus and cytoplasm and for preventing the transport of material not destined to cross the nuclear envelope. The large size of the NPC makes it difficult to study experimentally. Computational efforts can go a long way toward revealing properties of the NPC which are inaccessible by experiments. Recent molecular dynamics simulations have revealed interactions between the transport receptor importin-β and key nuclear pore proteins, bringing forth a better understanding of the selectivity of entry and exit from the nucleus.

In an optimistic future, cars and appliances will be powered by renewable energy produced by burning hydrogen gas, with water being the only waste product. To supply this hydrogen gas, scientists are turning their attention to an enzyme called hydrogenase that is found in certain microorganisms, which produce hydrogen gas from sunlight and water. This enzyme, however, is sensitive to oxygen gas, which irreversibly deactivates its hydrogen-producing active site. Understanding how oxygen reaches the active site will provide insight into how hydrogenase's oxygen tolerance can be increased through protein engineering, and in turn make hydrogenase an economical source of hydrogen fuel. In a recent paper (also described in this webpage), the programs NAMD and VMD are used to analyze the gas diffusion process inside hydrogenase, and how it correlates with the protein's internal fluctuations, thereby creating a map of the oxygen pathways. The calculations revealed two distinct pathways for oxygen to reach the active site. Gases participate in physiological processes of many organisms and the new computational strategy developed promises to image gas diffusion pathways for many relevant proteins. In fact, the researchers are currently inspecting hundreds of proteins for their ability to internally transport gas molecules.

Nuclear hormone receptors are cellular regulators which activate the transcription of specific genes in response to the binding of nuclear hormones. We studied the specificity of DNA recognition by the estrogen receptor protein. The role of water molecules at the protein-DNA interface and changes in the DNA structure between specific and non-specific binding were monitored and analyzed.

Through modelling and quantum chemical studies, the group is supporting the design of novel proteins in collaboration with M. Sarikaya, U. Washington, Seattle. It is hypothesized that electrostatic interactions between the polar residues of this genetically engineered polypeptide and the gold surface allow stronger adsorbtion onto the {111} surface than to other Au crystal faces, thus influencing the crystalization of gold in the presence of this polypeptide.

Many proteins interact with gas molecules such as oxygen to perform their functions. In most cases, the gas molecules must reach sites buried deep inside the proteins that bind the molecules, with no obvious way in. Understanding how, for example, oxygen enters the protein, and mapping out which pathways it takes has been a long-standing challenge. As reported recently, computational biologists, inspired by previous work on the hydrogenase enzyme (see the September 2005 highlight), have developed a method, called implicit ligand sampling, that maps the pathways taken by gas molecules inside proteins. The mapping is determined by monitoring fluctuations of the protein, surprisingly, in the absence of the gas molecules. The mapping method is available in the most recent version of the program VMD used for structure and sequence analysis of proteins. The researchers applied the method to myoglobin, an oxygen-storing protein present in muscle cells, and determined detailed three-dimensional maps of oxygen and carbon monoxide pathways inside the protein (for more information see our web page). While some details of these pathways were already known from experiment, the implicit ligand maps revealed a large number of new pathways and suggest that oxygen enters myoglobin using many different entrance doors.

Most classical molecular dynamics (MD) simulations employ potential functions that do not account for the effects of induced electronic polarization between atoms, instead treating atoms as simple fixed point charges. Incorporating the influence of polarization in large-scale simulations is a critical challenge in the progress toward computations of increased fidelity, providing a more realistic and accurate representation of microscopic and thermodynamic properties. The Drude oscillator model represents induced electronic polarization by introducing an auxiliary particle attached to each polarizable atom via a harmonic spring. The advantage with the Drude model is that it preserves the simple particle-particle Coulomb electrostatic interaction employed in nonpolarizable force fields, thus its implementation in NAMD is more straightforward than alternative models for polarization. Performance results, reported in a recent paper, show that the implementation of the Drude model maintains good parallel scalability, with an increase in computational cost by not more than twice that of using a nonpolarizable force field. More details are available on the research webpage.