We explore ultraviolet switch models for the dispersal of circumstellar discs in T Tauri stars that involve both photoevaporation by the central star and viscous evolution. We show that in combination these processes generate the observed ‘two‐time‐scale’ behaviour for the dispersal of such discs, whereby the disc is rapidly dispersed at the end of its life on a time‐scale that is a small fraction of the disc lifetime. This switch is activated when the accretion rate through the disc declines to a low level (a few times 10−10 M⊙ yr−1) such that it roughly matches the rate of photoevaporative mass loss from the disc outside . At this point, the inner disc is deprived of further replenishment from larger radii and empties on to the central star on its own short viscous time‐scale. This causes the rapid (∼105‐yr) decline in accretion rate on to the central star and in all disc‐related emission shortward of 10 μm. We discuss the implications of this model for the detection of millimetre emission around weak‐line T Tauri stars, and also point out the consequences of such a sudden draining for planet formation in the inner regions of circumstellar discs.

Single-molecule force experiments in vitro enable the characterization of the mechanical response of biological matter at the nanometer scale. However, they do not reveal the molecular mechanisms underlying mechanical function. These can only be readily studied through molecular dynamics simulations of atomic structural models: "in silico" (by computer analysis) single-molecule experiments. Steered molecular dynamics simulations, in which external forces are used to explore the response and function of macromolecules, have become a powerful tool complementing and guiding in vitro single-molecule experiments. The insights provided by in silico experiments are illustrated here through a review of recent research in three areas of protein mechanics: elasticity of the muscle protein titin and the extracellular matrix protein fibronectin; linker-mediated elasticity of the cytoskeleton protein spectrin; and elasticity of ankyrin repeats, a protein module found ubiquitously in cells but with an as-yet unclear function.

The proteins that form the permeation pathway of mechanosensory transduction channels in inner-ear hair cells have not been definitively identified. Genetic, anatomical, and physiological evidence support a role for transmembrane channel-like protein (TMC) 1 in hair cell sensory transduction, yet the molecular function of TMC proteins remains unclear. Here, we provide biochemical evidence suggesting TMC1 assembles as a dimer, along with structural and sequence analyses suggesting similarity to dimeric TMEM16 channels. To identify the pore region of TMC1, we used cysteine mutagenesis and expressed mutant TMC1 in hair cells of Tmc1/2-null mice. Cysteine-modification reagents rapidly and irreversibly altered permeation properties of mechanosensory transduction. We propose that TMC1 is structurally similar to TMEM16 channels and includes ten transmembrane domains with four domains, S4-S7, that line the channel pore. The data provide compelling evidence that TMC1 is a pore-forming component of sensory transduction channels in auditory and vestibular hair cells.

Hearing and balance use hair cells in the inner ear to transform mechanical stimuli into electrical signals1. Mechanical force from sound waves or head movements is conveyed to hair-cell transduction channels by tip links2,3, fine filaments formed by two atypical cadherins: protocadherin-15 and cadherin-234,5. These two proteins are products of deafness genes6–10 and feature long extracellular domains that interact tip-to-tip5,11 in a Ca2+-dependent manner. However, the molecular architecture of the complex is unknown. Here we combine crystallography, molecular dynamics simulations, and binding experiments to characterize the cadherin-23 and protocadherin-15 bond. We find a unique cadherin interaction mechanism, with the two most N-terminal cadherin repeats (EC1+2) of each protein interacting to form an overlapped, antiparallel heterodimer. Simulations predict that this tip-link bond is mechanically strong enough to resist forces in hair cells. In addition, the complex becomes unstable upon Ca2+ removal due to increased flexure of Ca2+-free cadherin repeats. Finally, we use structures and biochemical measurements to understand molecular mechanisms by which deafness mutations disrupt tip-link function. Overall, our results shed light on the molecular mechanics of hair-cell sensory transduction and on new interaction mechanisms for cadherins, a large protein family implicated in tissue and organ morphogenesis12,13, neural connectivity14, and cancer15.

Summary All-atom molecular dynamics simulations have become increasingly popular as a tool to investigate protein function and dynamics. However, researchers are concerned about the short time scales covered by simulations, the apparent impossibility to model large and integral biomolecular systems, and the actual predictive power of the molecular dynamics methodology. Here we review simulations that were in the past both hotly disputed and considered key successes, namely of proteins with mainly mechanical functions (titin, fibrinogen, ankyrin, and cadherin). The simulation work covered shows how state-of-the-art modeling alleviates some of the prior concerns, and how unrefuted discoveries are made through the “computational microscope".

Polyethylene terephthalate (PET) is one of the most-consumed synthetic polymers, with an annual production of 50 million tons. Unfortunately, PET accumulates as waste and is highly resistant to biodegradation. Recently, fungal and bacterial thermophilic hydrolases were found to catalyze PET hydrolysis with optimal activities at high temperatures. Strikingly, an enzyme from Ideonella sakaiensis, termed PETase, was described to efficiently degrade PET at room temperature, but the molecular basis of its activity is not currently understood. Here, a crystal structure of PETase was determined at 2.02 Å resolution and employed in molecular dynamics simulations showing that the active site of PETase has higher flexibility at room temperature than its thermophilic counterparts. This flexibility is controlled by a novel disulfide bond in its active site, with its removal leading to destabilization of the catalytic triad and reduction of the hydrolase activity. Molecular docking of a model substrate predicts that PET binds to PETase in a unique and energetically favorable conformation facilitated by several residue substitutions within its active site when compared to other enzymes. These computational predictions are in excellent agreement with recent mutagenesis and PET film degradation analyses. Finally, we rationalize the increased catalytic activity of PETase at room temperature through molecular dynamics simulations of enzyme-ligand complexes for PETase and other thermophilic PET-degrading enzymes at 298, 323, and 353 K. Our results reveal that both the binding pose and residue substitutions within PETase favor proximity between the catalytic residues and the labile carbonyl of the substrate at room temperature, suggesting a more favorable hydrolytic reaction. These results are valuable for enabling detailed evolutionary analysis of PET-degrading enzymes and for rational design endeavors aiming at increasing the efficiency of PETase and similar enzymes toward plastic degradation.

Advances in modern computational methods and technology make it possible to carry out extensive molecular dynamics simulations of complex membrane proteins based on detailed atomic models. The ultimate goal of such detailed simulations is to produce trajectories in which the behavior of the system is as realistic as possible. A critical aspect that requires consideration in the case of biological membrane systems is the existence of a net electric potential difference across the membrane. For meaningful computations, it is important to have well validated methodologies for incorporating the latter in molecular dynamics simulations. A widely used treatment of the membrane potential in molecular dynamics consists of applying an external uniform electric field E perpendicular to the membrane. The field acts on all charged particles throughout the simulated system, and the resulting applied membrane potential V is equal to the applied electric field times the length of the periodic cell in the direction perpendicular to the membrane. A series of test simulations based on simple membrane-slab models are carried out to clarify the consequences of the applied field. These illustrative tests demonstrate that the constant-field method is a simple and valid approach for accounting for the membrane potential in molecular dynamics studies of biomolecular systems.

Mechanosensitive channels are a class of ubiquitous membrane proteins gated by mechanical strain in the cellular membrane. MscS, the mechanosensitive channel of small conductance, is found in the inner membrane of Escherichia coli and its crystallographic structure in an open form has been recently solved. By means of molecular dynamics simulations we studied the stability of the channel conformation suggested by crystallography in a fully solvated lipid (POPC) bilayer, the combined system encompassing 224,340 atoms. When restraining the backbone of the protein, the channel remained in the open form and the simulation revealed intermittent permeation of water molecules through the channel. Abolishing the restraints under constant pressure conditions led to spontaneous closure of the transmembrane channel, whereas abolishing the restraints when surface tension (20 dyn/cm) was applied led to channel widening. The large balloon-shaped cytoplasmic domain of MscS exhibited spontaneous diffusion of ions through its side openings. Interaction between the transmembrane domain and the cytoplasmic domain of MscS was observed and involved formation of salt bridges between residues Asp62 and Arg128; this interaction may be essential for the gating of MscS. K+ and Cl- ions showed distinctively different distributions in and around the channel.